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