Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Cms hi/hydjet v1.9 implementation #4962

Closed
wants to merge 24 commits into from
Closed

Cms hi/hydjet v1.9 implementation #4962

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
24 commits
Select commit Hold shift + click to select a range
f459a20
Configuration updates for Online DQM / MWGR.
borrello Jul 1, 2014
11e3097
Fix for ROOT<->PB file conversion.
Jul 8, 2014
0ef7f7c
Merge pull request #4564 from diguida/fix_fasthadd_encoding
davidlange6 Jul 9, 2014
76410e2
Merge pull request #4545 from dmitrijus/SecondFixForJulyMWGR
davidlange6 Jul 9, 2014
36c15a3
Merge pull request #4561 from alkaloge/Skims_711a
davidlange6 Jul 9, 2014
b485e9d
Fix compilation issue on gcc490.
ktf Jul 9, 2014
2e6b83c
Merge pull request #4575 from ktf/fix-gcc490
ktf Jul 9, 2014
7add039
Fix workflow 510.0
ktf Jul 9, 2014
46a2f14
Merge pull request #4577 from ktf/fix-510-71X
ktf Jul 9, 2014
8388615
New GTs for Run1 simulations, Run2 simulations, and data.
Jul 9, 2014
d1c7f42
Merge pull request #4589 from diguida/autoCond_for_712
ktf Jul 9, 2014
ee1e231
Revert "Skims under 711-Modified WMu name path"
davidlange6 Jul 11, 2014
e426b3e
Merge pull request #4613 from cms-sw/revert-4561-Skims_711a
davidlange6 Jul 11, 2014
d4d143f
additional changes for hydjet 2015 PbPb settings
Jul 14, 2014
2da5f9e
updated hydjet with fix for dimension warning
Jul 30, 2014
e02a780
Configuration updates for Online DQM / MWGR.
borrello Jul 1, 2014
7dbfa31
Fix for ROOT<->PB file conversion.
Jul 8, 2014
9864e88
New GTs for Run1 simulations, Run2 simulations, and data.
Jul 9, 2014
bde1bc0
additional changes for hydjet 2015 PbPb settings
Jul 14, 2014
b05973a
updated hydjet with fix for dimension warning
Jul 30, 2014
3676067
Merge branch 'hydjet-v1.9-implementation' of https://github.com/CmsHI…
Aug 5, 2014
f4d5ae4
Merged hydjet-v1.9-implementation from repository CmsHI
Aug 14, 2014
7430ddf
some files were already modified from 7_1. Now checkout from 7_2_x
Aug 14, 2014
ac04c9d
updated energy for pyquen settings
Aug 15, 2014
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
245 changes: 245 additions & 0 deletions Configuration/Generator/python/Pyquen2013Settings_cff.py
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,245 @@
#Default Pythia Paramters for Hydjet & Pyquen

import FWCore.ParameterSet.Config as cms

from Configuration.Generator.PythiaUESettings_cfi import *

collisionParameters4TeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
comEnergy = cms.double(4000.0)
)

collisionParameters2760GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
comEnergy = cms.double(2760.0)
)

collisionParameters5020GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
comEnergy = cms.double(5020.0)
)

qgpParameters = cms.PSet(qgpInitialTemperature = cms.double(1.0), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
hadronFreezoutTemperature = cms.double(0.125),
doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
doCollisionalEnLoss = cms.bool(True),
qgpNumQuarkFlavor = cms.int32(0), ## number of active quark flavors in qgp; allowed values: 0,1,2,3
numQuarkFlavor = cms.int32(0) ## to be removed
)

pyquenParameters = cms.PSet(doIsospin = cms.bool(True),
angularSpectrumSelector = cms.int32(0), ## angular emitted gluon spectrum :
embeddingMode = cms.bool(False),
backgroundLabel = cms.InputTag("generator") ## ineffective in no mixing
)

hydjetParameters = cms.PSet(sigmaInelNN = cms.double(64),
shadowingSwitch = cms.int32(1),
nMultiplicity = cms.int32(18500),
fracSoftMultiplicity = cms.double(1.),
maxLongitudinalRapidity = cms.double(4.5),
maxTransverseRapidity = cms.double(1.15),
rotateEventPlane = cms.bool(True),
allowEmptyEvents = cms.bool(False),
embeddingMode = cms.bool(False)
)

pyquenPythiaDefaultBlock = cms.PSet(
pythiaUESettingsBlock,
hydjetPythiaDefault = cms.vstring('MSEL=0 ! user processes',
'CKIN(3)=7.5',# ! ptMin
'MSTP(81)=1'
),
ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes (only jets)',
'CKIN(3)=7.5',# ! ptMin
'MSTP(81)=1'
),
pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
'MSTU(21)=1', # ! to avoid stopping run',
'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
'MSTP(81)=0', # ! pp multiple scattering off',
'PMAS(5,1)=4.8', # ! b quark mass',
'PMAS(6,1)=175.0', # ! t quark mass'
'CKIN(3)=7.',# ! ptMin
'MSTJ(22)=2',
'PARJ(71)=10.', # Decays only if life time < 10mm
'PARP(67)=1.',
'PARP(82)=1.9',
'PARP(85)=0.33',
'PARP(86)=0.66',
'PARP(89)=1000.',
'PARP(91)=1.0',
'MSTJ(11)=3',
'MSTJ(22)=2'
),
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'
),

pythiaMuonCandidates = cms.vstring(
'CKIN(3)=20',
'MSTJ(22)=2',
'PARJ(71)=40.'
),
myParameters = cms.vstring('MDCY(310,1)=0')

)