Too many files, resources unavailable in parallel

Our university performance computing center just changed systems, so I lead with that this may be an error on their part and not with Madgraph. However, running version 2.6.0, I have been having problems doing long multi_runs. Doing a manual parameter scan, I first do
multi_run 50
then set the parameters and launch. When that one finishes (without exiting Madevent) I then start the next with another multi_run 50. It seems to get through about 3-4 parameter points and then fails with the error OSError: [Errno 11] Resource temporarily unavailable. When I contacted the computing center, they said that I am hitting the limit for the number of files being opened. Would it be better if I exited in between parameters? Is there something about doing a multi_run that keeps many files open?

Below is one of the debug logs.

Thanks for your help,
Bryan

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
multi_run 50
multi_run 50
multi_run 50
Traceback (most recent call last):
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1136, in do_import
    self.import_command_file(args[1])
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1580, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/madevent_interface.py", line 2555, in do_multi_run
    self.exec_cmd('generate_events %s_%s -f' % (main_name, i), postcmd=False)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/madevent_interface.py", line 2187, in do_generate_events
    self.exec_cmd('shower --no_default', postcmd=False, printcmd=False)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/madevent_interface.py", line 3502, in do_shower
    postcmd=False, printcmd=False)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/madevent_interface.py", line 4166, in do_pythia8
    self.cluster.wait(parallelization_dir,wait_monitoring)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/cluster.py", line 636, in worker
    proc = misc.Popen([exe] + arg, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/misc.py", line 788, in deco_f
    return f(arg, *args, **opt)
  File "/gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/StopsPythia8/bin/internal/misc.py", line 827, in Popen
    return subprocess.Popen(arg, *args, **opt)
  File "/packages/python/2.7.13/lib/python2.7/subprocess.py", line 390, in __init__
    errread, errwrite)
  File "/packages/python/2.7.13/lib/python2.7/subprocess.py", line 916, in _execute_child
    self.pid = os.fork()
OSError: [Errno 11] Resource temporarily unavailable
Related File: None
                              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 : 28 (user set)
               cluster_time : None
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /gpfs/projects/het/bostdiek/Tools/MG5_aMC_v2_6_0/Delphes (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 : None (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 2017-08-16 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model MSSM_SLHA2
generate p p > t1 t1~
output StopsPythia8
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DCINFO
###################################
BLOCK DCINFO # decay program information
      1 sdecay/hdecay # decay calculator
      2 1.5 /3.4 # version number
###################################
## INFORMATION FOR SPINFO
###################################
BLOCK SPINFO # spectrum calculator information
      1 suspect # rge +spectrum calculator
      2 2.430000e+00 # version number
###################################
## INFORMATION FOR MODSEL
###################################
BLOCK MODSEL # model selection
      1 0 # #general mssm low scale
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS # standard model inputs
      1 1.279340e+02 # alpha_em^-1(m_z)^msbar
      2 1.166390e-05 # g_f [gev^-2]
      3 1.172000e-01 # alpha_s(m_z)^msbar
      4 9.118700e+01 # m_z pole mass
      5 4.250000e+00 # mb(mb)^msbar
      6 1.731000e+02 # mt pole mass
      7 1.777000e+00 # mtau pole mass
###################################
## INFORMATION FOR MINPAR
###################################
BLOCK MINPAR # input parameters - minimal models

###################################
## INFORMATION FOR EXTPAR
###################################
BLOCK EXTPAR # input parameters - non-minimal models
      0 1.690729e+03 # ewsb
      1 1.000000e+02 # m_1
      2 4.500000e+03 # m_2
      3 3.000000e+03 # m_3
      11 0.000000e+00 # a_t
      12 0.000000e+00 # a_b
      13 0.000000e+00 # a_tau
      14 0.000000e+00 # a_u
      15 0.000000e+00 # a_d
      16 0.000000e+00 # a_e
      23 4.500000e+03 # mu(ewsb)
      25 1.000000e+01 # tanbeta(in)
      26 5.000000e+03 # ma_pole
      31 4.000000e+03 # m_el
      32 4.000000e+03 # m_mul
      33 4.000000e+03 # m_taul
      34 4.000000e+03 # m_er
      35 4.000000e+03 # m_mur
      36 4.000000e+03 # m_taur
      41 5.000000e+03 # m_q1l
      42 5.000000e+03 # m_q2l
      43 5.500000e+03 # m_q3l
      44 5.000000e+03 # m_ur
      45 5.000000e+03 # m_cr
      46 5.000000e+02 # m_tr
      47 5.000000e+03 # m_dr
      48 5.000000e+03 # m_sr
      49 5.000000e+03 # m_br
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS # mass spectrum
      24 8.052218e+01 # w+
      25 1.183660e+02 # h
      35 5.000043e+03 # h
      36 5.000000e+03 # a
      37 5.000501e+03 # h+
      5 4.878778e+00 # b-quark pole mass calculated from mb(mb)_msbar
      1000001 5.020456e+03 # ~d_l
      2000001 5.020197e+03 # ~d_r
      1000002 5.019872e+03 # ~u_l
      2000002 5.020003e+03 # ~u_r
      1000003 5.020456e+03 # ~s_l
      2000003 5.020197e+03 # ~s_r
      1000004 5.019872e+03 # ~c_l
      2000004 5.020003e+03 # ~c_r
      1000005 5.020005e+03 # ~b_1
      2000005 5.520424e+03 # ~b_2
      1000006 1.700000e+02 # ~t_1
      2000006 5.507735e+03 # ~t_2
      1000011 4.000244e+03 # ~e_l
      2000011 4.000243e+03 # ~e_r
      1000012 3.999513e+03 # ~nu_el
      1000013 4.000244e+03 # ~mu_l
      2000013 4.000243e+03 # ~mu_r
      1000014 3.999513e+03 # ~nu_mul
      1000015 3.990642e+03 # ~tau_1
      2000015 4.009823e+03 # ~tau_2
      1000016 3.999513e+03 # ~nu_taul
      1000021 3.270707e+03 # ~g
      1000022 1.000000e+00 # ~chi_10
      1000023 4.473892e+03 # ~chi_20
      1000025 -4.519247e+03 # ~chi_30
      1000035 4.598885e+03 # ~chi_40
      1000024 4.473885e+03 # ~chi_1+
      1000037 4.598883e+03 # ~chi_2+
      6 1.731000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118760e+01 # mz
      1 0.000000e+00 # d : 0.0
      2 0.000000e+00 # u : 0.0
      3 0.000000e+00 # s : 0.0
      4 0.000000e+00 # c : 0.0
      11 0.000000e+00 # e- : 0.0
      12 0.000000e+00 # ve : 0.0
      13 0.000000e+00 # mu- : 0.0
      14 0.000000e+00 # vm : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
###################################
## INFORMATION FOR NMIX
###################################
BLOCK NMIX # neutralino mixing matrix
      1 1 9.999509e-01 # n_11
      1 2 -3.909008e-05 # n_12
      1 3 9.831832e-03 # n_13
      1 4 -1.239245e-03 # n_14
      2 1 6.298358e-03 # n_21
      2 2 6.012981e-01 # n_22
      2 3 -5.672636e-01 # n_23
      2 4 5.626836e-01 # n_24
      3 1 6.075303e-03 # n_31
      3 2 -5.371076e-03 # n_32
      3 3 -7.070419e-01 # n_33
      3 4 -7.071252e-01 # n_34
      4 1 4.650113e-03 # n_41
      4 2 -7.990067e-01 # n_42
      4 3 -4.221458e-01 # n_43
      4 4 4.282050e-01 # n_44
###################################
## INFORMATION FOR UMIX
###################################
BLOCK UMIX # chargino mixing matrix u
      1 1 -5.969816e-01 # u_11
      1 2 8.022549e-01 # u_12
      2 1 8.022549e-01 # u_21
      2 2 5.969816e-01 # u_22
###################################
## INFORMATION FOR VMIX
###################################
BLOCK VMIX # chargino mixing matrix v
      1 1 -6.055647e-01 # v_11
      1 2 7.957961e-01 # v_12
      2 1 7.957961e-01 # v_21
      2 2 6.055647e-01 # v_22
###################################
## INFORMATION FOR STOPMIX
###################################
BLOCK STOPMIX # stop mixing matrix
      1 1 1.941337e-03 # cos(theta_t)
      1 2 9.999981e-01 # sin(theta_t)
      2 1 -9.999981e-01 # -sin(theta_t)
      2 2 1.941337e-03 # cos(theta_t)
###################################
## INFORMATION FOR SBOTMIX
###################################
BLOCK SBOTMIX # sbottom mixing matrix
      1 1 1.915475e-02 # cos(theta_b)
      1 2 9.998165e-01 # sin(theta_b)
      2 1 -9.998165e-01 # -sin(theta_b)
      2 2 1.915475e-02 # cos(theta_b)
###################################
## INFORMATION FOR STAUMIX
###################################
BLOCK STAUMIX # stau mixing matrix
      1 1 7.070742e-01 # cos(theta_tau)
      1 2 7.071394e-01 # sin(theta_tau)
      2 1 -7.071394e-01 # -sin(theta_tau)
      2 2 7.070742e-01 # cos(theta_tau)
###################################
## INFORMATION FOR ALPHA
###################################
BLOCK ALPHA # higgs mixing
       -1.038362e-01 # mixing angle in the neutral higgs boson sector
###################################
## INFORMATION FOR HMIX
###################################
BLOCK HMIX Q= 1.690729e+03 # drbar higgs parameters
      1 4.500000e+03 # mu(q)
      2 9.603993e+00 # tanbeta(q)
      3 2.460028e+02 # vev(q)
      4 2.198974e+07 # ma^2(q)
###################################
## INFORMATION FOR GAUGE
###################################
BLOCK GAUGE Q= 1.690729e+03 # the gauge couplings
      1 3.621302e-01 # gprime(q) drbar
      2 6.287497e-01 # g(q) drbar
      3 1.010183e+00 # g3(q) drbar
###################################
## INFORMATION FOR AU
###################################
BLOCK AU Q= 1.690729e+03 # the trilinear couplings
      1 1 0.000000e+00 # a_u(q) drbar
      2 2 0.000000e+00 # a_c(q) drbar
      3 3 0.000000e+00 # a_t(q) drbar
###################################
## INFORMATION FOR AD
###################################
BLOCK AD Q= 1.690729e+03 # the trilinear couplings
      1 1 0.000000e+00 # a_d(q) drbar
      2 2 0.000000e+00 # a_s(q) drbar
      3 3 0.000000e+00 # a_b(q) drbar
###################################
## INFORMATION FOR AE
###################################
BLOCK AE Q= 1.690729e+03 # the trilinear couplings
      1 1 0.000000e+00 # a_e(q) drbar
      2 2 0.000000e+00 # a_mu(q) drbar
      3 3 0.000000e+00 # a_tau(q) drbar
###################################
## INFORMATION FOR YU
###################################
BLOCK YU Q= 1.690729e+03 # the yukawa couplings
      1 1 0.000000e+00 # y_u(q) drbar
      2 2 0.000000e+00 # y_c(q) drbar
      3 3 8.457204e-01 # y_t(q) drbar
###################################
## INFORMATION FOR YD
###################################
BLOCK YD Q= 1.690729e+03 # the yukawa couplings
      1 1 0.000000e+00 # y_d(q) drbar
      2 2 0.000000e+00 # y_s(q) drbar
      3 3 1.292987e-01 # y_b(q) drbar
###################################
## INFORMATION FOR YE
###################################
BLOCK YE Q= 1.690729e+03 # the yukawa couplings
      1 1 0.000000e+00 # y_e(q) drbar
      2 2 0.000000e+00 # y_mu(q) drbar
      3 3 9.854853e-02 # y_tau(q) drbar
###################################
## INFORMATION FOR MSOFT
###################################
BLOCK MSOFT Q= 1.690729e+03 # the soft susy breaking masses at the scale q
      1 1.000000e+02 # m_1
      2 4.500000e+03 # m_2
      3 3.000000e+03 # m_3
      14 0.000000e+00 # a_u
      15 0.000000e+00 # a_d
      16 0.000000e+00 # a_e
      21 4.312761e+06 # m^2_hd
      22 -2.009445e+07 # m^2_hu
      31 4.000000e+03 # m_el
      32 4.000000e+03 # m_mul
      33 4.000000e+03 # m_taul
      34 4.000000e+03 # m_er
      35 4.000000e+03 # m_mur
      36 4.000000e+03 # m_taur
      41 5.000000e+03 # m_q1l
      42 5.000000e+03 # m_q2l
      43 5.500000e+03 # m_q3l
      44 5.000000e+03 # m_ur
      45 5.000000e+03 # m_cr
      46 5.000000e+02 # m_tr
      47 5.000000e+03 # m_dr
      48 5.000000e+03 # m_sr
      49 5.000000e+03 # m_br
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.384015e+00 # top decays
      9.998266e-01 2 5 24 # br(t -> b w+)
      1.734001e-04 2 1000006 1000022 # br(t -> ~t_1 ~chi_10)

DECAY 1000021 6.299104e+01 # gluino decays
      5.000000e-01 2 -6 1000006 # br(~g -> ~t_1 tb)
      5.000000e-01 2 -1000006 6 # br(~g -> ~t_1* t )

DECAY 1000006 2.383589e-04 # stop1 decays
      1.522873e-10 2 4 1000022 # br(~t_1 -> ~chi_10 c )
      1.174089e-12 2 2 1000022 # br(~t_1 -> ~chi_10 u )
      1.000000e+00 3 5 24 1000022 # br(~t_1 -> ~chi_10 b w+)

DECAY 2000006 1.726944e+02 # stop2 decays
      4.430768e-03 2 6 1000022 # br(~t_2 -> ~chi_10 t )
      1.786205e-02 2 6 1000023 # br(~t_2 -> ~chi_20 t )
      2.208935e-02 2 6 1000025 # br(~t_2 -> ~chi_30 t )
      1.632983e-02 2 6 1000035 # br(~t_2 -> ~chi_40 t )
      1.116970e-02 2 5 1000024 # br(~t_2 -> ~chi_1+ b )
      1.481126e-02 2 5 1000037 # br(~t_2 -> ~chi_2+ b )
      8.556523e-01 2 6 1000021 # br(~t_2 -> ~g t )
      1.533199e-03 2 25 1000006 # br(~t_2 -> ~t_1 h )
      2.699958e-02 2 35 1000006 # br(~t_2 -> ~t_1 h )
      2.700311e-02 2 36 1000006 # br(~t_2 -> ~t_1 a )
      9.274286e-04 2 23 1000006 # br(~t_2 -> ~t_1 z )
      1.191388e-03 2 24 1000005 # br(~t_2 -> ~b_1 w+)

DECAY 1000005 1.121689e+02 # sbottom1 decays
      2.484223e-02 2 5 1000022 # br(~b_1 -> ~chi_10 b )
      2.193717e-04 2 5 1000023 # br(~b_1 -> ~chi_20 b )
      2.536812e-04 2 5 1000025 # br(~b_1 -> ~chi_30 b )
      5.031790e-05 2 5 1000035 # br(~b_1 -> ~chi_40 b )
      3.961991e-04 2 -1000024 6 # br(~b_1 -> ~chi_1- t )
      8.426855e-05 2 -1000037 6 # br(~b_1 -> ~chi_2- t )
      9.741542e-01 2 5 1000021 # br(~b_1 -> ~g b )
      -2.893554e-07 2 -24 1000006 # br(~b_1 -> ~t_1 w-)

DECAY 2000005 1.739174e+02 # sbottom2 decays
      4.416236e-03 2 5 1000022 # br(~b_2 -> ~chi_10 b )
      5.537652e-03 2 5 1000023 # br(~b_2 -> ~chi_20 b )
      5.411492e-04 2 5 1000025 # br(~b_2 -> ~chi_30 b )
      7.615011e-03 2 5 1000035 # br(~b_2 -> ~chi_40 b )
      3.607299e-02 2 -1000024 6 # br(~b_2 -> ~chi_1- t )
      3.332891e-02 2 -1000037 6 # br(~b_2 -> ~chi_2- t )
      8.578812e-01 2 5 1000021 # br(~b_2 -> ~g b )
      5.368320e-04 2 25 1000005 # br(~b_2 -> ~b_1 h )
      5.462553e-02 2 -37 1000006 # br(~b_2 -> ~t_1 h-)
      7.131238e-04 2 23 1000005 # br(~b_2 -> ~b_1 z )
      -1.268599e-03 2 -24 1000006 # br(~b_2 -> ~t_1 w-)

DECAY 1000002 1.118491e+02 # sup_l decays
      6.224337e-03 2 2 1000022 # br(~u_l -> ~chi_10 u)
      2.720034e-03 2 2 1000023 # br(~u_l -> ~chi_20 u)
      1.127827e-07 2 2 1000025 # br(~u_l -> ~chi_30 u)
      2.939910e-03 2 2 1000035 # br(~u_l -> ~chi_40 u)
      5.495555e-03 2 1 1000024 # br(~u_l -> ~chi_1+ d)
      5.845764e-03 2 1 1000037 # br(~u_l -> ~chi_2+ d)
      9.767743e-01 2 2 1000021 # br(~u_l -> ~g u)

DECAY 2000002 1.204064e+02 # sup_r decays
      9.255171e-02 2 2 1000022 # br(~u_r -> ~chi_10 u)
      1.629001e-07 2 2 1000023 # br(~u_r -> ~chi_20 u)
      1.290450e-07 2 2 1000025 # br(~u_r -> ~chi_30 u)
      5.470263e-08 2 2 1000035 # br(~u_r -> ~chi_40 u)
      9.074479e-01 2 2 1000021 # br(~u_r -> ~g u)

DECAY 1000001 1.118966e+02 # sdown_l decays
      6.227513e-03 2 1 1000022 # br(~d_l -> ~chi_10 d)
      2.702173e-03 2 1 1000023 # br(~d_l -> ~chi_20 d)
      2.730859e-07 2 1 1000025 # br(~d_l -> ~chi_30 d)
      2.959211e-03 2 1 1000035 # br(~d_l -> ~chi_40 d)
      5.348654e-03 2 -1000024 2 # br(~d_l -> ~chi_1- u)
      5.953378e-03 2 -1000037 2 # br(~d_l -> ~chi_2- u)
      9.768088e-01 2 1 1000021 # br(~d_l -> ~g d)

DECAY 2000001 1.120653e+02 # sdown_r decays
      2.486109e-02 2 1 1000022 # br(~d_r -> ~chi_10 d)
      4.378352e-08 2 1 1000023 # br(~d_r -> ~chi_20 d)
      3.468627e-08 2 1 1000025 # br(~d_r -> ~chi_30 d)
      1.470574e-08 2 1 1000035 # br(~d_r -> ~chi_40 d)
      9.751388e-01 2 1 1000021 # br(~d_r -> ~g d)

DECAY 1000004 1.118491e+02 # scharm_l decays
      6.224337e-03 2 4 1000022 # br(~c_l -> ~chi_10 c)
      2.720034e-03 2 4 1000023 # br(~c_l -> ~chi_20 c)
      1.127827e-07 2 4 1000025 # br(~c_l -> ~chi_30 c)
      2.939910e-03 2 4 1000035 # br(~c_l -> ~chi_40 c)
      5.495555e-03 2 3 1000024 # br(~c_l -> ~chi_1+ s)
      5.845764e-03 2 3 1000037 # br(~c_l -> ~chi_2+ s)
      9.767743e-01 2 4 1000021 # br(~c_l -> ~g c)

DECAY 2000004 1.204064e+02 # scharm_r decays
      9.255171e-02 2 4 1000022 # br(~c_r -> ~chi_10 c)
      1.629001e-07 2 4 1000023 # br(~c_r -> ~chi_20 c)
      1.290450e-07 2 4 1000025 # br(~c_r -> ~chi_30 c)
      5.470263e-08 2 4 1000035 # br(~c_r -> ~chi_40 c)
      9.074479e-01 2 4 1000021 # br(~c_r -> ~g c)

DECAY 1000003 1.118966e+02 # sstrange_l decays
      6.227513e-03 2 3 1000022 # br(~s_l -> ~chi_10 s)
      2.702173e-03 2 3 1000023 # br(~s_l -> ~chi_20 s)
      2.730859e-07 2 3 1000025 # br(~s_l -> ~chi_30 s)
      2.959211e-03 2 3 1000035 # br(~s_l -> ~chi_40 s)
      5.348654e-03 2 -1000024 4 # br(~s_l -> ~chi_1- c)
      5.953378e-03 2 -1000037 4 # br(~s_l -> ~chi_2- c)
      9.768088e-01 2 3 1000021 # br(~s_l -> ~g s)

DECAY 2000003 1.120653e+02 # sstrange_r decays
      2.486109e-02 2 3 1000022 # br(~s_r -> ~chi_10 s)
      4.378352e-08 2 3 1000023 # br(~s_r -> ~chi_20 s)
      3.468627e-08 2 3 1000025 # br(~s_r -> ~chi_30 s)
      1.470574e-08 2 3 1000035 # br(~s_r -> ~chi_40 s)
      9.751388e-01 2 3 1000021 # br(~s_r -> ~g s)

DECAY 1000011 5.216924e+00 # selectron_l decays
      1.000000e+00 2 11 1000022 # br(~e_l -> ~chi_10 e-)

DECAY 2000011 2.087052e+01 # selectron_r decays
      1.000000e+00 2 11 1000022 # br(~e_r -> ~chi_10 e-)

DECAY 1000013 5.216924e+00 # smuon_l decays
      1.000000e+00 2 13 1000022 # br(~mu_l -> ~chi_10 mu-)

DECAY 2000013 2.087052e+01 # smuon_r decays
      1.000000e+00 2 13 1000022 # br(~mu_r -> ~chi_10 mu-)

DECAY 1000015 1.303291e+01 # stau_1 decays
      1.000000e+00 2 15 1000022 # br(~tau_1 -> ~chi_10 tau-)

DECAY 2000015 1.305451e+01 # stau_2 decays
      1.000000e+00 2 15 1000022 # br(~tau_2 -> ~chi_10 tau-)

DECAY 1000012 5.217386e+00 # snu_el decays
      1.000000e+00 2 12 1000022 # br(~nu_el -> ~chi_10 nu_e)

DECAY 1000014 5.217386e+00 # snu_mul decays
      1.000000e+00 2 14 1000022 # br(~nu_mul -> ~chi_10 nu_mu)

DECAY 1000016 5.217386e+00 # snu_taul decays
      1.000000e+00 2 16 1000022 # br(~nu_taul -> ~chi_10 nu_tau)

DECAY 1000024 6.626610e+01 # chargino1+ decays
      9.502651e-01 2 -5 1000006 # br(~chi_1+ -> ~t_1 bb)
      3.926311e-03 2 -11 1000012 # br(~chi_1+ -> ~nu_el e+ )
      3.926311e-03 2 -13 1000014 # br(~chi_1+ -> ~nu_mul mu+ )
      4.102026e-03 2 -15 1000016 # br(~chi_1+ -> ~nu_tau1 tau+)
      3.804695e-03 2 -1000011 12 # br(~chi_1+ -> ~e_l+ nu_e)
      3.804695e-03 2 -1000013 14 # br(~chi_1+ -> ~mu_l+ nu_mu)
      2.895584e-03 2 -1000015 16 # br(~chi_1+ -> ~tau_1+ nu_tau)
      1.140635e-03 2 -2000015 16 # br(~chi_1+ -> ~tau_2+ nu_tau)
      2.613462e-02 2 24 1000022 # br(~chi_1+ -> ~chi_10 w+)

DECAY 1000037 4.318634e+01 # chargino2+ decays
      8.614224e-01 2 -5 1000006 # br(~chi_2+ -> ~t_1 bb)
      1.574640e-02 2 -11 1000012 # br(~chi_2+ -> ~nu_el e+ )
      1.574640e-02 2 -13 1000014 # br(~chi_2+ -> ~nu_mul mu+ )
      1.595229e-02 2 -15 1000016 # br(~chi_2+ -> ~nu_tau1 tau+)
      1.596671e-02 2 -1000011 12 # br(~chi_2+ -> ~e_l+ nu_e)
      1.596671e-02 2 -1000013 14 # br(~chi_2+ -> ~mu_l+ nu_mu)
      6.416065e-03 2 -1000015 16 # br(~chi_2+ -> ~tau_1+ nu_tau)
      9.660400e-03 2 -2000015 16 # br(~chi_2+ -> ~tau_2+ nu_tau)
      8.076885e-03 2 23 1000024 # br(~chi_2+ -> ~chi_1+ z )
      2.351522e-02 2 24 1000022 # br(~chi_2+ -> ~chi_10 w+)
      1.088108e-02 2 24 1000023 # br(~chi_2+ -> ~chi_20 w+)
      6.494651e-04 2 25 1000024 # br(~chi_2+ -> ~chi_1+ h )

DECAY 1000022 0.000000e+00 # neutralino1 decays
DECAY 1000023 6.632060e+01 # neutralino2 decays
      1.012535e-02 2 23 1000022 # br(~chi_20 -> ~chi_10 z )
      1.694038e-02 2 25 1000022 # br(~chi_20 -> ~chi_10 h )
      4.746754e-01 2 -6 1000006 # br(~chi_20 -> ~t_1 tb)
      4.746754e-01 2 -1000006 6 # br(~chi_20 -> ~t_1* t )
      1.951762e-03 2 -11 1000011 # br(~chi_20 -> ~e_l- e+)
      1.951762e-03 2 -1000011 11 # br(~chi_20 -> ~e_l+ e-)
      2.807466e-07 2 -11 2000011 # br(~chi_20 -> ~e_r- e+)
      2.807466e-07 2 -2000011 11 # br(~chi_20 -> ~e_r+ e-)
      1.951762e-03 2 -13 1000013 # br(~chi_20 -> ~mu_l- mu+)
      1.951762e-03 2 -1000013 13 # br(~chi_20 -> ~mu_l+ mu-)
      2.807466e-07 2 -13 2000013 # br(~chi_20 -> ~mu_r- mu+)
      2.807466e-07 2 -2000013 13 # br(~chi_20 -> ~mu_r+ mu-)
      1.519337e-03 2 -15 1000015 # br(~chi_20 -> ~tau_1- tau+)
      1.519337e-03 2 -1000015 15 # br(~chi_20 -> ~tau_1+ tau-)
      6.359525e-04 2 -15 2000015 # br(~chi_20 -> ~tau_2- tau+)
      6.359525e-04 2 -2000015 15 # br(~chi_20 -> ~tau_2+ tau-)
      1.910787e-03 2 -12 1000012 # br(~chi_20 -> ~nu_el nu_eb)
      1.910787e-03 2 -1000012 12 # br(~chi_20 -> ~nu_el* nu_e )
      1.910787e-03 2 -14 1000014 # br(~chi_20 -> ~nu_mul nu_mub)
      1.910787e-03 2 -1000014 14 # br(~chi_20 -> ~nu_mul* nu_mu )
      1.910787e-03 2 -16 1000016 # br(~chi_20 -> ~nu_tau1 nu_taub)
      1.910787e-03 2 -1000016 16 # br(~chi_20 -> ~nu_tau1* nu_tau )

DECAY 1000025 -1.031242e+02 # neutralino3 decays
      1.726999e-02 2 23 1000022 # br(~chi_30 -> ~chi_10 z )
      1.132677e-02 2 25 1000022 # br(~chi_30 -> ~chi_10 h )
      4.855019e-01 2 -6 1000006 # br(~chi_30 -> ~t_1 tb)
      4.855019e-01 2 -1000006 6 # br(~chi_30 -> ~t_1* t )
      1.415272e-08 2 -11 1000011 # br(~chi_30 -> ~e_l- e+)
      1.415272e-08 2 -1000011 11 # br(~chi_30 -> ~e_l+ e-)
      1.977911e-07 2 -11 2000011 # br(~chi_30 -> ~e_r- e+)
      1.977911e-07 2 -2000011 11 # br(~chi_30 -> ~e_r+ e-)
      1.415272e-08 2 -13 1000013 # br(~chi_30 -> ~mu_l- mu+)
      1.415272e-08 2 -1000013 13 # br(~chi_30 -> ~mu_l+ mu-)
      1.977911e-07 2 -13 2000013 # br(~chi_30 -> ~mu_r- mu+)
      1.977911e-07 2 -2000013 13 # br(~chi_30 -> ~mu_r+ mu-)
      9.662209e-05 2 -15 1000015 # br(~chi_30 -> ~tau_1- tau+)
      9.662209e-05 2 -1000015 15 # br(~chi_30 -> ~tau_1+ tau-)
      1.016818e-04 2 -15 2000015 # br(~chi_30 -> ~tau_2- tau+)
      1.016818e-04 2 -2000015 15 # br(~chi_30 -> ~tau_2+ tau-)
      3.186004e-07 2 -12 1000012 # br(~chi_30 -> ~nu_el nu_eb)
      3.186004e-07 2 -1000012 12 # br(~chi_30 -> ~nu_el* nu_e )
      3.186004e-07 2 -14 1000014 # br(~chi_30 -> ~nu_mul nu_mub)
      3.186004e-07 2 -1000014 14 # br(~chi_30 -> ~nu_mul* nu_mu )
      3.186004e-07 2 -16 1000016 # br(~chi_30 -> ~nu_tau1 nu_taub)
      3.186004e-07 2 -1000016 16 # br(~chi_30 -> ~nu_tau1* nu_tau )

DECAY 1000035 4.334558e+01 # neutralino4 decays
      9.258896e-03 2 23 1000022 # br(~chi_40 -> ~chi_10 z )
      3.350083e-07 2 23 1000023 # br(~chi_40 -> ~chi_20 z )
      1.084087e-02 2 -24 1000024 # br(~chi_40 -> ~chi_1+ w-)
      1.084087e-02 2 -1000024 24 # br(~chi_40 -> ~chi_1- w+)
      1.504795e-02 2 25 1000022 # br(~chi_40 -> ~chi_10 h )
      6.470002e-04 2 25 1000023 # br(~chi_40 -> ~chi_20 h )
      4.291304e-01 2 -6 1000006 # br(~chi_40 -> ~t_1 tb)
      4.291304e-01 2 -1000006 6 # br(~chi_40 -> ~t_1* t )
      7.836989e-03 2 -11 1000011 # br(~chi_40 -> ~e_l- e+)
      7.836989e-03 2 -1000011 11 # br(~chi_40 -> ~e_l+ e-)
      3.545913e-07 2 -11 2000011 # br(~chi_40 -> ~e_r- e+)
      3.545913e-07 2 -2000011 11 # br(~chi_40 -> ~e_r+ e-)
      7.836989e-03 2 -13 1000013 # br(~chi_40 -> ~mu_l- mu+)
      7.836989e-03 2 -1000013 13 # br(~chi_40 -> ~mu_l+ mu-)
      3.545913e-07 2 -13 2000013 # br(~chi_40 -> ~mu_r- mu+)
      3.545913e-07 2 -2000013 13 # br(~chi_40 -> ~mu_r+ mu-)
      3.190043e-03 2 -15 1000015 # br(~chi_40 -> ~tau_1- tau+)
      3.190043e-03 2 -1000015 15 # br(~chi_40 -> ~tau_1+ tau-)
      4.804325e-03 2 -15 2000015 # br(~chi_40 -> ~tau_2- tau+)
      4.804325e-03 2 -2000015 15 # br(~chi_40 -> ~tau_2+ tau-)
      7.960846e-03 2 -12 1000012 # br(~chi_40 -> ~nu_el nu_eb)
      7.960846e-03 2 -1000012 12 # br(~chi_40 -> ~nu_el* nu_e )
      7.960846e-03 2 -14 1000014 # br(~chi_40 -> ~nu_mul nu_mub)
      7.960846e-03 2 -1000014 14 # br(~chi_40 -> ~nu_mul* nu_mu )
      7.960846e-03 2 -16 1000016 # br(~chi_40 -> ~nu_tau1 nu_taub)
      7.960846e-03 2 -1000016 16 # br(~chi_40 -> ~nu_tau1* nu_tau )

DECAY 25 3.681770e-03 # h decays
      6.867258e-01 2 -5 5 # br(h -> b bb )
      6.664984e-02 2 -15 15 # br(h -> tau+ tau- )
      2.359488e-04 2 -13 13 # br(h -> mu+ mu- )
      5.054153e-04 2 -3 3 # br(h -> s sb )
      2.197472e-02 2 -4 4 # br(h -> c cb )
      1.047574e-01 2 21 21 # br(h -> g g )
      1.724751e-03 2 22 22 # br(h -> gam gam )
      8.886820e-04 2 22 23 # br(h -> z gam )
      1.044273e-01 2 -24 24 # br(h -> w+ w- )
      1.171845e-02 2 23 23 # br(h -> z z )
      3.916498e-04 2 1000022 1000022 # br(h -> ~chi_10 ~chi_10)

DECAY 35 8.071660e+00 # h decays
      5.964804e-01 2 -5 5 # br(h -> b bb )
      1.184132e-01 2 -15 15 # br(h -> tau+ tau- )
      4.186331e-04 2 -13 13 # br(h -> mu+ mu- )
      4.805678e-04 2 -3 3 # br(h -> s sb )
      2.533196e-06 2 -4 4 # br(h -> c cb )
      2.566876e-01 2 -6 6 # br(h -> t tb )
      5.442847e-05 2 21 21 # br(h -> g g )
      2.487527e-07 2 22 22 # br(h -> gam gam )
      1.116878e-07 2 22 23 # br(h -> z gam )
      3.811707e-05 2 -24 24 # br(h -> w+ w- )
      1.905016e-05 2 23 23 # br(h -> z z )
      1.375160e-04 2 25 25 # br(h -> h h )
      -2.348497e-23 2 36 36 # br(h -> a a )
      1.981762e-19 2 23 36 # br(h -> z a )
      1.386824e-04 2 1000022 1000022 # br(h -> ~chi_10 ~chi_10)
      7.889579e-03 2 1000022 1000023 # br(h -> ~chi_10 ~chi_20)
      1.507047e-02 2 1000022 1000025 # br(h -> ~chi_10 ~chi_30)
      2.295605e-03 2 1000022 1000035 # br(h -> ~chi_10 ~chi_40)
      1.873258e-03 2 -1000006 1000006 # br(h -> ~t_1 ~t_1* )

DECAY 36 8.061114e+00 # a decays
      5.974659e-01 2 -5 5 # br(a -> b bb )
      1.185693e-01 2 -15 15 # br(a -> tau+ tau- )
      4.191848e-04 2 -13 13 # br(a -> mu+ mu- )
      4.815413e-04 2 -3 3 # br(a -> s sb )
      2.555718e-06 2 -4 4 # br(a -> c cb )
      2.573432e-01 2 -6 6 # br(a -> t tb )
      1.450918e-04 2 21 21 # br(a -> g g )
      4.280382e-07 2 22 22 # br(a -> gam gam )
      2.226292e-07 2 22 23 # br(a -> z gam )
      3.812304e-05 2 23 25 # br(a -> z h )
      1.463535e-04 2 1000022 1000022 # br(a -> ~chi_10 ~chi_10)
      1.191908e-02 2 1000022 1000023 # br(a -> ~chi_10 ~chi_20)
      9.891029e-03 2 1000022 1000025 # br(a -> ~chi_10 ~chi_30)
      3.577937e-03 2 1000022 1000035 # br(a -> ~chi_10 ~chi_40)

DECAY 37 7.711499e+00 # h+ decays
      9.319699e-04 2 -5 4 # br(h+ -> c bb )
      1.239573e-01 2 -15 16 # br(h+ -> tau+ nu_tau )
      4.382331e-04 2 -13 14 # br(h+ -> mu+ nu_mu )
      5.964578e-06 2 -5 2 # br(h+ -> u bb )
      2.411435e-05 2 -3 2 # br(h+ -> u sb )
      4.985375e-04 2 -3 4 # br(h+ -> c sb )
      8.472017e-01 2 -5 6 # br(h+ -> t bb )
      3.987218e-05 2 24 25 # br(h+ -> w+ h )
      5.360521e-15 2 24 36 # br(h+ -> w+ a )
      2.074788e-02 2 1000022 1000024 # br(h+ -> ~chi_1+ ~chi_10)
      6.154432e-03 2 1000022 1000037 # br(h+ -> ~chi_2+ ~chi_10)

DECAY 23 2.411433e+00 # wz
DECAY 24 2.002822e+00 # ww
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
###################################
## INFORMATION FOR SNUMIX
###################################
BLOCK SNUMIX #
      1 1 1.000000e+00 # rrn1x1
      2 2 1.000000e+00 # rrn2x2
      3 3 1.000000e+00 # rrn3x3
###################################
## INFORMATION FOR MSU2
###################################
BLOCK MSU2 #
      1 1 2.803821e+05 # set of param :1*rmu21x1, 1*rmu22x2
      2 2 2.803821e+05 # mg5 will not use this value use instead 1*mdl_rmu21x1
      3 3 1.791371e+05 # rmu23x3
###################################
## INFORMATION FOR VCKM
###################################
BLOCK VCKM #
      1 1 1.000000e+00 # rckm1x1
      2 2 1.000000e+00 # rckm2x2
      3 3 1.000000e+00 # rckm3x3
###################################
## INFORMATION FOR FRALPHA
###################################
BLOCK FRALPHA #
      1 -1.138252e-01 # alp
###################################
## INFORMATION FOR UPMNS
###################################
BLOCK UPMNS #
      1 1 1.000000e+00 # rmns1x1
      2 2 1.000000e+00 # rmns2x2
      3 3 1.000000e+00 # rmns3x3
###################################
## INFORMATION FOR MSD2
###################################
BLOCK MSD2 #
      1 1 2.736847e+05 # set of param :1*rmd21x1, 1*rmd22x2
      2 2 2.736847e+05 # mg5 will not use this value use instead 1*mdl_rmd21x1
      3 3 2.702620e+05 # rmd23x3
###################################
## INFORMATION FOR MSL2
###################################
BLOCK MSL2 #
      1 1 3.815567e+04 # set of param :1*rml21x1, 1*rml22x2
      2 2 3.815567e+04 # mg5 will not use this value use instead 1*mdl_rml21x1
      3 3 3.782868e+04 # rml23x3
###################################
## INFORMATION FOR TU
###################################
BLOCK TU #
      3 3 -4.447525e+02 # rtu3x3
###################################
## INFORMATION FOR SELMIX
###################################
BLOCK SELMIX #
      1 1 1.000000e+00 # rrl1x1
      2 2 1.000000e+00 # rrl2x2
      3 3 2.824872e-01 # rrl3x3
      3 6 9.592711e-01 # rrl3x6
      4 4 1.000000e+00 # rrl4x4
      5 5 1.000000e+00 # rrl5x5
      6 3 9.592711e-01 # rrl6x3
      6 6 -2.824872e-01 # rrl6x6
###################################
## INFORMATION FOR TD
###################################
BLOCK TD #
      3 3 -1.106937e+02 # rtd3x3
###################################
## INFORMATION FOR TE
###################################
BLOCK TE #
      3 3 -2.540197e+01 # rte3x3
###################################
## INFORMATION FOR USQMIX
###################################
BLOCK USQMIX #
      1 1 1.000000e+00 # rru1x1
      2 2 1.000000e+00 # rru2x2
      3 3 5.536450e-01 # rru3x3
      3 6 8.327528e-01 # rru3x6
      4 4 1.000000e+00 # rru4x4
      5 5 1.000000e+00 # rru5x5
      6 3 8.327528e-01 # rru6x3
      6 6 -5.536450e-01 # rru6x6
###################################
## INFORMATION FOR DSQMIX
###################################
BLOCK DSQMIX #
      1 1 1.000000e+00 # rrd1x1
      2 2 1.000000e+00 # rrd2x2
      3 3 9.387379e-01 # rrd3x3
      3 6 3.446319e-01 # rrd3x6
      4 4 1.000000e+00 # rrd4x4
      5 5 1.000000e+00 # rrd5x5
      6 3 -3.446319e-01 # rrd6x3
      6 6 9.387379e-01 # rrd6x6
###################################
## INFORMATION FOR MSE2
###################################
BLOCK MSE2 #
      1 1 1.863063e+04 # set of param :1*rme21x1, 1*rme22x2
      2 2 1.863063e+04 # mg5 will not use this value use instead 1*mdl_rme21x1
      3 3 1.796764e+04 # rme23x3
###################################
## INFORMATION FOR MSQ2
###################################
BLOCK MSQ2 #
      1 1 2.998367e+05 # set of param :1*rmq21x1, 1*rmq22x2
      2 2 2.998367e+05 # mg5 will not use this value use instead 1*mdl_rmq21x1
      3 3 2.487654e+05 # rmq23x3
#*********************************************************************
# 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 *
#*********************************************************************
#
#*******************
# 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. *
#*********************************************************************
  50000 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     4000.0 = ebeam1 ! beam 1 total energy in GeV
     4000.0 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
     0.0 = polbeam1 ! beam polarization for beam 1
     0.0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# 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)
#*********************************************************************
# 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
#*********************************************************************
# 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
#*********************************************************************
# 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
#*********************************************************************
# 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, 1000001, 1000002, 1000003, 1000004, 1000005, 1000006, 1000021, 2000001, 2000002, 2000003, 2000004, 2000005, 2000006 = 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
#
#**************************************
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#**************************************
#
0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_lo_as_0130_qed = sys_pdf # list of pdf sets
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf
#