reweighting interrupted with KeyError : 'm1_'

Hi,

indeed interference term are not compatible with the new reweighting module.
I have to fix that. If I'm correct this should be working if you use 2.5.x
Would that work for you?

Cheers,

Olivier

> On 14 May 2018, at 14:32, Ilaria Brivio <email address hidden> wrote:
>
> New question #669004 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/669004
>
> Hi,
>
> I am trying to use the reweighting to change the coefficient of an effective operator in a simple process but I keep getting an error.
>
> I have the latest version of MadGraph (2.6.2) and I am generating a pure interference term p p > w+ w- NP^2==1 NPcW^2==1
>
> the event generation goes smooth but when I launch the reweighting I get the following output
>
> -------------------------------------------
> INFO: process: p p > w+ w- QCD=0 NP^2==1 NPcW^2==1
> INFO: options:
> INFO: Running Reweighting
> launch
> INFO: detected model: /home/ilaria/Dropbox/SMEFT_Feynrules/SMEFTsim_public/SMEFTsim_A_MwScheme_UFO_int_orders-massless. Loading...
> INFO: generating the square matrix element for reweighting
> INFO: generate p p > w+ w- QCD=0 NP^2==1 NPcW^2==1 ;
> INFO: Done 1.689
> Command "reweight run_01 -f" interrupted with error:
> KeyError : 'm1_'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/home/ilaria/Work/SMEFTsim_tests/Fsymmetric_MwScheme_intOrders/pp_ww_cW_int/run_01_tag_2_debug.log'.
> Please attach this file to your report.
> INFO: storing files of previous run
> gzipping output file: unweighted_events.lhe
> INFO: Done
> INFO:
> INFO: Original cross-section: -0.2196334 +- 0.002832886 pb
> INFO: Computed cross-section:
> quit
> INFO:
>
> -------------------------------------------
>
> the reweight card is just
>
> launch
> set frblock 21 1.5
>
> launch
> set frblock 21 0.5
>
>
> and the log file is below.
> Any idea of where this could be coming from?
>
> Thanks a lot!
> Ilaria
>
>
>
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.2 20xx-xx-xx *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #************************************************************
> reweight run_01 -f
> Traceback (most recent call last):
> File "/home/ilaria/Work/SMEFTsim_tests/Fsymmetric_MwScheme_intOrders/pp_ww_cW_int/bin/internal/extended_cmd.py", line 1501, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/ilaria/Work/SMEFTsim_tests/Fsymmetric_MwScheme_intOrders/pp_ww_cW_int/bin/internal/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/ilaria/Work/SMEFTsim_tests/Fsymmetric_MwScheme_intOrders/pp_ww_cW_int/bin/internal/common_run_interface.py", line 2138, in do_reweight
> reweight_cmd.import_command_file(path)
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/interface/extended_cmd.py", line 1643, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/various/misc.py", line 100, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/interface/reweight_interface.py", line 487, in do_launch
> self.load_module()
> File "/home/ilaria/Physics_Software/MG5_aMC_v2_6_2/madgraph/interface/reweight_interface.py", line 1781, in load_module
> hel = hel_dict[prefix]
> KeyError: 'm1_'
> Run Options
> -----------
> stdout_level : None
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None
> cluster_size : 100
> cluster_queue : None
> nb_core : 4 (user set)
> cluster_time : None
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> thepeg_path : None (user set)
> cluster_type : condor
> madanalysis5_path : /home/ilaria/Physics_Software/MG5_aMC_v2_6_2/HEPTools/madanalysis5/madanalysis5 (user set)
> cluster_nb_retry : 1
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : lhapdf-config
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : /home/ilaria/Physics_Software/MG5_aMC_v2_6_2/ExRootAnalysis (user set)
> timeout : 60
> cpp_compiler : None
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.2 2018-04-29 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set default_unset_couplings 99
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model sm
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model /home/ilaria/Dropbox/SMEFT_Feynrules/SMEFTsim_public/SMEF\
> Tsim_A_MwScheme_UFO_int_orders-massless
> generate p p > w+ w- QCD=0 NP^2==1 NPcW^2==1
> output pp_ww_cW_int
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR FRBLOCK
> ###################################
> Block frblock
> 1 1.000000e-99 # set of param :1*ceWPh, 1*ceBPh, 1*cuGPh, 1*cuWPh, 1*cuBPh, 1*cdGPh, 1*cdWPh, 1*cdBPh, 1*cHudPh, 1*ceHPh, 1*cuHPh, 1*cdHPh, 1*cledqPh, 1*cquqd1Ph, 1*cquqd8Ph, 1*clequ1Ph, 1*clequ3Ph
> 18 1.000000e+03 # LambdaSMEFT
> 19 1.000000e+00 # cG
> 20 1.000000e+00 # cGtil
> 21 1.000000e+00 # cW
> 22 1.000000e+00 # cWtil
> 23 1.000000e+00 # cH
> 24 1.000000e+00 # cHbox
> 25 1.000000e+00 # cHDD
> 26 1.000000e+00 # cHG
> 27 1.000000e+00 # cHGtil
> 28 1.000000e+00 # cHW
> 29 1.000000e+00 # cHWtil
> 30 1.000000e+00 # cHB
> 31 1.000000e+00 # cHBtil
> 32 1.000000e+00 # cHWB
> 33 1.000000e+00 # cHWBtil
> 34 1.000000e+00 # ceHAbs
> 35 1.000000e+00 # cuHAbs
> 36 1.000000e+00 # cdHAbs
> 37 1.000000e+00 # ceWAbs
> 38 1.000000e+00 # ceBAbs
> 39 1.000000e+00 # cuGAbs
> 40 1.000000e+00 # cuWAbs
> 41 1.000000e+00 # cuBAbs
> 42 1.000000e+00 # cdGAbs
> 43 1.000000e+00 # cdWAbs
> 44 1.000000e+00 # cdBAbs
> 45 1.000000e+00 # cHl1
> 46 1.000000e+00 # cHl3
> 47 1.000000e+00 # cHe
> 48 1.000000e+00 # cHq1
> 49 1.000000e+00 # cHq3
> 50 1.000000e+00 # cHu
> 51 1.000000e+00 # cHd
> 52 1.000000e+00 # cHudAbs
> 53 1.000000e+00 # cll
> 54 1.000000e+00 # cll1
> 55 1.000000e+00 # cqq1
> 56 1.000000e+00 # cqq11
> 57 1.000000e+00 # cqq3
> 58 1.000000e+00 # cqq31
> 59 1.000000e+00 # clq1
> 60 1.000000e+00 # clq3
> 61 1.000000e+00 # cee
> 62 1.000000e+00 # cuu
> 63 1.000000e+00 # cuu1
> 64 1.000000e+00 # cdd
> 65 1.000000e+00 # cdd1
> 66 1.000000e+00 # ceu
> 67 1.000000e+00 # ced
> 68 1.000000e+00 # cud1
> 69 1.000000e+00 # cud8
> 70 1.000000e+00 # cle
> 71 1.000000e+00 # clu
> 72 1.000000e+00 # cld
> 73 1.000000e+00 # cqe
> 74 1.000000e+00 # cqu1
> 75 1.000000e+00 # cqu8
> 76 1.000000e+00 # cqd1
> 77 1.000000e+00 # cqd8
> 78 1.000000e+00 # cledqAbs
> 79 1.000000e+00 # cquqd1Abs
> 80 1.000000e+00 # cquqd8Abs
> 81 1.000000e+00 # clequ1Abs
> 82 1.000000e+00 # clequ3Abs
> 83 8.038700e+01 # MW0
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.180000e+00 # MB
> 6 1.732000e+02 # MT
> 23 9.118760e+01 # MZ
> 25 1.250900e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.000000 # d : 0.0
> 2 0.000000 # u : 0.0
> 3 0.000000 # s : 0.0
> 4 0.000000 # c : 0.0
> 11 0.000000 # e- : 0.0
> 12 0.000000 # ve : 0.0
> 13 0.000000 # mu- : 0.0
> 14 0.000000 # vm : 0.0
> 15 0.000000 # ta- : 0.0
> 16 0.000000 # vt : 0.0
> 21 0.000000 # g : 0.0
> 22 0.000000 # a : 0.0
> 24 80.387000 # w+ : MW0
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 2 1.166379e-05 # Gf
> 3 1.181000e-01 # aS
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.180000e+00 # ymb
> 6 1.732000e+02 # ymt
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.508336e+00 # WT
> DECAY 23 2.495200e+00 # WZ
> DECAY 24 2.085000e+00 # WW
> DECAY 25 4.070000e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.000000 # d : 0.0
> DECAY 2 0.000000 # u : 0.0
> DECAY 3 0.000000 # s : 0.0
> DECAY 4 0.000000 # c : 0.0
> DECAY 5 0.000000 # b : 0.0
> DECAY 11 0.000000 # e- : 0.0
> DECAY 12 0.000000 # ve : 0.0
> DECAY 13 0.000000 # mu- : 0.0
> DECAY 14 0.000000 # vm : 0.0
> DECAY 15 0.000000 # ta- : 0.0
> DECAY 16 0.000000 # vt : 0.0
> DECAY 21 0.000000 # g : 0.0
> DECAY 22 0.000000 # a : 0.0
> #*********************************************************************
> # 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 *
> #*********************************************************************
> #
> #*******************
> # Running parameters
> #*******************
> #
> #*********************************************************************
> # 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 *
> # If you want to run Pythia, avoid more than 50k events in a run. *
> #*********************************************************************
> 10000 = 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
> # To see polarised beam options: type "update beam_pol"
> #*********************************************************************
> # 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
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = 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)
> 3.0 = lhe_version ! Change the way clustering information pass to shower.
> True = clusinfo ! include clustering tag in output
> 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)
> #*********************************************************************
> # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #*********************************************************************
> False = cut_decays ! Cut decay products
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20.0 = ptj ! minimum pt for the jets
> 0.0 = ptb ! minimum pt for the b
> 10.0 = pta ! minimum pt for the photons
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptbmax ! maximum pt for the b
> -1.0 = ptamax ! maximum pt for the photons
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = 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) *
> #*********************************************************************
> 5.0 = etaj ! max rap for the jets
> -1.0 = etab ! max rap for the b
> 2.5 = etaa ! max rap for the photons
> 2.5 = etal ! max rap for the charged leptons
> 0.0 = etajmin ! min rap for the jets
> 0.0 = etabmin ! min rap for the b
> 0.0 = etaamin ! min rap for the photons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = 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 *
> #*********************************************************************
> 0.4 = drjj ! min distance between jets
> 0.0 = drbb ! min distance between b's
> 0.4 = drll ! min distance between leptons
> 0.4 = draa ! min distance between gammas
> 0.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.0 = drab ! min distance between gamma and b
> 0.0 = drbl ! min distance between b and lepton
> 0.4 = dral ! min distance between gamma and lepton
> -1.0 = drjjmax ! max distance between jets
> -1.0 = drbbmax ! max distance between b's
> -1.0 = drllmax ! max distance between leptons
> -1.0 = draamax ! max distance between gammas
> -1.0 = drbjmax ! max distance between b and jet
> -1.0 = drajmax ! max distance between gamma and jet
> -1.0 = drjlmax ! max distance between jet and lepton
> -1.0 = drabmax ! max distance between gamma and b
> -1.0 = drblmax ! max distance between b and lepton
> -1.0 = dralmax ! maxdistance between gamma and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> # WARNING: for four lepton final state mmll cut require to have *
> # different lepton masses for each flavor! *
> #*********************************************************************
> 0.0 = mmjj ! min invariant mass of a jet pair
> 0.0 = mmbb ! min invariant mass of a b pair
> 0.0 = mmaa ! min invariant mass of gamma gamma pair
> 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1.0 = mmjjmax ! max invariant mass of a jet pair
> -1.0 = mmbbmax ! max invariant mass of a b pair
> -1.0 = mmaamax ! max invariant mass of gamma gamma pair
> -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> {} = 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.
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons *
> #*********************************************************************
> 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0.0 = ptheavy ! minimum pt for at least one heavy final state
> 0.0 = xptj ! minimum pt for at least one jet
> 0.0 = xptb ! minimum pt for at least one b
> 0.0 = xpta ! minimum pt for at least one photon
> 0.0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0.0 = ptj1min ! minimum pt for the leading jet in pt
> 0.0 = ptj2min ! minimum pt for the second jet in pt
> 0.0 = ptj3min ! minimum pt for the third jet in pt
> 0.0 = ptj4min ! minimum pt for the fourth jet in pt
> -1.0 = ptj1max ! maximum pt for the leading jet in pt
> -1.0 = ptj2max ! maximum pt for the second jet in pt
> -1.0 = ptj3max ! maximum pt for the third jet in pt
> -1.0 = 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.0 = ptl1min ! minimum pt for the leading lepton in pt
> 0.0 = ptl2min ! minimum pt for the second lepton in pt
> 0.0 = ptl3min ! minimum pt for the third lepton in pt
> 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1.0 = ptl1max ! maximum pt for the leading lepton in pt
> -1.0 = ptl2max ! maximum pt for the second lepton in pt
> -1.0 = ptl3max ! maximum pt for the third lepton in pt
> -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
> 0.0 = ihtmin !inclusive Ht for all partons (including b)
> -1.0 = ihtmax !inclusive Ht for all partons (including b)
> 0.0 = ht2min ! minimum Ht for the two leading jets
> 0.0 = ht3min ! minimum Ht for the three leading jets
> 0.0 = ht4min ! minimum Ht for the four leading jets
> -1.0 = ht2max ! maximum Ht for the two leading jets
> -1.0 = ht3max ! maximum Ht for the three leading jets
> -1.0 = ht4max ! maximum Ht for the four leading jets
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442 *
> # When ptgmin=0, all the other parameters are ignored *
> # When ptgmin>0, pta and draj are not going to be used *
> #***********************************************************************
> 0.0 = ptgmin ! Min photon transverse momentum
> 0.4 = R0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> True = isoEM ! isolate photons from EM energy (photons and leptons)
> #*********************************************************************
> # WBF cuts *
> #*********************************************************************
> 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
> 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #***********************************************************************
> # Turn on either the ktdurham or ptlund cut to activate *
> # CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
> #***********************************************************************
> -1.0 = ktdurham
> 0.4 = dparameter
> -1.0 = ptlund
> 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above
> #*********************************************************************
> # 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
> #
> none = 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'
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier!

thanks for the answer. Indeed 2.5 works smooth.
Good to know it's just a temporary issue. ;)

Cheers,
Ilaria

and here is a diff which shoud work for 2.6.x version:

Cheers,

Olivier

=== modified file 'madgraph/iolibs/template_files/matrix_standalone_splitOrders_v4.inc'
--- madgraph/iolibs/template_files/matrix_standalone_splitOrders_v4.inc 2017-07-31 22:08:40 +0000
+++ madgraph/iolibs/template_files/matrix_standalone_splitOrders_v4.inc 2018-05-17 11:30:17 +0000
@@ -109,7 +109,7 @@
 C
       INTEGER NHEL(NEXTERNAL,NCOMB)
 %(helicity_lines)s
- COMMON/%(proc_prefix)sBORN_HEL_CONFIGS/NHEL
+ COMMON/%(proc_prefix)sPROCESS_NHEL/NHEL

       INTEGER USERHEL
    DATA USERHEL/-1/
@@ -357,7 +357,7 @@
 C GLOBALS
 C
       INTEGER HELC(NEXTERNAL,NCOMB)
- COMMON/%(proc_prefix)sBORN_HEL_CONFIGS/HELC
+ COMMON/%(proc_prefix)sPROCESS_NHEL/HELC

       INTEGER POLARIZATIONS(0:NEXTERNAL,0:5)
       COMMON/%(proc_prefix)sBORN_BEAM_POL/POLARIZATIONS

=== modified file 'madgraph/iolibs/template_files/matrix_standalone_v4.inc'
--- madgraph/iolibs/template_files/matrix_standalone_v4.inc 2017-07-31 22:08:40 +0000
+++ madgraph/iolibs/template_files/matrix_standalone_v4.inc 2018-05-17 11:31:04 +0000
@@ -275,7 +275,7 @@
 C GLOBALS
 C
       INTEGER HELC(NEXTERNAL,NCOMB)
- COMMON/%(proc_prefix)sBORN_HEL_CONFIGS/HELC
+ COMMON/%(proc_prefix)sPROCESS_NHEL/HELC

       INTEGER POLARIZATIONS(0:NEXTERNAL,0:5)

> On 17 May 2018, at 10:46, Ilaria Brivio <email address hidden> wrote:
>
> Question #669004 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/669004
>
> Status: Answered => Solved
>
> Ilaria Brivio confirmed that the question is solved:
> Hi Olivier!
>
> thanks for the answer. Indeed 2.5 works smooth.
> Good to know it's just a temporary issue. ;)
>
> Cheers,
> Ilaria
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.