CMS 3D CMS Logo

pyquenPythiaDefault_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 
6 pyquenPythiaDefaultBlock = cms.PSet(
7  pythiaUESettingsBlock,
8  pyquenMain = cms.vstring(
9  'MSTP(122)=0', # ! no printout of Pythia initialization information hereinafter
10  'MSTP(128)=2' # ! instructs pythia not to put multiple copies of resonances in the event record
11  ),
12 
13  ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes',
14  'CKIN(3)=7.',# ! ptMin
15  'MSTP(81)=0'
16  ),
17  pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
18  'MSTU(21)=1', # ! to avoid stopping run',
19  'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
20  'MSTP(81)=0', # ! pp multiple scattering off',
21  'PMAS(5,1)=4.8', # ! b quark mass',
22  'PMAS(6,1)=175.0', # ! t quark mass'
23  'CKIN(3)=7.',# ! ptMin
24  'MSTJ(22)=2',
25  'PARJ(71)=10.', # Decays only if life time < 10mm
26  'PARP(67)=1.',
27  'PARP(82)=1.9',
28  'PARP(85)=0.33',
29  'PARP(86)=0.66',
30  'PARP(89)=1000.',
31  'PARP(91)=1.0',
32  'MSTJ(11)=3',
33  'MSTJ(22)=2'
34  ),
35 
36  ppJets = cms.vstring('MSEL=1 ! QCD hight pT processes',
37  'MSTP(122)=0' # ! no printout of Pythia initialization information hereinafter
38  ),
39  pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
40  'MSUB(12)=1', # q+qbar->q+qbar
41  'MSUB(13)=1', # q+qbar->g+g
42  'MSUB(28)=1', # q+g->q+g
43  'MSUB(53)=1', # g+g->q+qbar
44  'MSUB(68)=1' # g+g->g+g
45  ),
46  pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
47  'MSUB(18)=1', # q+qbar->gamma+gamma
48  'MSUB(29)=1', # q+g->q+gamma
49  'MSUB(114)=1', # g+g->gamma+gamma
50  'MSUB(115)=1' # g+g->g+gamma
51  ),
52 
53  pythiaWeakBosons = cms.vstring('MSUB(1)=1',
54  'MSUB(2)=1'),
55 
56  pythiaZjets = cms.vstring('MSUB(15)=1',
57  'MSUB(30)=1'),
58 
59  pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
60  'MSUB(422) = 1',
61  'MSUB(423) = 1',
62  'MSUB(424) = 1',
63  'MSUB(425) = 1',
64  'MSUB(426) = 1',
65  'MSUB(427) = 1',
66  'MSUB(428) = 1',
67  'MSUB(429) = 1',
68  'MSUB(430) = 1',
69  'MSUB(431) = 1',
70  'MSUB(432) = 1',
71  'MSUB(433) = 1',
72  'MSUB(434) = 1',
73  'MSUB(435) = 1',
74  'MSUB(436) = 1',
75  'MSUB(437) = 1',
76  'MSUB(438) = 1',
77  'MSUB(439) = 1'
78  ),
79 
80  pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
81  'MSUB(462) = 1',
82  'MSUB(463) = 1',
83  'MSUB(464) = 1',
84  'MSUB(465) = 1',
85  'MSUB(466) = 1',
86  'MSUB(467) = 1',
87  'MSUB(468) = 1',
88  'MSUB(469) = 1',
89  'MSUB(470) = 1',
90  'MSUB(471) = 1',
91  'MSUB(472) = 1',
92  'MSUB(473) = 1',
93  'MSUB(474) = 1',
94  'MSUB(475) = 1',
95  'MSUB(476) = 1',
96  'MSUB(477) = 1',
97  'MSUB(478) = 1',
98  'MSUB(479) = 1',
99  ),
100  pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16', # Matrix Elements
101  'PARP(142)=0.0119',
102  'PARP(143)=0.01',
103  'PARP(144)=0.01',
104  'PARP(145)=0.05',
105  'PARP(146)=9.28',
106  'PARP(147)=0.15',
107  'PARP(148)=0.02',
108  'PARP(149)=0.02',
109  'PARP(150)=0.085',
110  # Meson spin
111  'PARJ(13)=0.60',
112  'PARJ(14)=0.162',
113  'PARJ(15)=0.018',
114  'PARJ(16)=0.054',
115  # Polarization
116  'MSTP(145)=0',
117  'MSTP(146)=0',
118  'MSTP(147)=0',
119  'MSTP(148)=1',
120  'MSTP(149)=1',
121  # Chi_c branching ratios
122  'BRAT(861)=0.202',
123  'BRAT(862)=0.798',
124  'BRAT(1501)=0.013',
125  'BRAT(1502)=0.987',
126  'BRAT(1555)=0.356',
127  'BRAT(1556)=0.644'
128  ),
129  pythiaZtoMuons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
130  "MDME(175,1)=0", # !Z decay into u ubar,
131  "MDME(176,1)=0", # !Z decay into s sbar,
132  "MDME(177,1)=0", # !Z decay into c cbar,
133  "MDME(178,1)=0", # !Z decay into b bbar,
134  "MDME(179,1)=0", # !Z decay into t tbar,
135  "MDME(182,1)=0", # !Z decay into e- e+,
136  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
137  "MDME(184,1)=1", # !Z decay into mu- mu+,
138  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
139  "MDME(186,1)=0", # !Z decay into tau- tau+,
140  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
141  ),
142 
143 
144  pythiaZtoElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
145  "MDME(175,1)=0", # !Z decay into u ubar,
146  "MDME(176,1)=0", # !Z decay into s sbar,
147  "MDME(177,1)=0", # !Z decay into c cbar,
148  "MDME(178,1)=0", # !Z decay into b bbar,
149  "MDME(179,1)=0", # !Z decay into t tbar,
150  "MDME(182,1)=1", # !Z decay into e- e+,
151  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
152  "MDME(184,1)=0", # !Z decay into mu- mu+,
153  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
154  "MDME(186,1)=0", # !Z decay into tau- tau+,
155  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
156  ),
157 
158  pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0", # !Z decay into d dbar,
159  "MDME(175,1)=0", # !Z decay into u ubar,
160  "MDME(176,1)=0", # !Z decay into s sbar,
161  "MDME(177,1)=0", # !Z decay into c cbar,
162  "MDME(178,1)=0", # !Z decay into b bbar,
163  "MDME(179,1)=0", # !Z decay into t tbar,
164  "MDME(182,1)=1", # !Z decay into e- e+,
165  "MDME(183,1)=0", # !Z decay into nu_e nu_ebar,
166  "MDME(184,1)=1", # !Z decay into mu- mu+,
167  "MDME(185,1)=0", # !Z decay into nu_mu nu_mubar,
168  "MDME(186,1)=0", # !Z decay into tau- tau+,
169  "MDME(187,1)=0" # !Z decay into nu_tau nu_taubar
170  ),
171 
172  pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ', # switch off',
173  'BRAT(1035) = 1 ', # switch on',
174  'BRAT(1036) = 0 ', # switch off',
175  'BRAT(1037) = 0 ', # switch off',
176  'BRAT(1038) = 0 ', # switch off',
177  'BRAT(1039) = 0 ', # switch off',
178  'BRAT(1040) = 0 ', # switch off',
179  'BRAT(1041) = 0 ', # switch off',
180  'BRAT(1042) = 0 ', # switch off',
181  'MDME(1034,1) = 0 ', # switch off',
182  'MDME(1035,1) = 1 ', # switch on',
183  'MDME(1036,1) = 0 ', # switch off',
184  'MDME(1037,1) = 0 ', # switch off',
185  'MDME(1038,1) = 0 ', # switch off',
186  'MDME(1039,1) = 0 ', # switch off',
187  'MDME(1040,1) = 0 ', # switch off',
188  'MDME(1041,1) = 0 ', # switch off',
189  'MDME(1042,1) = 0 ', # switch off'
190  ),
191 
192  pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ', # switch off',
193  'BRAT(859) = 1 ', # switch on',
194  'BRAT(860) = 0 ', # switch off',
195  'MDME(858,1) = 0 ', # switch off',
196  'MDME(859,1) = 1 ', # switch on',
197  'MDME(860,1) = 0 ', # switch off'
198  ),
199 
200  pythiaMuonCandidates = cms.vstring(
201  'CKIN(3)=20',
202  'MSTJ(22)=2',
203  'PARJ(71)=40.'
204  ),
205  myParameters = cms.vstring(
206  )
207 
208  )