| 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 | 'fermi' = 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 | #*********************************************************************
|
|---|
| 31 | 10000 = nevents ! Number of unweighted events requested
|
|---|
| 32 | 0 = iseed ! rnd seed (0=assigned automatically=default))
|
|---|
| 33 | #*********************************************************************
|
|---|
| 34 | # Collider type and energy *
|
|---|
| 35 | #*********************************************************************
|
|---|
| 36 | 1 = lpp1 ! beam 1 type (0=NO PDF)
|
|---|
| 37 | 1 = lpp2 ! beam 2 type (0=NO PDF)
|
|---|
| 38 | 7000 = ebeam1 ! beam 1 energy in GeV
|
|---|
| 39 | 7000 = ebeam2 ! beam 2 energy in GeV
|
|---|
| 40 | #*********************************************************************
|
|---|
| 41 | # Beam polarization from -100 (left-handed) to 100 (right-handed) *
|
|---|
| 42 | #*********************************************************************
|
|---|
| 43 | 0 = polbeam1 ! beam polarization for beam 1
|
|---|
| 44 | 0 = polbeam2 ! beam polarization for beam 2
|
|---|
| 45 | #*********************************************************************
|
|---|
| 46 | # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
|
|---|
| 47 | #*********************************************************************
|
|---|
| 48 | 'cteq6l1' = pdlabel ! PDF set
|
|---|
| 49 | #*********************************************************************
|
|---|
| 50 | # Renormalization and factorization scales *
|
|---|
| 51 | #*********************************************************************
|
|---|
| 52 | T = fixed_ren_scale ! if .true. use fixed ren scale
|
|---|
| 53 | T = fixed_fac_scale ! if .true. use fixed fac scale
|
|---|
| 54 | 91.1880 = scale ! fixed ren scale
|
|---|
| 55 | 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
|
|---|
| 56 | 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
|
|---|
| 57 | 1 = scalefact ! scale factor for event-by-event scales
|
|---|
| 58 | #*********************************************************************
|
|---|
| 59 | # Matching - Warning! ickkw > 0 is still beta
|
|---|
| 60 | #*********************************************************************
|
|---|
| 61 | 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
|
|---|
| 62 | #*********************************************************************
|
|---|
| 63 | #
|
|---|
| 64 | #**********************************
|
|---|
| 65 | # BW cutoff (M+/-bwcutoff*Gamma)
|
|---|
| 66 | #**********************************
|
|---|
| 67 | 5 = bwcutoff
|
|---|
| 68 | #*******************
|
|---|
| 69 | # Standard Cuts
|
|---|
| 70 | #*******************
|
|---|
| 71 | #
|
|---|
| 72 | #*********************************************************************
|
|---|
| 73 | # Minimum and maximum pt's *
|
|---|
| 74 | #*********************************************************************
|
|---|
| 75 | 20 = ptj ! minimum pt for the jets
|
|---|
| 76 | 0 = ptb ! minimum pt for the b
|
|---|
| 77 | 10 = pta ! minimum pt for the photons
|
|---|
| 78 | 10 = ptl ! minimum pt for the charged leptons
|
|---|
| 79 | 0 = misset ! minimum missing Et (sum of neutrino's momenta)
|
|---|
| 80 | 0 = ptheavy ! minimum pt for one heavy final state
|
|---|
| 81 | 1.0 = ptonium ! minimum pt for the quarkonium states
|
|---|
| 82 | 1d5 = ptjmax ! maximum pt for the jets
|
|---|
| 83 | 1d5 = ptbmax ! maximum pt for the b
|
|---|
| 84 | 1d5 = ptamax ! maximum pt for the photons
|
|---|
| 85 | 1d5 = ptlmax ! maximum pt for the charged leptons
|
|---|
| 86 | 1d5 = missetmax ! maximum missing Et (sum of neutrino's momenta)
|
|---|
| 87 | #*********************************************************************
|
|---|
| 88 | # Minimum and maximum E's (in the lab frame) *
|
|---|
| 89 | #*********************************************************************
|
|---|
| 90 | 0 = ej ! minimum E for the jets
|
|---|
| 91 | 0 = eb ! minimum E for the b
|
|---|
| 92 | 0 = ea ! minimum E for the photons
|
|---|
| 93 | 0 = el ! minimum E for the charged leptons
|
|---|
| 94 | 1d5 = ejmax ! maximum E for the jets
|
|---|
| 95 | 1d5 = ebmax ! maximum E for the b
|
|---|
| 96 | 1d5 = eamax ! maximum E for the photons
|
|---|
| 97 | 1d5 = elmax ! maximum E for the charged leptons
|
|---|
| 98 | #*********************************************************************
|
|---|
| 99 | # Maximum and minimum rapidity *
|
|---|
| 100 | #*********************************************************************
|
|---|
| 101 | 2.5 = etaj ! max rap for the jets
|
|---|
| 102 | 1d2 = etab ! max rap for the b
|
|---|
| 103 | 2.5 = etaa ! max rap for the photons
|
|---|
| 104 | 2.5 = etal ! max rap for the charged leptons
|
|---|
| 105 | 0.6 = etaonium ! max rap for the quarkonium states
|
|---|
| 106 | 0d0 = etajmin ! min rap for the jets
|
|---|
| 107 | 0d0 = etabmin ! min rap for the b
|
|---|
| 108 | 0d0 = etaamin ! min rap for the photons
|
|---|
| 109 | 0d0 = etalmin ! main rap for the charged leptons
|
|---|
| 110 | #*********************************************************************
|
|---|
| 111 | # Minimum and maximum DeltaR distance *
|
|---|
| 112 | #*********************************************************************
|
|---|
| 113 | 0.4 = drjj ! min distance between jets
|
|---|
| 114 | 0 = drbb ! min distance between b's
|
|---|
| 115 | 0.4 = drll ! min distance between leptons
|
|---|
| 116 | 0.4 = draa ! min distance between gammas
|
|---|
| 117 | 0 = drbj ! min distance between b and jet
|
|---|
| 118 | 0.4 = draj ! min distance between gamma and jet
|
|---|
| 119 | 0.4 = drjl ! min distance between jet and lepton
|
|---|
| 120 | 0 = drab ! min distance between gamma and b
|
|---|
| 121 | 0 = drbl ! min distance between b and lepton
|
|---|
| 122 | 0.4 = dral ! min distance between gamma and lepton
|
|---|
| 123 | 1d2 = drjjmax ! max distance between jets
|
|---|
| 124 | 1d2 = drbbmax ! max distance between b's
|
|---|
| 125 | 1d2 = drllmax ! max distance between leptons
|
|---|
| 126 | 1d2 = draamax ! max distance between gammas
|
|---|
| 127 | 1d2 = drbjmax ! max distance between b and jet
|
|---|
| 128 | 1d2 = drajmax ! max distance between gamma and jet
|
|---|
| 129 | 1d2 = drjlmax ! max distance between jet and lepton
|
|---|
| 130 | 1d2 = drabmax ! max distance between gamma and b
|
|---|
| 131 | 1d2 = drblmax ! max distance between b and lepton
|
|---|
| 132 | 1d2 = dralmax ! maxdistance between gamma and lepton
|
|---|
| 133 | #*********************************************************************
|
|---|
| 134 | # Minimum and maximum invariant mass for pairs *
|
|---|
| 135 | #*********************************************************************
|
|---|
| 136 | 0 = mmjj ! min invariant mass of a jet pair
|
|---|
| 137 | 0 = mmbb ! min invariant mass of a b pair
|
|---|
| 138 | 0 = mmaa ! min invariant mass of gamma gamma pair
|
|---|
| 139 | 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
|
|---|
| 140 | 1d5 = mmjjmax ! max invariant mass of a jet pair
|
|---|
| 141 | 1d5 = mmbbmax ! max invariant mass of a b pair
|
|---|
| 142 | 1d5 = mmaamax ! max invariant mass of gamma gamma pair
|
|---|
| 143 | 1d5 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
|
|---|
| 144 | #*********************************************************************
|
|---|
| 145 | # Minimum and maximum invariant mass for all letpons *
|
|---|
| 146 | #*********************************************************************
|
|---|
| 147 | 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
|
|---|
| 148 | 1d5 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
|
|---|
| 149 | #*********************************************************************
|
|---|
| 150 | # Inclusive cuts *
|
|---|
| 151 | #*********************************************************************
|
|---|
| 152 | 0 = xptj ! minimum pt for at least one jet
|
|---|
| 153 | 0 = xptb ! minimum pt for at least one b
|
|---|
| 154 | 0 = xpta ! minimum pt for at least one photon
|
|---|
| 155 | 0 = xptl ! minimum pt for at least one charged lepton
|
|---|
| 156 | #*********************************************************************
|
|---|
| 157 | # Control the pt's of the jets sorted by pt *
|
|---|
| 158 | #*********************************************************************
|
|---|
| 159 | 0 = ptj1min ! minimum pt for the leading jet in pt
|
|---|
| 160 | 0 = ptj2min ! minimum pt for the second jet in pt
|
|---|
| 161 | 0 = ptj3min ! minimum pt for the third jet in pt
|
|---|
| 162 | 0 = ptj4min ! minimum pt for the fourth jet in pt
|
|---|
| 163 | 1d5 = ptj1max ! maximum pt for the leading jet in pt
|
|---|
| 164 | 1d5 = ptj2max ! maximum pt for the second jet in pt
|
|---|
| 165 | 1d5 = ptj3max ! maximum pt for the third jet in pt
|
|---|
| 166 | 1d5 = ptj4max ! maximum pt for the fourth jet in pt
|
|---|
| 167 | 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
|
|---|
| 168 | #*********************************************************************
|
|---|
| 169 | # Control the Ht(k)=Sum of k leading jets *
|
|---|
| 170 | #*********************************************************************
|
|---|
| 171 | 0 = htjmin ! minimum jet HT=Sum(jet pt)
|
|---|
| 172 | 1d5 = htjmax ! maximum jet HT=Sum(jet pt)
|
|---|
| 173 | 0 = ht2min ! minimum Ht for the two leading jets
|
|---|
| 174 | 0 = ht3min ! minimum Ht for the three leading jets
|
|---|
| 175 | 0 = ht4min ! minimum Ht for the four leading jets
|
|---|
| 176 | 1d5 = ht2max ! maximum Ht for the two leading jets
|
|---|
| 177 | 1d5 = ht3max ! maximum Ht for the three leading jets
|
|---|
| 178 | 1d5 = ht4max ! maximum Ht for the four leading jets
|
|---|
| 179 | #*********************************************************************
|
|---|
| 180 | # WBF cuts *
|
|---|
| 181 | #*********************************************************************
|
|---|
| 182 | 0 = xetamin ! minimum rapidity for two jets in the WBF case
|
|---|
| 183 | 0 = deltaeta ! minimum rapidity for two jets in the WBF case
|
|---|
| 184 | #*********************************************************************
|
|---|
| 185 | # maximal pdg code for quark to be considered as a jet *
|
|---|
| 186 | # otherwise b cuts are applied *
|
|---|
| 187 | #*********************************************************************
|
|---|
| 188 | 4 = maxjetflavor
|
|---|
| 189 | #*********************************************************************
|
|---|
| 190 | # Jet measure cuts *
|
|---|
| 191 | #*********************************************************************
|
|---|
| 192 | 0 = xqcut ! minimum kt jet measure between partons
|
|---|
| 193 | #*********************************************************************
|
|---|
