| 1 | #*********************************************************************
|
|---|
| 2 | # MadGraph/MadEvent *
|
|---|
| 3 | # http://madgraph.hep.uiuc.edu *
|
|---|
| 4 | # *
|
|---|
| 5 | # run_card.dat *
|
|---|
| 6 | # *
|
|---|
| 7 | # This file is used to set the parameters of the run. *
|
|---|
| 8 | # *
|
|---|
| 9 | # Some notation/conventions: *
|
|---|
| 10 | # *
|
|---|
| 11 | # Lines starting with a '# ' are info or comments *
|
|---|
| 12 | # *
|
|---|
| 13 | # mind the format: value = variable ! comment *
|
|---|
| 14 | #*********************************************************************
|
|---|
| 15 | #
|
|---|
| 16 | #*******************
|
|---|
| 17 | # Running parameters
|
|---|
| 18 | #*******************
|
|---|
| 19 | #
|
|---|
| 20 | #*********************************************************************
|
|---|
| 21 | # Tag name for the run (one word) *
|
|---|
| 22 | #*********************************************************************
|
|---|
| 23 | tag_1 = run_tag ! name of the run
|
|---|
| 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 | # Warning: Do not generate more than 1M events in a single run *
|
|---|
| 31 | # If you want to run Pythia, avoid more than 50k events in a run. *
|
|---|
| 32 | #*********************************************************************
|
|---|
| 33 | 1000000 = nevents ! Number of unweighted events requested
|
|---|
| 34 | 0 = iseed ! rnd seed (0=assigned automatically=default))
|
|---|
| 35 | #*********************************************************************
|
|---|
| 36 | # Collider type and energy *
|
|---|
| 37 | # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
|
|---|
| 38 | # 3=photon from electron *
|
|---|
| 39 | #*********************************************************************
|
|---|
| 40 | 1 = lpp1 ! beam 1 type
|
|---|
| 41 | 1 = lpp2 ! beam 2 type
|
|---|
| 42 | 4000 = ebeam1 ! beam 1 total energy in GeV
|
|---|
| 43 | 4000 = ebeam2 ! beam 2 total energy in GeV
|
|---|
| 44 | #*********************************************************************
|
|---|
| 45 | # Beam polarization from -100 (left-handed) to 100 (right-handed) *
|
|---|
| 46 | #*********************************************************************
|
|---|
| 47 | 0 = polbeam1 ! beam polarization for beam 1
|
|---|
| 48 | 0 = polbeam2 ! beam polarization for beam 2
|
|---|
| 49 | #*********************************************************************
|
|---|
| 50 | # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
|
|---|
| 51 | #*********************************************************************
|
|---|
| 52 | 'cteq6l1' = pdlabel ! PDF set
|
|---|
| 53 | #*********************************************************************
|
|---|
| 54 | # Renormalization and factorization scales *
|
|---|
| 55 | #*********************************************************************
|
|---|
| 56 | T = fixed_ren_scale ! if .true. use fixed ren scale
|
|---|
| 57 | T = fixed_fac_scale ! if .true. use fixed fac scale
|
|---|
| 58 | 500.00 = scale ! fixed ren scale
|
|---|
| 59 | 500.00 = dsqrt_q2fact1 ! fixed fact scale for pdf1
|
|---|
| 60 | 500.00 = dsqrt_q2fact2 ! fixed fact scale for pdf2
|
|---|
| 61 | 1 = scalefact ! scale factor for event-by-event scales
|
|---|
| 62 | #*********************************************************************
|
|---|
| 63 | # Matching - Warning! ickkw > 1 is still beta
|
|---|
| 64 | #*********************************************************************
|
|---|
| 65 | 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
|
|---|
| 66 | 1 = highestmult ! for ickkw=2, highest mult group
|
|---|
| 67 | 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
|
|---|
| 68 | 1 = alpsfact ! scale factor for QCD emission vx
|
|---|
| 69 | F = chcluster ! cluster only according to channel diag
|
|---|
| 70 | T = pdfwgt ! for ickkw=1, perform pdf reweighting
|
|---|
| 71 | 5 = asrwgtflavor ! highest quark flavor for a_s reweight
|
|---|
| 72 | #*********************************************************************
|
|---|
| 73 | # Automatic ptj and mjj cuts if xqcut > 0
|
|---|
| 74 | # (turn off for VBF and single top processes)
|
|---|
| 75 | #**********************************************************
|
|---|
| 76 | F = auto_ptj_mjj ! Automatic setting of ptj and mjj
|
|---|
| 77 | #**********************************************************
|
|---|
| 78 | #
|
|---|
| 79 | #**********************************
|
|---|
| 80 | # BW cutoff (M+/-bwcutoff*Gamma)
|
|---|
| 81 | #**********************************
|
|---|
| 82 | 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
|
|---|
| 83 | #**********************************************************
|
|---|
| 84 | # Apply pt/E/eta/dr/mij cuts on decay products or not
|
|---|
| 85 | # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
|
|---|
| 86 | #**********************************************************
|
|---|
| 87 | F = cut_decays ! Cut decay products
|
|---|
| 88 | #*************************************************************
|
|---|
| 89 | # Number of helicities to sum per event (0 = all helicities)
|
|---|
| 90 | # 0 gives more stable result, but longer run time (needed for
|
|---|
| 91 | # long decay chains e.g.).
|
|---|
| 92 | # Use >=2 if most helicities contribute, e.g. pure QCD.
|
|---|
| 93 | #*************************************************************
|
|---|
| 94 | 0 = nhel ! Number of helicities used per event
|
|---|
| 95 | #*******************
|
|---|
| 96 | # Standard Cuts
|
|---|
| 97 | #*******************
|
|---|
| 98 | #
|
|---|
| 99 | #*********************************************************************
|
|---|
| 100 | # Minimum and maximum pt's (for max, -1 means no cut) *
|
|---|
| 101 | #*********************************************************************
|
|---|
| 102 | -1 = ptj ! minimum pt for the jets
|
|---|
| 103 | 0 = ptb ! minimum pt for the b
|
|---|
| 104 | -1 = pta ! minimum pt for the photons
|
|---|
| 105 | -1 = ptl ! minimum pt for the charged leptons
|
|---|
| 106 | 0 = misset ! minimum missing Et (sum of neutrino's momenta)
|
|---|
| 107 | 0 = ptheavy ! minimum pt for one heavy final state
|
|---|
| 108 | -1 = ptonium ! minimum pt for the quarkonium states
|
|---|
| 109 | -1 = ptjmax ! maximum pt for the jets
|
|---|
| 110 | -1 = ptbmax ! maximum pt for the b
|
|---|
| 111 | -1 = ptamax ! maximum pt for the photons
|
|---|
| 112 | -1 = ptlmax ! maximum pt for the charged leptons
|
|---|
| 113 | -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
|
|---|
| 114 | #*********************************************************************
|
|---|
| 115 | # Minimum and maximum E's (in the lab frame) *
|
|---|
| 116 | #*********************************************************************
|
|---|
| 117 | 0 = ej ! minimum E for the jets
|
|---|
| 118 | 0 = eb ! minimum E for the b
|
|---|
| 119 | 0 = ea ! minimum E for the photons
|
|---|
| 120 | 0 = el ! minimum E for the charged leptons
|
|---|
| 121 | -1 = ejmax ! maximum E for the jets
|
|---|
| 122 | -1 = ebmax ! maximum E for the b
|
|---|
| 123 | -1 = eamax ! maximum E for the photons
|
|---|
| 124 | -1 = elmax ! maximum E for the charged leptons
|
|---|
| 125 | #*********************************************************************
|
|---|
| 126 | # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
|
|---|
| 127 | #*********************************************************************
|
|---|
| 128 | -1 = etaj ! max rap for the jets
|
|---|
| 129 | -1 = etab ! max rap for the b
|
|---|
| 130 | -1 = etaa ! max rap for the photons
|
|---|
| 131 | -1 = etal ! max rap for the charged leptons
|
|---|
| 132 | -1 = etaonium ! max rap for the quarkonium states
|
|---|
| 133 | 0 = etajmin ! min rap for the jets
|
|---|
| 134 | 0 = etabmin ! min rap for the b
|
|---|
| 135 | 0 = etaamin ! min rap for the photons
|
|---|
| 136 | 0 = etalmin ! main rap for the charged leptons
|
|---|
| 137 | #*********************************************************************
|
|---|
| 138 | # Minimum and maximum DeltaR distance *
|
|---|
| 139 | #*********************************************************************
|
|---|
| 140 | -1 = drjj ! min distance between jets
|
|---|
| 141 | 0 = drbb ! min distance between b's
|
|---|
| 142 | -1 = drll ! min distance between leptons
|
|---|
| 143 | -1 = draa ! min distance between gammas
|
|---|
| 144 | 0 = drbj ! min distance between b and jet
|
|---|
| 145 | -1 = draj ! min distance between gamma and jet
|
|---|
| 146 | -1 = drjl ! min distance between jet and lepton
|
|---|
| 147 | 0 = drab ! min distance between gamma and b
|
|---|
| 148 | 0 = drbl ! min distance between b and lepton
|
|---|
| 149 | -1 = dral ! min distance between gamma and lepton
|
|---|
| 150 | -1 = drjjmax ! max distance between jets
|
|---|
| 151 | -1 = drbbmax ! max distance between b's
|
|---|
| 152 | -1 = drllmax ! max distance between leptons
|
|---|
| 153 | -1 = draamax ! max distance between gammas
|
|---|
| 154 | -1 = drbjmax ! max distance between b and jet
|
|---|
| 155 | -1 = drajmax ! max distance between gamma and jet
|
|---|
| 156 | -1 = drjlmax ! max distance between jet and lepton
|
|---|
| 157 | -1 = drabmax ! max distance between gamma and b
|
|---|
| 158 | -1 = drblmax ! max distance between b and lepton
|
|---|
| 159 | -1 = dralmax ! maxdistance between gamma and lepton
|
|---|
| 160 | #*********************************************************************
|
|---|
| 161 | # Minimum and maximum invariant mass for pairs *
|
|---|
| 162 | #*********************************************************************
|
|---|
| 163 | 0 = mmjj ! min invariant mass of a jet pair
|
|---|
| 164 | 0 = mmbb ! min invariant mass of a b pair
|
|---|
| 165 | 0 = mmaa ! min invariant mass of gamma gamma pair
|
|---|
| 166 | 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
|
|---|
| 167 | -1 = mmjjmax ! max invariant mass of a jet pair
|
|---|
| 168 | -1 = mmbbmax ! max invariant mass of a b pair
|
|---|
| 169 | -1 = mmaamax ! max invariant mass of gamma gamma pair
|
|---|
| 170 | -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
|
|---|
| 171 | #*********************************************************************
|
|---|
| 172 | # Minimum and maximum invariant mass for all letpons *
|
|---|
| 173 | #*********************************************************************
|
|---|
| 174 | 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
|
|---|
| 175 | -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
|
|---|
| 176 | #*********************************************************************
|
|---|
| 177 | # Minimum and maximum pt for 4-momenta sum of leptons *
|
|---|
| 178 | #*********************************************************************
|
|---|
| 179 | 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
|
|---|
| 180 | -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
|
|---|
| 181 | #*********************************************************************
|
|---|
| 182 | # Inclusive cuts *
|
|---|
| 183 | #*********************************************************************
|
|---|
| 184 | 0 = xptj ! minimum pt for at least one jet
|
|---|
| 185 | 0 = xptb ! minimum pt for at least one b
|
|---|
| 186 | 0 = xpta ! minimum pt for at least one photon
|
|---|
| 187 | 0 = xptl ! minimum pt for at least one charged lepton
|
|---|
| 188 | #*********************************************************************
|
|---|
| 189 | # Control the pt's of the jets sorted by pt *
|
|---|
| 190 | #*********************************************************************
|
|---|
| 191 | 0 = ptj1min ! minimum pt for the leading jet in pt
|
|---|
| 192 | 0 = ptj2min ! minimum pt for the second jet in pt
|
|---|
| 193 | 0 = ptj3min ! minimum pt for the third jet in pt
|
|---|
| 194 | 0 = ptj4min ! minimum pt for the fourth jet in pt
|
|---|
| 195 | -1 = ptj1max ! maximum pt for the leading jet in pt
|
|---|
| 196 | -1 = ptj2max ! maximum pt for the second jet in pt
|
|---|
| 197 | -1 = ptj3max ! maximum pt for the third jet in pt
|
|---|
| 198 | -1 = ptj4max ! maximum pt for the fourth jet in pt
|
|---|
| 199 | 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
|
|---|
| 200 | #*********************************************************************
|
|---|
| 201 | # Control the pt's of leptons sorted by pt *
|
|---|
| 202 | #*********************************************************************
|
|---|
| 203 | 0 = ptl1min ! minimum pt for the leading lepton in pt
|
|---|
| 204 | 0 = ptl2min ! minimum pt for the second lepton in pt
|
|---|
| 205 | 0 = ptl3min ! minimum pt for the third lepton in pt
|
|---|
| 206 | 0 = ptl4min ! minimum pt for the fourth lepton in pt
|
|---|
| 207 | -1 = ptl1max ! maximum pt for the leading lepton in pt
|
|---|
| 208 | -1 = ptl2max ! maximum pt for the second lepton in pt
|
|---|
| 209 | -1 = ptl3max ! maximum pt for the third lepton in pt
|
|---|
| 210 | -1 = ptl4max ! maximum pt for the fourth lepton in pt
|
|---|
| 211 | #*********************************************************************
|
|---|
| 212 | # Control the Ht(k)=Sum of k leading jets *
|
|---|
| 213 | #*********************************************************************
|
|---|
| 214 | 0 = htjmin ! minimum jet HT=Sum(jet pt)
|
|---|
| 215 | -1 = htjmax ! maximum jet HT=Sum(jet pt)
|
|---|
| 216 | 0 = ihtmin !inclusive Ht for all partons (including b)
|
|---|
| 217 | -1 = ihtmax !inclusive Ht for all partons (including b)
|
|---|
| 218 | 0 = ht2min ! minimum Ht for the two leading jets
|
|---|
| 219 | 0 = ht3min ! minimum Ht for the three leading jets
|
|---|
| 220 | 0 = ht4min ! minimum Ht for the four leading jets
|
|---|
| 221 | -1 = ht2max ! maximum Ht for the two leading jets
|
|---|
| 222 | -1 = ht3max ! maximum Ht for the three leading jets
|
|---|
| 223 | -1 = ht4max ! maximum Ht for the four leading jets
|
|---|
| 224 | #*********************************************************************
|
|---|
| 225 | # WBF cuts *
|
|---|
| 226 | #*********************************************************************
|
|---|
| 227 | 0 = xetamin ! minimum rapidity for two jets in the WBF case
|
|---|
| 228 | 0 = deltaeta ! minimum rapidity for two jets in the WBF case
|
|---|
| 229 | #*********************************************************************
|
|---|
| 230 | # maximal pdg code for quark to be considered as a light jet *
|
|---|
| 231 | # (otherwise b cuts are applied) *
|
|---|
| 232 | #*********************************************************************
|
|---|
| 233 | 4 = maxjetflavor ! Maximum jet pdg code
|
|---|
| 234 | #*********************************************************************
|
|---|
| 235 | # Jet measure cuts *
|
|---|
| 236 | #*********************************************************************
|
|---|
| 237 | 0 = xqcut ! minimum kt jet measure between partons
|
|---|
| 238 | #*********************************************************************
|
|---|
