16
15
#*******************
17
16
# Running parameters
20
#*********************************************************************
21
# Tag name for the run (one word) *
22
#*********************************************************************
23
'' = run_tag ! name of the run. overwritten by the MW card
24
#*********************************************************************
25
# Run to generate the grid pack *
26
#*********************************************************************
27
.false. = gridpack !True = setting up the grid pack
28
#*********************************************************************
29
# Number of events and rnd seed *
30
#*********************************************************************
18
#*********************************************************************
20
# Warning: Do not generate more than 1M events in a single run *
21
# If you want to run Pythia, avoid more than 50k events in a run. *
22
#*********************************************************************
23
%(iseed)s = iseed ! rnd seed (0=assigned automatically=default))
32
24
#*********************************************************************
33
25
# Collider type and energy *
26
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
27
# 3=photon from electron *
34
28
#*********************************************************************
35
1 = lpp1 ! beam 1 type (0=NO PDF)
36
1 = lpp2 ! beam 2 type (0=NO PDF)
37
7000 = ebeam1 ! beam 1 energy in GeV
38
7000 = ebeam2 ! beam 2 energy in GeV
29
%(lpp1)s = lpp1 ! beam 1 type
30
%(lpp2)s = lpp2 ! beam 2 type
31
%(ebeam1)s = ebeam1 ! beam 1 total energy in GeV
32
%(ebeam2)s = ebeam2 ! beam 2 total energy in GeV
39
33
#*********************************************************************
40
34
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
41
35
#*********************************************************************
42
0 = polbeam1 ! beam polarization for beam 1
43
0 = polbeam2 ! beam polarization for beam 2
36
%(polbeam1)s = polbeam1 ! beam polarization for beam 1
37
%(polbeam2)s = polbeam2 ! beam polarization for beam 2
44
38
#*********************************************************************
45
39
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
46
40
#*********************************************************************
47
'cteq6l1' = pdlabel ! PDF set
41
%(pdlabel)s = pdlabel ! PDF set
42
%(lhaid)s = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
48
43
#*********************************************************************
49
44
# Renormalization and factorization scales *
50
45
#*********************************************************************
51
T = fixed_ren_scale ! if .true. use fixed ren scale (false is beta)
52
T = fixed_fac_scale ! if .true. use fixed fac scale (false is beta)
53
91.1880 = scale ! fixed ren scale
54
91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
55
91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
56
1 = scalefact ! scale factor for event-by-event scales
46
%(fixed_ren_scale)s = fixed_ren_scale ! if .true. use fixed ren scale
47
%(fixed_fac_scale)s = fixed_fac_scale ! if .true. use fixed fac scale
48
%(scale)s = scale ! fixed ren scale
49
%(dsqrt_q2fact1)s = dsqrt_q2fact1 ! fixed fact scale for pdf1
50
%(dsqrt_q2fact2)s = dsqrt_q2fact2 ! fixed fact scale for pdf2
51
%(scalefact)s = scalefact ! scale factor for event-by-event scales
55
##################################################################################
56
##################################################################################
57
## INFORMATION USEFULL FOR MADWEIGHT STOPS HERE ##
58
##################################################################################
59
##################################################################################
60
#**********************************
61
# BW cutoff (M+/-bwcutoff*Gamma)
62
#**********************************
63
%(bwcutoff)s = bwcutoff ! (M+/-bwcutoff*Gamma)
64
#*********************************************************************
65
# Tag name for the run (one word) *
66
#*********************************************************************
67
' ' = run_tag ! name of the run
68
%(nevents)s = nevents ! Number of unweighted events requested
69
#*********************************************************************
70
# Run to generate the grid pack *
71
#*********************************************************************
72
%(gridpack)s = gridpack !True = setting up the grid pack
57
73
#*********************************************************************
58
74
# Matching - Warning! ickkw > 1 is still beta
59
75
#*********************************************************************
60
0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
61
1 = highestmult ! for ickkw=2, highest mult group
62
1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
63
1 = alpsfact ! scale factor for QCD emission vx
64
F = chcluster ! cluster only according to channel diag
65
F = pdfwgt ! for ickkw=1, perform pdf reweighting
66
#*********************************************************************
68
#**********************************
69
# BW cutoff (M+/-bwcutoff*Gamma)
70
#**********************************
71
40 = bwcutoff ! desactivate by default in the MadWeight_card
73
# Standard Cuts ! desactivate by default in the MadWeight_card
76
#*********************************************************************
77
# Minimum and maximum pt's *
78
#*********************************************************************
79
0 = ptj ! minimum pt for the jets
80
0 = ptb ! minimum pt for the b
81
0 = pta ! minimum pt for the photons
82
0 = ptl ! minimum pt for the charged leptons
83
0 = misset ! minimum missing Et (sum of neutrino's momenta)
84
0 = ptheavy ! minimum pt for one heavy final state
85
1.0 = ptonium ! minimum pt for the quarkonium states
86
1d5 = ptjmax ! maximum pt for the jets
87
1d5 = ptbmax ! maximum pt for the b
88
1d5 = ptamax ! maximum pt for the photons
89
1d5 = ptlmax ! maximum pt for the charged leptons
90
1d5 = missetmax ! maximum missing Et (sum of neutrino's momenta)
91
#*********************************************************************
92
# Minimum and maximum E's (in the lab frame) *
93
#*********************************************************************
94
0 = ej ! minimum E for the jets
95
0 = eb ! minimum E for the b
96
0 = ea ! minimum E for the photons
97
0 = el ! minimum E for the charged leptons
98
1d5 = ejmax ! maximum E for the jets
99
1d5 = ebmax ! maximum E for the b
100
1d5 = eamax ! maximum E for the photons
101
1d5 = elmax ! maximum E for the charged leptons
102
#*********************************************************************
103
# Maximum and minimum rapidity *
104
#*********************************************************************
105
1d2 = etaj ! max rap for the jets
106
1d2 = etab ! max rap for the b
107
1d2 = etaa ! max rap for the photons
108
1d2 = etal ! max rap for the charged leptons
109
1d2 = etaonium ! max rap for the quarkonium states
110
0d0 = etajmin ! min rap for the jets
111
0d0 = etabmin ! min rap for the b
112
0d0 = etaamin ! min rap for the photons
113
0d0 = etalmin ! main rap for the charged leptons
76
%(ickkw)s = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
77
%(highestmult)s = highestmult ! for ickkw=2, highest mult group
78
%(ktscheme)s = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
79
%(alpsfact)s = alpsfact ! scale factor for QCD emission vx
80
%(chcluster)s = chcluster ! cluster only according to channel diag
81
%(pdfwgt)s = pdfwgt ! for ickkw=1, perform pdf reweighting
82
%(asrwgtflavor)s = asrwgtflavor ! highest quark flavor for a_s reweight
83
%(clusinfo)s = clusinfo ! include clustering tag in output
84
%(lhe_version)s = lhe_version ! Change the way clustering information pass to shower.
85
#*********************************************************************
86
#**********************************************************
88
#**********************************************************
89
# Automatic ptj and mjj cuts if xqcut > 0
90
# (turn off for VBF and single top processes)
91
#**********************************************************
92
%(auto_ptj_mjj)s = auto_ptj_mjj ! Automatic setting of ptj and mjj
93
#**********************************************************
95
#**********************************************************
96
# Apply pt/E/eta/dr/mij cuts on decay products or not
97
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
98
#**********************************************************
99
%(cut_decays)s = cut_decays ! Cut decay products
100
#*************************************************************
101
# Number of helicities to sum per event (0 = all helicities)
102
# 0 gives more stable result, but longer run time (needed for
103
# long decay chains e.g.).
104
# Use >=2 if most helicities contribute, e.g. pure QCD.
105
#*************************************************************
106
%(nhel)s = nhel ! Number of helicities used per event
111
#*********************************************************************
112
# Minimum and maximum pt's (for max, -1 means no cut) *
113
#*********************************************************************
114
%(ptj)s = ptj ! minimum pt for the jets
115
%(ptb)s = ptb ! minimum pt for the b
116
%(pta)s = pta ! minimum pt for the photons
117
%(ptl)s = ptl ! minimum pt for the charged leptons
118
%(misset)s = misset ! minimum missing Et (sum of neutrino's momenta)
119
%(ptheavy)s = ptheavy ! minimum pt for one heavy final state
120
%(ptjmax)s = ptjmax ! maximum pt for the jets
121
%(ptbmax)s = ptbmax ! maximum pt for the b
122
%(ptamax)s = ptamax ! maximum pt for the photons
123
%(ptlmax)s = ptlmax ! maximum pt for the charged leptons
124
%(missetmax)s = missetmax ! maximum missing Et (sum of neutrino's momenta)
125
#*********************************************************************
126
# Minimum and maximum E's (in the center of mass frame) *
127
#*********************************************************************
128
%(ej)s = ej ! minimum E for the jets
129
%(eb)s = eb ! minimum E for the b
130
%(ea)s = ea ! minimum E for the photons
131
%(el)s = el ! minimum E for the charged leptons
132
%(ejmax)s = ejmax ! maximum E for the jets
133
%(ebmax)s = ebmax ! maximum E for the b
134
%(eamax)s = eamax ! maximum E for the photons
135
%(elmax)s = elmax ! maximum E for the charged leptons
136
#*********************************************************************
137
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
138
#*********************************************************************
139
%(etaj)s = etaj ! max rap for the jets
140
%(etab)s = etab ! max rap for the b
141
%(etaa)s = etaa ! max rap for the photons
142
%(etal)s = etal ! max rap for the charged leptons
143
%(etajmin)s = etajmin ! min rap for the jets
144
%(etabmin)s = etabmin ! min rap for the b
145
%(etaamin)s = etaamin ! min rap for the photons
146
%(etalmin)s = etalmin ! main rap for the charged leptons
114
147
#*********************************************************************
115
148
# Minimum and maximum DeltaR distance *
116
149
#*********************************************************************
117
0. = drjj ! min distance between jets
118
0 = drbb ! min distance between b's
119
0. = drll ! min distance between leptons
120
0. = draa ! min distance between gammas
121
0 = drbj ! min distance between b and jet
122
0. = draj ! min distance between gamma and jet
123
0. = drjl ! min distance between jet and lepton
124
0 = drab ! min distance between gamma and b
125
0 = drbl ! min distance between b and lepton
126
0. = dral ! min distance between gamma and lepton
127
1d2 = drjjmax ! max distance between jets
128
1d2 = drbbmax ! max distance between b's
129
1d2 = drllmax ! max distance between leptons
130
1d2 = draamax ! max distance between gammas
131
1d2 = drbjmax ! max distance between b and jet
132
1d2 = drajmax ! max distance between gamma and jet
133
1d2 = drjlmax ! max distance between jet and lepton
134
1d2 = drabmax ! max distance between gamma and b
135
1d2 = drblmax ! max distance between b and lepton
136
1d2 = dralmax ! maxdistance between gamma and lepton
150
%(drjj)s = drjj ! min distance between jets
151
%(drbb)s = drbb ! min distance between b's
152
%(drll)s = drll ! min distance between leptons
153
%(draa)s = draa ! min distance between gammas
154
%(drbj)s = drbj ! min distance between b and jet
155
%(draj)s = draj ! min distance between gamma and jet
156
%(drjl)s = drjl ! min distance between jet and lepton
157
%(drab)s = drab ! min distance between gamma and b
158
%(drbl)s = drbl ! min distance between b and lepton
159
%(dral)s = dral ! min distance between gamma and lepton
160
%(drjjmax)s = drjjmax ! max distance between jets
161
%(drbbmax)s = drbbmax ! max distance between b's
162
%(drllmax)s = drllmax ! max distance between leptons
163
%(draamax)s = draamax ! max distance between gammas
164
%(drbjmax)s = drbjmax ! max distance between b and jet
165
%(drajmax)s = drajmax ! max distance between gamma and jet
166
%(drjlmax)s = drjlmax ! max distance between jet and lepton
167
%(drabmax)s = drabmax ! max distance between gamma and b
168
%(drblmax)s = drblmax ! max distance between b and lepton
169
%(dralmax)s = dralmax ! maxdistance between gamma and lepton
137
170
#*********************************************************************
138
171
# Minimum and maximum invariant mass for pairs *
172
# WARNING: for four lepton final state mmll cut require to have *
173
# different lepton masses for each flavor! *
139
174
#*********************************************************************
140
0 = mmjj ! min invariant mass of a jet pair
141
0 = mmbb ! min invariant mass of a b pair
142
0 = mmaa ! min invariant mass of gamma gamma pair
143
0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
144
1d5 = mmjjmax ! max invariant mass of a jet pair
145
1d5 = mmbbmax ! max invariant mass of a b pair
146
1d5 = mmaamax ! max invariant mass of gamma gamma pair
147
1d5 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
175
%(mmjj)s = mmjj ! min invariant mass of a jet pair
176
%(mmbb)s = mmbb ! min invariant mass of a b pair
177
%(mmaa)s = mmaa ! min invariant mass of gamma gamma pair
178
%(mmll)s = mmll ! min invariant mass of l+l- (same flavour) lepton pair
179
%(mmjjmax)s = mmjjmax ! max invariant mass of a jet pair
180
%(mmbbmax)s = mmbbmax ! max invariant mass of a b pair
181
%(mmaamax)s = mmaamax ! max invariant mass of gamma gamma pair
182
%(mmllmax)s = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
148
183
#*********************************************************************
149
184
# Minimum and maximum invariant mass for all letpons *
150
185
#*********************************************************************
151
0 = mmnl ! min invariant mass for all letpons (l+- and vl)
152
1d5 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
186
%(mmnl)s = mmnl ! min invariant mass for all letpons (l+- and vl)
187
%(mmnlmax)s = mmnlmax ! max invariant mass for all letpons (l+- and vl)
188
#*********************************************************************
189
# Minimum and maximum pt for 4-momenta sum of leptons *
190
#*********************************************************************
191
%(ptllmin)s = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
192
%(ptllmax)s = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
153
193
#*********************************************************************
154
194
# Inclusive cuts *
155
195
#*********************************************************************
156
0 = xptj ! minimum pt for at least one jet
157
0 = xptb ! minimum pt for at least one b
158
0 = xpta ! minimum pt for at least one photon
159
0 = xptl ! minimum pt for at least one charged lepton
196
%(xptj)s = xptj ! minimum pt for at least one jet
197
%(xptb)s = xptb ! minimum pt for at least one b
198
%(xpta)s = xpta ! minimum pt for at least one photon
199
%(xptl)s = xptl ! minimum pt for at least one charged lepton
160
200
#*********************************************************************
161
201
# Control the pt's of the jets sorted by pt *
162
202
#*********************************************************************
163
0 = ptj1min ! minimum pt for the leading jet in pt
164
0 = ptj2min ! minimum pt for the second jet in pt
165
0 = ptj3min ! minimum pt for the third jet in pt
166
0 = ptj4min ! minimum pt for the fourth jet in pt
167
1d5 = ptj1max ! maximum pt for the leading jet in pt
168
1d5 = ptj2max ! maximum pt for the second jet in pt
169
1d5 = ptj3max ! maximum pt for the third jet in pt
170
1d5 = ptj4max ! maximum pt for the fourth jet in pt
171
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
203
%(ptj1min)s = ptj1min ! minimum pt for the leading jet in pt
204
%(ptj2min)s = ptj2min ! minimum pt for the second jet in pt
205
%(ptj3min)s = ptj3min ! minimum pt for the third jet in pt
206
%(ptj4min)s = ptj4min ! minimum pt for the fourth jet in pt
207
%(ptj1max)s = ptj1max ! maximum pt for the leading jet in pt
208
%(ptj2max)s = ptj2max ! maximum pt for the second jet in pt
209
%(ptj3max)s = ptj3max ! maximum pt for the third jet in pt
210
%(ptj4max)s = ptj4max ! maximum pt for the fourth jet in pt
211
%(cutuse)s = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
212
#*********************************************************************
213
# Control the pt's of leptons sorted by pt *
214
#*********************************************************************
215
%(ptl1min)s = ptl1min ! minimum pt for the leading lepton in pt
216
%(ptl2min)s = ptl2min ! minimum pt for the second lepton in pt
217
%(ptl3min)s = ptl3min ! minimum pt for the third lepton in pt
218
%(ptl4min)s = ptl4min ! minimum pt for the fourth lepton in pt
219
%(ptl1max)s = ptl1max ! maximum pt for the leading lepton in pt
220
%(ptl2max)s = ptl2max ! maximum pt for the second lepton in pt
221
%(ptl3max)s = ptl3max ! maximum pt for the third lepton in pt
222
%(ptl4max)s = ptl4max ! maximum pt for the fourth lepton in pt
172
223
#*********************************************************************
173
224
# Control the Ht(k)=Sum of k leading jets *
174
225
#*********************************************************************
175
0 = htjmin ! minimum jet HT=Sum(jet pt)
176
1d5 = htjmax ! maximum jet HT=Sum(jet pt)
177
0 = ht2min ! minimum Ht for the two leading jets
178
0 = ht3min ! minimum Ht for the three leading jets
179
0 = ht4min ! minimum Ht for the four leading jets
180
1d5 = ht2max ! maximum Ht for the two leading jets
181
1d5 = ht3max ! maximum Ht for the three leading jets
182
1d5 = ht4max ! maximum Ht for the four leading jets
226
%(htjmin)s = htjmin ! minimum jet HT=Sum(jet pt)
227
%(htjmax)s = htjmax ! maximum jet HT=Sum(jet pt)
228
%(ihtmin)s = ihtmin !inclusive Ht for all partons (including b)
229
%(ihtmax)s = ihtmax !inclusive Ht for all partons (including b)
230
%(ht2min)s = ht2min ! minimum Ht for the two leading jets
231
%(ht3min)s = ht3min ! minimum Ht for the three leading jets
232
%(ht4min)s = ht4min ! minimum Ht for the four leading jets
233
%(ht2max)s = ht2max ! maximum Ht for the two leading jets
234
%(ht3max)s = ht3max ! maximum Ht for the three leading jets
235
%(ht4max)s = ht4max ! maximum Ht for the four leading jets
236
#***********************************************************************
237
# Photon-isolation cuts, according to hep-ph/9801442 *
238
# When ptgmin=0, all the other parameters are ignored *
239
# When ptgmin>0, pta and draj are not going to be used *
240
#***********************************************************************
241
%(ptgmin)s = ptgmin ! Min photon transverse momentum
242
%(r0gamma)s = R0gamma ! Radius of isolation code
243
%(xn)s = xn ! n parameter of eq.(3.4) in hep-ph/9801442
244
%(epsgamma)s = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
245
%(isoem)s = isoEM ! isolate photons from EM energy (photons and leptons)
183
246
#*********************************************************************
185
248
#*********************************************************************
186
0 = xetamin ! minimum rapidity for two jets in the WBF case
187
0 = deltaeta ! minimum rapidity for two jets in the WBF case
188
#*********************************************************************
189
# maximal pdg code for quark to be considered as a jet *
190
# otherwise b cuts are applied *
191
#*********************************************************************
249
%(xetamin)s = xetamin ! minimum rapidity for two jets in the WBF case
250
%(deltaeta)s = deltaeta ! minimum rapidity for two jets in the WBF case
251
#*********************************************************************
253
#*********************************************************************
254
%(ktdurham)s = ktdurham
255
%(dparameter)s = dparameter
256
#*********************************************************************
257
# maximal pdg code for quark to be considered as a light jet *
258
# (otherwise b cuts are applied) *
259
#*********************************************************************
260
%(maxjetflavor)s = maxjetflavor ! Maximum jet pdg code
193
261
#*********************************************************************
194
262
# Jet measure cuts *
195
263
#*********************************************************************
196
0 = xqcut ! minimum kt jet measure between partons
197
#*********************************************************************
264
%(xqcut)s = xqcut ! minimum kt jet measure between partons
265
#*********************************************************************
267
#*********************************************************************
268
# Store info for systematics studies *
269
# WARNING: If use_syst is T, matched Pythia output is *
270
# meaningful ONLY if plotted taking matchscale *
271
# reweighting into account! *
272
#*********************************************************************
273
F = use_syst ! Enable systematics studies
275
#**************************************
276
# Parameter of the systematics study
277
# will be used by SysCalc (if installed)
278
#**************************************
280
%(sys_scalefact)s = sys_scalefact # factorization/renormalization scale factor
281
%(sys_alpsfact)s = sys_alpsfact # \alpha_s emission scale factors
282
%(sys_matchscale)s = sys_matchscale # variation of merging scale
283
# PDF sets and number of members (0 or none for all members).
284
%(sys_pdf)s = sys_pdf # matching scales
285
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf