MadGraph5_aMC@NLO

ValueError: invalid literal for float()

Asked by yuanfang

Hi
    I have generated 50000 events using madgraph and it worked well at first . Then it ran into error when madgraph tried to combine events ..
     I generated my spectrum using SUSY_HIT then I found there's two 'NAN' s in BLOCK MSOFT
                     41 NAN # M_q1L
                    42 NaN # M_q2L
    then Madgraph had some compiling error ,so I modified my spectrum these two lines like this(use # to coment these two lines):
                           41 5.00000000E+04 # NaN, M_q1L
                           42 5.00000000E+04 # NaN i, M_q2L
     then madgraph worked well and began compile and creating jobs ... All worked well .

       then it ran into error when combining runs ...
       Here's my log file . Hope it can help .

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.4.3 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 *
#* *
#************************************************************
generate_events
Traceback (most recent call last):
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/extended_cmd.py", line 1011, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/extended_cmd.py", line 966, in onecmd_orig
    return func(arg, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/madevent_interface.py", line 2015, in do_generate_events
    self.exec_cmd('combine_events', postcmd=False)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/extended_cmd.py", line 1037, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/extended_cmd.py", line 966, in onecmd_orig
    return func(arg, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/madevent_interface.py", line 3008, in do_combine_events
    get_wgt, log_level=logging.DEBUG, write_init=True)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 882, in unweight
    return super(MultiEventFile, self).unweight(outputpath, get_wgt_multi, **opts)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 315, in unweight
    all_wgt, cross, nb_event = self.initialize_unweighting(get_wgt, trunc_error)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 794, in initialize_unweighting
    for event in f:
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 248, in next
    return Event(text)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 910, in __init__
    self.parse(text)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 938, in parse
    self.assign_scale_line(line)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/proess_2/bin/internal/lhe_parser.py", line 1327, in assign_scale_line
    self.wgt=float(inputs[2])
ValueError: invalid literal for float(): 0.8317976-100
                              Run Options
                              -----------
               stdout_level : None

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : /cvmfs/cp3.uclouvain.be/madgraph/
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : /home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/pythia-pgs (user set)
                    td_path : /home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/td (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
        exrootanalysis_path : /home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/ExRootAnalysis (user set)
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : /home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/MadAnalysis (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)
           cluster_nb_retry : 1
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.4.3 2016-08-01 *
#* *
#* 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 low_mem_multicore_nlo_generation False
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
generate p p > x1+ x1- a , x1+ > n1 l+ vl , x1- > n1 l- vl~
output proess_2
#
# ======================
# | THE SUSYHIT OUTPUT |
# ======================
#
#
# ------------------------------------------------------
# | This is the output of the SUSY-HIT package |
# | created by A.Djouadi, M.Muehlleitner and M.Spira. |
# | In case of problems with SUSY-HIT email to |
# | <email address hidden> |
# | <email address hidden> |
# | <email address hidden> |
# ------------------------------------------------------
#
# ------------------------------------------------------
# | SUSY Les Houches Accord - MSSM Spectrum + Decays |
# | based on the decay programs |
# | |
# | SDECAY 1.5 |
# | |
# | Authors: M.Muhlleitner, A.Djouadi and Y.Mambrini |
# | Ref.: Comput.Phys.Commun.168(2005)46 |
# | [hep-ph/0311167] |
# | |
# | HDECAY 3.4 |
# | |
# | By: A.Djouadi,J.Kalinowski,M.Muhlleitner,M.Spira |
# | Ref.: Comput.Phys.Commun.108(1998)56 |
# | [hep-ph/9704448] |
# | |
# | |
# | If not stated otherwise all DRbar couplings and |
# | soft SUSY breaking masses are given at the scale |
# | Q= 0.30031502E+04
# | |
# ------------------------------------------------------
#
#
BLOCK DCINFO # Decay Program information
     1 SDECAY/HDECAY # decay calculator
     2 1.5 /3.4 # version number
#
BLOCK SPINFO # Spectrum calculator information
     1 SuSpect # RGE +Spectrum calculator
     2 2.41 # version number
#
BLOCK MODSEL # Model selection
     1 -1 # #bottom-up MSSM
#
BLOCK SMINPUTS # Standard Model inputs
         1 1.27934000E+02 # alpha_em^-1(M_Z)^MSbar
         2 1.16639000E-05 # G_F [GeV^-2]
         3 1.17200000E-01 # alpha_S(M_Z)^MSbar
         4 9.11870000E+01 # M_Z pole mass
         5 4.25000000E+00 # mb(mb)^MSbar
         6 1.72500000E+02 # mt pole mass
         7 1.77710000E+00 # mtau pole mass
#
BLOCK MINPAR # Input parameters - minimal models
#
BLOCK EXTPAR # Input parameters - non-minimal models
         0 3.00315018E+03 # EWSB
         1 1.60000000E+02 # M_1
         2 3.00000000E+03 # M_2
         3 3.00000000E+03 # M_3
        11 -1.00000000E+02 # A_t
        12 -1.00000000E+02 # A_b
        13 -2.51737263E+02 # A_tau
        14 -6.77437438E+02 # A_u
        15 -8.59633345E+02 # A_d
        16 -2.53493493E+02 # A_e
        23 -1.80000000E+02 # mu(EWSB)
        25 2.00000000E+00 # tanbeta(in)
        26 3.00000000E+03 # MA_pole
        31 3.00000000E+03 # M_eL
        32 3.00000000E+03 # M_muL
        33 3.00000000E+03 # M_tauL
        34 3.00000000E+03 # M_eR
        35 3.00000000E+03 # M_muR
        36 3.00000000E+03 # M_tauR
        41 3.00000000E+03 # M_q1L
        42 3.00000000E+03 # M_q2L
        43 3.00000000E+03 # M_q3L
        44 3.00000000E+03 # M_uR
        45 3.00000000E+03 # M_cR
        46 3.00000000E+03 # M_tR
        47 3.00000000E+03 # M_dR
        48 3.00000000E+03 # M_sR
        49 3.00000000E+03 # M_bR
#
BLOCK MASS # Mass Spectrum
# PDG code mass particle
        24 8.05138150E+01 # W+
        25 9.64328488E+01 # h
        35 3.00138358E+03 # H
        36 3.00000000E+03 # A
        37 3.00149756E+03 # H+
         5 4.87877839E+00 # b-quark pole mass calculated from mb(mb)_Msbar
   1000001 3.10116707E+03 # ~d_L
   2000001 3.10092490E+03 # ~d_R
   1000002 3.10062316E+03 # ~u_L
   2000002 3.10074630E+03 # ~u_R
   1000003 3.10116707E+03 # ~s_L
   2000003 3.10092490E+03 # ~s_R
   1000004 3.10062316E+03 # ~c_L
   2000004 3.10074630E+03 # ~c_R
   1000005 3.10089484E+03 # ~b_1
   2000005 3.10119898E+03 # ~b_2
   1000006 3.05300709E+03 # ~t_1
   2000006 3.11351467E+03 # ~t_2
   1000011 3.00018880E+03 # ~e_L
   2000011 3.00018459E+03 # ~e_R
   1000012 2.99962657E+03 # ~nu_eL
   1000013 3.00018880E+03 # ~mu_L
   2000013 3.00018459E+03 # ~mu_R
   1000014 2.99962657E+03 # ~nu_muL
   1000015 3.00016190E+03 # ~tau_1
   2000015 3.00021251E+03 # ~tau_2
   1000016 2.99962657E+03 # ~nu_tauL
   1000021 3.18403127E+03 # ~g
   1000022 1.57759136E+02 # ~chi_10
   1000023 -1.86971446E+02 # ~chi_20
   1000025 1.87293685E+02 # ~chi_30
   1000035 3.00191863E+03 # ~chi_40
   1000024 1.81553123E+02 # ~chi_1+
   1000037 3.00191841E+03 # ~chi_2+
#
BLOCK NMIX # Neutralino Mixing Matrix
  1 1 8.57643846E-01 # N_11
  1 2 -1.18004782E-03 # N_12
  1 3 -4.32097729E-01 # N_13
  1 4 -2.78813906E-01 # N_14
  2 1 1.21992538E-01 # N_21
  2 2 -2.31785793E-02 # N_22
  2 3 6.97551610E-01 # N_23
  2 4 -7.05692798E-01 # N_24
  3 1 4.99564515E-01 # N_31
  3 2 8.45994142E-03 # N_32
  3 3 5.71485854E-01 # N_33
  3 4 6.50974380E-01 # N_34
  4 1 3.86730082E-04 # N_41
  4 2 -9.99694849E-01 # N_42
  4 3 -1.08269264E-02 # N_43
  4 4 2.21999497E-02 # N_44
#
BLOCK UMIX # Chargino Mixing Matrix U
  1 1 -1.53045873E-02 # U_11
  1 2 9.99882878E-01 # U_12
  2 1 9.99882878E-01 # U_21
  2 2 1.53045873E-02 # U_22
#
BLOCK VMIX # Chargino Mixing Matrix V
  1 1 3.13948070E-02 # V_11
  1 2 -9.99507062E-01 # V_12
  2 1 9.99507062E-01 # V_21
  2 2 3.13948070E-02 # V_22
#
BLOCK STOPMIX # Stop Mixing Matrix
  1 1 8.65923055E-02 # cos(theta_t)
  1 2 9.96243832E-01 # sin(theta_t)
  2 1 -9.96243832E-01 # -sin(theta_t)
  2 2 8.65923055E-02 # cos(theta_t)
#
BLOCK SBOTMIX # Sbottom Mixing Matrix
  1 1 -3.19175691E-01 # cos(theta_b)
  1 2 9.47695562E-01 # sin(theta_b)
  2 1 -9.47695562E-01 # -sin(theta_b)
  2 2 -3.19175691E-01 # cos(theta_b)
#
BLOCK STAUMIX # Stau Mixing Matrix
  1 1 -6.77051606E-01 # cos(theta_tau)
  1 2 7.35935543E-01 # sin(theta_tau)
  2 1 -7.35935543E-01 # -sin(theta_tau)
  2 2 -6.77051606E-01 # cos(theta_tau)
#
BLOCK ALPHA # Higgs mixing
          -4.89836643E-01 # Mixing angle in the neutral Higgs boson sector
#
BLOCK HMIX Q= 3.00315018E+03 # DRbar Higgs Parameters
         1 -2.02082050E+02 # mu(Q)
         2 1.87933992E+00 # tanbeta(Q)
         3 2.43953322E+02 # vev(Q)
         4 8.94693319E+06 # MA^2(Q)
#
BLOCK GAUGE Q= 3.00315018E+03 # The gauge couplings
     1 3.65027786E-01 # gprime(Q) DRbar
     2 6.37035901E-01 # g(Q) DRbar
     3 1.00557214E+00 # g3(Q) DRbar
#
BLOCK AU Q= 3.00315018E+03 # The trilinear couplings
  1 1 7.01003557E+03 # A_u(Q) DRbar
  2 2 7.01003557E+03 # A_c(Q) DRbar
  3 3 1.06946843E+04 # A_t(Q) DRbar
#
BLOCK AD Q= 3.00315018E+03 # The trilinear couplings
  1 1 3.65113809E+03 # A_d(Q) DRbar
  2 2 3.65113809E+03 # A_s(Q) DRbar
  3 3 5.43866475E+03 # A_b(Q) DRbar
#
BLOCK AE Q= 3.00315018E+03 # The trilinear couplings
  1 1 1.42947879E+03 # A_e(Q) DRbar
  2 2 1.42947879E+03 # A_mu(Q) DRbar
  3 3 1.43139275E+03 # A_tau(Q) DRbar
#
BLOCK Yu Q= 3.00315018E+03 # The Yukawa couplings
  1 1 0.00000000E+00 # y_u(Q) DRbar
  2 2 0.00000000E+00 # y_c(Q) DRbar
  3 3 9.29389411E-01 # y_t(Q) DRbar
#
BLOCK Yd Q= 3.00315018E+03 # The Yukawa couplings
  1 1 0.00000000E+00 # y_d(Q) DRbar
  2 2 0.00000000E+00 # y_s(Q) DRbar
  3 3 2.95519314E-02 # y_b(Q) DRbar
#
BLOCK Ye Q= 3.00315018E+03 # The Yukawa couplings
  1 1 0.00000000E+00 # y_e(Q) DRbar
  2 2 0.00000000E+00 # y_mu(Q) DRbar
  3 3 2.15751919E-02 # y_tau(Q) DRbar
#
BLOCK MSOFT Q= 3.00315018E+03 # The soft SUSY breaking masses at the scale Q
         1 4.10183440E+02 # M_1
         2 3.72724831E+03 # M_2
         3 1.48671796E+03 # M_3
        14 -6.77437438E+02 # A_u
        15 -8.59633345E+02 # A_d
        16 -2.53493493E+02 # A_e
        21 1.74531711E+06 # M^2_Hd
        22 4.34479995E+07 # M^2_Hu
        31 1.95413381E+03 # M_eL
        32 1.95413381E+03 # M_muL
        33 1.95481571E+03 # M_tauL
        34 2.92316702E+03 # M_eR
        35 2.92316702E+03 # M_muR
        36 2.92409003E+03 # M_tauR
        41 5.00000000E+04 # NaN, M_q1L
        42 5.00000000E+04 # NaN i, M_q2L
        43 3.51110209E+03 # M_q3L
        44 1.71058133E+03 # M_uR
        45 1.71058133E+03 # M_cR
        46 5.78133610E+03 # M_tR
        47 1.56382433E+03 # M_dR
        48 1.56382433E+03 # M_sR
        49 1.56611984E+03 # M_bR
#
#
#
# =================
# |The decay table|
# =================
#
# - The QCD corrections to the decays gluino -> squark + quark
# squark -> gaugino + quark_prime
# squark -> squark_prime + Higgs
# squark -> gluino + quark
# are included.
#
# - The multi-body decays for the inos, stops and sbottoms are included.
#
# - The loop induced decays for the gluino, neutralinos and stops
# are included.
#
# - The SUSY decays of the top quark are included.
#
#
# PDG Width
DECAY 6 1.40367633E+00 # top decays
# BR NDA ID1 ID2
     1.00000000E+00 2 5 24 # BR(t -> b W+)
#
# PDG Width
DECAY 1000021 1.65773042E+00 # gluino decays
# BR NDA ID1 ID2
     4.96966483E-02 2 1000001 -1 # BR(~g -> ~d_L db)
     4.96966483E-02 2 -1000001 1 # BR(~g -> ~d_L* d )
     4.99836970E-02 2 2000001 -1 # BR(~g -> ~d_R db)
     4.99836970E-02 2 -2000001 1 # BR(~g -> ~d_R* d )
     5.03424801E-02 2 1000002 -2 # BR(~g -> ~u_L ub)
     5.03424801E-02 2 -1000002 2 # BR(~g -> ~u_L* u )
     5.01959101E-02 2 2000002 -2 # BR(~g -> ~u_R ub)
     5.01959101E-02 2 -2000002 2 # BR(~g -> ~u_R* u )
     4.96966483E-02 2 1000003 -3 # BR(~g -> ~s_L sb)
     4.96966483E-02 2 -1000003 3 # BR(~g -> ~s_L* s )
     4.99836970E-02 2 2000003 -3 # BR(~g -> ~s_R sb)
     4.99836970E-02 2 -2000003 3 # BR(~g -> ~s_R* s )
     5.03424801E-02 2 1000004 -4 # BR(~g -> ~c_L cb)
     5.03424801E-02 2 -1000004 4 # BR(~g -> ~c_L* c )
     5.01959101E-02 2 2000004 -4 # BR(~g -> ~c_R cb)
     5.01959101E-02 2 -2000004 4 # BR(~g -> ~c_R* c )
     4.96053716E-02 2 1000005 -5 # BR(~g -> ~b_1 bb)
     4.96053716E-02 2 -1000005 5 # BR(~g -> ~b_1* b )
     4.99571574E-02 2 2000005 -5 # BR(~g -> ~b_2 bb)
     4.99571574E-02 2 -2000005 5 # BR(~g -> ~b_2* b )
#
# PDG Width
DECAY 1000006 1.07675707E+02 # stop1 decays
# BR NDA ID1 ID2
     9.27998475E-02 2 1000022 6 # BR(~t_1 -> ~chi_10 t )
     2.40628617E-01 2 1000023 6 # BR(~t_1 -> ~chi_20 t )
     2.01346869E-01 2 1000025 6 # BR(~t_1 -> ~chi_30 t )
     4.65222566E-01 2 1000024 5 # BR(~t_1 -> ~chi_1+ b )
     2.10055565E-06 2 1000037 5 # BR(~t_1 -> ~chi_2+ b )
#
# PDG Width
DECAY 2000006 5.28843700E+01 # stop2 decays
# BR NDA ID1 ID2
     6.23780605E-02 2 1000022 6 # BR(~t_2 -> ~chi_10 t )
     4.74077943E-01 2 1000023 6 # BR(~t_2 -> ~chi_20 t )
     4.45501088E-01 2 1000025 6 # BR(~t_2 -> ~chi_30 t )
     1.55673087E-02 2 1000024 5 # BR(~t_2 -> ~chi_1+ b )
     2.47559993E-03 2 1000037 5 # BR(~t_2 -> ~chi_2+ b )
#
# PDG Width
DECAY 1000005 -4.99401042E+00 # sbottom1 decays
# BR NDA ID1 ID2
    -2.86796319E-01 2 1000022 5 # BR(~b_1 -> ~chi_10 b )
    -1.05911380E-02 2 1000023 5 # BR(~b_1 -> ~chi_20 b )
    -1.01439254E-01 2 1000025 5 # BR(~b_1 -> ~chi_30 b )
     1.23406677E-03 2 1000035 5 # BR(~b_1 -> ~chi_40 b )
     1.39759264E+00 2 -1000024 6 # BR(~b_1 -> ~chi_1- t )
#
# PDG Width
DECAY 2000005 5.87087321E+01 # sbottom2 decays
# BR NDA ID1 ID2
     3.53816556E-03 2 1000022 5 # BR(~b_2 -> ~chi_10 b )
     8.47377398E-04 2 1000023 5 # BR(~b_2 -> ~chi_20 b )
     9.53578556E-04 2 1000025 5 # BR(~b_2 -> ~chi_30 b )
     1.00477521E-03 2 1000035 5 # BR(~b_2 -> ~chi_40 b )
     9.93656103E-01 2 -1000024 6 # BR(~b_2 -> ~chi_1- t )
#
# PDG Width
DECAY 1000002 6.28556200E-01 # sup_L decays
# BR NDA ID1 ID2
     5.06691551E-01 2 1000022 2 # BR(~u_L -> ~chi_10 u)
     2.86419157E-07 2 1000023 2 # BR(~u_L -> ~chi_20 u)
     2.06317072E-01 2 1000025 2 # BR(~u_L -> ~chi_30 u)
     8.31075243E-02 2 1000035 2 # BR(~u_L -> ~chi_40 u)
     3.77056873E-02 2 1000024 1 # BR(~u_L -> ~chi_1+ d)
     1.66177879E-01 2 1000037 1 # BR(~u_L -> ~chi_2+ d)
#
# PDG Width
DECAY 2000002 7.02530785E+00 # sup_R decays
# BR NDA ID1 ID2
     7.35930219E-01 2 1000022 2 # BR(~u_R -> ~chi_10 u)
     1.48612823E-02 2 1000023 2 # BR(~u_R -> ~chi_20 u)
     2.49208498E-01 2 1000025 2 # BR(~u_R -> ~chi_30 u)
     6.51163983E-10 2 1000035 2 # BR(~u_R -> ~chi_40 u)
#
# PDG Width
DECAY 1000001 6.08756706E-01 # sdown_L decays
# BR NDA ID1 ID2
     5.38545830E-01 2 1000022 1 # BR(~d_L -> ~chi_10 d)
     4.26563588E-02 2 1000023 1 # BR(~d_L -> ~chi_20 d)
     1.49304882E-01 2 1000025 1 # BR(~d_L -> ~chi_30 d)
     8.67330698E-02 2 1000035 1 # BR(~d_L -> ~chi_40 d)
     9.25336129E-03 2 -1000024 2 # BR(~d_L -> ~chi_1- u)
     1.73506498E-01 2 -1000037 2 # BR(~d_L -> ~chi_2- u)
#
# PDG Width
DECAY 2000001 1.75641265E+00 # sdown_R decays
# BR NDA ID1 ID2
     7.35930171E-01 2 1000022 1 # BR(~d_R -> ~chi_10 d)
     1.48612849E-02 2 1000023 1 # BR(~d_R -> ~chi_20 d)
     2.49208544E-01 2 1000025 1 # BR(~d_R -> ~chi_30 d)
     6.53360504E-10 2 1000035 1 # BR(~d_R -> ~chi_40 d)
#
# PDG Width
DECAY 1000004 6.28556200E-01 # scharm_L decays
# BR NDA ID1 ID2
     5.06691551E-01 2 1000022 4 # BR(~c_L -> ~chi_10 c)
     2.86419157E-07 2 1000023 4 # BR(~c_L -> ~chi_20 c)
     2.06317072E-01 2 1000025 4 # BR(~c_L -> ~chi_30 c)
     8.31075243E-02 2 1000035 4 # BR(~c_L -> ~chi_40 c)
     3.77056873E-02 2 1000024 3 # BR(~c_L -> ~chi_1+ s)
     1.66177879E-01 2 1000037 3 # BR(~c_L -> ~chi_2+ s)
#
# PDG Width
DECAY 2000004 7.02530785E+00 # scharm_R decays
# BR NDA ID1 ID2
     7.35930219E-01 2 1000022 4 # BR(~c_R -> ~chi_10 c)
     1.48612823E-02 2 1000023 4 # BR(~c_R -> ~chi_20 c)
     2.49208498E-01 2 1000025 4 # BR(~c_R -> ~chi_30 c)
     6.51163983E-10 2 1000035 4 # BR(~c_R -> ~chi_40 c)
#
# PDG Width
DECAY 1000003 6.08756706E-01 # sstrange_L decays
# BR NDA ID1 ID2
     5.38545830E-01 2 1000022 3 # BR(~s_L -> ~chi_10 s)
     4.26563588E-02 2 1000023 3 # BR(~s_L -> ~chi_20 s)
     1.49304882E-01 2 1000025 3 # BR(~s_L -> ~chi_30 s)
     8.67330698E-02 2 1000035 3 # BR(~s_L -> ~chi_40 s)
     9.25336129E-03 2 -1000024 4 # BR(~s_L -> ~chi_1- c)
     1.73506498E-01 2 -1000037 4 # BR(~s_L -> ~chi_2- c)
#
# PDG Width
DECAY 2000003 1.75641265E+00 # sstrange_R decays
# BR NDA ID1 ID2
     7.35930171E-01 2 1000022 3 # BR(~s_R -> ~chi_10 s)
     1.48612849E-02 2 1000023 3 # BR(~s_R -> ~chi_20 s)
     2.49208544E-01 2 1000025 3 # BR(~s_R -> ~chi_30 s)
     6.53360504E-10 2 1000035 3 # BR(~s_R -> ~chi_40 s)
#
# PDG Width
DECAY 1000011 3.97042096E+00 # selectron_L decays
# BR NDA ID1 ID2
     7.29096223E-01 2 1000022 11 # BR(~e_L -> ~chi_10 e-)
     6.60763872E-03 2 1000023 11 # BR(~e_L -> ~chi_20 e-)
     2.62877647E-01 2 1000025 11 # BR(~e_L -> ~chi_30 e-)
     1.41849133E-03 2 -1000024 12 # BR(~e_L -> ~chi_1- nu_e)
#
# PDG Width
DECAY 2000011 1.58086360E+01 # selectron_R decays
# BR NDA ID1 ID2
     7.35994474E-01 2 1000022 11 # BR(~e_R -> ~chi_10 e-)
     1.48577184E-02 2 1000023 11 # BR(~e_R -> ~chi_20 e-)
     2.49147808E-01 2 1000025 11 # BR(~e_R -> ~chi_30 e-)
#
# PDG Width
DECAY 1000013 3.97042096E+00 # smuon_L decays
# BR NDA ID1 ID2
     7.29096223E-01 2 1000022 13 # BR(~mu_L -> ~chi_10 mu-)
     6.60763872E-03 2 1000023 13 # BR(~mu_L -> ~chi_20 mu-)
     2.62877647E-01 2 1000025 13 # BR(~mu_L -> ~chi_30 mu-)
     1.41849133E-03 2 -1000024 14 # BR(~mu_L -> ~chi_1- nu_mu)
#
# PDG Width
DECAY 2000013 1.58086360E+01 # smuon_R decays
# BR NDA ID1 ID2
     7.35994474E-01 2 1000022 13 # BR(~mu_R -> ~chi_10 mu-)
     1.48577184E-02 2 1000023 13 # BR(~mu_R -> ~chi_20 mu-)
     2.49147808E-01 2 1000025 13 # BR(~mu_R -> ~chi_30 mu-)
#
# PDG Width
DECAY 1000015 1.04055797E+01 # stau_1 decays
# BR NDA ID1 ID2
     7.45336157E-01 2 1000022 15 # BR(~tau_1 -> ~chi_10 tau-)
     1.10918060E-02 2 1000023 15 # BR(~tau_1 -> ~chi_20 tau-)
     2.43082188E-01 2 1000025 15 # BR(~tau_1 -> ~chi_30 tau-)
     4.89848467E-04 2 -1000024 16 # BR(~tau_1 -> ~chi_1- nu_tau)
#
# PDG Width
DECAY 2000015 9.45620076E+00 # stau_2 decays
# BR NDA ID1 ID2
     7.17471643E-01 2 1000022 15 # BR(~tau_2 -> ~chi_10 tau-)
     1.82454730E-02 2 1000023 15 # BR(~tau_2 -> ~chi_20 tau-)
     2.61310334E-01 2 1000025 15 # BR(~tau_2 -> ~chi_30 tau-)
     2.97255002E-03 2 -1000024 16 # BR(~tau_2 -> ~chi_1- nu_tau)
#
# PDG Width
DECAY 1000012 3.97715079E+00 # snu_eL decays
# BR NDA ID1 ID2
     7.34747964E-01 2 1000022 12 # BR(~nu_eL -> ~chi_10 nu_e)
     2.61739068E-02 2 1000023 12 # BR(~nu_eL -> ~chi_20 nu_e)
     2.33120397E-01 2 1000025 12 # BR(~nu_eL -> ~chi_30 nu_e)
     5.95773227E-03 2 1000024 11 # BR(~nu_eL -> ~chi_1+ e-)
#
# PDG Width
DECAY 1000014 3.97715079E+00 # snu_muL decays
# BR NDA ID1 ID2
     7.34747964E-01 2 1000022 14 # BR(~nu_muL -> ~chi_10 nu_mu)
     2.61739068E-02 2 1000023 14 # BR(~nu_muL -> ~chi_20 nu_mu)
     2.33120397E-01 2 1000025 14 # BR(~nu_muL -> ~chi_30 nu_mu)
     5.95773227E-03 2 1000024 13 # BR(~nu_muL -> ~chi_1+ mu-)
#
# PDG Width
DECAY 1000016 4.00472319E+00 # snu_tauL decays
# BR NDA ID1 ID2
     7.29689246E-01 2 1000022 16 # BR(~nu_tauL -> ~chi_10 nu_tau)
     2.59937003E-02 2 1000023 16 # BR(~nu_tauL -> ~chi_20 nu_tau)
     2.31515370E-01 2 1000025 16 # BR(~nu_tauL -> ~chi_30 nu_tau)
     1.28016837E-02 2 1000024 15 # BR(~nu_tauL -> ~chi_1+ tau-)
#
# PDG Width
DECAY 1000024 2.24131263E-06 # chargino1+ decays
# BR NDA ID1 ID2 ID3
     3.34340277E-01 3 1000022 2 -1 # BR(~chi_1+ -> ~chi_10 u db)
     3.34340277E-01 3 1000022 4 -3 # BR(~chi_1+ -> ~chi_10 c sb)
     1.11449942E-01 3 1000022 -11 12 # BR(~chi_1+ -> ~chi_10 e+ nu_e)
     1.11449942E-01 3 1000022 -13 14 # BR(~chi_1+ -> ~chi_10 mu+ nu_mu)
     1.08419562E-01 3 1000022 -15 16 # BR(~chi_1+ -> ~chi_10 tau+ nu_tau)
#
# PDG Width
DECAY 1000037 2.08548176E+01 # chargino2+ decays
# BR NDA ID1 ID2
     1.35235706E-06 2 1000012 -11 # BR(~chi_2+ -> ~nu_eL e+ )
     1.35235706E-06 2 1000014 -13 # BR(~chi_2+ -> ~nu_muL mu+ )
     8.53477371E-07 2 1000016 -15 # BR(~chi_2+ -> ~nu_tau1 tau+)
     7.70951490E-07 2 -1000011 12 # BR(~chi_2+ -> ~e_L+ nu_e)
     7.70951490E-07 2 -1000013 14 # BR(~chi_2+ -> ~mu_L+ nu_mu)
     3.64889092E-07 2 -1000015 16 # BR(~chi_2+ -> ~tau_1+ nu_tau)
     4.05794781E-07 2 -2000015 16 # BR(~chi_2+ -> ~tau_2+ nu_tau)
     2.52166575E-01 2 1000024 23 # BR(~chi_2+ -> ~chi_1+ Z )
     4.75899014E-02 2 1000022 24 # BR(~chi_2+ -> ~chi_10 W+)
     2.39941161E-01 2 1000023 24 # BR(~chi_2+ -> ~chi_20 W+)
     1.99929769E-01 2 1000025 24 # BR(~chi_2+ -> ~chi_30 W+)
     2.60366723E-01 2 1000024 25 # BR(~chi_2+ -> ~chi_1+ h )
#
# PDG Width
DECAY 1000022 0.00000000E+00 # neutralino1 decays
#
# PDG Width
DECAY 1000023 4.34496172E-06 # neutralino2 decays
# BR NDA ID1 ID2
     2.82404940E-03 2 1000022 22 # BR(~chi_20 -> ~chi_10 gam)
# BR NDA ID1 ID2 ID3
     1.19999008E-01 3 1000022 -2 2 # BR(~chi_20 -> ~chi_10 ub u)
     1.55811510E-01 3 1000022 -1 1 # BR(~chi_20 -> ~chi_10 db d)
     1.19999008E-01 3 1000022 -4 4 # BR(~chi_20 -> ~chi_10 cb c)
     1.55811510E-01 3 1000022 -3 3 # BR(~chi_20 -> ~chi_10 sb s)
     1.21569174E-01 3 1000022 -5 5 # BR(~chi_20 -> ~chi_10 bb b)
     3.58220251E-02 3 1000022 -11 11 # BR(~chi_20 -> ~chi_10 e+ e-)
     3.58220251E-02 3 1000022 -13 13 # BR(~chi_20 -> ~chi_10 mu+ mu-)
     3.50140279E-02 3 1000022 -15 15 # BR(~chi_20 -> ~chi_10 tau+ tau-)
     7.16748171E-02 3 1000022 -12 12 # BR(~chi_20 -> ~chi_10 nu_eb nu_e)
     7.16748171E-02 3 1000022 -14 14 # BR(~chi_20 -> ~chi_10 nu_mub nu_mu)
     7.16748171E-02 3 1000022 -16 16 # BR(~chi_20 -> ~chi_10 nu_taub nu_tau)
     4.01816245E-04 3 1000024 -2 1 # BR(~chi_20 -> ~chi_1+ ub d)
     4.01816245E-04 3 -1000024 -1 2 # BR(~chi_20 -> ~chi_1- db u)
     4.01816245E-04 3 1000024 -4 3 # BR(~chi_20 -> ~chi_1+ cb s)
     4.01816245E-04 3 -1000024 -3 4 # BR(~chi_20 -> ~chi_1- sb c)
     1.33938906E-04 3 1000024 -12 11 # BR(~chi_20 -> ~chi_1+ nu_eb e-)
     1.33938906E-04 3 -1000024 12 -11 # BR(~chi_20 -> ~chi_1- nu_e e+)
     1.33938906E-04 3 1000024 -14 13 # BR(~chi_20 -> ~chi_1+ nu_mub mu-)
     1.33938906E-04 3 -1000024 14 -13 # BR(~chi_20 -> ~chi_1- nu_mu mu+)
     8.00957581E-05 3 1000024 -16 15 # BR(~chi_20 -> ~chi_1+ nu_taub tau-)
     8.00957581E-05 3 -1000024 16 -15 # BR(~chi_20 -> ~chi_1- nu_tau tau+)
#
# PDG Width
DECAY 1000025 2.76133829E-07 # neutralino3 decays
# BR NDA ID1 ID2
     8.62699174E-02 2 1000022 22 # BR(~chi_30 -> ~chi_10 gam)
     1.15050345E-07 2 1000023 22 # BR(~chi_30 -> ~chi_20 gam)
# BR NDA ID1 ID2 ID3
     1.07045214E-01 3 1000022 -2 2 # BR(~chi_30 -> ~chi_10 ub u)
     1.38486191E-01 3 1000022 -1 1 # BR(~chi_30 -> ~chi_10 db d)
     1.07045214E-01 3 1000022 -4 4 # BR(~chi_30 -> ~chi_10 cb c)
     1.38486191E-01 3 1000022 -3 3 # BR(~chi_30 -> ~chi_10 sb s)
     1.01081567E-01 3 1000022 -5 5 # BR(~chi_30 -> ~chi_10 bb b)
     3.14195180E-02 3 1000022 -11 11 # BR(~chi_30 -> ~chi_10 e+ e-)
     3.14195180E-02 3 1000022 -13 13 # BR(~chi_30 -> ~chi_10 mu+ mu-)
     2.98419912E-02 3 1000022 -15 15 # BR(~chi_30 -> ~chi_10 tau+ tau-)
     6.39087444E-02 3 1000022 -12 12 # BR(~chi_30 -> ~chi_10 nu_eb nu_e)
     6.39087444E-02 3 1000022 -14 14 # BR(~chi_30 -> ~chi_10 nu_mub nu_mu)
     6.39087444E-02 3 1000022 -16 16 # BR(~chi_30 -> ~chi_10 nu_taub nu_tau)
     1.09209316E-09 3 1000023 -2 2 # BR(~chi_30 -> ~chi_20 ub u)
     1.41825175E-09 3 1000023 -1 1 # BR(~chi_30 -> ~chi_20 db d)
     1.09209316E-09 3 1000023 -4 4 # BR(~chi_30 -> ~chi_20 cb c)
     1.41825175E-09 3 1000023 -3 3 # BR(~chi_30 -> ~chi_20 sb s)
     3.26241544E-10 3 1000023 -11 11 # BR(~chi_30 -> ~chi_20 e+ e-)
     3.26241544E-10 3 1000023 -13 13 # BR(~chi_30 -> ~chi_20 mu+ mu-)
     6.52273885E-10 3 1000023 -12 12 # BR(~chi_30 -> ~chi_20 nu_eb nu_e)
     6.52273885E-10 3 1000023 -14 14 # BR(~chi_30 -> ~chi_20 nu_mub nu_mu)
     6.52273885E-10 3 1000023 -16 16 # BR(~chi_30 -> ~chi_20 nu_taub nu_tau)
     6.46016754E-03 3 1000024 -2 1 # BR(~chi_30 -> ~chi_1+ ub d)
     6.46016754E-03 3 -1000024 -1 2 # BR(~chi_30 -> ~chi_1- db u)
     6.46016754E-03 3 1000024 -4 3 # BR(~chi_30 -> ~chi_1+ cb s)
     6.46016754E-03 3 -1000024 -3 4 # BR(~chi_30 -> ~chi_1- sb c)
     2.15335595E-03 3 1000024 -12 11 # BR(~chi_30 -> ~chi_1+ nu_eb e-)
     2.15335595E-03 3 -1000024 12 -11 # BR(~chi_30 -> ~chi_1- nu_e e+)
     2.15335595E-03 3 1000024 -14 13 # BR(~chi_30 -> ~chi_1+ nu_mub mu-)
     2.15335595E-03 3 -1000024 14 -13 # BR(~chi_30 -> ~chi_1- nu_mu mu+)
     1.36211410E-03 3 1000024 -16 15 # BR(~chi_30 -> ~chi_1+ nu_taub tau-)
     1.36211410E-03 3 -1000024 16 -15 # BR(~chi_30 -> ~chi_1- nu_tau tau+)
#
# PDG Width
DECAY 1000035 2.08552223E+01 # neutralino4 decays
# BR NDA ID1 ID2
     4.90050391E-02 2 1000022 23 # BR(~chi_40 -> ~chi_10 Z )
     2.72068226E-02 2 1000023 23 # BR(~chi_40 -> ~chi_20 Z )
     1.76145018E-01 2 1000025 23 # BR(~chi_40 -> ~chi_30 Z )
     2.43545311E-01 2 1000024 -24 # BR(~chi_40 -> ~chi_1+ W-)
     2.43545311E-01 2 -1000024 24 # BR(~chi_40 -> ~chi_1- W+)
     2.99385398E-04 2 1000022 25 # BR(~chi_40 -> ~chi_10 h )
     2.28436062E-01 2 1000023 25 # BR(~chi_40 -> ~chi_20 h )
     3.18114482E-02 2 1000025 25 # BR(~chi_40 -> ~chi_30 h )
     3.85250425E-07 2 1000011 -11 # BR(~chi_40 -> ~e_L- e+)
     3.85250425E-07 2 -1000011 11 # BR(~chi_40 -> ~e_L+ e-)
     7.61222252E-14 2 2000011 -11 # BR(~chi_40 -> ~e_R- e+)
     7.61222252E-14 2 -2000011 11 # BR(~chi_40 -> ~e_R+ e-)
     3.85250425E-07 2 1000013 -13 # BR(~chi_40 -> ~mu_L- mu+)
     3.85250425E-07 2 -1000013 13 # BR(~chi_40 -> ~mu_L+ mu-)
     7.61222252E-14 2 2000013 -13 # BR(~chi_40 -> ~mu_R- mu+)
     7.61222252E-14 2 -2000013 13 # BR(~chi_40 -> ~mu_R+ mu-)
     6.76849249E-07 2 1000012 -12 # BR(~chi_40 -> ~nu_eL nu_eb)
     6.76849249E-07 2 -1000012 12 # BR(~chi_40 -> ~nu_eL* nu_e )
     6.76849249E-07 2 1000014 -14 # BR(~chi_40 -> ~nu_muL nu_mub)
     6.76849249E-07 2 -1000014 14 # BR(~chi_40 -> ~nu_muL* nu_mu )
     6.76849249E-07 2 1000016 -16 # BR(~chi_40 -> ~nu_tau1 nu_taub)
     6.76849249E-07 2 -1000016 16 # BR(~chi_40 -> ~nu_tau1* nu_tau )
#
# PDG Width
DECAY 25 2.39614095E-03 # h decays
# BR NDA ID1 ID2
     8.24879747E-01 2 5 -5 # BR(h -> b bb )
     8.35097945E-02 2 -15 15 # BR(h -> tau+ tau- )
     2.95804008E-04 2 -13 13 # BR(h -> mu+ mu- )
     6.60208986E-04 2 3 -3 # BR(h -> s sb )
     2.86435751E-02 2 4 -4 # BR(h -> c cb )
     5.41812023E-02 2 21 21 # BR(h -> g g )
     1.46620721E-03 2 22 22 # BR(h -> gam gam )
     1.10636421E-05 2 22 23 # BR(h -> Z gam )
     5.60830902E-03 2 24 -24 # BR(h -> W+ W- )
     7.44088090E-04 2 23 23 # BR(h -> Z Z )
#
# PDG Width
DECAY 35 3.79901720E+01 # H decays
# BR NDA ID1 ID2
     3.20326848E-03 2 5 -5 # BR(H -> b bb )
     5.77849387E-04 2 -15 15 # BR(H -> tau+ tau- )
     2.04267749E-06 2 -13 13 # BR(H -> mu+ mu- )
     2.48631266E-06 2 3 -3 # BR(H -> s sb )
     9.02012113E-06 2 4 -4 # BR(H -> c cb )
     9.01183981E-01 2 6 -6 # BR(H -> t tb )
     4.86338992E-04 2 21 21 # BR(H -> g g )
     1.59413603E-06 2 22 22 # BR(H -> gam gam )
     7.90409194E-07 2 23 22 # BR(H -> Z gam )
     1.63567345E-04 2 24 -24 # BR(H -> W+ W- )
     8.16839087E-05 2 23 23 # BR(H -> Z Z )
     6.69328860E-04 2 25 25 # BR(H -> h h )
     2.54518355E-25 2 36 36 # BR(H -> A A )
     1.44117724E-13 2 23 36 # BR(H -> Z A )
     2.76416403E-04 2 1000024 -1000024 # BR(H -> ~chi_1+ ~chi_1-)
     3.62920906E-02 2 1000022 1000022 # BR(H -> ~chi_10 ~chi_10)
     4.09130438E-04 2 1000023 1000023 # BR(H -> ~chi_20 ~chi_20)
     2.85281239E-02 2 1000025 1000025 # BR(H -> ~chi_30 ~chi_30)
     2.38377982E-03 2 1000022 1000023 # BR(H -> ~chi_10 ~chi_20)
     1.87870068E-02 2 1000022 1000025 # BR(H -> ~chi_10 ~chi_30)
     6.94150032E-03 2 1000023 1000025 # BR(H -> ~chi_20 ~chi_30)
#
# PDG Width
DECAY 36 3.81607355E+01 # A decays
# BR NDA ID1 ID2
     3.21559483E-03 2 5 -5 # BR(A -> b bb )
     5.75516198E-04 2 -15 15 # BR(A -> tau+ tau- )
     2.03442690E-06 2 -13 13 # BR(A -> mu+ mu- )
     2.52245902E-06 2 3 -3 # BR(A -> s sb )
     9.03008493E-06 2 4 -4 # BR(A -> c cb )
     9.01673078E-01 2 6 -6 # BR(A -> t tb )
     5.78357669E-04 2 21 21 # BR(A -> g g )
     2.02335467E-06 2 22 22 # BR(A -> gam gam )
     9.03039602E-07 2 23 22 # BR(A -> Z gam )
     1.62356839E-04 2 23 25 # BR(A -> Z h )
     8.08700072E-04 2 1000024 -1000024 # BR(A -> ~chi_1+ ~chi_1-)
     8.72437135E-03 2 1000022 1000022 # BR(A -> ~chi_10 ~chi_10)
     4.61677006E-03 2 1000023 1000023 # BR(A -> ~chi_20 ~chi_20)
     1.73324677E-03 2 1000025 1000025 # BR(A -> ~chi_30 ~chi_30)
     5.53915927E-02 2 1000022 1000023 # BR(A -> ~chi_10 ~chi_20)
     2.32196626E-04 2 1000022 1000025 # BR(A -> ~chi_10 ~chi_30)
     2.22717059E-02 2 1000023 1000025 # BR(A -> ~chi_20 ~chi_30)
#
# PDG Width
DECAY 37 3.75313360E+01 # H+ decays
# BR NDA ID1 ID2
     5.31355446E-06 2 4 -5 # BR(H+ -> c bb )
     5.85459697E-04 2 -15 16 # BR(H+ -> tau+ nu_tau )
     2.06957678E-06 2 -13 14 # BR(H+ -> mu+ nu_mu )
     3.39165209E-08 2 2 -5 # BR(H+ -> u bb )
     1.21458282E-07 2 2 -3 # BR(H+ -> u sb )
     1.13055838E-05 2 4 -3 # BR(H+ -> c sb )
     9.04541183E-01 2 6 -5 # BR(H+ -> t bb )
     1.65429190E-04 2 24 25 # BR(H+ -> W+ h )
     2.64026628E-13 2 24 36 # BR(H+ -> W+ A )
     6.77276277E-02 2 1000024 1000022 # BR(H+ -> ~chi_1+ ~chi_10)
     1.24474368E-03 2 1000024 1000023 # BR(H+ -> ~chi_1+ ~chi_20)
     2.57167126E-02 2 1000024 1000025 # BR(H+ -> ~chi_1+ ~chi_30)
#*********************************************************************
# 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
#*********************************************************************
# Run to generate the grid pack *
#*********************************************************************
  False = gridpack !True = setting up the grid pack
#*********************************************************************
# 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
     7000.0 = ebeam1 ! beam 1 total energy in GeV
     7000.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
#*********************************************************************
# Time of flight information. (-1 means not run)
#*********************************************************************
 -1.0 = time_of_flight ! threshold below which info is not written
#*********************************************************************
# Matching - Warning! ickkw > 1 is still beta
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
 1 = highestmult ! for ickkw=2, highest mult group
 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
 1.0 = alpsfact ! scale factor for QCD emission vx
 False = chcluster ! cluster only according to channel diag
 True = pdfwgt ! for ickkw=1, perform pdf reweighting
 4 = asrwgtflavor ! highest quark flavor for a_s reweight
 True = clusinfo ! include clustering tag in output
 3.0 = lhe_version ! Change the way clustering information pass to shower.
#*********************************************************************
#**********************************************************
#
#**********************************************************
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#**********************************************************
   True = auto_ptj_mjj ! Automatic setting of ptj and mjj
#**********************************************************
#
#**********************************
# BW cutoff (M+/-bwcutoff*Gamma)
#**********************************
  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
#*************************************************************
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#*************************************************************
   0 = nhel ! Number of helicities used per event
#*******************
# 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)
 0.0 = ptheavy ! minimum pt for one heavy final state
 -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 = 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
#*********************************************************************
# KT DURHAM CUT *
#*********************************************************************
 -1.0 = ktdurham
 0.4 = dparameter
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
# Jet measure cuts *
#*********************************************************************
 0.0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: If use_syst is T, matched Pythia output is *
# meaningful ONLY if plotted taking matchscale *
# reweighting into account! *
#*********************************************************************
   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
30 50 = 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 # matching scales
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
yuanfang (yuanfang405) said : 		#1

I modfied my susy_hit cards like this:
              41 5.00000000E+04 # M_q1L
              42 5.00000000E+04 # M_q2L
and tried again , still the same problem :
     generate_events
Traceback (most recent call last):
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/extended_cmd.py", line 1011, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/extended_cmd.py", line 966, in onecmd_orig
    return func(arg, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/madevent_interface.py", line 2015, in do_generate_events
    self.exec_cmd('combine_events', postcmd=False)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/extended_cmd.py", line 1037, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/extended_cmd.py", line 966, in onecmd_orig
    return func(arg, **opt)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/madevent_interface.py", line 3008, in do_combine_events
    get_wgt, log_level=logging.DEBUG, write_init=True)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 882, in unweight
    return super(MultiEventFile, self).unweight(outputpath, get_wgt_multi, **opts)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 315, in unweight
    all_wgt, cross, nb_event = self.initialize_unweighting(get_wgt, trunc_error)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 794, in initialize_unweighting
    for event in f:
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 248, in next
    return Event(text)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 910, in __init__
    self.parse(text)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 938, in parse
    self.assign_scale_line(line)
  File "/home/yuanfang/Desktop/working/well_forged_neutralinos/MG5_aMC_v2_4_2/bin/process_2_debug2/bin/internal/lhe_parser.py", line 1327, in assign_scale_line
    self.wgt=float(inputs[2])
ValueError: invalid literal for float(): 0.8386327-100

          I don't know what does this message "ValueError: invalid literal for float(): 0.8386327-100 " mean? Is is a bug ?

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#2

Hi,

This is a typical fortran formatting issue. Some of the computed number are so small that the number written are ill formatted.
Therefore when the python code tries to read it, that fails due to the wrong formatting.
This bug only trigger when one of your channel of integration is strongly suppressed and goes below 1e-98 pb.

Thanks a lot for letting us know this problem. In principle the following patch should solve the problem.

Cheers,

Olivier

=== modified file 'Template/LO/Source/rw_events.f'
--- Template/LO/Source/rw_events.f 2016-08-13 21:44:19 +0000
+++ Template/LO/Source/rw_events.f 2016-09-02 09:21:15 +0000
@@ -179,7 +179,7 @@
 c aqcd = g*g/4d0/pi
       write(largeBuff,'(a)') '<event>'
       evt_record=trim(evt_record)//trim(largeBuff)
- write(largeBuff,'(i2,i4,4e15.7)') nexternal,ievent,wgt,scale,
+ write(largeBuff,'(i2,i4,e15.7e3,3e15.7)') nexternal,ievent,wgt,scale,
      $ aqed,aqcd
       evt_record=trim(evt_record)//CHAR(13)//CHAR(10)//trim(largeBuff)
       do i=1,nexternal
@@ -277,7 +277,8 @@
 c aqcd = g*g/4d0/pi

       write(lun,'(a)') '<event>'
- write(lun,'(i2,i4,4e15.7)') nexternal,ievent,wgt,scale,aqed,aqcd
+ write(lun,'(i2,i4,e15.7e3,3e15.7)') nexternal,ievent,wgt,scale,aqed,aqcd
+ write(*,'(i2,i4,e15.7e3,3e15.7)') nexternal,ievent,wgt,scale,aqed,aqcd
       do i=1,nexternal
          write(lun,51) ic(1,i),ic(6,i),(ic(j,i),j=2,5),
      $ (p(j,i),j=1,3),p(0,i),p(4,i),0.,real(ic(7,i))

Revision history for this message
yuanfang (yuanfang405) said : 		#3

Hi ，
    Thank you Olivier. I am kind of confused . I have found rw_events.f in the directory

          MG5_aMC_v2_4_2/Template/LO/Source/

     but this file only has 221 lines . So I cannot find the line 277 to modify my file. My Madgraph version is 2.4.3 . I downloaded it from https://launchpad.net/mg5amcnlo

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#4

Sorry the patch was indeed created from the 2.5.0 version. which should be released as golden master next week.
Indeed in 2.4.3, you only need to modify one line instead of the twin 2.5.0.
If you want I can create a dedicated patch for 2.4.3.

Cheers,

Olivier

> On Sep 2, 2016, at 15:52, yuanfang <email address hidden> wrote:
>
> Question #376751 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/376751
>
> yuanfang posted a new comment:
> Hi ，
> Thank you Olivier. I am kind of confused . I have found rw_events.f in the directory
>
> MG5_aMC_v2_4_2/Template/LO/Source/
>
> but this file only has 221 lines . So I cannot find the line 277 to modify my file. My Madgraph version is 2.4.3 . I downloaded it from https://launchpad.net/mg5amcnlo
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
yuanfang (yuanfang405) said : 		#5

Thank you Olivier. I have put my susy_hit card in madgraph 2.3.3 and it worked well.I have no idea why it worked in 2.3.3 .
 Since you will release the version 2.5.0 next week ,I can wait for the new version.

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said : 		#6

Actually the patch was not fully correct. (was wrong in presence of non definite positive matrix-element)
Here is the correct one (and for 2.4.3). Note that I expect the bug to be there in 2.3.3 so I would be cautious to use the 2.3.3 results.

Cheers,

Olivier

=== modified file 'Template/LO/Source/rw_events.f'
--- Template/LO/Source/rw_events.f 2014-10-07 13:31:57 +0000
+++ Template/LO/Source/rw_events.f 2016-09-03 13:30:44 +0000
@@ -162,7 +162,7 @@
 c aqcd = g*g/4d0/pi

       write(lun,'(a)') '<event>'
- write(lun,'(i2,i4,4e15.7)') nexternal,ievent,wgt,scale,aqed,aqcd
+ write(lun,'(i2,i4,e16.7e3,3e15.7)') nexternal,ievent,wgt,scale,aqed,aqcd
       do i=1,nexternal
          write(lun,51) ic(1,i),ic(6,i),(ic(j,i),j=2,5),
      $ (p(j,i),j=1,3),p(0,i),p(4,i),0.,real(ic(7,i))

Cheers,

Olivier

Revision history for this message
yuanfang (yuanfang405) said : 		#7

Thanks Olivier Mattelaer, that solved my question.