HEFT with matching (repeat arxiv:1110.1728)

Asked by vaborla

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=default))
#*********************************************************************
# 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*Gamma)
#**********************************
  15 = bwcutoff ! (M+/-bwcutoff*Gamma)
#**********************************************************
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#**********************************************************
   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

!...Fragmentation/hadronization on or off
      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
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Johan Alwall
Solved:
Last query:
Last reply:
Revision history for this message
Johan Alwall (johan-alwall) said :
#1

Dear Violin,

The figure in the paper is after matching, not at parton level. Distributions at parton level have no meaning when you combine different multiplicity samples, since you have double counting and also don't regulate the collinear QCD divergencies. The Higgs pT below the xqcut scale is given purely by the parton shower, so at parton level it must be 0. Above the cutoff, the distribution is given by the matrix element but needs the additional Sudakov suppression from the matching procedure. Therefore, also there you need to look at the distributions after matching (i.e., after Pythia).

All the best,
Johan

Revision history for this message
vaborla (vaborla-tt) said :
#2

Thanks for your fast answer!
After Pythia, the events file could not give the distribution of higgs PT.
May I have other param errors?

Thanks again

Best wishes
violin

Revision history for this message
Best Johan Alwall (johan-alwall) said :
#3

Hello Violin,

In the pythia_events.lhe file, the Higgs is included as an intermediate propagator. You can plot this with MadAnalysis by setting plot_decayed to yes. However, I would suggest to use some more advanced event analysis package, such as MadAnalysis 5, for your analysis. Also note that the Pythia .lhe file doesn't include the full Pythia information, so for more proper analysis you want to use the full .hep file.

All the best,
Johan

Revision history for this message
vaborla (vaborla-tt) said :
#4

Thanks Johan Alwall, that solved my question.