CMS 3D CMS Logo
Pyquen2015Settings_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),
8  comEnergy = cms.double(4000.0)
9  )
10 
11 collisionParameters2760GeV = cms.PSet(aBeamTarget = cms.double(208.0),
12  comEnergy = cms.double(2760.0)
13  )
14 
15 collisionParameters5020GeV = cms.PSet(aBeamTarget = cms.double(208.0),
16  comEnergy = cms.double(5020.0)
17  )
18 
19 collisionParameters5362GeV = cms.PSet(aBeamTarget = cms.double(208.0),
20  comEnergy = cms.double(5362.0)
21  )
22 
23 collisionParameters5519GeV = cms.PSet(aBeamTarget = cms.double(208.0),
24  comEnergy = cms.double(5519.0)
25  )
26 
27 qgpParameters = cms.PSet(qgpInitialTemperature = cms.double(1.1),
28  qgpProperTimeFormation = cms.double(0.1),
29  hadronFreezoutTemperature = cms.double(0.125),
30  doRadiativeEnLoss = cms.bool(True),
31  doCollisionalEnLoss = cms.bool(True),
32  qgpNumQuarkFlavor = cms.int32(0),
33  )
34 
35 qgpParameters2023 = cms.PSet(qgpInitialTemperature = cms.double(1.0),
36  qgpProperTimeFormation = cms.double(0.1),
37  hadronFreezoutTemperature = cms.double(0.125),
38  doRadiativeEnLoss = cms.bool(True),
39  doCollisionalEnLoss = cms.bool(True),
40  qgpNumQuarkFlavor = cms.int32(0),
41  )
42 
43 pyquenParameters = cms.PSet(doIsospin = cms.bool(True),
44  angularSpectrumSelector = cms.int32(0),
45  embeddingMode = cms.int32(0),
46  )
47 
48 hydjetParameters = cms.PSet(sigmaInelNN = cms.double(70),
49  shadowingSwitch = cms.int32(1),
50  nMultiplicity = cms.int32(25000),
51  fracSoftMultiplicity = cms.double(1.),
52  maxLongitudinalRapidity = cms.double(4.5),
53  maxTransverseRapidity = cms.double(1.25),
54  rotateEventPlane = cms.bool(True),
55  allowEmptyEvents = cms.bool(False),
56  angularSpectrumSelector = cms.int32(1),
57  embeddingMode = cms.int32(0)
58  )
59 
60 hydjetParameters2023 = cms.PSet(sigmaInelNN = cms.double(70),
61  shadowingSwitch = cms.int32(1),
62  nMultiplicity = cms.int32(9000),
63  fracSoftMultiplicity = cms.double(1.),
64  maxLongitudinalRapidity = cms.double(2.3),
65  maxTransverseRapidity = cms.double(1.15),
66  rotateEventPlane = cms.bool(True),
67  allowEmptyEvents = cms.bool(False),
68  angularSpectrumSelector = cms.int32(0),
69  embeddingMode = cms.int32(0)
70  )
71 
72 pyquenPythiaDefaultBlock = cms.PSet(
73  pythiaUESettingsBlock,
74  hydjetPythiaDefault = cms.vstring('MSEL=0 ! user processes',
75  'CKIN(3)=10.',# ! ptMin
76  'MSTP(81)=1'
77  ),
78  ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes (only jets)',
79  'CKIN(3)=10.',# ! ptMin
80  'MSTP(81)=1'
81  ),
82  pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
83  'MSTU(21)=1', # ! to avoid stopping run',
84  'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
85  'MSTP(81)=0', # ! pp multiple scattering off',
86  'PMAS(5,1)=4.8', # ! b quark mass',
87  'PMAS(6,1)=175.0', # ! t quark mass'
88  'CKIN(3)=7.',# ! ptMin
89  'MSTJ(22)=2',
90  'PARJ(71)=10.', # Decays only if life time < 10mm
91  'PARP(67)=1.',
92  'PARP(82)=1.9',
93  'PARP(85)=0.33',
94  'PARP(86)=0.66',
95  'PARP(89)=1000.',
96  'PARP(91)=1.0',
97  'MSTJ(11)=3',
98  'MSTJ(22)=2'
99  ),
100  ppJets = cms.vstring('MSEL=1 ! QCD hight pT processes'),
101  customProcesses = cms.vstring('MSEL=0 ! User processes'),
102  pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
103  'MSUB(12)=1', # q+qbar->q+qbar
104  'MSUB(13)=1', # q+qbar->g+g
105  'MSUB(28)=1', # q+g->q+g
106  'MSUB(53)=1', # g+g->q+qbar
107  'MSUB(68)=1' # g+g->g+g
108  ),
109  pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
110  'MSUB(18)=1', # q+qbar->gamma+gamma
111  'MSUB(29)=1', # q+g->q+gamma
112  'MSUB(114)=1', # g+g->gamma+gamma
113  'MSUB(115)=1' # g+g->g+gamma
114  ),
115 
116  pythiaWeakBosons = cms.vstring('MSUB(1)=1',
117  'MSUB(2)=1'),
118 
119  pythiaZjets = cms.vstring('MSUB(15)=1',
120  'MSUB(30)=1'),
121 
122  pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
123  'MSUB(422) = 1',
124  'MSUB(423) = 1',
125  'MSUB(424) = 1',
126  'MSUB(425) = 1',
127  'MSUB(426) = 1',
128  'MSUB(427) = 1',
129  'MSUB(428) = 1',
130  'MSUB(429) = 1',
131  'MSUB(430) = 1',
132  'MSUB(431) = 1',
133  'MSUB(432) = 1',
134  'MSUB(433) = 1',
135  'MSUB(434) = 1',
136  'MSUB(435) = 1',
137  'MSUB(436) = 1',
138  'MSUB(437) = 1',
139  'MSUB(438) = 1',
140  'MSUB(439) = 1'
141  ),
142 
143  pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
144  'MSUB(462) = 1',
145  'MSUB(463) = 1',
146  'MSUB(464) = 1',
147  'MSUB(465) = 1',
148  'MSUB(466) = 1',
149  'MSUB(467) = 1',
150  'MSUB(468) = 1',
151  'MSUB(469) = 1',
152  'MSUB(470) = 1',
153  'MSUB(471) = 1',
154  'MSUB(472) = 1',
155  'MSUB(473) = 1',
156  'MSUB(474) = 1',
157  'MSUB(475) = 1',
158  'MSUB(476) = 1',
159  'MSUB(477) = 1',
160  'MSUB(478) = 1',
161  'MSUB(479) = 1',
162  ),
163 
164  pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16', # Matrix Elements
165  'PARP(142)=0.0119',
166  'PARP(143)=0.01',
167  'PARP(144)=0.01',
168  'PARP(145)=0.05',
169  'PARP(146)=9.28',
170  'PARP(147)=0.15',
171  'PARP(148)=0.02',
172  'PARP(149)=0.02',
173  'PARP(150)=0.085',
174  # Meson spin
175  'PARJ(13)=0.60',
176  'PARJ(14)=0.162',
177  'PARJ(15)=0.018',
178  'PARJ(16)=0.054',
179  # Polarization
180  'MSTP(145)=0',
181  'MSTP(146)=0',
182  'MSTP(147)=0',
183  'MSTP(148)=1',
184  'MSTP(149)=1',
185  # Chi_c branching ratios
186  'BRAT(861)=0.202',
187  'BRAT(862)=0.798',
188  'BRAT(1501)=0.013',
189  'BRAT(1502)=0.987',
190  'BRAT(1555)=0.356',
191  'BRAT(1556)=0.644'
192  ),
193 
194  pythiaZtoMuons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
195  "MDME(175,1)=0", # !Z decay into u ubar,
196  "MDME(176,1)=0", # !Z decay into s sbar,
197  "MDME(177,1)=0", # !Z decay into c cbar,
198  "MDME(178,1)=0", # !Z decay into b bbar,
199  "MDME(179,1)=0", # !Z decay into t tbar,
200  "MDME(182,1)=0", # !Z decay into e- e+,
201  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
202  "MDME(184,1)=1", # !Z decay into mu- mu+,
203  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
204  "MDME(186,1)=0", # !Z decay into tau- tau+,
205  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
206  ),
207 
208 
209  pythiaZtoElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
210  "MDME(175,1)=0", # !Z decay into u ubar,
211  "MDME(176,1)=0", # !Z decay into s sbar,
212  "MDME(177,1)=0", # !Z decay into c cbar,
213  "MDME(178,1)=0", # !Z decay into b bbar,
214  "MDME(179,1)=0", # !Z decay into t tbar,
215  "MDME(182,1)=1", # !Z decay into e- e+,
216  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
217  "MDME(184,1)=0", # !Z decay into mu- mu+,
218  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
219  "MDME(186,1)=0", # !Z decay into tau- tau+,
220  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
221  ),
222 
223  pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
224  "MDME(175,1)=0", # !Z decay into u ubar,
225  "MDME(176,1)=0", # !Z decay into s sbar,
226  "MDME(177,1)=0", # !Z decay into c cbar,
227  "MDME(178,1)=0", # !Z decay into b bbar,
228  "MDME(179,1)=0", # !Z decay into t tbar,
229  "MDME(182,1)=1", # !Z decay into e- e+,
230  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
231  "MDME(184,1)=1", # !Z decay into mu- mu+,
232  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
233  "MDME(186,1)=0", # !Z decay into tau- tau+,
234  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
235  ),
236 
237  pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ', # switch off',
238  'BRAT(1035) = 1 ', # switch on',
239  'BRAT(1036) = 0 ', # switch off',
240  'BRAT(1037) = 0 ', # switch off',
241  'BRAT(1038) = 0 ', # switch off',
242  'BRAT(1039) = 0 ', # switch off',
243  'BRAT(1040) = 0 ', # switch off',
244  'BRAT(1041) = 0 ', # switch off',
245  'BRAT(1042) = 0 ', # switch off',
246  'MDME(1034,1) = 0 ', # switch off',
247  'MDME(1035,1) = 1 ', # switch on',
248  'MDME(1036,1) = 0 ', # switch off',
249  'MDME(1037,1) = 0 ', # switch off',
250  'MDME(1038,1) = 0 ', # switch off',
251  'MDME(1039,1) = 0 ', # switch off',
252  'MDME(1040,1) = 0 ', # switch off',
253  'MDME(1041,1) = 0 ', # switch off',
254  'MDME(1042,1) = 0 ', # switch off'
255  ),
256 
257  pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ', # switch off',
258  'BRAT(859) = 1 ', # switch on',
259  'BRAT(860) = 0 ', # switch off',
260  'MDME(858,1) = 0 ', # switch off',
261  'MDME(859,1) = 1 ', # switch on',
262  'MDME(860,1) = 0 ', # switch off'
263  ),
264 
265  pythiaMuonCandidates = cms.vstring(
266  'CKIN(3)=20',
267  'MSTJ(22)=2',
268  'PARJ(71)=40.'
269  ),
270  myParameters = cms.vstring('MDCY(310,1)=0')
271 
272 )
PythiaUESettings_cfi
Generated for CMSSW Reference Manual by doxygen 1.8.14