HEFT with matching (repeat arxiv:1110.1728)
Hi Johan
I want to repeat HEFT with matching in fig1 of the paper 1110.1728. I run madevent+pythia level.
My fig in the parton level has some question:
1, the horizontal ordinate is pt(h1), not pt(h), in the plot_card.dat, I set the param as follow
jet 1 -1 2 -2 3 -3 4 -4 21 # Class number 1
b 5 -5 # Class number 2
h 25 # Class number 3
mET 12 -12 14 -14 16 -16 1000022 # Missing ET class, name is reserved
pt 3 1
pt 5 0 100 # bin size, min value, max value
2, In my fig, the cross section is 0 in pt(h1)<30GeV. I set the run_card.dat
#**
# Tag name for the run (one word) *
#******
scalar = run_tag ! name of the run
#******
# Run to generate the grid pack *
#******
.false. = gridpack !True = setting up the grid pack
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
30000 = nevents ! Number of unweighted events requested
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
#******
1 = lpp1 ! beam 1 type (0=NO PDF)
-1 = lpp2 ! beam 2 type (0=NO PDF)
980 = ebeam1 ! beam 1 energy in GeV
980 = ebeam2 ! beam 2 energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0 = polbeam1 ! beam polarization for beam 1
0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
'cteq6l1' = pdlabel ! PDF set
#******
# Renormalization and factorization scales *
#******
F = fixed_ren_scale ! if .true. use fixed ren scale
F = fixed_fac_scale ! if .true. use fixed fac scale
91.1880 = scale ! fixed ren scale
91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
1 = scalefact ! scale factor for event-by-event scales
#******
# Matching - Warning! ickkw > 1 is still beta
#******
1 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
1 = highestmult ! for ickkw=2, highest mult group
1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
1 = alpsfact ! scale factor for QCD emission vx
F = chcluster ! cluster only according to channel diag
T = pdfwgt ! for ickkw=1, perform pdf reweighting
#******
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#******
T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#******
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/
#******
T = cut_decays ! Cut decay products
#******
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#******
0 = nhel ! Number of helicities used per event
#******
# Standard Cuts
#******
#
#******
# Minimum and maximum pt's *
#******
30 = ptj ! minimum pt for the jets
30 = ptb ! minimum pt for the b
10 = pta ! minimum pt for the photons
10 = ptl ! minimum pt for the charged leptons
0 = misset ! minimum missing Et (sum of neutrino's momenta)
0 = ptheavy ! minimum pt for one heavy final state
1.0 = ptonium ! minimum pt for the quarkonium states
1d5 = ptjmax ! maximum pt for the jets
1d5 = ptbmax ! maximum pt for the b
1d5 = ptamax ! maximum pt for the photons
1d5 = ptlmax ! maximum pt for the charged leptons
1d5 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#******
# Minimum and maximum E's (in the lab frame) *
#******
0 = ej ! minimum E for the jets
0 = eb ! minimum E for the b
0 = ea ! minimum E for the photons
0 = el ! minimum E for the charged leptons
1d5 = ejmax ! maximum E for the jets
1d5 = ebmax ! maximum E for the b
1d5 = eamax ! maximum E for the photons
1d5 = elmax ! maximum E for the charged leptons
#******
# Maximum and minimum rapidity *
#******
4.5 = etaj ! max rap for the jets
1d2 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.6 = etaonium ! max rap for the quarkonium states
0d0 = etajmin ! min rap for the jets
0d0 = etabmin ! min rap for the b
0d0 = etaamin ! min rap for the photons
0d0 = etalmin ! main rap for the charged leptons
#******
# Minimum and maximum DeltaR distance *
#******
0 = drjj ! min distance between jets
0 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.4 = drjl ! min distance between jet and lepton
0 = drab ! min distance between gamma and b
0 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
1d2 = drjjmax ! max distance between jets
1d2 = drbbmax ! max distance between b's
1d2 = drllmax ! max distance between leptons
1d2 = draamax ! max distance between gammas
1d2 = drbjmax ! max distance between b and jet
1d2 = drajmax ! max distance between gamma and jet
1d2 = drjlmax ! max distance between jet and lepton
1d2 = drabmax ! max distance between gamma and b
1d2 = drblmax ! max distance between b and lepton
1d2 = dralmax ! maxdistance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
#******
0 = mmjj ! min invariant mass of a jet pair
0 = mmbb ! min invariant mass of a b pair
0 = mmaa ! min invariant mass of gamma gamma pair
0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
1d5 = mmjjmax ! max invariant mass of a jet pair
1d5 = mmbbmax ! max invariant mass of a b pair
1d5 = mmaamax ! max invariant mass of gamma gamma pair
1d5 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#******
# Minimum and maximum invariant mass for all letpons *
#******
0 = mmnl ! min invariant mass for all letpons (l+- and vl)
1d5 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
1d5 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#******
# Inclusive cuts *
#******
0 = xptj ! minimum pt for at least one jet
0 = xptb ! minimum pt for at least one b
0 = xpta ! minimum pt for at least one photon
0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0 = ptj1min ! minimum pt for the leading jet in pt
0 = ptj2min ! minimum pt for the second jet in pt
0 = ptj3min ! minimum pt for the third jet in pt
0 = ptj4min ! minimum pt for the fourth jet in pt
1d5 = ptj1max ! maximum pt for the leading jet in pt
1d5 = ptj2max ! maximum pt for the second jet in pt
1d5 = ptj3max ! maximum pt for the third jet in pt
1d5 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0 = ptl1min ! minimum pt for the leading lepton in pt
0 = ptl2min ! minimum pt for the second lepton in pt
0 = ptl3min ! minimum pt for the third lepton in pt
0 = ptl4min ! minimum pt for the fourth lepton in pt
1d5 = ptl1max ! maximum pt for the leading lepton in pt
1d5 = ptl2max ! maximum pt for the second lepton in pt
1d5 = ptl3max ! maximum pt for the third lepton in pt
1d5 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0 = htjmin ! minimum jet HT=Sum(jet pt)
1d5 = htjmax ! maximum jet HT=Sum(jet pt)
0 = ihtmin !inclusive Ht for all partons (including b)
1d5 = ihtmax !inclusive Ht for all partons (including b)
0 = ht2min ! minimum Ht for the two leading jets
0 = ht3min ! minimum Ht for the three leading jets
0 = ht4min ! minimum Ht for the four leading jets
1d5 = ht2max ! maximum Ht for the two leading jets
1d5 = ht3max ! maximum Ht for the three leading jets
1d5 = ht4max ! maximum Ht for the four leading jets
#******
# WBF cuts *
#******
0 = xetamin ! minimum rapidity for two jets in the WBF case
0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
# Jet measure cuts *
#******
30 = xqcut ! minimum kt jet measure between partons
#******
The pythia_card.dat is
#******
!...Parton showering on or off
MSTP(61)=1
MSTP(71)=1
!...Fragmentati
MSTJ(1)=1
!...Multiple interactions on or off
MSTP(81)=20
!...Don't stop execution after 10 errors
MSTU(21)=1
!...Cutoff in jet measure for matching
QCUT = 50
#******
I dont know where the error exist.
Thanks for you help
violin
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Johan Alwall
- Solved:
- Last query:
- Last reply: