MadGraph5_aMC@NLO

Need a newer LHE version.

Asked by tom ingebretsen carlson

Hi,

I am generating and using MG 3.1 or 3.2, however I am in need of a newer LHE version. The LHE file says:
<LesHouchesEvents version="1.0">
While I need version 3. Is there some way to get an output with LHE version 3? Or is there some way to convert a .root file to LHE version 3?

In case you want some more backround info. I am generating with MG and using MGs internal reweight function. Then I want to shower the samples, but separately from MG, however the weights from the reweighting are not stored. According to some experts it is because the lhe version is too old.

Thanks!

All best,

Tom

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
tom ingebretsen carlson (makronymen) said : 		#1

Some additional information.

I use SMEFT@NLO (https://feynrules.irmp.ucl.ac.be/wiki/SMEFTatNLO) and then I get LHE version 1, while if I don't use it I get LHE version 1.

Thanks again.

All best,

Tom

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#2

Hi,

The default format for our lhef output is version 3.0.
Are you sure that you do not use an old run_card.dat where you do request the format to be version 1?

If not maybe some tool that you use are not compatible with version 3.0 and therefore we do output version 1 for compatibility with such tools (like the old pythia-pgs package)

Cheers,

Olivier

> On 5 Aug 2022, at 16:25, tom ingebretsen carlson <email address hidden> wrote:
>
> Question #702670 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/702670
>
> tom ingebretsen carlson posted a new comment:
> Some additional information.
>
> I use SMEFT@NLO (https://feynrules.irmp.ucl.ac.be/wiki/SMEFTatNLO) and
> then I get LHE version 1, while if I don't use it I get LHE version 1.
>
>
> Thanks again.
>
>
> All best,
>
> Tom
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
tom ingebretsen carlson (makronymen) said : 		#3

Hi,

yes I actully just realized I use an old run_card (which you can see below).

Where would I need to change it to get LHE version 3?

Thanks!

All best,

Tom

#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
  100000 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
  1 = lpp1 ! beam 1 type
  1 = lpp2 ! beam 2 type
  6500.0 = ebeam1 ! beam 1 total energy in GeV
  6500.0 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
     0.0 = polbeam1 ! beam polarization for beam 1
     0.0 = polbeam2 ! beam polarization for beam 2

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
  nn23lo1 = pdlabel ! PDF set
  230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
  False = fixed_ren_scale ! if .true. use fixed ren scale
  False = fixed_fac_scale ! if .true. use fixed fac scale
  125.0 = scale ! fixed ren scale
  125.0 = dsqrt_q2fact1 ! fixed fact scale for pdf1
  125.0 = dsqrt_q2fact2 ! fixed fact scale for pdf2
  -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
  1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM
 1.0 = alpsfact ! scale factor for QCD emission vx
 False = chcluster ! cluster only according to channel diag
 4 = asrwgtflavor ! highest quark flavor for a_s reweight
 False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
                                   ! (turn off for VBF and single top processes)
 0.0 = xqcut ! minimum kt jet measure between partons

#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
  0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
  None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
  {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
  {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
  {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0.0 = ej ! minimum E for the jets
  0.0 = eb ! minimum E for the b
  0.0 = ea ! minimum E for the photons
  0.0 = el ! minimum E for the charged leptons
  -1.0 = ejmax ! maximum E for the jets
 -1.0 = ebmax ! maximum E for the b
 -1.0 = eamax ! maximum E for the photons
 -1.0 = elmax ! maximum E for the charged leptons
 {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
 {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})

#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
  {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
  {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
  {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
  {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
  0.0 = ptheavy ! minimum pt for at least one heavy final state
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
  4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
  False = use_syst ! Enable systematics studies
#
  systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'#*********************************************************************
# Additional hidden parameters
#*********************************************************************
  10.0 = ptl # hidden_parameter
  20.0 = ptj # hidden_parameter
  0.4 = drll # hidden_parameter
  -1.0 = ptj4max # hidden_parameter
  460.0 = ptb # hidden_parameter
  10.0 = pta # hidden_parameter
  0.0 = mmnl # hidden_parameter
  0.0 = ht3min # hidden_parameter
  -1.0 = draamax # hidden_parameter
  True = isoem # hidden_parameter
  True = clusinfo # hidden_parameter
  -1.0 = ptl1max # hidden_parameter
  0.0 = ptj4min # hidden_parameter
  0.0 = mmaa # hidden_parameter
  1.0 = epsgamma # hidden_parameter
  0.0 = xetamin # hidden_parameter
  0.4 = r0gamma # hidden_parameter
  0.0 = misset # hidden_parameter
  -1.0 = drjjmax # hidden_parameter
  -1.0 = ptl3max # hidden_parameter
  0.0 = etajmin # hidden_parameter
  -1.0 = ptlmax # hidden_parameter
  0.0 = ptl1min # hidden_parameter
  0.4 = drjl # hidden_parameter
  False = cut_decays # hidden_parameter
  0.0 = ptl2min # hidden_parameter
  2.0 = lhe_version # hidden_parameter
  -1.0 = drabmax # hidden_parameter
  0.4 = dparameter # hidden_parameter
  0.0 = mmjj # hidden_parameter
  0.4 = drjj # hidden_parameter
  0.0 = ptl3min # hidden_parameter
  -1.0 = drblmax # hidden_parameter
  -1.0 = mmjjmax # hidden_parameter
  0.0 = drbl # hidden_parameter
  0.0 = drbj # hidden_parameter
  -1.0 = drbbmax # hidden_parameter
  0.0 = drbb # hidden_parameter
  0.0 = htjmin # hidden_parameter
  0.0 = mmll # hidden_parameter
  0.0 = ptllmin # hidden_parameter
  -1.0 = htjmax # hidden_parameter
  -1.0 = missetmax # hidden_parameter
  -1.0 = ptllmax # hidden_parameter
  2.5 = etaa # hidden_parameter
  -1.0 = etab # hidden_parameter
  2.5 = etal # hidden_parameter
  5.0 = etaj # hidden_parameter
  -1.0 = drajmax # hidden_parameter
  -1.0 = ihtmax # hidden_parameter
  0.0 = ihtmin # hidden_parameter
  0.0 = xpta # hidden_parameter
  0.0 = xptb # hidden_parameter
  0.0 = xptj # hidden_parameter
  -1.0 = mmbbmax # hidden_parameter
  0.0 = xptl # hidden_parameter
  -1.0 = ktdurham # hidden_parameter
  -1.0 = ptjmax # hidden_parameter
  -1.0 = ht3max # hidden_parameter
  -1.0 = mmaamax # hidden_parameter
  1.0 = xn # hidden_parameter
  -1.0 = drjlmax # hidden_parameter
  0.0 = ht4min # hidden_parameter
  -1.0 = ht4max # hidden_parameter
  -1.0 = ptbmax # hidden_parameter
  -1.0 = drllmax # hidden_parameter
  -1.0 = ht2max # hidden_parameter
  0.0 = ptj2min # hidden_parameter
  0.0 = ht2min # hidden_parameter
  0.0 = etabmin # hidden_parameter
  0.0 = deltaeta # hidden_parameter
  0.0 = etalmin # hidden_parameter
  -1.0 = drbjmax # hidden_parameter
  -1.0 = mmnlmax # hidden_parameter
  -1.0 = ptj2max # hidden_parameter
  -1.0 = ptlund # hidden_parameter
  ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments # Choose the argment to pass to the systematics command. like --mur=0.25,1,4. Look at the help of the systematics function for more details.
  0.0 = ptl4min # hidden_parameter
  -1.0 = ptl2max # hidden_parameter
  -1.0 = ptj3max # hidden_parameter
  0.0 = ptj3min # hidden_parameter
  -1.0 = ptl4max # hidden_parameter
  0 = cutuse # hidden_parameter
  -1.0 = ptj1max # hidden_parameter
  -1.0 = ptamax # hidden_parameter
  0.0 = ptj1min # hidden_parameter
  0.0 = ptgmin # hidden_parameter
  0.0 = etaamin # hidden_parameter
  0.4 = draa # hidden_parameter
  0.0 = drab # hidden_parameter
  0.4 = dral # hidden_parameter
  0.4 = draj # hidden_parameter
  -1.0 = mmllmax # hidden_parameter
  -1.0 = dralmax # hidden_parameter
  0.0 = mmbb # hidden_parameter
  [1, 2, 3, 4, 5, 6, 21] = pdgs_for_merging_cut # hidden_parameter

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said : 		#4

You should change the line
 2.0 = lhe_version # hidden_parameter

to
3.0 = lhe_version

Cheers,

Olivier

> On 5 Aug 2022, at 16:45, tom ingebretsen carlson <email address hidden> wrote:
>
> 2.0 = lhe_version # hidden_parameter

Revision history for this message
tom ingebretsen carlson (makronymen) said : 		#5

Thanks Olivier Mattelaer, that solved my question.

Revision history for this message
tom ingebretsen carlson (makronymen) said : 		#6

Oh, right!

Thanks for your help!

All best,

Tom