CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
PyquenDefaultSettings_cff.py
Go to the documentation of this file.
1 #Default Pythia Paramters for Hydjet & Pyquen
2 
3 import FWCore.ParameterSet.Config as cms
4 
5 from Configuration.Generator.PythiaUESettings_cfi import *
6 
7 collisionParameters4TeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
8  comEnergy = cms.double(4000.0)
9  )
10 
11 collisionParameters2760GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
12  comEnergy = cms.double(2760.0)
13  )
14 
15 collisionParameters = collisionParameters2760GeV.clone()
16 
17 qgpParameters = cms.PSet(qgpInitialTemperature = cms.double(1.0), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
18  qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
19  hadronFreezoutTemperature = cms.double(0.14),
20  doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
21  doCollisionalEnLoss = cms.bool(False),
22  qgpNumQuarkFlavor = cms.int32(0), ## number of active quark flavors in qgp; allowed values: 0,1,2,3
23  numQuarkFlavor = cms.int32(0) ## to be removed
24  )
25 
26 pyquenParameters = cms.PSet(doIsospin = cms.bool(True),
27  angularSpectrumSelector = cms.int32(0), ## angular emitted gluon spectrum :
28  embeddingMode = cms.bool(False),
29  )
30 
31 hydjetParameters = cms.PSet(sigmaInelNN = cms.double(58),
32  shadowingSwitch = cms.int32(0),
33  nMultiplicity = cms.int32(21500),
34  fracSoftMultiplicity = cms.double(1.),
35  maxLongitudinalRapidity = cms.double(4.5),
36  maxTransverseRapidity = cms.double(1.),
37  rotateEventPlane = cms.bool(True),
38  allowEmptyEvents = cms.bool(False),
39  embeddingMode = cms.bool(False)
40  )
41 
42 pyquenPythiaDefaultBlock = cms.PSet(
43  pythiaUESettingsBlock,
44  hydjetPythiaDefault = cms.vstring('MSEL=0 ! user processes',
45  'CKIN(3)=6.',# ! ptMin
46  'MSTP(81)=0'
47  ),
48  ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes (only jets)',
49  'CKIN(3)=6.',# ! ptMin
50  'MSTP(81)=0'
51  ),
52  pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
53  'MSTU(21)=1', # ! to avoid stopping run',
54  'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
55  'MSTP(81)=0', # ! pp multiple scattering off',
56  'PMAS(5,1)=4.8', # ! b quark mass',
57  'PMAS(6,1)=175.0', # ! t quark mass'
58  'CKIN(3)=7.',# ! ptMin
59  'MSTJ(22)=2',
60  'PARJ(71)=10.', # Decays only if life time < 10mm
61  'PARP(67)=1.',
62  'PARP(82)=1.9',
63  'PARP(85)=0.33',
64  'PARP(86)=0.66',
65  'PARP(89)=1000.',
66  'PARP(91)=1.0',
67  'MSTJ(11)=3',
68  'MSTJ(22)=2'
69  ),
70  ppJets = cms.vstring('MSEL=1 ! QCD hight pT processes'),
71  customProcesses = cms.vstring('MSEL=0 ! User processes'),
72  pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
73  'MSUB(12)=1', # q+qbar->q+qbar
74  'MSUB(13)=1', # q+qbar->g+g
75  'MSUB(28)=1', # q+g->q+g
76  'MSUB(53)=1', # g+g->q+qbar
77  'MSUB(68)=1' # g+g->g+g
78  ),
79  pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
80  'MSUB(18)=1', # q+qbar->gamma+gamma
81  'MSUB(29)=1', # q+g->q+gamma
82  'MSUB(114)=1', # g+g->gamma+gamma
83  'MSUB(115)=1' # g+g->g+gamma
84  ),
85 
86  pythiaWeakBosons = cms.vstring('MSUB(1)=1',
87  'MSUB(2)=1'),
88 
89  pythiaZjets = cms.vstring('MSUB(15)=1',
90  'MSUB(30)=1'),
91 
92  pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
93  'MSUB(422) = 1',
94  'MSUB(423) = 1',
95  'MSUB(424) = 1',
96  'MSUB(425) = 1',
97  'MSUB(426) = 1',
98  'MSUB(427) = 1',
99  'MSUB(428) = 1',
100  'MSUB(429) = 1',
101  'MSUB(430) = 1',
102  'MSUB(431) = 1',
103  'MSUB(432) = 1',
104  'MSUB(433) = 1',
105  'MSUB(434) = 1',
106  'MSUB(435) = 1',
107  'MSUB(436) = 1',
108  'MSUB(437) = 1',
109  'MSUB(438) = 1',
110  'MSUB(439) = 1'
111  ),
112 
113  pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
114  'MSUB(462) = 1',
115  'MSUB(463) = 1',
116  'MSUB(464) = 1',
117  'MSUB(465) = 1',
118  'MSUB(466) = 1',
119  'MSUB(467) = 1',
120  'MSUB(468) = 1',
121  'MSUB(469) = 1',
122  'MSUB(470) = 1',
123  'MSUB(471) = 1',
124  'MSUB(472) = 1',
125  'MSUB(473) = 1',
126  'MSUB(474) = 1',
127  'MSUB(475) = 1',
128  'MSUB(476) = 1',
129  'MSUB(477) = 1',
130  'MSUB(478) = 1',
131  'MSUB(479) = 1',
132  ),
133 
134  pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16', # Matrix Elements
135  'PARP(142)=0.0119',
136  'PARP(143)=0.01',
137  'PARP(144)=0.01',
138  'PARP(145)=0.05',
139  'PARP(146)=9.28',
140  'PARP(147)=0.15',
141  'PARP(148)=0.02',
142  'PARP(149)=0.02',
143  'PARP(150)=0.085',
144  # Meson spin
145  'PARJ(13)=0.60',
146  'PARJ(14)=0.162',
147  'PARJ(15)=0.018',
148  'PARJ(16)=0.054',
149  # Polarization
150  'MSTP(145)=0',
151  'MSTP(146)=0',
152  'MSTP(147)=0',
153  'MSTP(148)=1',
154  'MSTP(149)=1',
155  # Chi_c branching ratios
156  'BRAT(861)=0.202',
157  'BRAT(862)=0.798',
158  'BRAT(1501)=0.013',
159  'BRAT(1502)=0.987',
160  'BRAT(1555)=0.356',
161  'BRAT(1556)=0.644'
162  ),
163 
164  pythiaZtoMuons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
165  "MDME(175,1)=0", # !Z decay into u ubar,
166  "MDME(176,1)=0", # !Z decay into s sbar,
167  "MDME(177,1)=0", # !Z decay into c cbar,
168  "MDME(178,1)=0", # !Z decay into b bbar,
169  "MDME(179,1)=0", # !Z decay into t tbar,
170  "MDME(182,1)=0", # !Z decay into e- e+,
171  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
172  "MDME(184,1)=1", # !Z decay into mu- mu+,
173  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
174  "MDME(186,1)=0", # !Z decay into tau- tau+,
175  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
176  ),
177 
178 
179  pythiaZtoElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
180  "MDME(175,1)=0", # !Z decay into u ubar,
181  "MDME(176,1)=0", # !Z decay into s sbar,
182  "MDME(177,1)=0", # !Z decay into c cbar,
183  "MDME(178,1)=0", # !Z decay into b bbar,
184  "MDME(179,1)=0", # !Z decay into t tbar,
185  "MDME(182,1)=1", # !Z decay into e- e+,
186  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
187  "MDME(184,1)=0", # !Z decay into mu- mu+,
188  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
189  "MDME(186,1)=0", # !Z decay into tau- tau+,
190  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
191  ),
192 
193  pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
194  "MDME(175,1)=0", # !Z decay into u ubar,
195  "MDME(176,1)=0", # !Z decay into s sbar,
196  "MDME(177,1)=0", # !Z decay into c cbar,
197  "MDME(178,1)=0", # !Z decay into b bbar,
198  "MDME(179,1)=0", # !Z decay into t tbar,
199  "MDME(182,1)=1", # !Z decay into e- e+,
200  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
201  "MDME(184,1)=1", # !Z decay into mu- mu+,
202  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
203  "MDME(186,1)=0", # !Z decay into tau- tau+,
204  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
205  ),
206 
207  pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ', # switch off',
208  'BRAT(1035) = 1 ', # switch on',
209  'BRAT(1036) = 0 ', # switch off',
210  'BRAT(1037) = 0 ', # switch off',
211  'BRAT(1038) = 0 ', # switch off',
212  'BRAT(1039) = 0 ', # switch off',
213  'BRAT(1040) = 0 ', # switch off',
214  'BRAT(1041) = 0 ', # switch off',
215  'BRAT(1042) = 0 ', # switch off',
216  'MDME(1034,1) = 0 ', # switch off',
217  'MDME(1035,1) = 1 ', # switch on',
218  'MDME(1036,1) = 0 ', # switch off',
219  'MDME(1037,1) = 0 ', # switch off',
220  'MDME(1038,1) = 0 ', # switch off',
221  'MDME(1039,1) = 0 ', # switch off',
222  'MDME(1040,1) = 0 ', # switch off',
223  'MDME(1041,1) = 0 ', # switch off',
224  'MDME(1042,1) = 0 ', # switch off'
225  ),
226 
227  pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ', # switch off',
228  'BRAT(859) = 1 ', # switch on',
229  'BRAT(860) = 0 ', # switch off',
230  'MDME(858,1) = 0 ', # switch off',
231  'MDME(859,1) = 1 ', # switch on',
232  'MDME(860,1) = 0 ', # switch off'
233  ),
234 
235  pythiaMuonCandidates = cms.vstring(
236  'CKIN(3)=20',
237  'MSTJ(22)=2',
238  'PARJ(71)=40.'
239  ),
240  myParameters = cms.vstring(
241  )
242 
243 )
244 
245 # This one is not to be used
246 impactParameters = cms.PSet(cFlag = cms.int32(1),
247  bFixed = cms.double(0),
248  bMin = cms.double(0),
249  bMax = cms.double(30)
250  )
251 
PythiaUESettings_cfi
Generated for CMSSW Reference Manual by doxygen 1.8.5