CMS 3D CMS Logo

hydjet2DefaultParameters_cff.py
Go to the documentation of this file.
1 import FWCore.ParameterSet.Config as cms
2 
4 
5 collisionParametersRHIC200GeV = cms.PSet(
6  fAw = cms.double(197.0), # beam/target atomic number
7  fSqrtS = cms.double(200.0), #
8  fMuB = cms.double(0.0285), # Chemical baryon potential per unit charge, fMuB [GeV]
9  fMuS = cms.double(0.007), # Chemical strangeness potential per unit charge, fMuS [GeV]
10  fMuI3 = cms.double(-0.001), # Chemical isospin potential per unit charge, fMuI3 [GeV]
11 
12 
13  fThFO = cms.double(0.1), # Temperature at thermal freeze-out, fTthFO [GeV]
14  fMu_th_pip = cms.double(0.06), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV]
15 
16 
17  fYlmax = cms.double(3.3), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
18  fUmax = cms.double(1.1), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
19  fPtmin = cms.double(3.55), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c]
20  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;
21 
22 
23  fTau = cms.double(8.), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
24  fR = cms.double(10.), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
25  fSigmaTau = cms.double(2.), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
26 
27  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
28 )
29 
30 collisionParameters2760GeV = cms.PSet(
31  fAw = cms.double(208.0), # beam/target atomic number
32  fSqrtS = cms.double(2760.0), #
33  fMuB = cms.double(0.), # Chemical baryon potential per unit charge, fMuB [GeV]
34  fMuS = cms.double(0.), # Chemical strangeness potential per unit charge, fMuS [GeV]
35  fMuI3 = cms.double(0.), # Chemical isospin potential per unit charge, fMuI3 [GeV]
36 
37 
38  fThFO = cms.double(0.105), # Temperature at thermal freeze-out, fTthFO [GeV]
39  fMu_th_pip = cms.double(0.), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV]
40 
41 
42  fYlmax = cms.double(4.5), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
43 
44  fUmax = cms.double(1.265), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
45  fPtmin = cms.double(8.2), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c]
46  fT0 = cms.double(1.), # Initial QGP temperature for central Pb+Pb collisions in mid-rapidity, fT0 [GeV]; allowed range [0.2,2.0]GeV;
47 
48 
49  fTau = cms.double(12.2), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
50  fR = cms.double(13.45), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
51  fSigmaTau = cms.double(3.5), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
52 
53  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
54 
55 )
56 
57 collisionParameters5020GeV = cms.PSet(
58  fAw = cms.double(208.0), # beam/target atomic number
59  fSqrtS = cms.double(5020.0),
60  fMuB = cms.double(0.), # Chemical baryon potential per unit charge, fMuB [GeV]
61  fMuS = cms.double(0.), # Chemical strangeness potential per unit charge, fMuS [GeV]
62  fMuI3 = cms.double(0.), # Chemical isospin potential per unit charge, fMuI3 [GeV]
63 
64 
65  fThFO = cms.double(0.105), # Temperature at thermal freeze-out, fTthFO [GeV]
66  fMu_th_pip = cms.double(0.), # Chemical potential of pi+ at thermal freeze-out, fMu_th_pip [GeV]
67 
68 
69  fYlmax = cms.double(3.99), # Maximal longitudinal flow rapidity at thermal freeze-out, fYlmax
70 
71 
72  fUmax = cms.double(1.280), # Maximal transverse flow rapidity at thermal freeze-out for central collisions, fUmax
73  fPtmin = cms.double(9.06), # Minimal pt of parton-parton scattering in PYTHIA event, fPtmin [GeV/c]
74  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;
75 
76 
77  fTau = cms.double(11.5), # Proper time proper at thermal freeze-out for central collisions, fTau [fm/c]
78  fR = cms.double(16.), # Maximal transverse radius at thermal freeze-out for central collisions, fR [fm]
79  fSigmaTau = cms.double(2.), # Duration of emission at thermal freeze-out for central collisions, fSigmaTau [fm/c]
80 
81  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
82 )
83 
84 qgpParametersLHC = cms.PSet(
85  fTau0 = cms.double(0.1), # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10)
86  fNf = cms.int32(0), # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3)
87 )
88 
89 qgpParametersRHIC = cms.PSet(
90  fTau0 = cms.double(0.4), # Proper QGP formation time in fm/c, fTau0 (0.01<fTau0<10)
91  fNf = cms.int32(2), # Number of active quark flavours in QGP, fNf (0, 1, 2 or 3)
92 )
93 
94 
95 hydjet2Parameters = cms.PSet(
96 
97  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)
98  fT = cms.double(0.165), # Temperature at chemical freeze-out, fT [GeV]
99  fMuC = cms.double(0.), # Chemical charm potential per unit charge, fMuC [GeV] (used if charm production is turned on)
100 
101  fCorrS = cms.double(1.), # Strangeness supression factor gamma_s with fCorrS value (0<fCorrS <=1, if fCorrS <= 0 then it is calculated)
102 
103 
104  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
105 
106  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
107 
108  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
109 
110 
111  fIfDeltaEpsilon = cms.double(1.), # Flag to specify fDelta and fEpsilon values, fIfDeltaEpsilon (=0 user's ones, >=1 calculated)
112  fDelta = cms.double(0.1), # Momentum azimuthal anizotropy parameter at thermal freeze-out, fDelta
113  fEpsilon = cms.double(0.05), # Spatial azimuthal anisotropy parameter at thermal freeze-out, fEpsilon
114 
115  fKeps2 = cms.double(0.34), # Parameter to specify/switch ON(>0)/OFF(0) elliptic flow fluctuations
116 
117  fKeps3 = cms.double(0.52), # Parameter to specify/switch ON(>0)/OFF(0) triangular flow fluctuations
118 
119 
120  fDecay = cms.int32(1), # Flag to switch on/off hadron decays, fDecay (=0 decays off, >=1 decays on)
121 
122 
123  fCharmProd = cms.int32(1), # Flag to include thermal charm production, fCharmProd (=0 no charm production, >=1 charm production)
124  fCorrC = cms.double(-1.), # Charmness enhancement factor gamma_c with fCorrC value (fCorrC >0, if fCorrC<0 then it is calculated)
125  fEtaType = cms.double(1.), # Flag to choose longitudinal flow rapidity distribution, fEtaType (=0 uniform, >0 Gaussian with the dispersion Ylmax)
126  fIshad = cms.int32(1), # Flag to switch on/off nuclear shadowing, fIshad (0 shadowing off, 1 shadowing on)
127  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)
128  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)
129  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).
130  embeddingMode = cms.int32(0),
131  rotateEventPlane= cms.bool(True)
132 )
133 
134 PythiaDefaultBlock = cms.PSet(
135  pythiaUESettingsBlock,
136  TDB = cms.vstring(
137  'PARJ(14)=0.' # ! replacing the same parametr from pythiaQuarkoniaSettings block to avoid producing h'_1 (pdg - 10333)
138  ),
139  hydjet2PythiaDefault = cms.vstring(
140  'MSEL=1', # ! type of hard QCD production process
141  'MSTU(21) = 1', # ! controle parameter to avoid stopping run
142  'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation'
143  'MSTP(81)=1', # ! pp multiple scattering on (UE model)
144  'MSTJ(21) = 1', # ! hadron decays on (if off - decays by FASTMC decayer)
145  'MSTP(2) = 1', # ! which order running alphaS
146  'MSTP(33) = 0', # ! inclusion of k factor in cross section (on/off)
147 
148  #Mod for proQ20
149  'parp(67)=1.', # ! ISR Q2max factor (amount of initial-state radiation)
150  'parp(82)=2.', # ! UE IR cutoff at reference ecm
151  'mstj(11)=3', # ! HAD choice of fragmentation function(s)
152 
153  #'MSTJ(22)=2', # ! particle decays if lifetime < parj(71)
154  #'PARJ(71)=10.',# ! ctau=10 mm
155  #'MSTP(52) = 1',# ! NO LAPDF
156  #'mstp(122)=0' # ! no printout of Pythia initialization information hereinafter
157  ),
158  ProQ2Otune = cms.vstring(
159  'mstp(51)=7', # ! PDF set: structure function chosen - CTEQ5M pdf
160  'mstp(3)=2', # ! QCD switch for choice of LambdaQCD
161  'parp(62)=2.9', # ! ISR IR cutoff
162  'parp(64)=0.14',# ! ISR renormalization scale prefactor
163  #'parp(67)=2.65',# ! ISR Q2max factor
164  'mstp(68)=3', # ! ISR phase space choice & ME corrections
165  'parp(71)=4.', # ! FSR Q2max factor for non-s-channel procs
166  'parj(81)=0.29',# ! FSR Lambda_QCD scale
167  'parj(82)=1.65',# ! FSR IR cutoff
168  'mstp(33)=0', # ! "K" switch for K-factor on/off & type
169  'mstp(81)=1', # ! UE model
170  #'parp(82)=1.9',# ! UE IR cutoff at reference ecm
171  'parp(89)=1800.',# ! UE IR cutoff reference ecm
172  'parp(90)=0.22',# ! UE IR cutoff ecm scaling power
173  'mstp(82)=4', # ! UE hadron transverse mass distribution
174  'parp(83)=0.83',# ! UE mass distribution parameter
175  'parp(84)=0.6', # ! UE mass distribution parameter
176  'parp(85)=0.86',# ! UE gg colour correlated fraction
177  'parp(86)=0.93',# ! UE total gg fraction
178  'mstp(91)=1', # ! BR primordial kT distribution
179  'parp(91)=2.1', # ! BR primordial kT width <|kT|>
180  'parp(93)=5.', # ! BR primordial kT UV cutoff
181  #'mstj(11)=5', # ! HAD choice of fragmentation function(s)
182  'parj(1)=0.073',# ! HAD diquark suppression
183  'parj(2)=0.2', # ! HAD strangeness suppression
184  'parj(3)=0.94', # ! HAD strange diquark suppression
185  'parj(4)=0.032',# ! HAD vector diquark suppression
186  'parj(11)=0.31',# ! HAD P(vector meson), u and d only
187  'parj(12)=0.4', # ! HAD P(vector meson), contains s
188  'parj(13)=0.54',# ! HAD P(vector meson), heavy quarks
189  'parj(21)=0.325',# ! HAD fragmentation pT
190  'parj(25)=0.63',# ! HAD eta0 suppression
191  'parj(26)=0.12',# ! HAD eta0' suppression
192  'parj(41)=0.5', # ! HAD string parameter a
193  'parj(42)=0.6', # ! HAD string parameter b
194  'parj(46)=1.', # ! HAD Lund(=0)-Bowler(=1) rQ (rc)
195  'parj(47)=0.67' # ! HAD Lund(=0)-Bowler(=1) rb
196 
197  ),
198  ppJets = cms.vstring('MSEL=1'), # ! QCD hight pT processes
199  customProcesses = cms.vstring('MSEL=0'),# ! User processes
200  pythiaJets = cms.vstring(
201  'MSUB(11)=1', # ! q+q->q+q
202  'MSUB(12)=1', # ! q+qbar->q+qbar
203  'MSUB(13)=1', # ! q+qbar->g+g
204  'MSUB(28)=1', # ! q+g->q+g
205  'MSUB(53)=1', # ! g+g->q+qbar
206  'MSUB(68)=1' # ! g+g->g+g
207  ),
208  pythiaPromptPhotons = cms.vstring(
209  'MSUB(14)=1', # ! q+qbar->g+gamma
210  'MSUB(18)=1', # ! q+qbar->gamma+gamma
211  'MSUB(29)=1', # ! q+g->q+gamma
212  'MSUB(114)=1', # ! g+g->gamma+gamma
213  'MSUB(115)=1' # ! g+g->g+gamma
214  ),
215  pythiaWeakBosons = cms.vstring(
216  'MSUB(1)=1',
217  'MSUB(2)=1'
218  ),
219  pythiaZjets = cms.vstring(
220  'MSUB(15)=1',
221  'MSUB(30)=1'
222  ),
223  pythiaCharmoniumNRQCD = cms.vstring(
224  'MSUB(421) = 1',
225  'MSUB(422) = 1',
226  'MSUB(423) = 1',
227  'MSUB(424) = 1',
228  'MSUB(425) = 1',
229  'MSUB(426) = 1',
230  'MSUB(427) = 1',
231  'MSUB(428) = 1',
232  'MSUB(429) = 1',
233  'MSUB(430) = 1',
234  'MSUB(431) = 1',
235  'MSUB(432) = 1',
236  'MSUB(433) = 1',
237  'MSUB(434) = 1',
238  'MSUB(435) = 1',
239  'MSUB(436) = 1',
240  'MSUB(437) = 1',
241  'MSUB(438) = 1',
242  'MSUB(439) = 1'
243  ),
244  pythiaBottomoniumNRQCD = cms.vstring(
245  'MSUB(461) = 1',
246  'MSUB(462) = 1',
247  'MSUB(463) = 1',
248  'MSUB(464) = 1',
249  'MSUB(465) = 1',
250  'MSUB(466) = 1',
251  'MSUB(467) = 1',
252  'MSUB(468) = 1',
253  'MSUB(469) = 1',
254  'MSUB(470) = 1',
255  'MSUB(471) = 1',
256  'MSUB(472) = 1',
257  'MSUB(473) = 1',
258  'MSUB(474) = 1',
259  'MSUB(475) = 1',
260  'MSUB(476) = 1',
261  'MSUB(477) = 1',
262  'MSUB(478) = 1',
263  'MSUB(479) = 1',
264  ),
265  pythiaQuarkoniaSettings = cms.vstring(
266  'PARP(141)=1.16', # Matrix Elements
267  'PARP(142)=0.0119',
268  'PARP(143)=0.01',
269  'PARP(144)=0.01',
270  'PARP(145)=0.05',
271  'PARP(146)=9.28',
272  'PARP(147)=0.15',
273  'PARP(148)=0.02',
274  'PARP(149)=0.02',
275  'PARP(150)=0.085',
276  # Meson spin
277  'PARJ(13)=0.60',
278  'PARJ(14)=0.162',
279  'PARJ(15)=0.018',
280  'PARJ(16)=0.054',
281  # Polarization
282  'MSTP(145)=0',
283  'MSTP(146)=0',
284  'MSTP(147)=0',
285  'MSTP(148)=1',
286  'MSTP(149)=1',
287  # Chi_c branching ratios
288  'BRAT(861)=0.202',
289  'BRAT(862)=0.798',
290  'BRAT(1501)=0.013',
291  'BRAT(1502)=0.987',
292  'BRAT(1555)=0.356',
293  'BRAT(1556)=0.644'
294  ),
295  pythiaZtoMuons = cms.vstring(
296  "MDME(174,1)=0", # !Z decay into d dbar,
297  "MDME(175,1)=0", # !Z decay into u ubar,
298  "MDME(176,1)=0", # !Z decay into s sbar,
299  "MDME(177,1)=0", # !Z decay into c cbar,
300  "MDME(178,1)=0", # !Z decay into b bbar,
301  "MDME(179,1)=0", # !Z decay into t tbar,
302  "MDME(182,1)=0", # !Z decay into e- e+,
303  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
304  "MDME(184,1)=1", # !Z decay into mu- mu+,
305  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
306  "MDME(186,1)=0", # !Z decay into tau- tau+,
307  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
308  ),
309  pythiaZtoElectrons = cms.vstring(
310  "MDME(174,1)=0", # !Z decay into d dbar,
311  "MDME(175,1)=0", # !Z decay into u ubar,
312  "MDME(176,1)=0", # !Z decay into s sbar,
313  "MDME(177,1)=0", # !Z decay into c cbar,
314  "MDME(178,1)=0", # !Z decay into b bbar,
315  "MDME(179,1)=0", # !Z decay into t tbar,
316  "MDME(182,1)=1", # !Z decay into e- e+,
317  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
318  "MDME(184,1)=0", # !Z decay into mu- mu+,
319  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
320  "MDME(186,1)=0", # !Z decay into tau- tau+,
321  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
322  ),
323  pythiaZtoMuonsAndElectrons = cms.vstring(
324  "MDME(174,1)=0", # !Z decay into d dbar,
325  "MDME(175,1)=0", # !Z decay into u ubar,
326  "MDME(176,1)=0", # !Z decay into s sbar,
327  "MDME(177,1)=0", # !Z decay into c cbar,
328  "MDME(178,1)=0", # !Z decay into b bbar,
329  "MDME(179,1)=0", # !Z decay into t tbar,
330  "MDME(182,1)=1", # !Z decay into e- e+,
331  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
332  "MDME(184,1)=1", # !Z decay into mu- mu+,
333  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
334  "MDME(186,1)=0", # !Z decay into tau- tau+,
335  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
336  ),
337  pythiaUpsilonToMuons = cms.vstring(
338  'BRAT(1034) = 0 ', # switch off',
339  'BRAT(1035) = 1 ', # switch on',
340  'BRAT(1036) = 0 ', # switch off',
341  'BRAT(1037) = 0 ', # switch off',
342  'BRAT(1038) = 0 ', # switch off',
343  'BRAT(1039) = 0 ', # switch off',
344  'BRAT(1040) = 0 ', # switch off',
345  'BRAT(1041) = 0 ', # switch off',
346  'BRAT(1042) = 0 ', # switch off',
347  'MDME(1034,1) = 0 ', # switch off',
348  'MDME(1035,1) = 1 ', # switch on',
349  'MDME(1036,1) = 0 ', # switch off',
350  'MDME(1037,1) = 0 ', # switch off',
351  'MDME(1038,1) = 0 ', # switch off',
352  'MDME(1039,1) = 0 ', # switch off',
353  'MDME(1040,1) = 0 ', # switch off',
354  'MDME(1041,1) = 0 ', # switch off',
355  'MDME(1042,1) = 0 ', # switch off'
356  ),
357  pythiaJpsiToMuons = cms.vstring(
358  'BRAT(858) = 0 ', # switch off',
359  'BRAT(859) = 1 ', # switch on',
360  'BRAT(860) = 0 ', # switch off',
361  'MDME(858,1) = 0 ', # switch off',
362  'MDME(859,1) = 1 ', # switch on',
363  'MDME(860,1) = 0 ', # switch off'
364  ),
365  pythiaBToJpsi = cms.vstring(
366  'BRAT(889) = 0 ', # switch off B0->J/Psi K0',
367  'BRAT(890) = 0 ', # switch off B0->J/Psi K*0',
368  'BRAT(934) = 0 ', # switch off B+->J/Psi K+',
369  'BRAT(935) = 0 ', # switch off B+->J/Psi K*+',
370  'BRAT(980) = 0 ', # switch off B_s0->J/Psi eta',
371  'BRAT(981) = 0 ', # switch off B_s0->J/Psi eta"',
372  'BRAT(982) = 0 ', # switch off B_s0->J/Psi phi',
373  'BRAT(1001) = 0 ', # switch off B_c+>J/Psi nu_e e+',
374  'BRAT(1003) = 0 ', # switch off B_c+>J/Psi nu_mu mu+',
375  'BRAT(1005) = 0 ', # switch off B_c+>J/Psi nu_tau tau+',
376 
377  'MDME(889,1) = 0 ', # switch off',
378  'MDME(890,1) = 0 ', # switch off',
379  'MDME(934,1) = 0 ', # switch off'
380  'MDME(935,1) = 0 ', # switch off'
381  'MDME(980,1) = 0 ', # switch off'
382  'MDME(981,1) = 0 ', # switch off'
383  'MDME(982,1) = 0 ', # switch off'
384  'MDME(1001,1) = 0 ', # switch off'
385  'MDME(1003,1) = 0 ', # switch off'
386  'MDME(1005,1) = 0 ', # switch off'
387  ),
388  pythiaXToJpsi = cms.vstring(
389  'BRAT(1228) = 0 ', # switch off Lambda_b0->J/Psi Lambda0',
390  'BRAT(1501) = 0 ', # switch off chi_0c->J/Psi gamma',
391  'BRAT(1555) = 0 ', # switch off chi_1c->J/Psi gamma',
392  'BRAT(1570) = 0 ', # switch off psi"->J/Psi pi+ pi-',
393  'BRAT(1571) = 0 ', # switch off psi"->J/Psi pi0 pi0',
394  'BRAT(1572) = 0 ', # switch off psi"->J/Psi eta',
395  'BRAT(1573) = 0 ', # switch off psi"->J/Psi pi0',
396 
397  'MDME(1228,1) = 0 ', # switch off'
398  'MDME(1501,1) = 0 ', # switch off'
399  'MDME(1555,1) = 0 ', # switch off'
400  'MDME(1570,1) = 0 ', # switch off'
401  'MDME(1571,1) = 0 ', # switch off'
402  'MDME(1572,1) = 0 ', # switch off'
403  'MDME(1573,1) = 0 ', # switch off'
404  ),
405  pythiaPromptJpsi = cms.vstring(
406  'BRAT(4285) = 1 ', # switch off cc~[3S18]->J/Psi g',
407  'BRAT(4286) = 1 ', # switch off cc~[1S08]->J/Psi g',
408  'BRAT(4287) = 1 ', # switch off cc~[3P08]->J/Psi g',
409 
410  'MDME(4285,1) = 1 ', # switch off'
411  'MDME(4286,1) = 1 ', # switch off'
412  'MDME(4287,1) = 1 ', # switch off'
413  ),
414  pythiaMuonCandidates = cms.vstring(
415  'CKIN(3)=20',
416  'MSTJ(22)=2',
417  'PARJ(71)=40.'
418  ),
419  myParameters = cms.vstring('MDCY(310,1)=0')
420 )