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