How to use the PyquenInterface in CMSSW
Goal of this page
The aim of this page is to document the CMSSW interface to the Pyquen Monte-Carlo generator—where the code lives, how to run it, and details of the implementation.
General Info
CMSSW Release Status
The PyquenInterface is fully integrated into CMSSW as of release CMSSW_2_1_0.
Ongoing developments
None
Code
The PyquenInterface code lives in the package
GeneratorInterface/PyquenInterface
.
Original Generator
The present version of Pyquen used in CMSSW is version 1.5, which was released on December 19th 2007 on the
Pyquen webpage
.
Running Instructions
There is an example configuration file,
testPyquen.py
, in the test directory of the PyquenInterface package, which can be run with the following commands:
cd $CMSSW_BASE/src
cvs co GeneratorInterface/PyquenInterface
cd GeneratorInterface/PyquenInterface/test
cmsRun testPyquen.py
To customise the level of output messages, you may need to check the MessageLogger documentation at
SWGuideMessageLogger.
Official Configurations (Used in Production and CMS.RelVal)
Some configuration files containing pre-defined Pyquen settings for various scenarios exist. These are used for official Monte-Carlo production, and for producing Release Validation (CMS.RelVal) samples. These cfi files are listed in the table below.
The cmsDriver command can be used to generate full configuration files from the above cfi's which mix/embed a "signal" Pyquen event into a "background" Hydjet event as in the following example:
cmsDriver.py Pyquen_DiJet_pt80to120_4TeV_cfi.py \
-s GEN:hiSignal,SIM,DIGI,L1,DIGI2RAW,HLT:HIon \
--himix --scenario CMS.HeavyIons \
--conditions FrontierConditions_CMS.GlobalTag,MC_31X_V9::All \
--datatier 'GEN-SIM-RAW-RECO' --eventcontent=RAWDEBUGHLT \
--processName 'HISIGNAL' \
--filein=inputfile.root --fileout=outputfile.root \
-n 1 \
--no_exec
See
SWGuideHeavyIonEvtMixing for more details.
In addition, common settings for Pythia, Pyquen and Hydjet are defined in
Configuration/Generator/python/PyquenDefaultSettings_cff.py
.
Physics settings
The following sections describe which parameters need to be modified in your configuration file for various physics scenarios.
Impact parameter
(Note that although the nucleon-nucleon collisions that Pyquen generates do not have an impact parameter, the heavy-ion Hydjet events into which they are often embedded do. The correct description of centrality-dependent effects (eg quenching) requires the Pyquen event to have an impact parameter specified as well.)
process.generator.embeddingMode = True
This mode will pick up the impact parameter from the input background event.
Parameters cFlag, bFixed, bMin and bMax are irrelevant when embeddingMode = True.
- Fixed impact parameter example
process.generator.embeddingMode = False
process.generator.cFlag = 0
process.generator.bFixed = 0.2
Parameters bMin and bMax are irrelevant when cFlag = 0.
- Varying impact parameter example
process.generator.embeddingMode = False
process.generator.cFlag = 1
process.generator.bMin = 4.2
process.generator.bMax = 7.3
Please remember that the variation of impact parameter in Pyquen is not the same as that of Hydjet, in fact the distribution is scaled by Ncoll. Therefore there is no such thing as "Min Bias Pyquen events".
Parameter bFixed is irrelevant when cFlag != 0.
Quenching
process.generator.doQuench = False
- Only radiative energy loss (sometimes abbreviated noColl)
process.generator.doQuench = True
process.generator.doRadiativeEnLoss = True
process.generator.doCollisionalEnLoss = False
- Only collisional energy loss
process.generator.doQuench = True
process.generator.doCollisionalEnLoss = True
process.generator.doRadiativeEnLoss = False
- Both collisional and radiative energy loss (denoted simply as Quenched)
process.generator.doQuench = True
process.generator.doCollisionalEnLoss = True
process.generator.doRadiativeEnLoss = True
Implementation details
PYQUEN (PYthia - QUENched) is a Fortran sub-routine that applies re-scattering and energy loss to Pythia generated partons in order to simulate Quark-Gluon-Plasma induced effects. For this to work one has to set hadronization off in a Pythia run. The interface for Pyquen in CMSSW is the same as the one for Pythia, making use of HepMC::IO_HEPEVT class.
Further information
Contact
Any questions, suggestions, complains you have, please send email to
the CMS-HI software list (
hn-cms-hiswDevelopment@cernNOSPAMPLEASE.ch).
Review status
Responsible:
YetkinYilmaz