Memory error

Asked by TMartin

I am getting a bizarre error. I am doing a large batch run of MG/Py/Delphes. The first bunch of runs worked fine, but now I am getting this error message, and it happens for all runs.

  === Results Summary for run: f_5700_700 tag: tag_1 ===

     Cross-section : 9.277e-06 +- 6.774e-09 pb
     Nb of events : 50000

INFO: storring files of Previous run
INFO: Storing Pythia files of Previous run
python(65833,0x7fff79309180) malloc: *** mmap(size=140730350571520) failed (error code=12)
*** error: can't allocate region
*** set a breakpoint in malloc_error_break to debug
Command "generate_events f_5700_700 -f" interrupted with error:
MemoryError :
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/Users/mactop/madgraph/incmsugra/f_5700_700_tag_1_debug.log'.
Please attach this file to your report.
quit
INFO: storring files of Previous run
INFO: Storing Pythia files of Previous run
python(65833,0x7fff79309180) malloc: *** mmap(size=140730350571520) failed (error code=12)
*** error: can't allocate region
*** set a breakpoint in malloc_error_break to debug
INFO:

INFO: storring files of Previous run

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

Do you use the latest version of the code?
Could you put the information of the file pointed in the bug message?

Cheers,

Olivier

Revision history for this message
TMartin (tmartin-4) said :
#2

I am not sure how to attach a file, so I have just copy-pasted it here. I apologize for the length.

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.2.1 *
#* *
#* 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 f_5700_300 -f
Traceback (most recent call last):
  File "/Users/mactop/madgraph/incmsugra/bin/internal/extended_cmd.py", line 879, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/mactop/madgraph/incmsugra/bin/internal/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/Users/mactop/madgraph/incmsugra/bin/internal/madevent_interface.py", line 1947, in do_generate_events
    self.store_result()
  File "/Users/mactop/madgraph/incmsugra/bin/internal/madevent_interface.py", line 3149, in store_result
    stdout=pjoin(p,'%s/%s_pythia_events.hep' % (n,t)))
  File "/Users/mactop/madgraph/incmsugra/bin/internal/misc.py", line 726, in gzip
    ziplib.open(stdout,"w").write(open(path).read())
MemoryError
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
          cluster_temp_path : None
              cluster_queue : madgraph
             cluster_memory : None
               cluster_time : None
                    nb_core : 4 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : /Users/mactop/madgraph/pythia-pgs (user set)
                    td_path : /Users/mactop/madgraph/td (user set)
               delphes_path : /Users/mactop/madgraph/Delphes (user set)
                thepeg_path : None (user set)
               cluster_type : condor
        exrootanalysis_path : /Users/mactop/madgraph/ExRootAnalysis (user set)
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : /Users/mactop/madgraph/MadAnalysis (user set)
                     lhapdf : lhapdf-config
                 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.2.1 2014-09-25 *
#* *
#* 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 complex_mass_scheme False
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
define su = go ul cl t1 ur cr t2 dl sl b1 dr sr b2 ul~ cl~ t1~ ur~ cr~\
 t2~ dl~ sl~ b1~ dr~ sr~ b2~ sve svm svt el- mul- ta1- er- mur- ta2- s\
ve~ svm~ svt~ el+ mul+ ta1+ er+ mur+ ta2+ n1 n2 n3 n4 x1+ x2+ x1- x2-
generate p p > su su
output incmsugra
#
# ======================
# | 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.3b |
# | |
# | 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.22148333E+04
# | |
# ------------------------------------------------------
#
#
BLOCK DCINFO # Decay Program information
     1 SDECAY/HDECAY # decay calculator
     2 1.3b /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 # #SUGRA
#
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
         1 5.70000000E+03 # m0
         2 3.00000000E+02 # m_1
         3 3.00000000E+01 # tanbeta(mZ)
         4 1.00000000E+00 # sign(mu)
         5 -1.14000000E+04 # A0
#
BLOCK EXTPAR # Input parameters - non-minimal models
         0 2.21483325E+03 # EWSB
#
BLOCK MASS # Mass Spectrum
# PDG code mass particle
        24 8.04800953E+01 # W+
        25 1.27961208E+02 # h
        35 4.87098610E+03 # H
        36 4.86983611E+03 # A
        37 4.87084059E+03 # H+
         5 4.87877839E+00 # b-quark pole mass calculated from mb(mb)_Msbar
   1000001 5.66047210E+03 # ~d_L
   2000001 5.68111477E+03 # ~d_R
   1000002 5.65993233E+03 # ~u_L
   2000002 5.67806281E+03 # ~u_R
   1000003 5.66047210E+03 # ~s_L
   2000003 5.68111477E+03 # ~s_R
   1000004 5.65993233E+03 # ~c_L
   2000004 5.67806281E+03 # ~c_R
   1000005 3.41725245E+03 # ~b_1
   2000005 4.60672867E+03 # ~b_2
   1000006 1.51301868E+03 # ~t_1
   2000006 3.45125179E+03 # ~t_2
   1000011 5.67714497E+03 # ~e_L
   2000011 5.69050932E+03 # ~e_R
   1000012 5.67661594E+03 # ~nu_eL
   1000013 5.67714497E+03 # ~mu_L
   2000013 5.69050932E+03 # ~mu_R
   1000014 5.67661594E+03 # ~nu_muL
   1000015 4.68903707E+03 # ~tau_1
   2000015 5.20435079E+03 # ~tau_2
   1000016 5.20299325E+03 # ~nu_tauL
   1000021 9.40365357E+02 # ~g
   1000022 1.42454262E+02 # ~chi_10
   1000023 2.87825946E+02 # ~chi_20
   1000025 -3.62778897E+03 # ~chi_30
   1000035 3.62809441E+03 # ~chi_40
   1000024 2.87825840E+02 # ~chi_1+
   1000037 3.62851117E+03 # ~chi_2+
#
BLOCK NMIX # Neutralino Mixing Matrix
  1 1 9.99924660E-01 # N_11
  1 2 -7.02593152E-04 # N_12
  1 3 1.22218140E-02 # N_13
  1 4 -8.99084694E-04 # N_14
  2 1 9.68186382E-04 # N_21
  2 2 9.99764123E-01 # N_22
  2 3 -2.15587663E-02 # N_23
  2 4 2.44554054E-03 # N_24
  3 1 7.99518188E-03 # N_31
  3 2 -1.35219555E-02 # N_32
  3 3 -7.06919605E-01 # N_33
  3 4 -7.07119443E-01 # N_34
  4 1 9.26360733E-03 # N_41
  4 2 -1.69812111E-02 # N_42
  4 3 -7.06859618E-01 # N_43
  4 4 7.07089319E-01 # N_44
#
BLOCK UMIX # Chargino Mixing Matrix U
  1 1 -9.99534811E-01 # U_11
  1 2 3.04985585E-02 # U_12
  2 1 3.04985585E-02 # U_21
  2 2 9.99534811E-01 # U_22
#
BLOCK VMIX # Chargino Mixing Matrix V
  1 1 -9.99994014E-01 # V_11
  1 2 3.46011023E-03 # V_12
  2 1 3.46011023E-03 # V_21
  2 2 9.99994014E-01 # V_22
#
BLOCK STOPMIX # Stop Mixing Matrix
  1 1 7.13258995E-02 # cos(theta_t)
  1 2 9.97453065E-01 # sin(theta_t)
  2 1 -9.97453065E-01 # -sin(theta_t)
  2 2 7.13258995E-02 # cos(theta_t)
#
BLOCK SBOTMIX # Sbottom Mixing Matrix
  1 1 9.99520446E-01 # cos(theta_b)
  1 2 3.09657557E-02 # sin(theta_b)
  2 1 -3.09657557E-02 # -sin(theta_b)
  2 2 9.99520446E-01 # cos(theta_b)
#
BLOCK STAUMIX # Stau Mixing Matrix
  1 1 3.99058334E-02 # cos(theta_tau)
  1 2 9.99203445E-01 # sin(theta_tau)
  2 1 -9.99203445E-01 # -sin(theta_tau)
  2 2 3.99058334E-02 # cos(theta_tau)
#
BLOCK ALPHA # Higgs mixing
          -3.43894775E-02 # Mixing angle in the neutral Higgs boson sector
#
BLOCK HMIX Q= 2.21483325E+03 # DRbar Higgs Parameters
         1 3.63468568E+03 # mu(Q)
         2 2.90821047E+01 # tanbeta(Q)
         3 2.43301873E+02 # vev(Q)
         4 3.12268517E+07 # MA^2(Q)
#
BLOCK GAUGE Q= 2.21483325E+03 # The gauge couplings
     1 3.63115338E-01 # gprime(Q) DRbar
     2 6.37833109E-01 # g(Q) DRbar
     3 1.02792197E+00 # g3(Q) DRbar
#
BLOCK AU Q= 2.21483325E+03 # The trilinear couplings
  1 1 -8.49183249E+03 # A_u(Q) DRbar
  2 2 -8.49183249E+03 # A_c(Q) DRbar
  3 3 -4.53531747E+03 # A_t(Q) DRbar
#
BLOCK AD Q= 2.21483325E+03 # The trilinear couplings
  1 1 -1.09324143E+04 # A_d(Q) DRbar
  2 2 -1.09324143E+04 # A_s(Q) DRbar
  3 3 -8.73344121E+03 # A_b(Q) DRbar
#
BLOCK AE Q= 2.21483325E+03 # The trilinear couplings
  1 1 -1.02781750E+04 # A_e(Q) DRbar
  2 2 -1.02781750E+04 # A_mu(Q) DRbar
  3 3 -9.41185563E+03 # A_tau(Q) DRbar
#
BLOCK Yu Q= 2.21483325E+03 # The Yukawa couplings
  1 1 0.00000000E+00 # y_u(Q) DRbar
  2 2 0.00000000E+00 # y_c(Q) DRbar
  3 3 8.28318423E-01 # y_t(Q) DRbar
#
BLOCK Yd Q= 2.21483325E+03 # The Yukawa couplings
  1 1 0.00000000E+00 # y_d(Q) DRbar
  2 2 0.00000000E+00 # y_s(Q) DRbar
  3 3 4.29432059E-01 # y_b(Q) DRbar
#
BLOCK Ye Q= 2.21483325E+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.98692486E-01 # y_tau(Q) DRbar
#
BLOCK MSOFT Q= 2.21483325E+03 # The soft SUSY breaking masses at the scale Q
         1 1.41119984E+02 # M_1
         2 2.61323546E+02 # M_2
         3 6.77390641E+02 # M_3
        21 1.03895273E+07 # M^2_Hd
        22 -1.32419076E+07 # M^2_Hu
        31 5.67696619E+03 # M_eL
        32 5.67696619E+03 # M_muL
        33 5.20337538E+03 # M_tauL
        34 5.69033827E+03 # M_eR
        35 5.69033827E+03 # M_muR
        36 4.68969483E+03 # M_tauR
        41 5.61151365E+03 # M_q1L
        42 5.61151365E+03 # M_q2L
        43 3.38742245E+03 # M_q3L
        44 5.62936680E+03 # M_uR
        45 5.62936680E+03 # M_cR
        46 1.45370715E+03 # M_tR
        47 5.63221997E+03 # M_dR
        48 5.63221997E+03 # M_sR
        49 4.56579922E+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.40867016E+00 # top decays
# BR NDA ID1 ID2
     1.00000000E+00 2 5 24 # BR(t -> b W+)
#
# PDG Width
DECAY 1000021 1.79374043E-04 # gluino decays
# BR NDA ID1 ID2
     2.26332830E-03 2 1000022 21 # BR(~g -> ~chi_10 g)
     1.97528935E-04 2 1000023 21 # BR(~g -> ~chi_20 g)
# BR NDA ID1 ID2 ID3
     1.71321308E-03 3 1000022 1 -1 # BR(~g -> ~chi_10 d db)
     6.41845195E-03 3 1000023 1 -1 # BR(~g -> ~chi_20 d db)
     5.81349774E-03 3 1000022 2 -2 # BR(~g -> ~chi_10 u ub)
     6.42565687E-03 3 1000023 2 -2 # BR(~g -> ~chi_20 u ub)
     1.71321308E-03 3 1000022 3 -3 # BR(~g -> ~chi_10 s sb)
     6.41845195E-03 3 1000023 3 -3 # BR(~g -> ~chi_20 s sb)
     5.81349774E-03 3 1000022 4 -4 # BR(~g -> ~chi_10 c cb)
     6.42565687E-03 3 1000023 4 -4 # BR(~g -> ~chi_20 c cb)
     5.94523677E-03 3 1000022 5 -5 # BR(~g -> ~chi_10 b bb)
     5.05334092E-02 3 1000023 5 -5 # BR(~g -> ~chi_20 b bb)
     7.16835211E-01 3 1000022 6 -6 # BR(~g -> ~chi_10 t tb)
     1.37172337E-02 3 1000023 6 -6 # BR(~g -> ~chi_20 t tb)
     1.28441173E-02 3 1000024 1 -2 # BR(~g -> ~chi_1+ d ub)
     1.28441173E-02 3 -1000024 2 -1 # BR(~g -> ~chi_1- u db)
     1.28441173E-02 3 1000024 3 -4 # BR(~g -> ~chi_1+ s cb)
     1.28441173E-02 3 -1000024 4 -3 # BR(~g -> ~chi_1- c sb)
     5.91949717E-02 3 1000024 5 -6 # BR(~g -> ~chi_1+ b tb)
     5.91949717E-02 3 -1000024 6 -5 # BR(~g -> ~chi_1- t bb)
#
# PDG Width
DECAY 1000006 4.78300861E+01 # stop1 decays
# BR NDA ID1 ID2
     6.76852365E-02 2 1000022 6 # BR(~t_1 -> ~chi_10 t )
     6.81028640E-04 2 1000023 6 # BR(~t_1 -> ~chi_20 t )
     1.39463920E-03 2 1000024 5 # BR(~t_1 -> ~chi_1+ b )
     9.30239096E-01 2 1000021 6 # BR(~t_1 -> ~g t )
#
# PDG Width
DECAY 2000006 3.11381596E+02 # stop2 decays
# BR NDA ID1 ID2
     1.65555159E-03 2 1000022 6 # BR(~t_2 -> ~chi_10 t )
     4.24772969E-02 2 1000023 6 # BR(~t_2 -> ~chi_20 t )
     0.00000000E+00 2 1000025 6 # BR(~t_2 -> ~chi_30 t )
     8.53121062E-02 2 1000024 5 # BR(~t_2 -> ~chi_1+ b )
     6.46155959E-01 2 1000021 6 # BR(~t_2 -> ~g t )
     1.10388837E-01 2 1000006 25 # BR(~t_2 -> ~t_1 h )
     1.14010249E-01 2 1000006 23 # BR(~t_2 -> ~t_1 Z )
#
# PDG Width
DECAY 1000005 3.08470070E+02 # sbottom1 decays
# BR NDA ID1 ID2
     1.59911345E-03 2 1000022 5 # BR(~b_1 -> ~chi_10 b )
     4.28634061E-02 2 1000023 5 # BR(~b_1 -> ~chi_20 b )
     8.53606271E-02 2 -1000024 6 # BR(~b_1 -> ~chi_1- t )
     6.49567826E-01 2 1000021 5 # BR(~b_1 -> ~g b )
     2.20609027E-01 2 1000006 -24 # BR(~b_1 -> ~t_1 W-)
#
# PDG Width
DECAY 2000005 2.96310797E+02 # sbottom2 decays
# BR NDA ID1 ID2
     8.75666998E-03 2 1000022 5 # BR(~b_2 -> ~chi_10 b )
     5.71564161E-06 2 1000023 5 # BR(~b_2 -> ~chi_20 b )
     4.04356799E-03 2 1000025 5 # BR(~b_2 -> ~chi_30 b )
     4.04402983E-03 2 1000035 5 # BR(~b_2 -> ~chi_40 b )
     1.12824269E-05 2 -1000024 6 # BR(~b_2 -> ~chi_1- t )
     8.10956847E-03 2 -1000037 6 # BR(~b_2 -> ~chi_2- t )
     9.36598114E-01 2 1000021 5 # BR(~b_2 -> ~g b )
     9.09701231E-03 2 1000005 25 # BR(~b_2 -> ~b_1 h )
     1.04683893E-02 2 1000005 23 # BR(~b_2 -> ~b_1 Z )
     7.42130953E-04 2 1000006 -24 # BR(~b_2 -> ~t_1 W-)
     1.81235195E-02 2 2000006 -24 # BR(~b_2 -> ~t_2 W-)
#
# PDG Width
DECAY 1000002 4.08028060E+02 # sup_L decays
# BR NDA ID1 ID2
     1.94577039E-03 2 1000022 2 # BR(~u_L -> ~chi_10 u)
     5.42288758E-02 2 1000023 2 # BR(~u_L -> ~chi_20 u)
     2.72958695E-06 2 1000025 2 # BR(~u_L -> ~chi_30 u)
     4.38840916E-06 2 1000035 2 # BR(~u_L -> ~chi_40 u)
     1.08467766E-01 2 1000024 1 # BR(~u_L -> ~chi_1+ d)
     4.53274114E-07 2 1000037 1 # BR(~u_L -> ~chi_2+ d)
     8.35350016E-01 2 1000021 2 # BR(~u_L -> ~g u)
#
# PDG Width
DECAY 2000002 3.54901178E+02 # sup_R decays
# BR NDA ID1 ID2
     3.61751968E-02 2 1000022 2 # BR(~u_R -> ~chi_10 u)
     3.37858025E-08 2 1000023 2 # BR(~u_R -> ~chi_20 u)
     8.11705348E-07 2 1000025 2 # BR(~u_R -> ~chi_30 u)
     1.08943597E-06 2 1000035 2 # BR(~u_R -> ~chi_40 u)
     9.63822868E-01 2 1000021 2 # BR(~u_R -> ~g u)
#
# PDG Width
DECAY 1000001 4.08041431E+02 # sdown_L decays
# BR NDA ID1 ID2
     1.97492792E-03 2 1000022 1 # BR(~d_L -> ~chi_10 d)
     5.41924805E-02 2 1000023 1 # BR(~d_L -> ~chi_20 d)
     4.28528686E-06 2 1000025 1 # BR(~d_L -> ~chi_30 d)
     6.65165483E-06 2 1000035 1 # BR(~d_L -> ~chi_40 d)
     1.08375076E-01 2 -1000024 2 # BR(~d_L -> ~chi_1- u)
     3.52274852E-05 2 -1000037 2 # BR(~d_L -> ~chi_2- u)
     8.35411351E-01 2 1000021 1 # BR(~d_L -> ~g d)
#
# PDG Width
DECAY 2000001 3.45478024E+02 # sdown_R decays
# BR NDA ID1 ID2
     9.29548911E-03 2 1000022 1 # BR(~d_R -> ~chi_10 d)
     8.68155184E-09 2 1000023 1 # BR(~d_R -> ~chi_20 d)
     2.08880436E-07 2 1000025 1 # BR(~d_R -> ~chi_30 d)
     2.80350444E-07 2 1000035 1 # BR(~d_R -> ~chi_40 d)
     9.90704013E-01 2 1000021 1 # BR(~d_R -> ~g d)
#
# PDG Width
DECAY 1000004 4.08028060E+02 # scharm_L decays
# BR NDA ID1 ID2
     1.94577039E-03 2 1000022 4 # BR(~c_L -> ~chi_10 c)
     5.42288758E-02 2 1000023 4 # BR(~c_L -> ~chi_20 c)
     2.72958695E-06 2 1000025 4 # BR(~c_L -> ~chi_30 c)
     4.38840916E-06 2 1000035 4 # BR(~c_L -> ~chi_40 c)
     1.08467766E-01 2 1000024 3 # BR(~c_L -> ~chi_1+ s)
     4.53274114E-07 2 1000037 3 # BR(~c_L -> ~chi_2+ s)
     8.35350016E-01 2 1000021 4 # BR(~c_L -> ~g c)
#
# PDG Width
DECAY 2000004 3.54901178E+02 # scharm_R decays
# BR NDA ID1 ID2
     3.61751968E-02 2 1000022 4 # BR(~c_R -> ~chi_10 c)
     3.37858025E-08 2 1000023 4 # BR(~c_R -> ~chi_20 c)
     8.11705348E-07 2 1000025 4 # BR(~c_R -> ~chi_30 c)
     1.08943597E-06 2 1000035 4 # BR(~c_R -> ~chi_40 c)
     9.63822868E-01 2 1000021 4 # BR(~c_R -> ~g c)
#
# PDG Width
DECAY 1000003 4.08041431E+02 # sstrange_L decays
# BR NDA ID1 ID2
     1.97492792E-03 2 1000022 3 # BR(~s_L -> ~chi_10 s)
     5.41924805E-02 2 1000023 3 # BR(~s_L -> ~chi_20 s)
     4.28528686E-06 2 1000025 3 # BR(~s_L -> ~chi_30 s)
     6.65165483E-06 2 1000035 3 # BR(~s_L -> ~chi_40 s)
     1.08375076E-01 2 -1000024 4 # BR(~s_L -> ~chi_1- c)
     3.52274852E-05 2 -1000037 4 # BR(~s_L -> ~chi_2- c)
     8.35411351E-01 2 1000021 3 # BR(~s_L -> ~g s)
#
# PDG Width
DECAY 2000003 3.45478024E+02 # sstrange_R decays
# BR NDA ID1 ID2
     9.29548911E-03 2 1000022 3 # BR(~s_R -> ~chi_10 s)
     8.68155184E-09 2 1000023 3 # BR(~s_R -> ~chi_20 s)
     2.08880436E-07 2 1000025 3 # BR(~s_R -> ~chi_30 s)
     2.80350444E-07 2 1000035 3 # BR(~s_R -> ~chi_40 s)
     9.90704013E-01 2 1000021 3 # BR(~s_R -> ~g s)
#
# PDG Width
DECAY 1000011 7.59750448E+01 # selectron_L decays
# BR NDA ID1 ID2
     9.76254699E-02 2 1000022 11 # BR(~e_L -> ~chi_10 e-)
     3.01031279E-01 2 1000023 11 # BR(~e_L -> ~chi_20 e-)
     8.51789421E-06 2 1000025 11 # BR(~e_L -> ~chi_30 e-)
     1.45057746E-05 2 1000035 11 # BR(~e_L -> ~chi_40 e-)
     6.01123409E-01 2 -1000024 12 # BR(~e_L -> ~chi_1- nu_e)
     1.96817921E-04 2 -1000037 12 # BR(~e_L -> ~chi_2- nu_e)
#
# PDG Width
DECAY 2000011 2.98135622E+01 # selectron_R decays
# BR NDA ID1 ID2
     9.99946245E-01 2 1000022 11 # BR(~e_R -> ~chi_10 e-)
     9.33855212E-07 2 1000023 11 # BR(~e_R -> ~chi_20 e-)
     2.25523945E-05 2 1000025 11 # BR(~e_R -> ~chi_30 e-)
     3.02688613E-05 2 1000035 11 # BR(~e_R -> ~chi_40 e-)
#
# PDG Width
DECAY 1000013 7.59750448E+01 # smuon_L decays
# BR NDA ID1 ID2
     9.76254699E-02 2 1000022 13 # BR(~mu_L -> ~chi_10 mu-)
     3.01031279E-01 2 1000023 13 # BR(~mu_L -> ~chi_20 mu-)
     8.51789421E-06 2 1000025 13 # BR(~mu_L -> ~chi_30 mu-)
     1.45057746E-05 2 1000035 13 # BR(~mu_L -> ~chi_40 mu-)
     6.01123409E-01 2 -1000024 14 # BR(~mu_L -> ~chi_1- nu_mu)
     1.96817921E-04 2 -1000037 14 # BR(~mu_L -> ~chi_2- nu_mu)
#
# PDG Width
DECAY 2000013 2.98135622E+01 # smuon_R decays
# BR NDA ID1 ID2
     9.99946245E-01 2 1000022 13 # BR(~mu_R -> ~chi_10 mu-)
     9.33855212E-07 2 1000023 13 # BR(~mu_R -> ~chi_20 mu-)
     2.25523945E-05 2 1000025 13 # BR(~mu_R -> ~chi_30 mu-)
     3.02688613E-05 2 1000035 13 # BR(~mu_R -> ~chi_40 mu-)
#
# PDG Width
DECAY 1000015 2.73606557E+01 # stau_1 decays
# BR NDA ID1 ID2
     8.96530494E-01 2 1000022 15 # BR(~tau_1 -> ~chi_10 tau-)
     2.02109889E-03 2 1000023 15 # BR(~tau_1 -> ~chi_20 tau-)
     2.44309962E-02 2 1000025 15 # BR(~tau_1 -> ~chi_30 tau-)
     2.43976136E-02 2 1000035 15 # BR(~tau_1 -> ~chi_40 tau-)
     4.03784862E-03 2 -1000024 16 # BR(~tau_1 -> ~chi_1- nu_tau)
     4.85819486E-02 2 -1000037 16 # BR(~tau_1 -> ~chi_2- nu_tau)
#
# PDG Width
DECAY 2000015 7.28468955E+01 # stau_2 decays
# BR NDA ID1 ID2
     9.36756424E-02 2 1000022 15 # BR(~tau_2 -> ~chi_10 tau-)
     2.86803182E-01 2 1000023 15 # BR(~tau_2 -> ~chi_20 tau-)
     1.68038245E-02 2 1000025 15 # BR(~tau_2 -> ~chi_30 tau-)
     1.68157177E-02 2 1000035 15 # BR(~tau_2 -> ~chi_40 tau-)
     5.72594134E-01 2 -1000024 16 # BR(~tau_2 -> ~chi_1- nu_tau)
     3.68946680E-04 2 -1000037 16 # BR(~tau_2 -> ~chi_2- nu_tau)
     6.63655363E-03 2 1000015 25 # BR(~tau_2 -> ~tau_1 h)
     6.30199868E-03 2 1000015 23 # BR(~tau_2 -> ~tau_1 Z)
#
# PDG Width
DECAY 1000012 7.59862944E+01 # snu_eL decays
# BR NDA ID1 ID2
     9.80849456E-02 2 1000022 12 # BR(~nu_eL -> ~chi_10 nu_e)
     3.00295417E-01 2 1000023 12 # BR(~nu_eL -> ~chi_20 nu_e)
     3.45610024E-05 2 1000025 12 # BR(~nu_eL -> ~chi_30 nu_e)
     5.23901982E-05 2 1000035 12 # BR(~nu_eL -> ~chi_40 nu_e)
     6.01530154E-01 2 1000024 11 # BR(~nu_eL -> ~chi_1+ e-)
     2.53203257E-06 2 1000037 11 # BR(~nu_eL -> ~chi_2+ e-)
#
# PDG Width
DECAY 1000014 7.59862944E+01 # snu_muL decays
# BR NDA ID1 ID2
     9.80849456E-02 2 1000022 14 # BR(~nu_muL -> ~chi_10 nu_mu)
     3.00295417E-01 2 1000023 14 # BR(~nu_muL -> ~chi_20 nu_mu)
     3.45610024E-05 2 1000025 14 # BR(~nu_muL -> ~chi_30 nu_mu)
     5.23901982E-05 2 1000035 14 # BR(~nu_muL -> ~chi_40 nu_mu)
     6.01530154E-01 2 1000024 13 # BR(~nu_muL -> ~chi_1+ mu-)
     2.53203257E-06 2 1000037 13 # BR(~nu_muL -> ~chi_2+ mu-)
#
# PDG Width
DECAY 1000016 7.29186747E+01 # snu_tauL decays
# BR NDA ID1 ID2
     9.36609190E-02 2 1000022 16 # BR(~nu_tauL -> ~chi_10 nu_tau)
     2.86538319E-01 2 1000023 16 # BR(~nu_tauL -> ~chi_20 nu_tau)
     2.49043910E-05 2 1000025 16 # BR(~nu_tauL -> ~chi_30 nu_tau)
     3.77487196E-05 2 1000035 16 # BR(~nu_tauL -> ~chi_40 nu_tau)
     5.74089253E-01 2 1000024 15 # BR(~nu_tauL -> ~chi_1+ tau-)
     3.33849062E-02 2 1000037 15 # BR(~nu_tauL -> ~chi_2+ tau-)
     1.22639503E-02 2 -1000015 -24 # BR(~nu_tauL -> ~tau_1+ W-)
#
# PDG Width
DECAY 1000024 5.71056243E-06 # chargino1+ decays
# BR NDA ID1 ID2
     1.00000000E+00 2 1000022 24 # BR(~chi_1+ -> ~chi_10 W+)
#
# PDG Width
DECAY 1000037 8.20069011E+01 # chargino2+ decays
# BR NDA ID1 ID2
     6.98107600E-01 2 1000006 -5 # BR(~chi_2+ -> ~t_1 bb)
     3.00868336E-03 2 2000006 -5 # BR(~chi_2+ -> ~t_2 bb)
     1.08135946E-02 2 -1000005 6 # BR(~chi_2+ -> ~b_1* t )
     8.66768212E-02 2 1000024 23 # BR(~chi_2+ -> ~chi_1+ Z )
     2.70442228E-02 2 1000022 24 # BR(~chi_2+ -> ~chi_10 W+)
     8.41684335E-02 2 1000023 24 # BR(~chi_2+ -> ~chi_20 W+)
     9.01806449E-02 2 1000024 25 # BR(~chi_2+ -> ~chi_1+ h )
#
# PDG Width
DECAY 1000022 0.00000000E+00 # neutralino1 decays
#
# PDG Width
DECAY 1000023 4.27768297E-06 # neutralino2 decays
# BR NDA ID1 ID2
     6.85605737E-02 2 1000022 23 # BR(~chi_20 -> ~chi_10 Z )
     9.31439426E-01 2 1000022 25 # BR(~chi_20 -> ~chi_10 h )
#
# PDG Width
DECAY 1000025 8.07910286E+01 # neutralino3 decays
# BR NDA ID1 ID2
     1.63045694E-02 2 1000022 23 # BR(~chi_30 -> ~chi_10 Z )
     5.37879020E-02 2 1000023 23 # BR(~chi_30 -> ~chi_20 Z )
     8.53591839E-02 2 1000024 -24 # BR(~chi_30 -> ~chi_1+ W-)
     8.53591839E-02 2 -1000024 24 # BR(~chi_30 -> ~chi_1- W+)
     1.26350094E-02 2 1000022 25 # BR(~chi_30 -> ~chi_10 h )
     3.54889300E-02 2 1000023 25 # BR(~chi_30 -> ~chi_20 h )
     3.50755485E-01 2 1000006 -6 # BR(~chi_30 -> ~t_1 tb)
     3.50755485E-01 2 -1000006 6 # BR(~chi_30 -> ~t_1* t )
     2.68145260E-03 2 2000006 -6 # BR(~chi_30 -> ~t_2 tb)
     2.68145260E-03 2 -2000006 6 # BR(~chi_30 -> ~t_2* t )
     2.09567285E-03 2 1000005 -5 # BR(~chi_30 -> ~b_1 bb)
     2.09567285E-03 2 -1000005 5 # BR(~chi_30 -> ~b_1* b )
#
# PDG Width
DECAY 1000035 8.21169019E+01 # neutralino4 decays
# BR NDA ID1 ID2
     1.19492784E-02 2 1000022 23 # BR(~chi_40 -> ~chi_10 Z )
     3.35695424E-02 2 1000023 23 # BR(~chi_40 -> ~chi_20 Z )
     8.39927878E-02 2 1000024 -24 # BR(~chi_40 -> ~chi_1+ W-)
     8.39927878E-02 2 -1000024 24 # BR(~chi_40 -> ~chi_1- W+)
     1.66914274E-02 2 1000022 25 # BR(~chi_40 -> ~chi_10 h )
     5.50707305E-02 2 1000023 25 # BR(~chi_40 -> ~chi_20 h )
     3.53286735E-01 2 1000006 -6 # BR(~chi_40 -> ~t_1 tb)
     3.53286735E-01 2 -1000006 6 # BR(~chi_40 -> ~t_1* t )
     2.01231367E-03 2 2000006 -6 # BR(~chi_40 -> ~t_2 tb)
     2.01231367E-03 2 -2000006 6 # BR(~chi_40 -> ~t_2* t )
     2.06767460E-03 2 1000005 -5 # BR(~chi_40 -> ~b_1 bb)
     2.06767460E-03 2 -1000005 5 # BR(~chi_40 -> ~b_1* b )
#
# PDG Width
DECAY 25 4.73450864E-03 # h decays
# BR NDA ID1 ID2
     5.87699791E-01 2 5 -5 # BR(h -> b bb )
     5.60129060E-02 2 -15 15 # BR(h -> tau+ tau- )
     1.98231805E-04 2 -13 13 # BR(h -> mu+ mu- )
     4.18208944E-04 2 3 -3 # BR(h -> s sb )
     1.81941908E-02 2 4 -4 # BR(h -> c cb )
     6.18922485E-02 2 21 21 # BR(h -> g g )
     2.19702829E-03 2 22 22 # BR(h -> gam gam )
     1.70943300E-03 2 22 23 # BR(h -> Z gam )
     2.39994079E-01 2 24 -24 # BR(h -> W+ W- )
     3.16838825E-02 2 23 23 # BR(h -> Z Z )
#
# PDG Width
DECAY 35 7.66492128E+01 # H decays
# BR NDA ID1 ID2
     6.96217748E-01 2 5 -5 # BR(H -> b bb )
     1.11404477E-01 2 -15 15 # BR(H -> tau+ tau- )
     3.93810441E-04 2 -13 13 # BR(H -> mu+ mu- )
     4.53643619E-04 2 3 -3 # BR(H -> s sb )
     2.84139744E-08 2 4 -4 # BR(H -> c cb )
     2.85710993E-03 2 6 -6 # BR(H -> t tb )
     3.78875419E-06 2 21 21 # BR(H -> g g )
     6.62502692E-08 2 22 22 # BR(H -> gam gam )
     1.93430993E-09 2 23 22 # BR(H -> Z gam )
     1.53246220E-07 2 24 -24 # BR(H -> W+ W- )
     7.65875038E-08 2 23 23 # BR(H -> Z Z )
     1.20418001E-06 2 25 25 # BR(H -> h h )
    -3.53965237E-25 2 36 36 # BR(H -> A A )
     2.83222631E-14 2 23 36 # BR(H -> Z A )
     2.28105596E-18 2 24 -37 # BR(H -> W+ H- )
     2.28105596E-18 2 -24 37 # BR(H -> W- H+ )
     4.87659951E-04 2 1000024 -1000024 # BR(H -> ~chi_1+ ~chi_1-)
     5.22849884E-02 2 1000024 -1000037 # BR(H -> ~chi_1+ ~chi_2-)
     5.22849884E-02 2 1000037 -1000024 # BR(H -> ~chi_2+ ~chi_1-)
     2.26993103E-05 2 1000022 1000022 # BR(H -> ~chi_10 ~chi_10)
     2.43308490E-04 2 1000023 1000023 # BR(H -> ~chi_20 ~chi_20)
     1.48911339E-04 2 1000022 1000023 # BR(H -> ~chi_10 ~chi_20)
     8.85027117E-03 2 1000022 1000025 # BR(H -> ~chi_10 ~chi_30)
     6.32873115E-03 2 1000022 1000035 # BR(H -> ~chi_10 ~chi_40)
     3.30144180E-02 2 1000023 1000025 # BR(H -> ~chi_20 ~chi_30)
     1.91900181E-02 2 1000023 1000035 # BR(H -> ~chi_20 ~chi_40)
     1.58118968E-02 2 1000006 -1000006 # BR(H -> ~t_1 ~t_1* )
#
# PDG Width
DECAY 36 7.54076785E+01 # A decays
# BR NDA ID1 ID2
     7.07555349E-01 2 5 -5 # BR(A -> b bb )
     1.13212136E-01 2 -15 15 # BR(A -> tau+ tau- )
     4.00200232E-04 2 -13 13 # BR(A -> mu+ mu- )
     4.61053357E-04 2 3 -3 # BR(A -> s sb )
     2.91131713E-08 2 4 -4 # BR(A -> c cb )
     2.90946099E-03 2 6 -6 # BR(A -> t tb )
     1.21704367E-05 2 21 21 # BR(A -> g g )
     3.52464677E-09 2 22 22 # BR(A -> gam gam )
     1.26392835E-08 2 23 22 # BR(A -> Z gam )
     1.55427508E-07 2 23 25 # BR(A -> Z h )
     5.10492760E-04 2 1000024 -1000024 # BR(A -> ~chi_1+ ~chi_1-)
     5.30619442E-02 2 1000024 -1000037 # BR(A -> ~chi_1+ ~chi_2-)
     5.30619442E-02 2 1000037 -1000024 # BR(A -> ~chi_2+ ~chi_1-)
     2.33822417E-05 2 1000022 1000022 # BR(A -> ~chi_10 ~chi_10)
     2.54701760E-04 2 1000023 1000023 # BR(A -> ~chi_20 ~chi_20)
     1.54341470E-04 2 1000022 1000023 # BR(A -> ~chi_10 ~chi_20)
     6.42903039E-03 2 1000022 1000025 # BR(A -> ~chi_10 ~chi_30)
     8.97625483E-03 2 1000022 1000035 # BR(A -> ~chi_10 ~chi_40)
     1.94909674E-02 2 1000023 1000025 # BR(A -> ~chi_20 ~chi_30)
     3.34863695E-02 2 1000023 1000035 # BR(A -> ~chi_20 ~chi_40)
#
# PDG Width
DECAY 37 7.04062406E+01 # H+ decays
# BR NDA ID1 ID2
     1.09796723E-03 2 4 -5 # BR(H+ -> c bb )
     1.21279381E-01 2 -15 16 # BR(H+ -> tau+ nu_tau )
     4.28717610E-04 2 -13 14 # BR(H+ -> mu+ nu_mu )
     7.02698677E-06 2 2 -5 # BR(H+ -> u bb )
     2.36655787E-05 2 2 -3 # BR(H+ -> u sb )
     4.86773176E-04 2 4 -3 # BR(H+ -> c sb )
     6.88987818E-01 2 6 -5 # BR(H+ -> t bb )
     1.66610647E-07 2 24 25 # BR(H+ -> W+ h )
     1.91143883E-14 2 24 36 # BR(H+ -> W+ A )
     3.26813759E-04 2 1000024 1000022 # BR(H+ -> ~chi_1+ ~chi_10)
     3.07425958E-10 2 1000024 1000023 # BR(H+ -> ~chi_1+ ~chi_20)
     5.69682162E-02 2 1000024 1000025 # BR(H+ -> ~chi_1+ ~chi_30)
     5.69436981E-02 2 1000024 1000035 # BR(H+ -> ~chi_1+ ~chi_40)
     1.64046200E-02 2 1000037 1000022 # BR(H+ -> ~chi_2+ ~chi_10)
     5.70451348E-02 2 1000037 1000023 # BR(H+ -> ~chi_2+ ~chi_20)
#*********************************************************************
# 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
     4000 = ebeam1 ! beam 1 total energy in GeV
     4000 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
        0 = polbeam1 ! beam polarization for beam 1
        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 *
#*********************************************************************
 F = fixed_ren_scale ! if .true. use fixed ren scale
 F = fixed_fac_scale ! if .true. use fixed fac scale
 91.1880 = scale ! fixed ren scale
 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 1 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# 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 = alpsfact ! scale factor for QCD emission vx
 F = chcluster ! cluster only according to channel diag
 T = pdfwgt ! for ickkw=1, perform pdf reweighting
 5 = asrwgtflavor ! highest quark flavor for a_s reweight
 T = 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)
#**********************************************************
   T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#**********************************************************
#
#**********************************
# BW cutoff (M+/-bwcutoff*Gamma)
#**********************************
  15 = bwcutoff ! (M+/-bwcutoff*Gamma)
#**********************************************************
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#**********************************************************
   T = 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 = ptj ! minimum pt for the jets
  0 = ptb ! minimum pt for the b
 10 = pta ! minimum pt for the photons
 10 = ptl ! minimum pt for the charged leptons
  0 = misset ! minimum missing Et (sum of neutrino's momenta)
  0 = ptheavy ! minimum pt for one heavy final state
 1.0 = ptonium ! minimum pt for the quarkonium states
 -1 = ptjmax ! maximum pt for the jets
 -1 = ptbmax ! maximum pt for the b
 -1 = ptamax ! maximum pt for the photons
 -1 = ptlmax ! maximum pt for the charged leptons
 -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0 = ej ! minimum E for the jets
  0 = eb ! minimum E for the b
  0 = ea ! minimum E for the photons
  0 = el ! minimum E for the charged leptons
 -1 = ejmax ! maximum E for the jets
 -1 = ebmax ! maximum E for the b
 -1 = eamax ! maximum E for the photons
 -1 = elmax ! maximum E for the charged leptons
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
   5 = etaj ! max rap for the jets
  -1 = etab ! max rap for the b
 2.5 = etaa ! max rap for the photons
 2.5 = etal ! max rap for the charged leptons
 0.6 = etaonium ! max rap for the quarkonium states
   0 = etajmin ! min rap for the jets
   0 = etabmin ! min rap for the b
   0 = etaamin ! min rap for the photons
   0 = etalmin ! main rap for the charged leptons
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.4 = drjj ! min distance between jets
 0 = drbb ! min distance between b's
 0.4 = drll ! min distance between leptons
 0.4 = draa ! min distance between gammas
 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 = drab ! min distance between gamma and b
 0 = drbl ! min distance between b and lepton
 0.4 = dral ! min distance between gamma and lepton
 -1 = drjjmax ! max distance between jets
 -1 = drbbmax ! max distance between b's
 -1 = drllmax ! max distance between leptons
 -1 = draamax ! max distance between gammas
 -1 = drbjmax ! max distance between b and jet
 -1 = drajmax ! max distance between gamma and jet
 -1 = drjlmax ! max distance between jet and lepton
 -1 = drabmax ! max distance between gamma and b
 -1 = drblmax ! max distance between b and lepton
 -1 = 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 = mmjj ! min invariant mass of a jet pair
 0 = mmbb ! min invariant mass of a b pair
 0 = mmaa ! min invariant mass of gamma gamma pair
 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1 = mmjjmax ! max invariant mass of a jet pair
 -1 = mmbbmax ! max invariant mass of a b pair
 -1 = mmaamax ! max invariant mass of gamma gamma pair
 -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0 = xptj ! minimum pt for at least one jet
 0 = xptb ! minimum pt for at least one b
 0 = xpta ! minimum pt for at least one photon
 0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0 = ptj1min ! minimum pt for the leading jet in pt
 0 = ptj2min ! minimum pt for the second jet in pt
 0 = ptj3min ! minimum pt for the third jet in pt
 0 = ptj4min ! minimum pt for the fourth jet in pt
 -1 = ptj1max ! maximum pt for the leading jet in pt
 -1 = ptj2max ! maximum pt for the second jet in pt
 -1 = ptj3max ! maximum pt for the third jet in pt
 -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
 0 = ptl2min ! minimum pt for the second lepton in pt
 0 = ptl3min ! minimum pt for the third lepton in pt
 0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1 = ptl1max ! maximum pt for the leading lepton in pt
 -1 = ptl2max ! maximum pt for the second lepton in pt
 -1 = ptl3max ! maximum pt for the third lepton in pt
 -1 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1 = htjmax ! maximum jet HT=Sum(jet pt)
 0 = ihtmin !inclusive Ht for all partons (including b)
 -1 = ihtmax !inclusive Ht for all partons (including b)
 0 = ht2min ! minimum Ht for the two leading jets
 0 = ht3min ! minimum Ht for the three leading jets
 0 = ht4min ! minimum Ht for the four leading jets
 -1 = ht2max ! maximum Ht for the two leading jets
 -1 = ht3max ! maximum Ht for the three leading jets
 -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
 0 = deltaeta ! minimum rapidity for two jets in the WBF case
#*********************************************************************
# KT DURHAM CUT *
#*********************************************************************
 -1 = 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 = 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! *
#*********************************************************************
   F = 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
0.5 1 2 = 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).
CT10nlo.LHgrid = sys_pdf # matching scales
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf

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

Hi,

This is already fix in 2.2.3.

Cheers,

Olivier

Can you help with this problem?

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

To post a message you must log in.