1 | <!--
|
---|
2 | #*********************************************************************
|
---|
3 | # *
|
---|
4 | # MadGraph/MadEvent *
|
---|
5 | # *
|
---|
6 | # The new web generation *
|
---|
7 | # *
|
---|
8 | # http://madgraph.hep.uiuc.edu *
|
---|
9 | # http://madgraph.phys.ucl.ac.be *
|
---|
10 | # http://madgraph.roma2.infn.it *
|
---|
11 | # *
|
---|
12 | # CP3 team Tim Stelzer *
|
---|
13 | # Universite de Louvain University of Illinois *
|
---|
14 | # *
|
---|
15 | #....................................................................*
|
---|
16 | # *
|
---|
17 | # This file contains all the information necessary to reproduce *
|
---|
18 | # the events generated: *
|
---|
19 | # *
|
---|
20 | # 1. software version *
|
---|
21 | # 2. proc_card.dat : code generation info including model *
|
---|
22 | # 3. param_card.dat : model primary parameters in the LH format *
|
---|
23 | # 4. run_card.dat : running parameters (collider and cuts) *
|
---|
24 | # 5. pythia_card.dat : present only in the pythia event file *
|
---|
25 | # 6. pgs_card.dat : present only in the pgs event file *
|
---|
26 | # *
|
---|
27 | # *
|
---|
28 | #*********************************************************************
|
---|
29 | -->
|
---|
30 | <MGVersion>
|
---|
31 | # MG/ME version : 4.2.0
|
---|
32 | # madgraph version : 4.2.0
|
---|
33 | # template version : 2.2.0
|
---|
34 | # helas version : 3.2
|
---|
35 | # model version : sm_1.1
|
---|
36 | </MGVersion>
|
---|
37 | <MGProcCard>
|
---|
38 | #*********************************************************************
|
---|
39 | # MadGraph/MadEvent *
|
---|
40 | # http://madgraph.hep.uiuc.edu *
|
---|
41 | # *
|
---|
42 | # proc_card.dat *
|
---|
43 | # *
|
---|
44 | # This file is used to generate the code for a specific process. *
|
---|
45 | # Some notation/conventions: *
|
---|
46 | # *
|
---|
47 | # 0. Do not modify the TAGS and their order. *
|
---|
48 | # 1. hash/pound is a comment. *
|
---|
49 | # 2. The charcter after the @ is used as an identifier for the class *
|
---|
50 | # of processes. It can be a single or a digit. *
|
---|
51 | # 3. The number of lines for the max couplings depends on how many *
|
---|
52 | # different classes of couplings are present in the model *
|
---|
53 | # In the SM these are just two: QED (which include EW) and QCD *
|
---|
54 | # 4. Write "end_coup" after the couplings list, *
|
---|
55 | # to tell MG that the couplings input is over. *
|
---|
56 | # 5. Write "done" after the proc list to *
|
---|
57 | # to tell MG that the proc input is over. *
|
---|
58 | # 6. Some model names available at present are: *
|
---|
59 | # sm = Standard Model *
|
---|
60 | # sm_ckm = Standard Model with Cabibbo matrix *
|
---|
61 | # mssm = Minimal Supersymmetric Standard Model *
|
---|
62 | # 2hdm = Generic Two Higgs Doublet model *
|
---|
63 | # heft = Higgs EFT (+Standard Model) *
|
---|
64 | # usrmod = User Model *
|
---|
65 | # 7. Don't leave spaces between the particles name in the *
|
---|
66 | # definition of the multiparticles. *
|
---|
67 | #*********************************************************************
|
---|
68 | #--------------------------------------------------------------------*
|
---|
69 | # Process(es) requested : mg2 input *
|
---|
70 | #--------------------------------------------------------------------*
|
---|
71 | # Begin PROCESS # This is TAG. Do not modify this line
|
---|
72 | pp>z>mu+mu- @0 # First Process
|
---|
73 | QCD=99 # Max QCD couplings
|
---|
74 | QED=2 # Max QED couplings
|
---|
75 | end_coup
|
---|
76 | pp>z>mu+mu-j @1 # First Process
|
---|
77 | QCD=99 # Max QCD couplings
|
---|
78 | QED=2 # Max QED couplings
|
---|
79 | end_coup
|
---|
80 | pp>z>mu+mu-jj@2 # First Process
|
---|
81 | QCD=99 # Max QCD couplings
|
---|
82 | QED=2 # Max QED couplings
|
---|
83 | end_coup
|
---|
84 | pp>z>mu+mu-jjj @3
|
---|
85 | QCD=99 # Max QCD couplings
|
---|
86 | QED=2 # Max QED couplings
|
---|
87 | end_coup
|
---|
88 | done
|
---|
89 | # End PROCESS # This is TAG. Do not modify this line
|
---|
90 | #--------------------------------------------------------------------*
|
---|
91 | # Model information *
|
---|
92 | #--------------------------------------------------------------------*
|
---|
93 | # Begin MODEL # This is TAG. Do not modify this line
|
---|
94 | smzerobmass
|
---|
95 | # End MODEL # This is TAG. Do not modify this line
|
---|
96 | #--------------------------------------------------------------------*
|
---|
97 | # Start multiparticle definitions *
|
---|
98 | #--------------------------------------------------------------------*
|
---|
99 | # Begin MULTIPARTICLES # This is TAG. Do not modify this line
|
---|
100 | P uu~dd~ss~cc~bb~g
|
---|
101 | J uu~dd~ss~cc~bb~g
|
---|
102 | L e+e-mu+mu-ta+ta-
|
---|
103 | # End MULTIPARTICLES # This is TAG. Do not modify this line
|
---|
104 | </MGProcCard>
|
---|
105 | <slha>
|
---|
106 | #******************************************************************
|
---|
107 | # MadGraph/MadEvent *
|
---|
108 | #******************************************************************
|
---|
109 | # Les Houches friendly file for the SM parameters of MadGraph *
|
---|
110 | # Spectrum and decay widths produced by SMCalc *
|
---|
111 | #******************************************************************
|
---|
112 | #*Please note the following IMPORTANT issues: *
|
---|
113 | # *
|
---|
114 | #0. REFRAIN from editing this file by hand! Some of the parame- *
|
---|
115 | # ters are not independent *
|
---|
116 | # (such as G_Fermi, alpha_em, sin(theta_W),MZ,MW) and serious *
|
---|
117 | # problems might be encountered (such as violation of unitarity *
|
---|
118 | # or gauge invariance). Always use a calculator. *
|
---|
119 | # *
|
---|
120 | #1. alpha_S(MZ) has been used in the calculation of the parameters*
|
---|
121 | # but, for consistency, it will be reset by madgraph to the *
|
---|
122 | # value expected IF the pdfs for collisions with hadrons are *
|
---|
123 | # used. This value is KEPT by madgraph when no pdf are used *
|
---|
124 | # lpp(i)=0 . *
|
---|
125 | # *
|
---|
126 | #2. Values of the charm and bottom kinematic (pole) masses are *
|
---|
127 | # those used in the matrix elements and phase space UNLESS they *
|
---|
128 | # are set to ZERO from the start in the model (particles.dat) *
|
---|
129 | # This happens, for example, when using 5-flavor QCD where *
|
---|
130 | # charm and bottom are treated as partons in the initial state *
|
---|
131 | # and a zero mass might be hardwired in the model definition. *
|
---|
132 | # *
|
---|
133 | #******************************************************************
|
---|
134 | Block MODSEL # Select Model
|
---|
135 | 0 1 # 0 1 = SM
|
---|
136 | Block SMINPUTS # Standard Model inputs
|
---|
137 | 1 1.32506980E+02 # alpha_em(MZ)(-1) SM MSbar
|
---|
138 | 2 1.16639000E-05 # G_Fermi
|
---|
139 | 3 1.18000000E-01 # alpha_s(MZ) SM MSbar
|
---|
140 | 4 9.11880000E+01 # Z mass (as input parameter)
|
---|
141 | Block MGYUKAWA # Yukawa masses m/v=y/sqrt(2)
|
---|
142 | # PDG YMASS
|
---|
143 | 5 4.70000000E+00 # mbottom for the Yukawa y_b
|
---|
144 | 4 1.42000000E+00 # mcharm for the Yukawa y_c
|
---|
145 | 6 1.74300000E+02 # mtop for the Yukawa y_t
|
---|
146 | 15 1.77700000E+00 # mtau for the Yukawa y_ta
|
---|
147 | Block MGCKM # CKM elements for MadGraph
|
---|
148 | 1 1 9.75000000E-01 # Vud for Cabibbo matrix
|
---|
149 | Block MASS # Mass spectrum (kinematic masses)
|
---|
150 | # PDG Mass
|
---|
151 | 6 1.74300000E+02 # top pole mass
|
---|
152 | 15 1.77700000E+00 # tau mass
|
---|
153 | 23 9.11880000E+01 # Z mass
|
---|
154 | 24 8.04190000E+01 # W mass
|
---|
155 | 25 1.20000000E+02 # H mass
|
---|
156 | # PDG Width
|
---|
157 | DECAY 6 1.50833649E+00 # top width
|
---|
158 | DECAY 23 2.44140351E+00 # Z width
|
---|
159 | DECAY 24 2.04759951E+00 # W width
|
---|
160 | DECAY 25 6.79485838E-03 # H width
|
---|
161 | # BR NDA ID1 ID2
|
---|
162 | 7.00588389E-02 2 4 -4 # BR( H -> c cbar )
|
---|
163 | 7.61074375E-01 2 5 -5 # BR( H -> b bbar )
|
---|
164 | 0.00000000E+00 2 6 -6 # BR( H -> t tbar )
|
---|
165 | 3.65529942E-02 2 15 -15 # BR( H -> tau- tau+)
|
---|
166 | 5.84716539E-03 2 23 23 # BR( H -> Z Z^(*))
|
---|
167 | 6.31545979E-02 2 24 -24 # BR( H -> W W^(*))
|
---|
168 | 2.55499652E-02 2 21 21 # BR( H -> g g )
|
---|
169 | 1.20906932E-03 2 22 22 # BR( H -> A A )
|
---|
170 | </slha>
|
---|
171 | <MGRunCard>
|
---|
172 | #*********************************************************************
|
---|
173 | # MadGraph/MadEvent *
|
---|
174 | # http://madgraph.hep.uiuc.edu *
|
---|
175 | # *
|
---|
176 | # run_card.dat *
|
---|
177 | # *
|
---|
178 | # This file is used to set the parameters of the run. *
|
---|
179 | # *
|
---|
180 | # Some notation/conventions: *
|
---|
181 | # *
|
---|
182 | # Lines starting with a '# ' are info or comments *
|
---|
183 | # *
|
---|
184 | # mind the format: value = variable ! comment *
|
---|
185 | #*********************************************************************
|
---|
186 | #
|
---|
187 | #*******************
|
---|
188 | # Running parameters
|
---|
189 | #*******************
|
---|
190 | #
|
---|
191 | #*********************************************************************
|
---|
192 | # Tag name for the run (one word) *
|
---|
193 | #*********************************************************************
|
---|
194 | 'AFB-tevatron' = run_tag ! name of the run
|
---|
195 | #*********************************************************************
|
---|
196 | # Number of events and rnd seed *
|
---|
197 | #*********************************************************************
|
---|
198 | 50000 = nevents ! Number of unweighted events requested
|
---|
199 | 201 = iseed ! rnd seed (0=assigned automatically=default))
|
---|
200 | #*********************************************************************
|
---|
201 | # Collider type and energy *
|
---|
202 | #*********************************************************************
|
---|
203 | 1 = lpp1 ! beam 1 type (0=NO PDF)
|
---|
204 | -1 = lpp2 ! beam 2 type (0=NO PDF)
|
---|
205 | 980 = ebeam1 ! beam 1 energy in GeV
|
---|
206 | 980 = ebeam2 ! beam 2 energy in GeV
|
---|
207 | #*********************************************************************
|
---|
208 | # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
|
---|
209 | #*********************************************************************
|
---|
210 | 'cteq6l1' = pdlabel ! PDF set
|
---|
211 | #*********************************************************************
|
---|
212 | # Renormalization and factorization scales *
|
---|
213 | #*********************************************************************
|
---|
214 | F = fixed_ren_scale ! if .true. use fixed ren scale
|
---|
215 | F = fixed_fac_scale ! if .true. use fixed fac scale
|
---|
216 | 174.0 = scale ! fixed ren scale
|
---|
217 | 174.0 = dsqrt_q2fact1 ! fixed fact scale for pdf1
|
---|
218 | 174.0 = dsqrt_q2fact2 ! fixed fact scale for pdf2
|
---|
219 | 1 = scalefact ! scale factor for event-by-event scales
|
---|
220 | #*********************************************************************
|
---|
221 | # Matching - Warning! ickkw > 0 is still beta
|
---|
222 | #*********************************************************************
|
---|
223 | 1 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
|
---|
224 | #*********************************************************************
|
---|
225 | #
|
---|
226 | #*******************
|
---|
227 | # Standard Cuts
|
---|
228 | #*******************
|
---|
229 | #
|
---|
230 | #*********************************************************************
|
---|
231 | # Minimum pt's *
|
---|
232 | #*********************************************************************
|
---|
233 | 10 = ptj ! minimum pt for the jets
|
---|
234 | 10 = ptb ! minimum pt for the b
|
---|
235 | 10 = pta ! minimum pt for the photons
|
---|
236 | 5 = ptl ! minimum pt for the charged leptons
|
---|
237 | #*********************************************************************
|
---|
238 | # Maximum rapidity *
|
---|
239 | #*********************************************************************
|
---|
240 | 5 = etaj ! max rap for the jets
|
---|
241 | 5 = etab ! max rap for the b
|
---|
242 | 2.5 = etaa ! max rap for the photons
|
---|
243 | 5 = etal ! max rap for the charged leptons
|
---|
244 | #*********************************************************************
|
---|
245 | # Minimum DeltaR distance *
|
---|
246 | #*********************************************************************
|
---|
247 | 0.001 = drjj ! distance between jets
|
---|
248 | 0 = drbb ! distance between b's
|
---|
249 | 0.3 = drll ! distance between leptons
|
---|
250 | 0.3 = draa ! distance between gammas
|
---|
251 | 0.3 = drbj ! distance between b and jet
|
---|
252 | 0.3 = draj ! distance between gamma and jet
|
---|
253 | 0.3 = drjl ! distance between jet and lepton
|
---|
254 | 0.3 = drab ! distance between gamma and b
|
---|
255 | 0.3 = drbl ! distance between b and lepton
|
---|
256 | 0.3 = dral ! distance between gamma and lepton
|
---|
257 | #*********************************************************************
|
---|
258 | # Minimum invariant mass for pairs *
|
---|
259 | #*********************************************************************
|
---|
260 | 10 = mmjj ! min invariant mass of a jet pair
|
---|
261 | 0 = mmbb ! min invariant mass of a b pair
|
---|
262 | 0 = mmaa ! min invariant mass of gamma gamma pair
|
---|
263 | 10 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
|
---|
264 | #*********************************************************************
|
---|
265 | # Inclusive cuts *
|
---|
266 | #*********************************************************************
|
---|
267 | 0 = xptj ! minimum pt for at least one jet
|
---|
268 | 0 = xptb ! minimum pt for at least one b
|
---|
269 | 0 = xpta ! minimum pt for at least one photon
|
---|
270 | 0 = xptl ! minimum pt for at least one charged lepton
|
---|
271 | #*********************************************************************
|
---|
272 | # WBF cuts *
|
---|
273 | #*********************************************************************
|
---|
274 | 0 = xetamin ! minimum rapidity for two jets in the WBF case
|
---|
275 | #*********************************************************************
|
---|
276 | # Jet measure cuts *
|
---|
277 | #*********************************************************************
|
---|
278 | 10 = xqcut ! minimum kt jet measure between partons
|
---|
279 | #*********************************************************************
|
---|
280 | </MGRunCard>
|
---|
281 | <MGGenerationInfo>
|
---|
282 | # Number of Events : 50000
|
---|
283 | # Integrated weight (pb) : .23399E+03
|
---|
284 | # Truncated wgt (pb) : .86211E-01
|
---|
285 | # Unit wgt : .46798E-02
|
---|
286 | </MGGenerationInfo>
|
---|
287 | <MGPythiaCard>
|
---|
288 | !...Fix (0) or event-based (1) PS scale
|
---|
289 | PSSCALE=1
|
---|
290 |
|
---|
291 | !...Parton showering on or off
|
---|
292 | MSTP(61)=1
|
---|
293 | MSTP(71)=1
|
---|
294 |
|
---|
295 | !...Fragmentation/hadronization on or off
|
---|
296 | MSTJ(1)=1
|
---|
297 |
|
---|
298 | !...Multiple interactions on or off
|
---|
299 | MSTP(81)=0
|
---|
300 |
|
---|
301 | !...Cutoff in jet measure for matching
|
---|
302 | QCUT = 20
|
---|
303 | !...ETCLUS: Relative ET in jet (KTCLUS)
|
---|
304 | ETCLUS = 20
|
---|
305 | !...ETMIN: min ET per jet (GETJET)
|
---|
306 | ETMIN = 20
|
---|
307 | !...RMAX: Jet cone radius (GETJET)
|
---|
308 | RMAX = 0.4
|
---|
309 | !...ETAMAX: Maximum Eta for calorimeter/tracker
|
---|
310 | ETAMAX = 5
|
---|
311 | !...ETAJETMAX: Maximum Eta for jet
|
---|
312 | ETAJETMAX = 4.5
|
---|
313 | !...NCY: Eta cells in calorimeter
|
---|
314 | NCY = 100
|
---|
315 | !...NCPHI: Phi cells in calorimeter
|
---|
316 | NCPHI = 60
|
---|
317 |
|
---|
318 | LHAPATH=/home/madgraph/MG_ME/WebBin/pythia-pgs/PDFsets
|
---|
319 | </MGPythiaCard>
|
---|
320 | <MGPGSCard>
|
---|
321 | CDF ! parameter set name
|
---|
322 | 80 ! eta cells in calorimeter
|
---|
323 | 24 ! phi cells in calorimeter
|
---|
324 | 0.1 ! eta width of calorimeter cells
|
---|
325 | 0.261799388 ! phi width of calorimeter cells
|
---|
326 | 0.01 ! electromagnetic calorimeter resolution const
|
---|
327 | 0.2 ! electromagnetic calorimeter resolution * sqrt(E)
|
---|
328 | 0.8 ! hadronic calolrimeter resolution * sqrt(E)
|
---|
329 | 0.2 ! MET resolution
|
---|
330 | 0.05 ! calorimeter cell edge crack fraction
|
---|
331 | cone ! jet finding algorithm (cone or ktjet)
|
---|
332 | 3.0 ! calorimeter cluster finding seed threshold (GeV)
|
---|
333 | 0.5 ! calorimeter cluster finding shoulder threshold (GeV)
|
---|
334 | 0.4 ! calorimeter cluster finder cone size (delta R)
|
---|
335 | 1.0 ! outer radius of tracker (m)
|
---|
336 | 1.4 ! magnetic field (T)
|
---|
337 | 0.000040 ! sagitta resolution (m)
|
---|
338 | 0.98 ! track finding efficiency
|
---|
339 | 0.30 ! minimum track pt (GeV/c)
|
---|
340 | 2.0 ! tracking eta coverage
|
---|
341 | 2.0 ! e/gamma eta coverage
|
---|
342 | 2.0 ! muon eta coverage
|
---|
343 | 2.0 ! tau eta coverage
|
---|
344 | </MGPGSCard>
|
---|