End code 2.0

Asked by Giacomo on 2020-10-19

I'm running the MadDump plugin and get the following error from MadGraph. I'd appreciate any advice you have on what to. Here is the contents of the error log file. Thank you.

generate_events run_01 -f
Traceback (most recent call last):
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/extended_cmd.py", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/madevent_interface.py", line 2468, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/common_run_interface.py", line 7011, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/madevent_interface.py", line 2511, in run_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/extended_cmd.py", line 1544, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/madevent_interface.py", line 3484, in do_refine
    x_improve.launch() # create the ajob for the refinment.
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/gen_ximprove.py", line 864, in launch
    main_dir=pjoin(self.cmd.me_dir,'SubProcesses')) #main_dir is for gridpack readonly mode
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/sum_html.py", line 747, in collect_result
    P_comb.add_results(os.path.basename(G), path, mfactors[G])
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/sum_html.py", line 425, in add_results
    oneresult.read_results(filepath)
  File "/home/maroccol/Code/MG5_aMC_v2_7_0/output/BEBC/mQ_phi/interaction_electron/bin/internal/sum_html.py", line 315, in read_results
    % (error_code, open(log).read())
Exception: Reported error: End code 2.0
 Full associated log:
 Process in group number 1200
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23_lo_as_0130_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23lo1: 0.13000000000000000
 Define smin to 1.2327121000000002E-005
 *****************************************************
 * MadGraph/MadEvent *
 * -------------------------------- *
 * http://madgraph.hep.uiuc.edu *
 * http://madgraph.phys.ucl.ac.be *
 * http://madgraph.roma2.infn.it *
 * -------------------------------- *
 * *
 * PARAMETER AND COUPLING VALUES *
 * *
 *****************************************************

  External Params
  ---------------------------------

 mdl_cabi = 0.22773599999999999
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.11839999999999999
 mdl_ymdo = 5.0400000000000002E-003
 mdl_ymup = 2.5500000000000002E-003
 mdl_yms = 0.10100000000000001
 mdl_ymc = 1.2700000000000000
 mdl_ymb = 4.7000000000000002
 mdl_ymt = 172.00000000000000
 mdl_yme = 5.1099999999999995E-004
 mdl_ymm = 0.10566000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_muBI = 1.0000000000000000
 mdl_muChi = 0.0000000000000000
 mdl_dChi = 0.0000000000000000
 mdl_eps = 1.0000000000000000
 mdl_MZ = 91.187600000000003
 mdl_Me = 5.1099999999999995E-004
 mdl_MMU = 0.10566000000000000
 mdl_MTA = 1.7769999999999999
 mdl_MU = 2.5500000000000002E-003
 mdl_MC = 1.2700000000000000
 mdl_MT = 172.00000000000000
 mdl_MD = 5.0400000000000002E-003
 mdl_MS = 0.10100000000000001
 mdl_MB = 4.7000000000000002
 mdl_MH = 125.00000000000000
 mdl_mChi = 3.0000000000000001E-003
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WT = 1.5083359999999999
 mdl_WH = 4.0699999999999998E-003
 mdl_Mpi0 = 0.13497999999999999
 mdl_Wpi0 = 7.8000000000000004E-009
 mdl_Meta = 0.54800000000000004
 mdl_Weta = 1.3100000000000000E-006
 mdl_Momega = 0.78300000000000003
 mdl_Womega = 8.5000000000000006E-003
 mdl_Mrho = 0.77500000000000002
 mdl_Wrho = 0.14699999999999999
 mdl_Mphi = 1.0189999999999999
 mdl_Wphi = 4.1999999999999997E-003
 mdl_Mjpsi = 3.0970000000000000
 mdl_Wjpsi = 9.2999999999999999E-002
  Internal Params
  ---------------------------------

 mdl_cos__cabi = 0.97418004031982097
 mdl_CKM1x1 = (0.97418004031982097,0.0000000000000000)
 mdl_sin__cabi = 0.22577256042856930
 mdl_CKM1x2 = (0.22577256042856930,0.0000000000000000)
 mdl_CKM1x3 = (0.0000000000000000,0.0000000000000000)
 mdl_CKM2x1 = (-0.22577256042856930,0.0000000000000000)
 mdl_CKM2x2 = (0.97418004031982097,0.0000000000000000)
 mdl_CKM2x3 = (0.0000000000000000,0.0000000000000000)
 mdl_CKM3x1 = (0.0000000000000000,0.0000000000000000)
 mdl_CKM3x2 = (0.0000000000000000,0.0000000000000000)
 mdl_CKM3x3 = (1.0000000000000000,0.0000000000000000)
 mdl_MZ__exp__2 = 8315.1783937600012
 mdl_MZ__exp__4 = 69142191.720053151
 mdl_sqrt__2 = 1.4142135623730951
 mdl_MH__exp__2 = 15625.000000000000
 mdl_complexi = (0.0000000000000000,1.0000000000000000)
 mdl_conjg__CKM1x1 = (0.97418004031982097,-0.0000000000000000)
 mdl_conjg__CKM1x2 = (0.22577256042856930,-0.0000000000000000)
 mdl_conjg__CKM1x3 = (0.0000000000000000,-0.0000000000000000)
 mdl_conjg__CKM2x1 = (-0.22577256042856930,-0.0000000000000000)
At line 68 of file splev.f
Fortran runtime error: Index '4' of dimension 1 of array 't' above upper bound of 0

Error termination. Backtrace:
 mdl_conjg__CKM2x2 = (0.97418004031982097,-0.0000000000000000)
 mdl_conjg__CKM2x3 = (0.0000000000000000,-0.0000000000000000)
 mdl_conjg__CKM3x1 = (0.0000000000000000,-0.0000000000000000)
 mdl_conjg__CKM3x2 = (0.0000000000000000,-0.0000000000000000)
 mdl_conjg__CKM3x3 = (1.0000000000000000,-0.0000000000000000)
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824359746197842
 mdl_sqrt__aEW = 8.8422894590285753E-002
 mdl_ee = 0.31345100004952897
 mdl_MW__exp__2 = 6371.9284088904105
 mdl_sw2 = 0.23369913342182447
 mdl_cw = 0.87538612427783857
 mdl_sqrt__sw2 = 0.48342438232036300
 mdl_sw = 0.48342438232036300
 mdl_g1 = 0.35807170271074895
 mdl_gw = 0.64839716719502682
 mdl_vev = 246.22056907348590
 mdl_vev__exp__2 = 60624.568634871241
 mdl_lam = 0.12886689630821144
 mdl_yb = 2.6995322804122722E-002
 mdl_yc = 7.2944808428161397E-003
 mdl_ydo = 2.8948175943144366E-005
 mdl_ye = 2.9350233942354699E-006
 mdl_ym = 6.0687783137949077E-004
 mdl_ys = 5.8011225600348830E-004
 mdl_yt = 0.98791394091683138
 mdl_ytau = 1.0206529494239589E-002
 mdl_yup = 1.4646398542662327E-005
 mdl_muH = 88.388347648318444
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_sw__exp__2 = 0.23369913342182450
 mdl_cw__exp__2 = 0.76630086657817542
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34409301068170506
 mdl_G__exp__2 = 1.4878582807401259
  Couplings of milliDipoles
  ---------------------------------

         GC_3 = -0.00000E+00 -0.31345E+00
         GC_6 = 0.00000E+00 0.31345E+00
        GC_11 = 0.00000E+00 0.00000E+00
        GC_12 = 0.00000E+00 0.00000E+00

 Collider parameters:
 --------------------

 Running at aUe machine @ 0.62330872802055315 GeV
 PDF set = nn23lo1
 alpha_s(Mz)= 0.1300 running at 2 loops.
 alpha_s(Mz)= 0.1300 running at 2 loops.
 Renormalization scale set on event-by-event basis
 Factorization scale set on event-by-event basis

 getting user params
Enter number of events and max and min iterations:
 Number of events and iterations 1000 5 3
Enter desired fractional accuracy:
 Desired fractional accuracy: 0.10000000000000001
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
 Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
 Explicitly summing over helicities
Enter Configuration Number:
Running Configuration Number: 1
 Not subdividing B.W.
 Attempting mappinvarients 1 4
 Completed mapping 4
 about to integrate 3 1000 5 3 3 1
 Using non-zero grid deformation.
  3 dimensions 1000 events 3 invarients 5 iterations 1 config(s), (0.99)
 Using h-tuple random number sequence.
 Error opening grid
 Using Uniform Grid! 16
 Using uniform alpha 1.0000000000000000
 Grid defined OK
 Set CM energy to 0.44
 Mapping Graph 1 to config 1
Setting grid 1 0.13441E-05 1
 Warning: No minimum found for integration
          Setting minimum to 5.1475772741132415
 Transforming s_hat 1/s 3 5.1475772741132415 1.2327121000000002E-005 0.19426614633429998
   1 1 2 3
 Masses: 0.300E-02 0.511E-03 0.300E-02 0.511E-03
Using random seed offsets 1 : 1
  with seed 24
 Ranmar initialization seeds 31059 9398
#0 0x7ff69d82c32a
#1 0x7ff69d82ced5
#2 0x7ff69d82d2a7
#3 0x556159e12b65
#4 0x556159dc5a8a
#5 0x556159dc5d73
#6 0x556159d9f0cb
#7 0x556159d9f9e7
#8 0x556159ddc52b
#9 0x556159d98bdc
#10 0x556159d98d4f
#11 0x7ff69ce8bb96
#12 0x556159d96a69
#13 0xffffffffffffffff

ls status:
input_app.txt
param.log
run1_app.log

                              Run Options
                              -----------
               stdout_level : None

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : /home/maroccol/Code/MG5_aMC_v2_7_0/pythia-pgs (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : lhapdf-config
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : /home/maroccol/Code/MG5_aMC_v2_7_0/ExRootAnalysis (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.7.3 2020-06-21 *
#* *
#* 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 milliDipoles
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_events decay /home/maroccol/Code/pythia8244/examples/phi400GeV.\
hepmc
decay phi > chi chi~
define darkmatter chi
add process interaction @electron
output output/BEBC/mQ_phi
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR CKMBLOCK
###################################
BLOCK CKMBLOCK #
      1 2.277360e-01 # cabi
###################################
## INFORMATION FOR FRBLOCK
###################################
BLOCK FRBLOCK #
      1 1.000000e+00 # mubi
      2 0.000000e+00 # muchi
      3 0.000000e+00 # dchi
      4 1.000000e+00 # eps
###################################
## 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
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh
      111 1.349800e-01 # mpi0
      113 7.750000e-01 # mrho
      221 5.480000e-01 # meta
      223 7.830000e-01 # momega
      333 1.019000e+00 # mphi
      443 3.097000e+00 # mjpsi
      9000005 0.003 # mchi
      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
      24 7.982436e+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.279000e+02 # aewm1
      2 1.166370e-05 # gf
      3 1.184000e-01 # as
###################################
## 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
DECAY 111 7.800000e-09 # wpi0
DECAY 113 1.470000e-01 # wrho
DECAY 221 1.310000e-06 # weta
DECAY 223 8.500000e-03 # womega
DECAY 333 4.200000e-03 # wphi
DECAY 443 9.300000e-02 # wjpsi
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
DECAY 9000005 0.000000e+00 # chi : 0.0
###################################
## INFORMATION FOR QNUMBERS 9000005
###################################
BLOCK QNUMBERS 9000005 # chi
      1 0 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 111
###################################
BLOCK QNUMBERS 111 # pi0
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 221
###################################
BLOCK QNUMBERS 221 # eta
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 223
###################################
BLOCK QNUMBERS 223 # omega
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 113
###################################
BLOCK QNUMBERS 113 # rho
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 333
###################################
BLOCK QNUMBERS 333 # phi
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 443
###################################
BLOCK QNUMBERS 443 # jpsi
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)

#*********************************************************************
# 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 *
#*********************************************************************
  312 = 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 *
#*********************************************************************
     9 = lpp1 ! beam 1 type
     0 = lpp2 ! beam 2 type
     190.07523015 = ebeam1 ! beam 1 total energy in GeV
     0.000511 = 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"#*********************************************************************
#
#*********************************************************************
# 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)
#IF(d)# #*********************************************************************
#IF(d)# # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
#IF(d)# # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#IF(d)# #*********************************************************************
   False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 0.0 = ptj ! minimum pt for the jets
 0.0 = ptb ! minimum pt for the b
 0.0 = pta ! minimum pt for the photons
 0.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})
#
# 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) *
#*********************************************************************
 -1.0 = etaj ! max rap for the jets
 -1.0 = etab ! max rap for the b
 -1.0 = etaa ! max rap for the photons
 -1.0 = 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.0 = drjj ! min distance between jets
 0.0 = drbb ! min distance between b's
 0.0 = drll ! min distance between leptons
 0.0 = draa ! min distance between gammas
 0.0 = drbj ! min distance between b and jet
 0.0 = draj ! min distance between gamma and jet
 0.0 = 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.0 = 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 *
#IF(llll)# # WARNING: for four lepton final state mmll cut require to have *
#IF(llll)# # 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.
#IF(LL)# #*********************************************************************
#IF(LL)# # Minimum and maximum invariant mass for all letpons *
#IF(LL)# #*********************************************************************
 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#IF(LL)# #*********************************************************************
#IF(LL)# # Minimum and maximum pt for 4-momenta sum of leptons / neutrino *
#IF(LL)# # for pair of lepton includes only same flavor, opposite charge
#IF(LL)# #*********************************************************************
 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
#IF(jj)# #*********************************************************************
#IF(jj)# # Control the pt's of the jets sorted by pt *
#IF(jj)# #*********************************************************************
 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
#IF(ll)# #*********************************************************************
#IF(ll)# # Control the pt's of leptons sorted by pt *
#IF(ll)# #*********************************************************************
 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
#IF(JJ)# #*********************************************************************
#IF(JJ)# # Control the Ht(k)=Sum of k leading jets *
#IF(JJ)# #*********************************************************************
 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
#IF(aj)# #***********************************************************************
#IF(aj)# # Photon-isolation cuts, according to hep-ph/9801442 *
#IF(aj)# # When ptgmin=0, all the other parameters are ignored *
#IF(aj)# # When ptgmin>0, pta and draj are not going to be used *
#IF(aj)# #***********************************************************************
#IF(aj)# %(ptgmin)s = ptgmin ! Min photon transverse momentum
#IF(aj)# %(r0gamma)s = R0gamma ! Radius of isolation code
#IF(aj)# %(xn)s = xn ! n parameter of eq.(3.4) in hep-ph/9801442
#IF(aj)# %(epsgamma)s = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
#IF(aj)# %(isoem)s = isoEM ! isolate photons from EM energy (photons and leptons)
#IF(jj)# #*********************************************************************
#IF(jj)# # WBF cuts *
#IF(jj)# #*********************************************************************
 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
 -1.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   False = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'#*********************************************************************
# Additional hidden parameters
#*********************************************************************
  1.0 = xn # hidden_parameter
  True = isoem # hidden_parameter
  1.0 = epsgamma # hidden_parameter
  0.4 = r0gamma # hidden_parameter
  0.0 = ptgmin # hidden_parameter

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Luca Edit question
Last query:
2020-10-21
Last reply:
2020-10-25

This question was reopened

The important part of the debug file is

>At line 68 of file splev.f
Fortran runtime error: Index '4' of dimension 1 of array 't' above upper bound of 0

Luca do you know if a more recent version of the file exists? or if this is an issue in the way you call that function?

Cheers,

Olivier

Giacomo (gmarocco) said : #2

Just to give slightly more information, everything runs fine when I have a larger number of decaying mesons in the input .hepmc file, with a corresponding larger n_passed. This could just be a coincidence, but this is the only difference I can find between the case when it runs vs when it doesn't.

Also thanks Olivier for the clarification.

Then it might be that it needs a minimum amount of event in order to be able to do the extrapolation use in maddump (this makes sense actually)

Cheers,

Olivier

> On 21 Oct 2020, at 11:50, Giacomo <email address hidden> wrote:
>
> Question #693528 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/693528
>
> Giacomo gave more information on the question:
> Just to give slightly more information, everything runs fine when I have
> a larger number of decaying mesons in the input .hepmc file, with a
> corresponding larger n_passed. This could just be a coincidence, but
> this is the only difference I can find between the case when it runs vs
> when it doesn't.
>
> Also thanks Olivier for the clarification.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Giacomo (gmarocco) said : #4

Yes, it might be that, and indeed MadDump will not let you simulate the interaction if below a threshold. This seems to be above MadDump's hardcoded threshold (which is perhaps too still low for the subsequent processes)

Luca (lbuono) said : #5

Dear Giacomo and Olivier,

sorry, I was not available the last few days.
I think you already understood the main problem,
but let me have some time to look more carefully at the issue.

I'll let you know tomorrow.

Cheers,
Luca

Luca (lbuono) said : #6

Dear Giacomo,

so far, I have been not available to reproduce your issue.
If you observe that the error does not occur increasing the
number of mesons in input, I would suggest to stay with
a high statistics.
The minimum value required by maddump
should be interpreted as a very poor statistics run that
should go through the end, but this it it not guaranteed,
and the result might be affected by a large uncertainty.
For some cases, I experienced stability problems when
generating the interaction events within the detector
because of a very poor dark matter energy fit.

In case, you can send me your model file and the input mesons file,
so that I can reproduce your error and investigate further.

Cheers,
Luca

Can you help with this problem?

Provide an answer of your own, or ask Giacomo for more information if necessary.

To post a message you must log in.