| 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.1
|
|---|
| 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>mu+mu- # Process
|
|---|
| 73 | QCD=99 # max qcd order
|
|---|
| 74 | QED=99 # max qed order
|
|---|
| 75 | end_coup # the coupling list is over
|
|---|
| 76 | done # the process list is over
|
|---|
| 77 | # End PROCESS # This is TAG. Do not modify this line
|
|---|
| 78 | #*********************************************************************
|
|---|
| 79 | # Model information *
|
|---|
| 80 | #*********************************************************************
|
|---|
| 81 | # Begin MODEL # This is TAG. Do not modify this line
|
|---|
| 82 | sm
|
|---|
| 83 | # End MODEL # This is TAG. Do not modify this line
|
|---|
| 84 | #*********************************************************************
|
|---|
| 85 | # Start multiparticle definitions *
|
|---|
| 86 | #*********************************************************************
|
|---|
| 87 | # Begin MULTIPARTICLES # This is TAG. Do not modify this line
|
|---|
| 88 | p uu~dd~ss~cc~g
|
|---|
| 89 | j uu~dd~ss~cc~g
|
|---|
| 90 | l+ e+mu+
|
|---|
| 91 | l- e-mu-
|
|---|
| 92 | vl vevm
|
|---|
| 93 | vl~ ve~vm~
|
|---|
| 94 | # End MULTIPARTICLES # This is TAG. Do not modify this line
|
|---|
| 95 | </MGProcCard>
|
|---|
| 96 | <slha>
|
|---|
| 97 | #******************************************************************
|
|---|
| 98 | # MadGraph/MadEvent *
|
|---|
| 99 | #******************************************************************
|
|---|
| 100 | # Les Houches friendly file for the SM parameters of MadGraph *
|
|---|
| 101 | # Spectrum and decay widths produced by SMCalc *
|
|---|
| 102 | #******************************************************************
|
|---|
| 103 | #*Please note the following IMPORTANT issues: *
|
|---|
| 104 | # *
|
|---|
| 105 | #0. REFRAIN from editing this file by hand! Some of the parame- *
|
|---|
| 106 | # ters are not independent *
|
|---|
| 107 | # (such as G_Fermi, alpha_em, sin(theta_W),MZ,MW) and serious *
|
|---|
| 108 | # problems might be encountered (such as violation of unitarity *
|
|---|
| 109 | # or gauge invariance). Always use a calculator. *
|
|---|
| 110 | # *
|
|---|
| 111 | #1. alpha_S(MZ) has been used in the calculation of the parameters*
|
|---|
| 112 | # but, for consistency, it will be reset by madgraph to the *
|
|---|
| 113 | # value expected IF the pdfs for collisions with hadrons are *
|
|---|
| 114 | # used. This value is KEPT by madgraph when no pdf are used *
|
|---|
| 115 | # lpp(i)=0 . *
|
|---|
| 116 | # *
|
|---|
| 117 | #2. Values of the charm and bottom kinematic (pole) masses are *
|
|---|
| 118 | # those used in the matrix elements and phase space UNLESS they *
|
|---|
| 119 | # are set to ZERO from the start in the model (particles.dat) *
|
|---|
| 120 | # This happens, for example, when using 5-flavor QCD where *
|
|---|
| 121 | # charm and bottom are treated as partons in the initial state *
|
|---|
| 122 | # and a zero mass might be hardwired in the model definition. *
|
|---|
| 123 | # *
|
|---|
| 124 | #******************************************************************
|
|---|
| 125 | Block SMINPUTS # Standard Model inputs
|
|---|
| 126 | 1 1.32506980E+02 # alpha_em(MZ)(-1) SM MSbar
|
|---|
| 127 | 2 1.16639000E-05 # G_Fermi
|
|---|
| 128 | 3 1.18000000E-01 # alpha_s(MZ) SM MSbar
|
|---|
| 129 | 4 9.11880000E+01 # Z mass (as input parameter)
|
|---|
| 130 | Block MGYUKAWA # Yukawa masses m/v=y/sqrt(2)
|
|---|
| 131 | # PDG YMASS
|
|---|
| 132 | 5 4.20000000E+00 # mbottom for the Yukawa y_b
|
|---|
| 133 | 4 1.42000000E+00 # mcharm for the Yukawa y_c
|
|---|
| 134 | 6 1.64500000E+02 # mtop for the Yukawa y_t
|
|---|
| 135 | 15 1.77700000E+00 # mtau for the Yukawa y_ta
|
|---|
| 136 | Block MGCKM # CKM elements for MadGraph
|
|---|
| 137 | 1 1 9.75000000E-01 # Vud for Cabibbo matrix
|
|---|
| 138 | Block MASS # Mass spectrum (kinematic masses)
|
|---|
| 139 | # PDG Mass
|
|---|
| 140 | 5 4.70000000E+00 # bottom pole mass
|
|---|
| 141 | 6 1.74300000E+02 # top pole mass
|
|---|
| 142 | 15 1.77700000E+00 # tau mass
|
|---|
| 143 | 23 9.11880000E+01 # Z mass
|
|---|
| 144 | 24 8.04190000E+01 # W mass
|
|---|
| 145 | 25 1.20000000E+02 # H mass
|
|---|
| 146 | # PDG Width
|
|---|
| 147 | DECAY 6 1.50833649E+00 # top width
|
|---|
| 148 | DECAY 23 2.44140351E+00 # Z width
|
|---|
| 149 | DECAY 24 2.04759951E+00 # W width
|
|---|
| 150 | DECAY 25 5.75308848E-03 # H width
|
|---|
| 151 | # BR NDA ID1 ID2
|
|---|
| 152 | 8.27451012E-02 2 4 -4 # BR( H -> c cbar )
|
|---|
| 153 | 7.17809696E-01 2 5 -5 # BR( H -> b bbar )
|
|---|
| 154 | 0.00000000E+00 2 6 -6 # BR( H -> t tbar )
|
|---|
| 155 | 4.31720144E-02 2 15 -15 # BR( H -> tau- tau+)
|
|---|
| 156 | 6.90597075E-03 2 23 23 # BR( H -> Z Z^(*))
|
|---|
| 157 | 7.45906395E-02 2 24 -24 # BR( H -> W W^(*))
|
|---|
| 158 | 3.01765558E-02 2 21 21 # BR( H -> g g )
|
|---|
| 159 | 1.42800773E-03 2 22 22 # BR( H -> A A )
|
|---|
| 160 | </slha>
|
|---|
| 161 | <MGRunCard>
|
|---|
| 162 | #*********************************************************************
|
|---|
| 163 | # MadGraph/MadEvent *
|
|---|
| 164 | # http://madgraph.hep.uiuc.edu *
|
|---|
| 165 | # *
|
|---|
| 166 | # run_card.dat *
|
|---|
| 167 | # *
|
|---|
| 168 | # This file is used to set the parameters of the run. *
|
|---|
| 169 | # *
|
|---|
| 170 | # Some notation/conventions: *
|
|---|
| 171 | # *
|
|---|
| 172 | # Lines starting with a '# ' are info or comments *
|
|---|
| 173 | # *
|
|---|
| 174 | # mind the format: value = variable ! comment *
|
|---|
| 175 | #*********************************************************************
|
|---|
| 176 | #
|
|---|
| 177 | #*******************
|
|---|
| 178 | # Running parameters
|
|---|
| 179 | #*******************
|
|---|
| 180 | #
|
|---|
| 181 | #*********************************************************************
|
|---|
| 182 | # Tag name for the run (one word) *
|
|---|
| 183 | #*********************************************************************
|
|---|
| 184 | 'fermi' = run_tag ! name of the run
|
|---|
| 185 | #*********************************************************************
|
|---|
| 186 | # Number of events and rnd seed *
|
|---|
| 187 | #*********************************************************************
|
|---|
| 188 | 50000 = nevents ! Number of unweighted events requested
|
|---|
| 189 | 25 = iseed ! rnd seed (0=assigned automatically=default))
|
|---|
| 190 | #*********************************************************************
|
|---|
| 191 | # Collider type and energy *
|
|---|
| 192 | #*********************************************************************
|
|---|
| 193 | 1 = lpp1 ! beam 1 type (0=NO PDF)
|
|---|
| 194 | -1 = lpp2 ! beam 2 type (0=NO PDF)
|
|---|
| 195 | 980 = ebeam1 ! beam 1 energy in GeV
|
|---|
| 196 | 980 = ebeam2 ! beam 2 energy in GeV
|
|---|
| 197 | #*********************************************************************
|
|---|
| 198 | # Beam polarization from -100 (left-handed) to 100 (right-handed) *
|
|---|
| 199 | #*********************************************************************
|
|---|
| 200 | 0 = polbeam1 ! beam polarization for beam 1
|
|---|
| 201 | 0 = polbeam2 ! beam polarization for beam 2
|
|---|
| 202 | #*********************************************************************
|
|---|
| 203 | # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
|
|---|
| 204 | #*********************************************************************
|
|---|
| 205 | 'cteq6l1' = pdlabel ! PDF set
|
|---|
| 206 | #*********************************************************************
|
|---|
| 207 | # Renormalization and factorization scales *
|
|---|
| 208 | #*********************************************************************
|
|---|
| 209 | F = fixed_ren_scale ! if .true. use fixed ren scale
|
|---|
| 210 | F = fixed_fac_scale ! if .true. use fixed fac scale
|
|---|
| 211 | 174.0 = scale ! fixed ren scale
|
|---|
| 212 | 174.0 = dsqrt_q2fact1 ! fixed fact scale for pdf1
|
|---|
| 213 | 174.0 = dsqrt_q2fact2 ! fixed fact scale for pdf2
|
|---|
| 214 | 1 = scalefact ! scale factor for event-by-event scales
|
|---|
| 215 | #*********************************************************************
|
|---|
| 216 | # Matching - Warning! ickkw > 0 is still beta
|
|---|
| 217 | #*********************************************************************
|
|---|
| 218 | 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
|
|---|
| 219 | #*********************************************************************
|
|---|
| 220 | #
|
|---|
| 221 | #*******************
|
|---|
| 222 | # Standard Cuts
|
|---|
| 223 | #*******************
|
|---|
| 224 | #
|
|---|
| 225 | #*********************************************************************
|
|---|
| 226 | # Minimum pt's *
|
|---|
| 227 | #*********************************************************************
|
|---|
| 228 | 10 = ptj ! minimum pt for the jets
|
|---|
| 229 | 10 = ptb ! minimum pt for the b
|
|---|
| 230 | 10 = pta ! minimum pt for the photons
|
|---|
| 231 | 5 = ptl ! minimum pt for the charged leptons
|
|---|
| 232 | #*********************************************************************
|
|---|
| 233 | # Minimum energy (in the lab frame) *
|
|---|
| 234 | #*********************************************************************
|
|---|
| 235 | 0 = ej ! minimum E for the jets
|
|---|
| 236 | 0 = eb ! minimum E for the b
|
|---|
| 237 | 0 = ea ! minimum E for the photons
|
|---|
| 238 | 0 = el ! minimum E for the charged leptons
|
|---|
| 239 | #*********************************************************************
|
|---|
| 240 | # Maximum rapidity *
|
|---|
| 241 | #*********************************************************************
|
|---|
| 242 | 5d0 = etaj ! max rap for the jets
|
|---|
| 243 | 5d0 = etab ! max rap for the b
|
|---|
| 244 | 2.5 = etaa ! max rap for the photons
|
|---|
| 245 | 5d0 = etal ! max rap for the charged leptons
|
|---|
| 246 | #*********************************************************************
|
|---|
| 247 | # Minimum DeltaR distance *
|
|---|
| 248 | #*********************************************************************
|
|---|
| 249 | 0.4 = drjj ! distance between jets
|
|---|
| 250 | 0 = drbb ! distance between b's
|
|---|
| 251 | 0.3 = drll ! distance between leptons
|
|---|
| 252 | 0.3 = draa ! distance between gammas
|
|---|
| 253 | 0.3 = drbj ! distance between b and jet
|
|---|
| 254 | 0.3 = draj ! distance between gamma and jet
|
|---|
| 255 | 0.3 = drjl ! distance between jet and lepton
|
|---|
| 256 | 0.3 = drab ! distance between gamma and b
|
|---|
| 257 | 0.3 = drbl ! distance between b and lepton
|
|---|
| 258 | 0.3 = dral ! distance between gamma and lepton
|
|---|
| 259 | #*********************************************************************
|
|---|
| 260 | # Minimum invariant mass for pairs *
|
|---|
| 261 | #*********************************************************************
|
|---|
| 262 | 10 = mmjj ! min invariant mass of a jet pair
|
|---|
| 263 | 0 = mmbb ! min invariant mass of a b pair
|
|---|
| 264 | 0 = mmaa ! min invariant mass of gamma gamma pair
|
|---|
| 265 | 10 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
|
|---|
| 266 | #*********************************************************************
|
|---|
| 267 | # Inclusive cuts *
|
|---|
| 268 | #*********************************************************************
|
|---|
| 269 | 0 = xptj ! minimum pt for at least one jet
|
|---|
| 270 | 0 = xptb ! minimum pt for at least one b
|
|---|
| 271 | 0 = xpta ! minimum pt for at least one photon
|
|---|
| 272 | 0 = xptl ! minimum pt for at least one charged lepton
|
|---|
| 273 | #*********************************************************************
|
|---|
| 274 | # WBF cuts *
|
|---|
| 275 | #*********************************************************************
|
|---|
| 276 | 0 = xetamin ! minimum rapidity for two jets in the WBF case
|
|---|
| 277 | 0 = deltaeta ! minimum rapidity for two jets in the WBF case
|
|---|
| 278 | #*********************************************************************
|
|---|
| 279 | # Jet measure cuts *
|
|---|
| 280 | #*********************************************************************
|
|---|
| 281 | 0 = xqcut ! minimum kt jet measure between partons
|
|---|
| 282 | #*********************************************************************
|
|---|
| 283 | </MGRunCard>
|
|---|
| 284 | <MGGenerationInfo>
|
|---|
| 285 | # Number of Events : 50000
|
|---|
| 286 | # Integrated weight (pb) : .67807E+03
|
|---|
| 287 | # Truncated wgt (pb) : .76413E+00
|
|---|
| 288 | # Unit wgt : .13561E-01
|
|---|
| 289 | </MGGenerationInfo>
|
|---|
| 290 | <MGPythiaCard>
|
|---|
| 291 | !...Fix (0) or event-based (1) PS scale
|
|---|
| 292 | PSSCALE=1
|
|---|
| 293 | !...Parton showering on or off
|
|---|
| 294 | MSTP(61)=1
|
|---|
| 295 | MSTP(71)=1
|
|---|
| 296 |
|
|---|
| 297 | !...Fragmentation/hadronization on or off
|
|---|
| 298 | MSTJ(1)=1
|
|---|
| 299 |
|
|---|
| 300 | !...Multiple interactions on or off
|
|---|
| 301 | MSTP(81)=0
|
|---|
| 302 |
|
|---|
| 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 | LHAPATH=/home/madgraph/MG_ME/WebBin/pythia-pgs/PDFsets
|
|---|
| 318 | </MGPythiaCard>
|
|---|
| 319 | <MGPGSCard>
|
|---|
| 320 | CDF ! parameter set name
|
|---|
| 321 | 80 ! eta cells in calorimeter
|
|---|
| 322 | 24 ! phi cells in calorimeter
|
|---|
| 323 | 0.1 ! eta width of calorimeter cells
|
|---|
| 324 | 0.261799388 ! phi width of calorimeter cells
|
|---|
| 325 | 0.01 ! electromagnetic calorimeter resolution const
|
|---|
| 326 | 0.2 ! electromagnetic calorimeter resolution * sqrt(E)
|
|---|
| 327 | 0.8 ! hadronic calolrimeter resolution * sqrt(E)
|
|---|
| 328 | 0.2 ! MET resolution
|
|---|
| 329 | 0.05 ! calorimeter cell edge crack fraction
|
|---|
| 330 | cone ! jet finding algorithm (cone or ktjet)
|
|---|
| 331 | 3.0 ! calorimeter cluster finding seed threshold (GeV)
|
|---|
| 332 | 0.5 ! calorimeter cluster finding shoulder threshold (GeV)
|
|---|
| 333 | 0.4 ! calorimeter cluster finder cone size (delta R)
|
|---|
| 334 | 1.0 ! outer radius of tracker (m)
|
|---|
| 335 | 1.4 ! magnetic field (T)
|
|---|
| 336 | 0.000040 ! sagitta resolution (m)
|
|---|
| 337 | 0.98 ! track finding efficiency
|
|---|
| 338 | 0.30 ! minimum track pt (GeV/c)
|
|---|
| 339 | 2.0 ! tracking eta coverage
|
|---|
| 340 | 2.0 ! e/gamma eta coverage
|
|---|
| 341 | 2.0 ! muon eta coverage
|
|---|
| 342 | 2.0 ! tau eta coverage
|
|---|
| 343 | </MGPGSCard>
|
|---|
