IndexError : list index out of range

Asked by Seungyoon Kim on 2018-06-11

I have installed boost, pythia8, and MadAnalysis5, and now I am trying to run the tutorial.

When I launch my process, I first get a warning when running Systematics Computation saying that:
         /PATH/TO/lhapdf6/bin/lhapdf-config launch ends with non zero status: 1.

Once the program reaches the end of the systematics computations, I get the following error:
        Command "generate_events run_01" interrupted with error:
        IndexError : list index out of range

What can I do to resolve these warnings/errors? I will attach the run_01_log_1_debug.log file below.

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 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 MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1501, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2467, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 6685, in new_fct
    original_fct(obj, *args, **opts)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2549, in run_generate_events
    postcmd=False, printcmd=False)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 1868, in do_systematics
    sys_obj.print_cross_sections(all_cross, nb_event, result_file)
  File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/various/systematics.py", line 417, in print_cross_sections
    to_print.append(all_cross[i])
IndexError: list index out of range
                              Run Options
                              -----------
               stdout_level : 20 (user set)

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/pythia8 (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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/lhapdf6/bin/lhapdf-config (user set)
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : None (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~
generate p p > t t~
output seungyoon
######################################################################
## 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.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
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 80.419002 # 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

###################################
## 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.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 *
#*********************************************************************
   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'

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2018-06-14
Last query:
2018-06-14
Last reply:
2018-06-13

Thanks, I will check what this is.
(maybe already solved in the development version)

In the mean time you can set in your run_card the parameter "use_syst" to False
to bypass the computation of systemtatic uncertainty which trigger that issue.

Cheers,

Olivier

> On 11 Jun 2018, at 19:22, Seungyoon Kim <email address hidden> wrote:
>
> New question #670139 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/670139
>
> I have installed boost, pythia8, and MadAnalysis5, and now I am trying to run the tutorial.
>
> When I launch my process, I first get a warning when running Systematics Computation saying that:
> /PATH/TO/lhapdf6/bin/lhapdf-config launch ends with non zero status: 1.
>
> Once the program reaches the end of the systematics computations, I get the following error:
> Command "generate_events run_01" interrupted with error:
> IndexError : list index out of range
>
> What can I do to resolve these warnings/errors? I will attach the run_01_log_1_debug.log file below.
>
>
>
>
>
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 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 MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #************************************************************
> generate_events run_01
> Traceback (most recent call last):
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1501, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2467, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 6685, in new_fct
> original_fct(obj, *args, **opts)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2549, in run_generate_events
> postcmd=False, printcmd=False)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 1868, in do_systematics
> sys_obj.print_cross_sections(all_cross, nb_event, result_file)
> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/various/systematics.py", line 417, in print_cross_sections
> to_print.append(all_cross[i])
> IndexError: list index out of range
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None (user set)
> cluster_size : 100
> cluster_queue : None (user set)
> nb_core : 40 (user set)
> cluster_time : 40 (user set)
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> pythia8_path : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/pythia8 (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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/lhapdf6/bin/lhapdf-config (user set)
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : None (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~
> generate p p > t t~
> output seungyoon
> ######################################################################
> ## 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.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
> 16 0.000000 # vt : 0.0
> 21 0.000000 # g : 0.0
> 22 0.000000 # a : 0.0
> 24 80.419002 # 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
>
> ###################################
> ## 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.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 *
> #*********************************************************************
> 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'
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

I have taken a deeper look,
and this seems to be an issue with lhapdf6 that was installed for pythia8.
What is printed to the screen when you do
 /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/lhapdf6/bin/lhapdf-config (user set)

Also what is your os?

Thanks,

Olivier

> On 12 Jun 2018, at 23:43, Olivier Mattelaer <email address hidden> wrote:
>
> Question #670139 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/670139
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Thanks, I will check what this is.
> (maybe already solved in the development version)
>
> In the mean time you can set in your run_card the parameter "use_syst" to False
> to bypass the computation of systemtatic uncertainty which trigger that issue.
>
> Cheers,
>
> Olivier
>
>> On 11 Jun 2018, at 19:22, Seungyoon Kim <email address hidden> wrote:
>>
>> New question #670139 on MadGraph5_aMC@NLO:
>> https://answers.launchpad.net/mg5amcnlo/+question/670139
>>
>> I have installed boost, pythia8, and MadAnalysis5, and now I am trying to run the tutorial.
>>
>> When I launch my process, I first get a warning when running Systematics Computation saying that:
>> /PATH/TO/lhapdf6/bin/lhapdf-config launch ends with non zero status: 1.
>>
>> Once the program reaches the end of the systematics computations, I get the following error:
>> Command "generate_events run_01" interrupted with error:
>> IndexError : list index out of range
>>
>> What can I do to resolve these warnings/errors? I will attach the run_01_log_1_debug.log file below.
>>
>>
>>
>>
>>
>>
>> #************************************************************
>> #* MadGraph5_aMC@NLO/MadEvent *
>> #* *
>> #* * * *
>> #* * * * * *
>> #* * * * * 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 MadEvent *
>> #* *
>> #* run as ./bin/madevent.py filename *
>> #* *
>> #************************************************************
>> generate_events run_01
>> Traceback (most recent call last):
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1501, in onecmd
>> return self.onecmd_orig(line, **opt)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
>> return func(arg, **opt)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2467, in do_generate_events
>> self.run_generate_events(switch_mode, args)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 6685, in new_fct
>> original_fct(obj, *args, **opts)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/madevent_interface.py", line 2549, in run_generate_events
>> postcmd=False, printcmd=False)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
>> stop = Cmd.onecmd_orig(current_interface, line, **opt)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
>> return func(arg, **opt)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/interface/common_run_interface.py", line 1868, in do_systematics
>> sys_obj.print_cross_sections(all_cross, nb_event, result_file)
>> File "/cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/madgraph/various/systematics.py", line 417, in print_cross_sections
>> to_print.append(all_cross[i])
>> IndexError: list index out of range
>> Run Options
>> -----------
>> stdout_level : 20 (user set)
>>
>> MadEvent Options
>> ----------------
>> automatic_html_opening : False (user set)
>> notification_center : True
>> cluster_temp_path : None
>> cluster_memory : None (user set)
>> cluster_size : 100
>> cluster_queue : None (user set)
>> nb_core : 40 (user set)
>> cluster_time : 40 (user set)
>> run_mode : 2
>>
>> Configuration Options
>> ---------------------
>> text_editor : None
>> cluster_local_path : None
>> cluster_status_update : (600, 30)
>> pythia8_path : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/pythia8 (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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/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 : /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/lhapdf6/bin/lhapdf-config (user set)
>> f2py_compiler : None
>> hepmc_path : None (user set)
>> cluster_retry_wait : 300
>> fortran_compiler : None
>> auto_update : 7 (user set)
>> exrootanalysis_path : None (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~
>> generate p p > t t~
>> output seungyoon
>> ######################################################################
>> ## 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.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
>> 16 0.000000 # vt : 0.0
>> 21 0.000000 # g : 0.0
>> 22 0.000000 # a : 0.0
>> 24 80.419002 # 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
>>
>> ###################################
>> ## 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.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 *
>> #*********************************************************************
>> 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'
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Seungyoon Kim (senyen) said : #3

Hi Olivier,

My OS is CentOS 6.9.

I don't know if the user set in parenthesis was significant or not, but this is what printed when I typed /cluster/kappa/90-days-archive/beaucheminlab/skim38/mg5amc/HEPTools/lhapdf6/bin/lhapdf-config:

lhapdf-config: configuration tool for the LHAPDF
               parton density function evolution library
               http://projects.hepforge.org/lhapdf/

Usage: lhapdf-config [options]
Options:
  --help | -h : show this help message
  --prefix : show the installation prefix (cf. autoconf)
  --incdir : show the path to the LHAPDF C++ header directory
  --libdir : show the path to the LHAPDF library directory
  --datadir : show the path to the LHAPDF data directory

  --cflags : get compiler flags (aka --cppflags|--cxxflags)
  --libs : get linker flags (aka --ldflags)

  --version : return LHAPDF release version number

Thanks,
Seungyoon

Seungyoon Kim (senyen) said : #4

Also, just out of curiosity, how do I output a hepmc file instead of a lhe file or MA5 report analysis?

Hi,

You can convert lhe file to a hepmc file via some converter that you can find online.
Otherwise if you shower your events with Pythia8 the default output format is then hepmc.

Cheers,

Olivier

> On 13 Jun 2018, at 19:12, Seungyoon Kim <email address hidden> wrote:
>
> Question #670139 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/670139
>
> Seungyoon Kim posted a new comment:
> Also, just out of curiosity, how do I output a hepmc file instead of a
> lhe file or MA5 report analysis?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Seungyoon Kim (senyen) said : #6

Thanks Olivier Mattelaer, that solved my question.