multi_run Merge Issue

i have tried a quick run:

generate p p > t t~
output
launch -i
multi_run 2

and it worked can you check the same?

Olivier

Hello,

Well, now I am very confused because that also worked for me.

I have also ran my script for 10k events and a multi_run 2 and that worked as well. However, once I do 100k events and multi_run 10 it again encounters the "FileNotFound Error" as before.

Thanks for the help!
Torrey

Hello,

Since then I have tried a variety of things.

I have tried with multi_run 2:
100k runs (failed)
90k runs (failed)
50k runs (failed)
10k runs (succeeded)

Below is the log of my latest attempt with 50k,

----------------------------------------------------------------------------------------------------------------------------------------------------

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.2.0 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 /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/EFTML/Events/run_10_0/unweighted_events.lhe.gz_1.lhe -from_cards --multicore=wait
Traceback (most recent call last):
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/EFTML/bin/internal/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/EFTML/bin/internal/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/EFTML/bin/internal/common_run_interface.py", line 2234, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1692, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/reweight_interface.py", line 515, in do_launch
    param_card_iterator, tag_name = self.handle_param_card(model_line, args, type_rwgt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/reweight_interface.py", line 759, in handle_param_card
    tmp_card = new_card.lower().split('block',1)[1]
IndexError: list index out of range
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : True
        notification_center : True
                   run_mode : 2
              cluster_queue : None
               cluster_time : None
               cluster_size : 100
             cluster_memory : None
                    nb_core : 24 (user set)
          cluster_temp_path : None

                      Configuration Options
                      ---------------------
               pythia8_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : None (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/Delphes (user set)
        exrootanalysis_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/ExRootAnalysis (user set)
               syscalc_path : None (user set)
                     lhapdf : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.2.0 2021-08-22 *
#* *
#* 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 group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
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 /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/models\
/HEL_UFO
generate p p > e+ e-
output EFTML
######################################################################
## 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 CKMBLOCK
###################################
Block ckmblock
    1 2.277360e-01 # cabi

###################################
## INFORMATION FOR MASS
###################################
Block mass
    1 5.040000e-03 # MD
    2 2.550000e-03 # MU
    3 1.010000e-01 # MS
    4 1.270000e+00 # MC
    5 4.700000e+00 # MB
    6 1.720000e+02 # MT
   11 5.110000e-04 # Me
   13 1.056600e-01 # MMU
   15 1.777000e+00 # MTA
   24 8.038500e+01 # MW
   25 1.250000e+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.
  12 0.000000e+00 # ve : 0.0
  14 0.000000e+00 # vm : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  23 9.163403e+01 # z : cmath.sqrt((gw__exp__2*vev__exp__2*(1 - cT + (ee__exp__2*(cw__exp__2*cWW + cB*sw__exp__2 + 4*cA*sw__exp__4)*vev__exp__2)/(2.*cw__exp__2*MW__exp__2*sw__exp__2)))/cw__exp__2)/2.

###################################
## INFORMATION FOR NEWCOUP
###################################
Block newcoup
    1 0.000000e+00 # cH
    2 0.000000e+00 # cT
    3 0.000000e+00 # c6
    4 0.000000e+00 # cu
    5 0.000000e+00 # cd
    6 0.000000e+00 # cl
    7 0.000000e+00 # cWW
    8 0.000000e+00 # cB
    9 0.000000e+00 # cHW
   10 0.000000e+00 # cHB
   11 0.000000e+00 # cA
   12 0.000000e+00 # cG
   13 0.000000e+00 # cHQ
   14 0.000000e+00 # cpHQ
   15 0.000000e+00 # cHu
   16 0.000000e+00 # cHd
   17 0.000000e+00 # cHud
   18 0.000000e+00 # cHL
   19 0.000000e+00 # cpHL
   20 0.000000e+00 # cHe
   21 0.000000e+00 # cuB
   22 0.000000e+00 # cuW
   23 0.000000e+00 # cuG
   24 0.000000e+00 # cdB
   25 0.000000e+00 # cdW
   26 0.000000e+00 # cdG
   27 0.000000e+00 # clB
   28 0.000000e+00 # clW
   29 0.000000e+00 # c3W
   30 0.000000e+00 # c3G
   31 0.000000e+00 # c2W
   32 0.000000e+00 # c2B
   33 0.000000e+00 # c2G
   34 0.000000e+00 # tcHW
   35 0.000000e+00 # tcHB
   36 0.000000e+00 # tcG
   37 0.000000e+00 # tcA
   38 0.000000e+00 # tc3W
   39 0.000000e+00 # tc3G

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.279000e+02 # aEWM1
    2 1.166370e-05 # Gf
    3 1.184000e-01 # aS (Note that Parameter not used if you use a PDF set)

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    1 5.040000e-03 # ymdo
    2 2.550000e-03 # ymup
    3 1.010000e-01 # yms
    4 1.270000e+00 # ymc
    5 4.700000e+00 # ymb
    6 1.720000e+02 # ymt
   11 5.110000e-04 # yme
   13 1.056600e-01 # ymm
   15 1.777000e+00 # ymtau

###################################
## 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.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # 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 *
#*********************************************************************
#
#*********************************************************************
# 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 *
#*********************************************************************
  50000 = 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, 4=photon from muon *
#*********************************************************************
     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)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   2 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# 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) *
#*********************************************************************
 10.0 = ptl ! minimum pt for the charged leptons
 -1.0 = ptlmax ! maximum pt for the charged leptons
 {} = 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})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
 2.5 = etal ! max rap for the charged leptons
 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 = drll ! min distance between leptons
 -1.0 = drllmax ! max distance between leptons
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton 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 / neutrino *
 # for pair of lepton includes only same flavor, opposite charge
 #*********************************************************************
 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 = xptl ! minimum pt for at least one charged lepton
 #*********************************************************************
 # 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
 -1.0 = ptl1max ! maximum pt for the leading lepton in pt
 -1.0 = ptl2max ! maximum pt for the second lepton in pt
#*********************************************************************
# 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'

----------------------------------------------------------------------------------------------------------------------------------------------------

It seems to be complaining about a list index being out of range.

Thanks for the help!
Torrey

Hello,

Since then I have also ran the following:

generate p p > t t~
output
launch -i
multi_run 2

And I set the number of events to 50k. I then encountered the same list index error as before.
--------------------------------------------------------------------------------------------------------------------------------------------
#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.2.0 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 /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/PROC_sm_0/Events/run_01_1/unweighted_events.lhe.gz_3.lhe -from_cards --multicore=wait
Traceback (most recent call last):
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/PROC_sm_0/bin/internal/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/PROC_sm_0/bin/internal/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/PROC_sm_0/bin/internal/common_run_interface.py", line 2234, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1692, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/reweight_interface.py", line 515, in do_launch
    param_card_iterator, tag_name = self.handle_param_card(model_line, args, type_rwgt)
  File "/afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/madgraph/interface/reweight_interface.py", line 759, in handle_param_card
    tmp_card = new_card.lower().split('block',1)[1]
IndexError: list index out of range
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : True
        notification_center : True
                   run_mode : 2
              cluster_queue : None
               cluster_time : None
               cluster_size : 100
             cluster_memory : None
                    nb_core : 24 (user set)
          cluster_temp_path : None

                      Configuration Options
                      ---------------------
               pythia8_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : None (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/Delphes (user set)
        exrootanalysis_path : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/ExRootAnalysis (user set)
               syscalc_path : None (user set)
                     lhapdf : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /afs/crc.nd.edu/user/t/tsaxton/Home/MG5_aMC_v3_2_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.2.0 2021-08-22 *
#* *
#* 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 group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
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~
generate p p > t t~
output
######################################################################
## 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 MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+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.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.700000e+00 # ymb
    6 1.730000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-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.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # 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 *
#*********************************************************************
#
#*********************************************************************
# 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 *
#*********************************************************************
  50000 = 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, 4=photon from muon *
#*********************************************************************
     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)
  average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
                             ! 0: sum over helicity, 1: importance sampling
   1 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx)
                             ! 1 is old strategy (using amp square)
        ! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#*********************************************************************
# 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})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# 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 *
#*********************************************************************
   True = 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'
-------------------------------------------------------------------------------------------------------------------------------------------

Thanks for the help!

Hi,

your first error was a filenotfound error,
so looks like another issue now.

If I'm doing this:
generate p p > t t~
output
launch -i
multi_run 2
set nevents 50k

I do not have any reweighting setup and therefore can potentially reproduce the first error of this thread ( I do not)
but certainly not the second since the crash is within the reweighting module that is not run with the above command.

So I did try the following:

generate p p > t t~
output
launch -i
multi_run 2
reweight=ON
set nevents 50k

and this one also works, but this might also depend of your reweighting_card which is actually not include in the debug file.
the above script is also fully working on my machine. (tested with latest LTS version, and 3.2.0)

Hello,

Below is how my reweight card looks for the HEL_UFO model I am using. I've used this card in 3.1.1, not doing multi_run and had no issues with it.

--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------#
#*************************************************************************
# Reweight Module *
# Matrix-Element reweighting at LO/NLO *
# Mattelaer Olivier arxiv:1607.00763 *
#*************************************************************************
#
# Note:
# 1) the value of alpha_s will be used from the event so the value in
# the param_card is not taken into account.
# 2) It is (in general) dangerous/wrong to change parameters by a large
# amount, if this changes the shape of the matrix elements a lot.
# (For example, changing a particle's mass by much more than its
# width leads to very inaccurate result). In such a case, separate
# event generation runs are needed.
#
#************************************************************************
# ENTER YOUR COMMANDS BELOW.
#************************************************************************

change mode LO

change model HEL_UFO

launch --rwgt_name=case1
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.0
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.0
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

launch --rwgt_name=case2
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.07184950398757557
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.02992736187170405
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

launch --rwgt_name=case3
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.17207662924114733
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.27906538067330533
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

launch --rwgt_name=case4
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.05876149135539773
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.14623164074985998
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

launch --rwgt_name=case5
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.07807894064778671
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.2934652981823629
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

launch --rwgt_name=case6
   set NEWCOUP 1 0.
   set NEWCOUP 2 0.
   set NEWCOUP 3 0.
   set NEWCOUP 4 0.
   set NEWCOUP 5 0.
   set NEWCOUP 6 0.
   set NEWCOUP 7 0.2211117415535981
   set NEWCOUP 8 0.
   set NEWCOUP 9 0.
   set NEWCOUP 10 0.
   set NEWCOUP 11 0.06019191704486092
   set NEWCOUP 12 0.
   set NEWCOUP 13 0.
   set NEWCOUP 14 0.
   set NEWCOUP 15 0.
   set NEWCOUP 16 0.
   set NEWCOUP 17 0.
   set NEWCOUP 18 0.
   set NEWCOUP 19 0.
   set NEWCOUP 20 0.
   set NEWCOUP 21 0.
   set NEWCOUP 22 0.
   set NEWCOUP 23 0.
   set NEWCOUP 24 0.
   set NEWCOUP 25 0.
   set NEWCOUP 26 0.
   set NEWCOUP 27 0.
   set NEWCOUP 28 0.
   set NEWCOUP 29 0.
   set NEWCOUP 30 0.
   set NEWCOUP 31 0.
   set NEWCOUP 32 0.
   set NEWCOUP 33 0.
   set NEWCOUP 34 0.
   set NEWCOUP 35 0.
   set NEWCOUP 36 0.
   set NEWCOUP 37 0.
   set NEWCOUP 38 0.
   set NEWCOUP 39 0.

#
# SPECIFY A PATH OR USE THE SET COMMAND LIKE THIS:
# set sminputs 1 130
# modify 1/alpha_EW
#************************************************************************
# Manual: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Reweight
#************************************************************************
#
# Example of (standard) code for the computation of two weights:
#
# launch ! tag to start the computation of the first weight
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# launch ! start to compute a second weight
# /home/Cards/param_card_2.dat ! you can also enter a path to a valid card
#
# Note: The command to specify the parameter are AFTER the associated "launch"
#
# Possible options:
# You can enter one of the following lines to customize the reweighting
# procedure. These need to be given before the 'launch' command.
#
# change model NAME : use another model for the matrix-elements to reweight
# with. In this case you need to provide the path to a correct
# param_card for the new model; you cannot modify the original one
# with the 'set' command.
# change process DEF [--add]: change the process by which you reweight.
# The initial and final state particles of the new process should
# be exactly identical to the ones in the original process.
# change helicity False: perform the reweighting by helicity summed
# matrix-elements even if the events have been written with a
# single helicity state.
# change mode XXX: change the type of reweighting performed.
# allowed values: LO, NLO, LO+NLO
# - This command has no effect for reweighting an .lhe event file with LO accuracy.
# In that case LO mode is always used (whatever entry is set).
# - When the .lhe file reweighted is at NLO accuracy, then all modes are allowed.
# * "LO" is an approximate leading order method
# * "NLO" is the NLO accurate method
# * "LO+NLO" runs both
# - "NLO" and "LO+NLO" modes requires 'store_rwgt_info' equals True (run_card.dat)
# If the reweighting is done at generation level this parameter will
# automatically be set on True.
#************************************************************************
---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------

For reasons I cannot fathom, when I do

launch EFTML -i
multi_run 10
set nevents 100k

it makes it to the end, and fails to combine the files.

However, if I do

launch EFTML -i
multi_run 2
set nevents 50k

it fails with the index out of range error.

Thanks for the help!

Still can not reproduce your issue.
Should be related to hardware specific issue I guess

Olivier