| 1 | ! This is a Pythia8 ini file (to set configuration flags)
|
|---|
| 2 | !
|
|---|
| 3 | ! This file should be used to implement Pythia8 jet matching from an LHE generated by MadGraph
|
|---|
| 4 | !
|
|---|
| 5 | ! Since the parameter (JetMatching:setmad) is not quite sophisticated enough to suck in
|
|---|
| 6 | ! all possible parameters from the LHE header, it is disabled
|
|---|
| 7 | !
|
|---|
| 8 | ! This means that this file will need to be reset before each new analysis is begun
|
|---|
| 9 | !
|
|---|
| 10 |
|
|---|
| 11 | ! <Standard Parameters> No need to change
|
|---|
| 12 |
|
|---|
| 13 | ! 1) Settings used in the main program.
|
|---|
| 14 | Main:timesAllowErrors = 10 ! how many aborts before run stops
|
|---|
| 15 | Main:spareMode1 = 0 ! skip n events at beginning of file
|
|---|
| 16 | Beams:frameType = 4 ! pp collision
|
|---|
| 17 |
|
|---|
| 18 | ! 2) Settings related to output in init(), next() and stat().
|
|---|
| 19 | Init:showChangedSettings = on ! list changed settings
|
|---|
| 20 | Init:showChangedParticleData = on ! list changed particle data
|
|---|
| 21 | Next:numberCount = 0 ! print message every n events
|
|---|
| 22 | Next:numberShowInfo = 1 ! print event information n times
|
|---|
| 23 | Next:numberShowProcess = 1 ! print process record n times
|
|---|
| 24 | Next:numberShowEvent = 1 ! print event record n times
|
|---|
| 25 |
|
|---|
| 26 | ! 3) Settings related to decay
|
|---|
| 27 | ParticleDecays:limitCylinder = on
|
|---|
| 28 | ParticleDecays:xyMax = 1290 ! in mm
|
|---|
| 29 | ParticleDecays:zMax = 3000 ! in mm
|
|---|
| 30 |
|
|---|
| 31 | ! 4) Settings related to matching
|
|---|
| 32 | JetMatching:merge = on
|
|---|
| 33 | JetMatching:setMad = off ! Do NOT set parameters (like qCut) from LHE header
|
|---|
| 34 | ! In both MadGraph and Pythia, qCut (xqcut) represents the pt_rel of the emission to the jet
|
|---|
| 35 |
|
|---|
| 36 | ! </StandardParameters>
|
|---|
| 37 |
|
|---|
| 38 |
|
|---|
| 39 | ! CHANGE BELOW
|
|---|
| 40 | Beams:LHEF = testing/unweighted_events.lhe
|
|---|
| 41 | Main:numberOfEvents = 10000 ! number of events to generate (0 for all)
|
|---|
| 42 | JetMatching:nJetMax = 2 ! The maximum jet multiplicity in the LHE file
|
|---|
| 43 | JetMatching:exclusive = 2 ! 2 = exclusive jet matching (1:1 parton:jet required) except for maximum jet multiplicity (Pythia fills in extra jet phase space), 1 = inclusive
|
|---|
| 44 | JetMatching:jetAllow = 1 ! 1 = Top quarks, leptons and photons are not passed to the jet clustering algorithm
|
|---|
| 45 |
|
|---|
| 46 | ! MadGraph style matching
|
|---|
| 47 | JetMatching:scheme = 1 ! 1 = MadGraph style matching
|
|---|
| 48 |
|
|---|
| 49 | !JetMatching:jetAlgorithm = 2 ! SlowJet
|
|---|
| 50 | !JetMatching:slowJetPower = 1 ! kt
|
|---|
| 51 |
|
|---|
| 52 | JetMatching:coneRadius = 1.0 ! It is not recommended to change this
|
|---|
| 53 | JetMatching:etaJetMax = 10.0 ! It is not recommended to change this
|
|---|
| 54 |
|
|---|
| 55 | JetMatching:qCut = 300 ! The qCut used by Pythia. ptMin is automatically set to qCut
|
|---|
| 56 | JetMatching:qCutME = 300 ! The qCut that was used by MadGraph (if doShowerKt, qCut=qCutME is OK)
|
|---|
| 57 |
|
|---|
| 58 | JetMatching:nQmatch = 4 ! The max absPDGid to be treated as a light quark (and thus participate in matching)
|
|---|
| 59 | JetMatching:clFact = 1 ! Jet emission vertex factor
|
|---|
| 60 | JetMatching:doShowerKt = on ! This checks for event veto after the first shower emission, so "the same value can be used for qCut and xqcut" <https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Matching>
|
|---|
| 61 |
|
|---|
| 62 | ! Settings to be Used for MadGraph
|
|---|
| 63 | ! 1 = ickkw ! 0 for no matching, 1 for MLM matching and 2 for CKKW matching
|
|---|
| 64 | ! ? = xqcut ! minimum jet measure (pT/kT) for QCD partons, if xqcut=0 use cone jet matching, if xqcut> 0 use kt jet matching. This value should be related to the hard scale
|
|---|
| 65 | ! 1 = highestmult ! for ickkw=2, highest mult group (leave as default)
|
|---|
| 66 | ! 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE (leave as default)
|
|---|
| 67 | ! 1 = alpsfact ! scale factor for QCD emission vx
|
|---|
| 68 | ! F = chcluster ! cluster only according to channel diag (leave as default)
|
|---|
| 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 |
|
|---|
| 73 | ! F = fixed_ren_scale
|
|---|
| 74 | ! 0 = drjj !Set to 0 when using kt-scheme
|
|---|
| 75 | ! 0 = drbb !Set to 0 when using kt-scheme
|
|---|
| 76 | ! 0 = drbj !Set to 0 when using kt-scheme
|
|---|
| 77 | ! 0 = ptj !Set to 0 when using kt-scheme (or xqcut to imporve speed)
|
|---|
| 78 | ! 0 = ptb !Set to 20% below min cut
|
|---|
| 79 | ! 0 = mjj !essentially drjj, (so may be xqcut to improve speed)
|
|---|
| 80 | ! etaj = -1 !MadGraph suggests 5. If your hard event is constraining the phase space, then perhaps remove any limit
|
|---|
| 81 | ! etab = 5 !Set to slightly below min cut
|
|---|
| 82 | ! 4 = maxjetflavor ! (default 4) Defines which partons are considered as "j" for the matching procedure. Only use 5 when bottom quarks are expected to emissions from some other jet (or if you expect a lot of them coming out of the proton in t-channel production)
|
|---|
