/
githubmirror
/
cmssw
Обзор
Документация
Войти
/
githubmirror
/
cmssw
Код
Запросы
0
Пакеты
0
Релизы
0
Аналитика
Безопасность
master
GeneratorInterface/Hydjet2Interface/python/hydjet2DefaultParameters_cff.py
420 строк
19 KB
Andrey
Hydjet2 bugfixes for status codes and units
20 июн 2022, 16:16
20 июн 2022, 16:16
b310a57
Код
Авторство
О чём код?
import FWCore.ParameterSet.Config as cms from Configuration.Generator.PythiaUESettings_cfi import * collisionParametersRHIC200GeV = cms.PSet( fAw = cms.double(197.0), # beam/target atomic number fSqrtS = cms.double(200.0), # fMuB = cms.double(0.0285), # Chemical baryon potential per unit charge, fMuB [GeV] fMuS = cms.double(0.007), # Chemical strangeness potential per unit charge, fMuS [GeV] fMuI3 = cms.double(-0.001), # Chemical isospin potential per unit charge, fMuI3 [GeV] ### Thermodinamic parameters at thermal freez-out ### fThFO = cms.double(0.1), # Temperature at thermal freeze-out, fTthFO [GeV] fMu_th_pip = cms.double(0.06), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] ### Maximal longitudinal flow rapidity at thermal freeze-out ### fYlmax = cms.double(3.3), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax fUmax = cms.double(1.1), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax fPtmin = cms.double(3.55), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] fT0 = cms.double(0.3), # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV; ### Volume parameters at thermal freeze-out ### fTau = cms.double(8.), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c] fR = cms.double(10.), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm] fSigmaTau = cms.double(2.), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c] fWeakDecay = cms.double(0.), # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays ) collisionParameters2760GeV = cms.PSet( fAw = cms.double(208.0), # beam/target atomic number fSqrtS = cms.double(2760.0), # fMuB = cms.double(0.), # Chemical baryon potential per unit charge, fMuB [GeV] fMuS = cms.double(0.), # Chemical strangeness potential per unit charge, fMuS [GeV] fMuI3 = cms.double(0.), # Chemical isospin potential per unit charge, fMuI3 [GeV] ### Thermodinamic parameters at thermal freez-out ### fThFO = cms.double(0.105), # Temperature at thermal freeze-out, fTthFO [GeV] fMu_th_pip = cms.double(0.), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] ### Maximal longitudinal flow rapidity at thermal freeze-out ### fYlmax = cms.double(4.5), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax fUmax = cms.double(1.265), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax fPtmin = cms.double(8.2), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] fT0 = cms.double(1.), # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV; ### Volume parameters at thermal freeze-out ### fTau = cms.double(12.2), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c] fR = cms.double(13.45), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm] fSigmaTau = cms.double(3.5), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c] fWeakDecay = cms.double(0.00000000000001), # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays ) collisionParameters5020GeV = cms.PSet( fAw = cms.double(208.0), # beam/target atomic number fSqrtS = cms.double(5020.0), fMuB = cms.double(0.), # Chemical baryon potential per unit charge, fMuB [GeV] fMuS = cms.double(0.), # Chemical strangeness potential per unit charge, fMuS [GeV] fMuI3 = cms.double(0.), # Chemical isospin potential per unit charge, fMuI3 [GeV] ### Thermodinamic parameters at thermal freez-out ### fThFO = cms.double(0.105), # Temperature at thermal freeze-out, fTthFO [GeV] fMu_th_pip = cms.double(0.), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV] ### Maximal longitudinal flow rapidity at thermal freeze-out ### fYlmax = cms.double(3.99), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax fUmax = cms.double(1.280), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax fPtmin = cms.double(9.06), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c] fT0 = cms.double(1.1), # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV; ### Volume parameters at thermal freeze-out ### fTau = cms.double(11.5), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c] fR = cms.double(16.), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm] fSigmaTau = cms.double(2.), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c] fWeakDecay = cms.double(0.00000000000001), # Low decay width threshold fWeakDecay[GeV]: width<fWeakDecay decay off, width>=fDecayWidth decay on; can be used to switch off weak decays ) qgpParametersLHC = cms.PSet( fTau0 = cms.double(0.1), # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10) fNf = cms.int32(0), # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3) ) qgpParametersRHIC = cms.PSet( fTau0 = cms.double(0.4), # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10) fNf = cms.int32(2), # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3) ) hydjet2Parameters = cms.PSet( ### Thermodinamic parameters at chemical freez-out ### fTMuType = cms.double(0.), # Flag to use calculated T_ch, mu_B and mu_S as a function of fSqrtS, fTMuType (=0 user's ones, >0 calculated) fT = cms.double(0.165), # Temperature at chemical freeze-out, fT [GeV] fMuC = cms.double(0.), # Chemical charm potential per unit charge, fMuC [GeV] (used if charm production is turned on) ### Strangeness suppression factor ### fCorrS = cms.double(1.), # Strangeness supression factor gamma_s with fCorrS value (0<fCorrS <=1, if fCorrS <= 0 then it is calculated) fRhou2 = cms.double(0.), # Parameter to specify/switch ON(>0)/OFF(0) elliptic modulations of collective velocity profile for the alternative treatment of anisotropic flow fRhou3 = cms.double(0.), # Parameter to specify/switch ON(>0)/OFF(0) triangular modulations of collective velocity profile for the alternative treatment of anisotropic flow fRhou4 = cms.double(0.), # Parameter to specify/switch ON(>0)/OFF(0) quadrangular modulations of collective velocity profile for the alternative treatment of anisotropic flow ### Anizotropy parameter at thermal freeze-out ### fIfDeltaEpsilon = cms.double(1.), # Flag to specify fDelta and fEpsilon values, fIfDeltaEpsilon (=0 user's ones, >=1 calculated) fDelta = cms.double(0.1), # Momentum azimuthal anizotropy parameter at thermal freeze-out, fDelta fEpsilon = cms.double(0.05), # Spatial azimuthal anisotropy parameter at thermal freeze-out, fEpsilon fKeps2 = cms.double(0.34), # Parameter to specify/switch ON(>0)/OFF(0) elliptic flow fluctuations fKeps3 = cms.double(0.52), # Parameter to specify/switch ON(>0)/OFF(0) triangular flow fluctuations ### Decays ### fDecay = cms.int32(1), # Flag to switch on/off hadron decays, fDecay (=0 decays off, >=1 decays on) ### Charm ### fCharmProd = cms.int32(1), # Flag to include thermal charm production, fCharmProd (=0 no charm production, >=1 charm production) fCorrC = cms.double(-1.), # Charmness enhancement factor gamma_c with fCorrC value (fCorrC >0, if fCorrC<0 then it is calculated) fEtaType = cms.double(1.), # Flag to choose longitudinal flow rapidity distribution, fEtaType (=0 uniform, >0 Gaussian with the dispersion Ylmax) fIshad = cms.int32(1), # Flag to switch on/off nuclear shadowing, fIshad (0 shadowing off, 1 shadowing on) fPyhist = cms.int32(0), # Flag to suppress the output of particle history from PYTHIA, fPyhist (=1 only final state particles; =0 full particle history from PYTHIA) fIenglu = cms.int32(0), # Flag to fix type of partonic energy loss, fIenglu (0 radiative and collisional loss, 1 radiative loss only, 2 collisional loss only) fIanglu = cms.int32(0), # Flag to fix type of angular distribution of in-medium emitted gluons, fIanglu (0 small-angular, 1 wide-angular, 2 collinear). embeddingMode = cms.int32(0), rotateEventPlane= cms.bool(True) ) PythiaDefaultBlock = cms.PSet( pythiaUESettingsBlock, TDB = cms.vstring( 'PARJ(14)=0.' # ! replacing the same parametr from pythiaQuarkoniaSettings block to avoid producing h'_1 (pdg - 10333) ), hydjet2PythiaDefault = cms.vstring( 'MSEL=1', # ! type of hard QCD production process 'MSTU(21) = 1', # ! controle parameter to avoid stopping run 'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation' 'MSTP(81)=1', # ! pp multiple scattering on (UE model) 'MSTJ(21) = 1', # ! hadron decays on (if off - decays by FASTMC decayer) 'MSTP(2) = 1', # ! which order running alphaS 'MSTP(33) = 0', # ! inclusion of k factor in cross section (on/off) #Mod for proQ20 'parp(67)=1.', # ! ISR Q2max factor (amount of initial-state radiation) 'parp(82)=2.', # ! UE IR cutoff at reference ecm 'mstj(11)=3', # ! HAD choice of fragmentation function(s) #'MSTJ(22)=2', # ! particle decays if lifetime < parj(71) #'PARJ(71)=10.',# ! ctau=10 mm #'MSTP(52) = 1',# ! NO LAPDF #'mstp(122)=0' # ! no printout of Pythia initialization information hereinafter ), ProQ2Otune = cms.vstring( 'mstp(51)=7', # ! PDF set: structure function chosen - CTEQ5M pdf 'mstp(3)=2', # ! QCD switch for choice of LambdaQCD 'parp(62)=2.9', # ! ISR IR cutoff 'parp(64)=0.14',# ! ISR renormalization scale prefactor #'parp(67)=2.65',# ! ISR Q2max factor 'mstp(68)=3', # ! ISR phase space choice & ME corrections 'parp(71)=4.', # ! FSR Q2max factor for non-s-channel procs 'parj(81)=0.29',# ! FSR Lambda_QCD scale 'parj(82)=1.65',# ! FSR IR cutoff 'mstp(33)=0', # ! "K" switch for K-factor on/off & type 'mstp(81)=1', # ! UE model #'parp(82)=1.9',# ! UE IR cutoff at reference ecm 'parp(89)=1800.',# ! UE IR cutoff reference ecm 'parp(90)=0.22',# ! UE IR cutoff ecm scaling power 'mstp(82)=4', # ! UE hadron transverse mass distribution 'parp(83)=0.83',# ! UE mass distribution parameter 'parp(84)=0.6', # ! UE mass distribution parameter 'parp(85)=0.86',# ! UE gg colour correlated fraction 'parp(86)=0.93',# ! UE total gg fraction 'mstp(91)=1', # ! BR primordial kT distribution 'parp(91)=2.1', # ! BR primordial kT width <|kT|> 'parp(93)=5.', # ! BR primordial kT UV cutoff #'mstj(11)=5', # ! HAD choice of fragmentation function(s) 'parj(1)=0.073',# ! HAD diquark suppression 'parj(2)=0.2', # ! HAD strangeness suppression 'parj(3)=0.94', # ! HAD strange diquark suppression 'parj(4)=0.032',# ! HAD vector diquark suppression 'parj(11)=0.31',# ! HAD P(vector meson), u and d only 'parj(12)=0.4', # ! HAD P(vector meson), contains s 'parj(13)=0.54',# ! HAD P(vector meson), heavy quarks 'parj(21)=0.325',# ! HAD fragmentation pT 'parj(25)=0.63',# ! HAD eta0 suppression 'parj(26)=0.12',# ! HAD eta0' suppression 'parj(41)=0.5', # ! HAD string parameter a 'parj(42)=0.6', # ! HAD string parameter b 'parj(46)=1.', # ! HAD Lund(=0)-Bowler(=1) rQ (rc) 'parj(47)=0.67' # ! HAD Lund(=0)-Bowler(=1) rb ), ppJets = cms.vstring('MSEL=1'), # ! QCD hight pT processes customProcesses = cms.vstring('MSEL=0'),# ! User processes pythiaJets = cms.vstring( 'MSUB(11)=1', # ! q+q->q+q 'MSUB(12)=1', # ! q+qbar->q+qbar 'MSUB(13)=1', # ! q+qbar->g+g 'MSUB(28)=1', # ! q+g->q+g 'MSUB(53)=1', # ! g+g->q+qbar 'MSUB(68)=1' # ! g+g->g+g ), pythiaPromptPhotons = cms.vstring( 'MSUB(14)=1', # ! q+qbar->g+gamma 'MSUB(18)=1', # ! q+qbar->gamma+gamma 'MSUB(29)=1', # ! q+g->q+gamma 'MSUB(114)=1', # ! g+g->gamma+gamma 'MSUB(115)=1' # ! g+g->g+gamma ), pythiaWeakBosons = cms.vstring( 'MSUB(1)=1', 'MSUB(2)=1' ), pythiaZjets = cms.vstring( 'MSUB(15)=1', 'MSUB(30)=1' ), pythiaCharmoniumNRQCD = cms.vstring( 'MSUB(421) = 1', 'MSUB(422) = 1', 'MSUB(423) = 1', 'MSUB(424) = 1', 'MSUB(425) = 1', 'MSUB(426) = 1', 'MSUB(427) = 1', 'MSUB(428) = 1', 'MSUB(429) = 1', 'MSUB(430) = 1', 'MSUB(431) = 1', 'MSUB(432) = 1', 'MSUB(433) = 1', 'MSUB(434) = 1', 'MSUB(435) = 1', 'MSUB(436) = 1', 'MSUB(437) = 1', 'MSUB(438) = 1', 'MSUB(439) = 1' ), pythiaBottomoniumNRQCD = cms.vstring( 'MSUB(461) = 1', 'MSUB(462) = 1', 'MSUB(463) = 1', 'MSUB(464) = 1', 'MSUB(465) = 1', 'MSUB(466) = 1', 'MSUB(467) = 1', 'MSUB(468) = 1', 'MSUB(469) = 1', 'MSUB(470) = 1', 'MSUB(471) = 1', 'MSUB(472) = 1', 'MSUB(473) = 1', 'MSUB(474) = 1', 'MSUB(475) = 1', 'MSUB(476) = 1', 'MSUB(477) = 1', 'MSUB(478) = 1', 'MSUB(479) = 1', ), pythiaQuarkoniaSettings = cms.vstring( 'PARP(141)=1.16', # Matrix Elements 'PARP(142)=0.0119', 'PARP(143)=0.01', 'PARP(144)=0.01', 'PARP(145)=0.05', 'PARP(146)=9.28', 'PARP(147)=0.15', 'PARP(148)=0.02', 'PARP(149)=0.02', 'PARP(150)=0.085', # Meson spin 'PARJ(13)=0.60', 'PARJ(14)=0.162', 'PARJ(15)=0.018', 'PARJ(16)=0.054', # Polarization 'MSTP(145)=0', 'MSTP(146)=0', 'MSTP(147)=0', 'MSTP(148)=1', 'MSTP(149)=1', # Chi_c branching ratios 'BRAT(861)=0.202', 'BRAT(862)=0.798', 'BRAT(1501)=0.013', 'BRAT(1502)=0.987', 'BRAT(1555)=0.356', 'BRAT(1556)=0.644' ), pythiaZtoMuons = cms.vstring( "MDME(174,1)=0", # !Z decay into d dbar, "MDME(175,1)=0", # !Z decay into u ubar, "MDME(176,1)=0", # !Z decay into s sbar, "MDME(177,1)=0", # !Z decay into c cbar, "MDME(178,1)=0", # !Z decay into b bbar, "MDME(179,1)=0", # !Z decay into t tbar, "MDME(182,1)=0", # !Z decay into e- e+, "MDME(183,1)=0", # !Z decay into nu_e nu_ebar, "MDME(184,1)=1", # !Z decay into mu- mu+, "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar, "MDME(186,1)=0", # !Z decay into tau- tau+, "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar ), pythiaZtoElectrons = cms.vstring( "MDME(174,1)=0", # !Z decay into d dbar, "MDME(175,1)=0", # !Z decay into u ubar, "MDME(176,1)=0", # !Z decay into s sbar, "MDME(177,1)=0", # !Z decay into c cbar, "MDME(178,1)=0", # !Z decay into b bbar, "MDME(179,1)=0", # !Z decay into t tbar, "MDME(182,1)=1", # !Z decay into e- e+, "MDME(183,1)=0", # !Z decay into nu_e nu_ebar, "MDME(184,1)=0", # !Z decay into mu- mu+, "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar, "MDME(186,1)=0", # !Z decay into tau- tau+, "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar ), pythiaZtoMuonsAndElectrons = cms.vstring( "MDME(174,1)=0", # !Z decay into d dbar, "MDME(175,1)=0", # !Z decay into u ubar, "MDME(176,1)=0", # !Z decay into s sbar, "MDME(177,1)=0", # !Z decay into c cbar, "MDME(178,1)=0", # !Z decay into b bbar, "MDME(179,1)=0", # !Z decay into t tbar, "MDME(182,1)=1", # !Z decay into e- e+, "MDME(183,1)=0", # !Z decay into nu_e nu_ebar, "MDME(184,1)=1", # !Z decay into mu- mu+, "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar, "MDME(186,1)=0", # !Z decay into tau- tau+, "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar ), pythiaUpsilonToMuons = cms.vstring( 'BRAT(1034) = 0 ', # switch off', 'BRAT(1035) = 1 ', # switch on', 'BRAT(1036) = 0 ', # switch off', 'BRAT(1037) = 0 ', # switch off', 'BRAT(1038) = 0 ', # switch off', 'BRAT(1039) = 0 ', # switch off', 'BRAT(1040) = 0 ', # switch off', 'BRAT(1041) = 0 ', # switch off', 'BRAT(1042) = 0 ', # switch off', 'MDME(1034,1) = 0 ', # switch off', 'MDME(1035,1) = 1 ', # switch on', 'MDME(1036,1) = 0 ', # switch off', 'MDME(1037,1) = 0 ', # switch off', 'MDME(1038,1) = 0 ', # switch off', 'MDME(1039,1) = 0 ', # switch off', 'MDME(1040,1) = 0 ', # switch off', 'MDME(1041,1) = 0 ', # switch off', 'MDME(1042,1) = 0 ', # switch off' ), pythiaJpsiToMuons = cms.vstring( 'BRAT(858) = 0 ', # switch off', 'BRAT(859) = 1 ', # switch on', 'BRAT(860) = 0 ', # switch off', 'MDME(858,1) = 0 ', # switch off', 'MDME(859,1) = 1 ', # switch on', 'MDME(860,1) = 0 ', # switch off' ), pythiaBToJpsi = cms.vstring( 'BRAT(889) = 0 ', # switch off B0->J/Psi K0', 'BRAT(890) = 0 ', # switch off B0->J/Psi K*0', 'BRAT(934) = 0 ', # switch off B+->J/Psi K+', 'BRAT(935) = 0 ', # switch off B+->J/Psi K*+', 'BRAT(980) = 0 ', # switch off B_s0->J/Psi eta', 'BRAT(981) = 0 ', # switch off B_s0->J/Psi eta"', 'BRAT(982) = 0 ', # switch off B_s0->J/Psi phi', 'BRAT(1001) = 0 ', # switch off B_c+>J/Psi nu_e e+', 'BRAT(1003) = 0 ', # switch off B_c+>J/Psi nu_mu mu+', 'BRAT(1005) = 0 ', # switch off B_c+>J/Psi nu_tau tau+', 'MDME(889,1) = 0 ', # switch off', 'MDME(890,1) = 0 ', # switch off', 'MDME(934,1) = 0 ', # switch off' 'MDME(935,1) = 0 ', # switch off' 'MDME(980,1) = 0 ', # switch off' 'MDME(981,1) = 0 ', # switch off' 'MDME(982,1) = 0 ', # switch off' 'MDME(1001,1) = 0 ', # switch off' 'MDME(1003,1) = 0 ', # switch off' 'MDME(1005,1) = 0 ', # switch off' ), pythiaXToJpsi = cms.vstring( 'BRAT(1228) = 0 ', # switch off Lambda_b0->J/Psi Lambda0', 'BRAT(1501) = 0 ', # switch off chi_0c->J/Psi gamma', 'BRAT(1555) = 0 ', # switch off chi_1c->J/Psi gamma', 'BRAT(1570) = 0 ', # switch off psi"->J/Psi pi+ pi-', 'BRAT(1571) = 0 ', # switch off psi"->J/Psi pi0 pi0', 'BRAT(1572) = 0 ', # switch off psi"->J/Psi eta', 'BRAT(1573) = 0 ', # switch off psi"->J/Psi pi0', 'MDME(1228,1) = 0 ', # switch off' 'MDME(1501,1) = 0 ', # switch off' 'MDME(1555,1) = 0 ', # switch off' 'MDME(1570,1) = 0 ', # switch off' 'MDME(1571,1) = 0 ', # switch off' 'MDME(1572,1) = 0 ', # switch off' 'MDME(1573,1) = 0 ', # switch off' ), pythiaPromptJpsi = cms.vstring( 'BRAT(4285) = 1 ', # switch off cc~[3S18]->J/Psi g', 'BRAT(4286) = 1 ', # switch off cc~[1S08]->J/Psi g', 'BRAT(4287) = 1 ', # switch off cc~[3P08]->J/Psi g', 'MDME(4285,1) = 1 ', # switch off' 'MDME(4286,1) = 1 ', # switch off' 'MDME(4287,1) = 1 ', # switch off' ), pythiaMuonCandidates = cms.vstring( 'CKIN(3)=20', 'MSTJ(22)=2', 'PARJ(71)=40.' ), myParameters = cms.vstring('MDCY(310,1)=0') )