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)
|
---|