Function not implemented error in generate?

Asked by Radha Mastandrea

I'm running a process using the SMEFTsim_U35_MwScheme_UFO scheme (version 3 from https://github.com/SMEFTsim/SMEFTsim/tree/main/UFO_models/SMEFTsim_U35_MwScheme_UFO, so perhaps the issue is with the model). I'm getting a strage "function not implemented" error with the generate events command. The bug report is given below:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 2023-07-11 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1543, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
    return func(arg, **opt)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2403, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 7656, in new_fct
    original_fct(obj, *args, **opts)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2597, in run_generate_events
    self.exec_cmd('survey %s %s' % (self.run_name,' '.join(args)),
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1572, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
    return func(arg, **opt)
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 3456, in do_survey
    jobs, P_zero_result = ajobcreator.get_helicity()
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/gen_ximprove.py", line 303, in get_helicity
    recycler.generate_output_file()
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 740, in generate_output_file
    self.read_orig()
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 669, in read_orig
    for line_num, line in tqdm(enumerate(input_file), total=get_num_lines(self.input_file)):
  File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 28, in get_num_lines
    buf = mmap.mmap(fp.fileno(),0)
OSError: [Errno 38] Function not implemented
Related File: None
                              Run Options
                              -----------
               stdout_level : 20 (user set)

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

                      Configuration Options
                      ---------------------
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
                thepeg_path : None (user set)
                 hepmc_path : None (user set)
           madanalysis_path : None (user set)
          madanalysis5_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
        exrootanalysis_path : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/ExRootAnalysis (user set)
               syscalc_path : None (user set)
                 rivet_path : None
                  yoda_path : None
                     lhapdf : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                    timeout : 60
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                web_browser : None
                 eps_viewer : None
                text_editor : None
           fortran_compiler : None
               cpp_compiler : None
                auto_update : 7 (user set)
               cluster_type : condor
      cluster_status_update : (600, 30)
           cluster_nb_retry : 1
         cluster_local_path : None
         cluster_retry_wait : 300
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 2023-07-11 *
#* *
#* 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 low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model SMEFTsim_U35_MwScheme_UFO-full
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > w- h, h > b b~, w- > l- vl~
output ./mg_processes/signal1
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.180000e+00 # MB
    6 1.727600e+02 # MT
   23 9.118760e+01 # MZ
   25 1.250900e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  15 0.000000e+00 # ta- : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  9000005 9.118760e+01 # z1 : MZ
  9000006 8.038700e+01 # w1+ : MWsm
  9000007 1.727600e+02 # t1 : MT
  9000008 1.250900e+02 # h1 : MH
  24 8.038700e+01 # w+ : MW

###################################
## INFORMATION FOR SMEFT
###################################
Block smeft
    1 0.000000e+00 # cG
    2 0.000000e+00 # cW
    3 0.000000e+00 # cH
    4 0.000000e+00 # cHbox
    5 0.0 # MadMiner
    6 0.000000e+00 # cHG
    7 0.0 # MadMiner
    8 0.000000e+00 # cHB
    9 0.000000e+00 # cHWB
   10 0.000000e+00 # ceHRe11
   11 0.000000e+00 # ceHRe22
   12 0.000000e+00 # ceHRe33
   19 0.000000e+00 # cuHRe11
   20 0.000000e+00 # cuHRe22
   21 0.000000e+00 # cuHRe33
   28 0.000000e+00 # cdHRe11
   29 0.000000e+00 # cdHRe22
   30 0.000000e+00 # cdHRe33
   37 0.000000e+00 # ceWRe11
   38 0.000000e+00 # ceWRe22
   39 0.000000e+00 # ceWRe33
   46 0.000000e+00 # ceBRe11
   47 0.000000e+00 # ceBRe22
   48 0.000000e+00 # ceBRe33
   55 0.000000e+00 # cuGRe11
   56 0.000000e+00 # cuGRe22
   57 0.000000e+00 # cuGRe33
   64 0.000000e+00 # cuWRe11
   65 0.000000e+00 # cuWRe22
   66 0.000000e+00 # cuWRe33
   73 0.000000e+00 # cuBRe11
   74 0.000000e+00 # cuBRe22
   75 0.000000e+00 # cuBRe33
   82 0.000000e+00 # cdGRe11
   83 0.000000e+00 # cdGRe22
   84 0.000000e+00 # cdGRe33
   91 0.000000e+00 # cdWRe11
   92 0.000000e+00 # cdWRe22
   93 0.000000e+00 # cdWRe33
  100 0.000000e+00 # cdBRe11
  101 0.000000e+00 # cdBRe22
  102 0.000000e+00 # cdBRe33
  109 0.000000e+00 # cHl1Re11
  112 0.000000e+00 # cHl1Re22
  114 0.000000e+00 # cHl1Re33
  115 0.000000e+00 # cHl3Re11
  118 0.000000e+00 # cHl3Re22
  120 0.000000e+00 # cHl3Re33
  121 0.000000e+00 # cHq1Re11
  124 0.000000e+00 # cHq1Re22
  126 0.000000e+00 # cHq1Re33
  127 0.000000e+00 # cHq3Re11
  130 0.000000e+00 # cHq3Re22
  132 0.000000e+00 # cHq3Re33
  133 0.000000e+00 # cHeRe11
  136 0.000000e+00 # cHeRe22
  138 0.000000e+00 # cHeRe33
  139 0.000000e+00 # cHuRe11
  142 0.000000e+00 # cHuRe22
  144 0.000000e+00 # cHuRe33
  145 0.000000e+00 # cHdRe11
  148 0.000000e+00 # cHdRe22
  150 0.000000e+00 # cHdRe33
  151 0.000000e+00 # cHudRe11
  152 0.000000e+00 # cHudRe22
  153 0.000000e+00 # cHudRe33
  160 0.000000e+00 # cllRe1111
  161 0.000000e+00 # cllRe1122
  162 0.000000e+00 # cllRe1221
  163 0.000000e+00 # cllRe1133
  164 0.000000e+00 # cllRe1331
  165 0.000000e+00 # cllRe2222
  166 0.000000e+00 # cllRe2233
  167 0.000000e+00 # cllRe2332
  168 0.000000e+00 # cllRe3333
  172 0.000000e+00 # cllRe1212
  179 0.000000e+00 # cllRe1313
  185 0.000000e+00 # cllRe2323
  187 0.000000e+00 # cqq1Re1111
  188 0.000000e+00 # cqq1Re1122
  189 0.000000e+00 # cqq1Re1221
  190 0.000000e+00 # cqq1Re1133
  191 0.000000e+00 # cqq1Re1331
  192 0.000000e+00 # cqq1Re2222
  193 0.000000e+00 # cqq1Re2233
  194 0.000000e+00 # cqq1Re2332
  195 0.000000e+00 # cqq1Re3333
  199 0.000000e+00 # cqq1Re1212
  206 0.000000e+00 # cqq1Re1313
  212 0.000000e+00 # cqq1Re2323
  214 0.000000e+00 # cqq3Re1111
  215 0.000000e+00 # cqq3Re1122
  216 0.000000e+00 # cqq3Re1221
  217 0.000000e+00 # cqq3Re1133
  218 0.000000e+00 # cqq3Re1331
  219 0.000000e+00 # cqq3Re2222
  220 0.000000e+00 # cqq3Re2233
  221 0.000000e+00 # cqq3Re2332
  222 0.000000e+00 # cqq3Re3333
  226 0.000000e+00 # cqq3Re1212
  233 0.000000e+00 # cqq3Re1313
  239 0.000000e+00 # cqq3Re2323
  241 0.000000e+00 # clq1Re1111
  245 0.000000e+00 # clq1Re1122
  246 0.000000e+00 # clq1Re1133
  265 0.000000e+00 # clq1Re2211
  268 0.000000e+00 # clq1Re2222
  270 0.000000e+00 # clq1Re2233
  280 0.000000e+00 # clq1Re3311
  283 0.000000e+00 # clq1Re3322
  284 0.000000e+00 # clq1Re3333
  286 0.000000e+00 # clq3Re1111
  290 0.000000e+00 # clq3Re1122
  291 0.000000e+00 # clq3Re1133
  310 0.000000e+00 # clq3Re2211
  313 0.000000e+00 # clq3Re2222
  315 0.000000e+00 # clq3Re2233
  325 0.000000e+00 # clq3Re3311
  328 0.000000e+00 # clq3Re3322
  329 0.000000e+00 # clq3Re3333
  331 0.000000e+00 # ceeRe1111
  332 0.000000e+00 # ceeRe1122
  333 0.000000e+00 # ceeRe1133
  334 0.000000e+00 # ceeRe2222
  335 0.000000e+00 # ceeRe2233
  336 0.000000e+00 # ceeRe3333
  341 0.000000e+00 # ceeRe1212
  344 0.000000e+00 # ceeRe1313
  351 0.000000e+00 # ceeRe2323
  352 0.000000e+00 # cuuRe1111
  353 0.000000e+00 # cuuRe1122
  354 0.000000e+00 # cuuRe1221
  355 0.000000e+00 # cuuRe1133
  356 0.000000e+00 # cuuRe1331
  357 0.000000e+00 # cuuRe2222
  358 0.000000e+00 # cuuRe2233
  359 0.000000e+00 # cuuRe2332
  360 0.000000e+00 # cuuRe3333
  364 0.000000e+00 # cuuRe1212
  371 0.000000e+00 # cuuRe1313
  377 0.000000e+00 # cuuRe2323
  379 0.000000e+00 # cddRe1111
  380 0.000000e+00 # cddRe1122
  381 0.000000e+00 # cddRe1221
  382 0.000000e+00 # cddRe1133
  383 0.000000e+00 # cddRe1331
  384 0.000000e+00 # cddRe2222
  385 0.000000e+00 # cddRe2233
  386 0.000000e+00 # cddRe2332
  387 0.000000e+00 # cddRe3333
  391 0.000000e+00 # cddRe1212
  398 0.000000e+00 # cddRe1313
  404 0.000000e+00 # cddRe2323
  406 0.000000e+00 # ceuRe1111
  410 0.000000e+00 # ceuRe1122
  411 0.000000e+00 # ceuRe1133
  430 0.000000e+00 # ceuRe2211
  433 0.000000e+00 # ceuRe2222
  434 0.000000e+00 # ceuRe2223
  435 0.000000e+00 # ceuRe2233
  445 0.000000e+00 # ceuRe3311
  448 0.000000e+00 # ceuRe3322
  449 0.000000e+00 # ceuRe3333
  451 0.000000e+00 # cedRe1111
  455 0.000000e+00 # cedRe1122
  456 0.000000e+00 # cedRe1133
  475 0.000000e+00 # cedRe2211
  478 0.000000e+00 # cedRe2222
  480 0.000000e+00 # cedRe2233
  490 0.000000e+00 # cedRe3311
  493 0.000000e+00 # cedRe3322
  494 0.000000e+00 # cedRe3333
  496 0.000000e+00 # cud1Re1111
  500 0.000000e+00 # cud1Re1122
  501 0.000000e+00 # cud1Re1133
  520 0.000000e+00 # cud1Re2211
  523 0.000000e+00 # cud1Re2222
  525 0.000000e+00 # cud1Re2233
  535 0.000000e+00 # cud1Re3311
  538 0.000000e+00 # cud1Re3322
  539 0.000000e+00 # cud1Re3333
  541 0.000000e+00 # cud8Re1111
  545 0.000000e+00 # cud8Re1122
  546 0.000000e+00 # cud8Re1133
  565 0.000000e+00 # cud8Re2211
  568 0.000000e+00 # cud8Re2222
  569 0.000000e+00 # cud8Re2223
  570 0.000000e+00 # cud8Re2233
  580 0.000000e+00 # cud8Re3311
  583 0.000000e+00 # cud8Re3322
  584 0.000000e+00 # cud8Re3333
  586 0.000000e+00 # cleRe1111
  590 0.000000e+00 # cleRe1122
  591 0.000000e+00 # cleRe1133
  610 0.000000e+00 # cleRe2211
  613 0.000000e+00 # cleRe2222
  615 0.000000e+00 # cleRe2233
  625 0.000000e+00 # cleRe3311
  628 0.000000e+00 # cleRe3322
  629 0.000000e+00 # cleRe3333
  631 0.000000e+00 # cluRe1111
  635 0.000000e+00 # cluRe1122
  636 0.000000e+00 # cluRe1133
  655 0.000000e+00 # cluRe2211
  658 0.000000e+00 # cluRe2222
  660 0.000000e+00 # cluRe2233
  670 0.000000e+00 # cluRe3311
  673 0.000000e+00 # cluRe3322
  674 0.000000e+00 # cluRe3333
  676 0.000000e+00 # cldRe1111
  680 0.000000e+00 # cldRe1122
  681 0.000000e+00 # cldRe1133
  700 0.000000e+00 # cldRe2211
  703 0.000000e+00 # cldRe2222
  705 0.000000e+00 # cldRe2233
  715 0.000000e+00 # cldRe3311
  718 0.000000e+00 # cldRe3322
  719 0.000000e+00 # cldRe3333
  721 0.000000e+00 # cqeRe1111
  725 0.000000e+00 # cqeRe1122
  726 0.000000e+00 # cqeRe1133
  745 0.000000e+00 # cqeRe2211
  748 0.000000e+00 # cqeRe2222
  750 0.000000e+00 # cqeRe2233
  760 0.000000e+00 # cqeRe3311
  763 0.000000e+00 # cqeRe3322
  764 0.000000e+00 # cqeRe3333
  766 0.000000e+00 # cqu1Re1111
  770 0.000000e+00 # cqu1Re1122
  771 0.000000e+00 # cqu1Re1133
  790 0.000000e+00 # cqu1Re2211
  793 0.000000e+00 # cqu1Re2222
  795 0.000000e+00 # cqu1Re2233
  805 0.000000e+00 # cqu1Re3311
  808 0.000000e+00 # cqu1Re3322
  809 0.000000e+00 # cqu1Re3333
  811 0.000000e+00 # cqu8Re1111
  815 0.000000e+00 # cqu8Re1122
  816 0.000000e+00 # cqu8Re1133
  835 0.000000e+00 # cqu8Re2211
  838 0.000000e+00 # cqu8Re2222
  840 0.000000e+00 # cqu8Re2233
  850 0.000000e+00 # cqu8Re3311
  853 0.000000e+00 # cqu8Re3322
  854 0.000000e+00 # cqu8Re3333
  856 0.000000e+00 # cqd1Re1111
  860 0.000000e+00 # cqd1Re1122
  861 0.000000e+00 # cqd1Re1133
  880 0.000000e+00 # cqd1Re2211
  883 0.000000e+00 # cqd1Re2222
  885 0.000000e+00 # cqd1Re2233
  895 0.000000e+00 # cqd1Re3311
  898 0.000000e+00 # cqd1Re3322
  899 0.000000e+00 # cqd1Re3333
  901 0.000000e+00 # cqd8Re1111
  905 0.000000e+00 # cqd8Re1122
  906 0.000000e+00 # cqd8Re1133
  925 0.000000e+00 # cqd8Re2211
  928 0.000000e+00 # cqd8Re2222
  930 0.000000e+00 # cqd8Re2233
  940 0.000000e+00 # cqd8Re3311
  943 0.000000e+00 # cqd8Re3322
  944 0.000000e+00 # cqd8Re3333
  946 0.000000e+00 # cledqRe1111
  950 0.000000e+00 # cledqRe1122
  954 0.000000e+00 # cledqRe1133
  955 0.000000e+00 # cledqRe2211
  959 0.000000e+00 # cledqRe2222
  963 0.000000e+00 # cledqRe2233
  964 0.000000e+00 # cledqRe3311
  968 0.000000e+00 # cledqRe3322
  972 0.000000e+00 # cledqRe3333
  1027 0.000000e+00 # cquqd1Re1111
  1031 0.000000e+00 # cquqd1Re1122
  1035 0.000000e+00 # cquqd1Re1133
  1036 0.000000e+00 # cquqd1Re2211
  1040 0.000000e+00 # cquqd1Re2222
  1044 0.000000e+00 # cquqd1Re2233
  1045 0.000000e+00 # cquqd1Re3311
  1049 0.000000e+00 # cquqd1Re3322
  1053 0.000000e+00 # cquqd1Re3333
  1108 0.000000e+00 # cquqd8Re1111
  1112 0.000000e+00 # cquqd8Re1122
  1116 0.000000e+00 # cquqd8Re1133
  1117 0.000000e+00 # cquqd8Re2211
  1121 0.000000e+00 # cquqd8Re2222
  1125 0.000000e+00 # cquqd8Re2233
  1126 0.000000e+00 # cquqd8Re3311
  1130 0.000000e+00 # cquqd8Re3322
  1134 0.000000e+00 # cquqd8Re3333
  1189 0.000000e+00 # clequ1Re1111
  1193 0.000000e+00 # clequ1Re1122
  1197 0.000000e+00 # clequ1Re1133
  1198 0.000000e+00 # clequ1Re2211
  1202 0.000000e+00 # clequ1Re2222
  1206 0.000000e+00 # clequ1Re2233
  1207 0.000000e+00 # clequ1Re3311
  1211 0.000000e+00 # clequ1Re3322
  1215 0.000000e+00 # clequ1Re3333
  1270 0.000000e+00 # clequ3Re1111
  1274 0.000000e+00 # clequ3Re1122
  1278 0.000000e+00 # clequ3Re1133
  1279 0.000000e+00 # clequ3Re2211
  1283 0.000000e+00 # clequ3Re2222
  1287 0.000000e+00 # clequ3Re2233
  1288 0.000000e+00 # clequ3Re3311
  1292 0.000000e+00 # clequ3Re3322
  1296 0.000000e+00 # clequ3Re3333

###################################
## INFORMATION FOR SMEFTCUTOFF
###################################
Block smeftcutoff
    1 1.000000e+03 # LambdaSMEFT

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 8.038700e+01 # MW
    2 1.166379e-05 # Gf
    3 1.179000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.180000e+00 # ymb
    6 1.727600e+02 # ymt

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.330000e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.085000e+00 # WW
DECAY 25 4.070000e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 9000005 2.495200e+00 # z1 : WZ
DECAY 9000006 2.085000e+00 # w1+ : WW
DECAY 9000007 1.330000e+00 # t1 : WT
DECAY 9000008 4.070000e-03 # h1 : WH
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 9000005 # z1
        1 0 # 3 times electric charge
        2 3 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000006 # w1+
        1 3 # 3 times electric charge
        2 3 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000007 # t1
        1 2 # 3 times electric charge
        2 2 # number of spin states (2S+1)
        3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000008 # h1
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#*********************************************************************

#*******************************
# Running parameters *
#*******************************

#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run

#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  10000 = nevents ! Number of unweighted events requested
      0 = iseed ! rnd seed (0=assigned automatically=default))

#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
  6500.0 = ebeam1 ! beam 1 total energy in GeV
  6500.0 = ebeam2 ! beam 2 total energy in GeV

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
  nn23lo1 = pdlabel ! PDF set
  230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number

#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
   False = fixed_ren_scale ! if .true. use fixed ren scale
   False = fixed_fac_scale ! if .true. use fixed fac scale
  91.188 = scale ! fixed ren scale
  91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
  91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
  -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
   1.0 = scalefact ! scale factor for event-by-event scales

#*********************************************************************
# Type and output format *
#*********************************************************************
  False = gridpack ! True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant lifetime is not written (-1 means not written)
   3.0 = lhe_version ! Change the way clustering information pass to shower.
  True = clusinfo ! include clustering tag in output
   sum = event_norm ! average/sum. Normalization of the weight in the LHEF

#*********************************************************************
# Matching parameter (MLM only) *
#*********************************************************************
  0 = ickkw ! 0 no matching, 1 MLM
  1.0 = alpsfact ! scale factor for QCD emission vx
  False = chcluster ! cluster only according to channel diag
  4 = asrwgtflavor ! highest quark flavor for a_s reweight
  False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 (turn off for VBF and single top processes)
  0.0 = xqcut ! minimum kt jet measure between partons

#*********************************************************************
# handling of the helicities: *
# 0: sum over all helicities *
# 1: importance sampling over helicities *
#*********************************************************************
  0 = nhel ! using helicities importance sampling or not.

#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
  None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
  {} = bias_parameters ! Specifies the parameters of the module.

#*******************************
# Parton level cuts definition *
#*******************************

#*********************************************************************
# BW cutoff. Define on/off-shell for "$" and decay *
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)

#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
   False = cut_decays ! Cut decay products

#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
  20.0 = ptj ! minimum pt for the jets
   0.0 = ptb ! minimum pt for the b
  10.0 = pta ! minimum pt for the photons
  10.0 = ptl ! minimum pt for the charged leptons
   0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
  -1.0 = ptjmax ! maximum pt for the jets
  -1.0 = ptbmax ! maximum pt for the b
  -1.0 = ptamax ! maximum pt for the photons
  -1.0 = ptlmax ! maximum pt for the charged leptons
  -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
   {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
   {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})

#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
   0.0 = ej ! minimum E for the jets
   0.0 = eb ! minimum E for the b
   0.0 = ea ! minimum E for the photons
   0.0 = el ! minimum E for the charged leptons
  -1.0 = ejmax ! maximum E for the jets
  -1.0 = ebmax ! maximum E for the b
  -1.0 = eamax ! maximum E for the photons
  -1.0 = elmax ! maximum E for the charged leptons
    {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
    {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})

#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
   5.0 = etaj ! max rap for the jets
  -1.0 = etab ! max rap for the b
   2.5 = etaa ! max rap for the photons
   2.5 = etal ! max rap for the charged leptons
   0.0 = etajmin ! min rap for the jets
   0.0 = etabmin ! min rap for the b
   0.0 = etaamin ! min rap for the photons
   0.0 = etalmin ! main rap for the charged leptons
    {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
    {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})

#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
   0.4 = drjj ! min distance between jets
   0.0 = drbb ! min distance between b's
   0.4 = drll ! min distance between leptons
   0.4 = draa ! min distance between gammas
   0.0 = drbj ! min distance between b and jet
   0.4 = draj ! min distance between gamma and jet
   0.4 = drjl ! min distance between jet and lepton
   0.0 = drab ! min distance between gamma and b
   0.0 = drbl ! min distance between b and lepton
   0.4 = dral ! min distance between gamma and lepton
  -1.0 = drjjmax ! max distance between jets
  -1.0 = drbbmax ! max distance between b's
  -1.0 = drllmax ! max distance between leptons
  -1.0 = draamax ! max distance between gammas
  -1.0 = drbjmax ! max distance between b and jet
  -1.0 = drajmax ! max distance between gamma and jet
  -1.0 = drjlmax ! max distance between jet and lepton
  -1.0 = drabmax ! max distance between gamma and b
  -1.0 = drblmax ! max distance between b and lepton
  -1.0 = dralmax ! max distance 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
  120.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
  130.0 = mmaamax ! max invariant mass of gamma gamma pair
   -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
     {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
  {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                                              ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.

#*********************************************************************
# Minimum and maximum invariant mass for all leptons *
#*********************************************************************
   0.0 = mmnl ! min invariant mass for all leptons (l+- and vl)
  -1.0 = mmnlmax ! max invariant mass for all leptons (l+- and vl)

#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
   0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
  -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)

#*********************************************************************
# Inclusive cuts *
#*********************************************************************
   0.0 = ptheavy ! minimum pt for at least one heavy final state
   0.0 = xptj ! minimum pt for at least one jet
   0.0 = xptb ! minimum pt for at least one b
   0.0 = xpta ! minimum pt for at least one photon
   0.0 = xptl ! minimum pt for at least one charged lepton

#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
   0.0 = ptj1min ! minimum pt for the leading jet in pt
   0.0 = ptj2min ! minimum pt for the second jet in pt
   0.0 = ptj3min ! minimum pt for the third jet in pt
   0.0 = ptj4min ! minimum pt for the fourth jet in pt
  -1.0 = ptj1max ! maximum pt for the leading jet in pt
  -1.0 = ptj2max ! maximum pt for the second jet in pt
  -1.0 = ptj3max ! maximum pt for the third jet in pt
  -1.0 = ptj4max ! maximum pt for the fourth jet in pt
   0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts

#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
   0.0 = ptl1min ! minimum pt for the leading lepton in pt
   0.0 = ptl2min ! minimum pt for the second lepton in pt
   0.0 = ptl3min ! minimum pt for the third lepton in pt
   0.0 = ptl4min ! minimum pt for the fourth lepton in pt
  -1.0 = ptl1max ! maximum pt for the leading lepton in pt
  -1.0 = ptl2max ! maximum pt for the second lepton in pt
  -1.0 = ptl3max ! maximum pt for the third lepton in pt
  -1.0 = ptl4max ! maximum pt for the fourth lepton in pt

#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
   0.0 = htjmin ! minimum jet HT=Sum(jet pt)
  -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
   0.0 = ihtmin ! inclusive Ht for all partons (including b)
  -1.0 = ihtmax ! inclusive Ht for all partons (including b)
   0.0 = ht2min ! minimum Ht for the two leading jets
   0.0 = ht3min ! minimum Ht for the three leading jets
   0.0 = ht4min ! minimum Ht for the four leading jets
  -1.0 = ht2max ! maximum Ht for the two leading jets
  -1.0 = ht3max ! maximum Ht for the three leading jets
  -1.0 = ht4max ! maximum Ht for the four leading jets

#*********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#*********************************************************************
   0.0 = ptgmin ! Min photon transverse momentum
   0.4 = R0gamma ! Radius of isolation code
   1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
   1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
  True = isoEM ! isolate photons from EM energy (photons and leptons)

#*********************************************************************
# WBF cuts *
#*********************************************************************
   0.0 = xetamin ! minimum rapidity for two jets in the WBF case
   0.0 = deltaeta ! minimum rapidity for two jets in the WBF case

#*********************************************************************
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#*********************************************************************
  -1.0 = ktdurham
   0.4 = dparameter
  -1.0 = ptlund
   1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above

#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
  4 = maxjetflavor ! Maximum jet pdg code

#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
# False = use_syst ! Enable systematics studies # Commented out by MadMiner
# systematics = systematics_program ! none, systematics [python], SysCalc [deprecated, C++] # Commented out by MadMiner
# ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule # Commented out by MadMiner

#*********************************************************************
# MadMiner systematics setup *
#*********************************************************************
False = use_syst

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

No this has nothing to do with the model, it is related to the os python library which does system call operation.
So looks like your OS does not support one of the function of this line:
  buf = mmap.mmap(fp.fileno(),0)

So which OS/kernel are you using?

Since this seems related to the helicity recycling optimization,
the quick work-around is to de-activate such optimization.
(one way is to set hel_recycling to false in the run_card -- you can add the corresponding line in the file if needed--)

Cheers,

Olivier

> On 17 Aug 2023, at 04:10, Radha Mastandrea <email address hidden> wrote:
>
> New question #707663 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/707663
>
> I'm running a process using the SMEFTsim_U35_MwScheme_UFO scheme (version 3 from https://github.com/SMEFTsim/SMEFTsim/tree/main/UFO_models/SMEFTsim_U35_MwScheme_UFO, so perhaps the issue is with the model). I'm getting a strage "function not implemented" error with the generate events command. The bug report is given below:
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.5.1 2023-07-11 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #************************************************************
> generate_events run_01
> Traceback (most recent call last):
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1543, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
> return func(arg, **opt)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2403, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 7656, in new_fct
> original_fct(obj, *args, **opts)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2597, in run_generate_events
> self.exec_cmd('survey %s %s' % (self.run_name,' '.join(args)),
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1572, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
> return func(arg, **opt)
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 3456, in do_survey
> jobs, P_zero_result = ajobcreator.get_helicity()
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/gen_ximprove.py", line 303, in get_helicity
> recycler.generate_output_file()
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 740, in generate_output_file
> self.read_orig()
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 669, in read_orig
> for line_num, line in tqdm(enumerate(input_file), total=get_num_lines(self.input_file)):
> File "/global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/madgraph/madevent/hel_recycle.py", line 28, in get_num_lines
> buf = mmap.mmap(fp.fileno(),0)
> OSError: [Errno 38] Function not implemented
> Related File: None
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> run_mode : 2
> cluster_queue : None (user set)
> cluster_time : None (user set)
> cluster_size : 100
> cluster_memory : 100 (user set)
> nb_core : 256 (user set)
> cluster_temp_path : None
>
> Configuration Options
> ---------------------
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> thepeg_path : None (user set)
> hepmc_path : None (user set)
> madanalysis_path : None (user set)
> madanalysis5_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> exrootanalysis_path : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/ExRootAnalysis (user set)
> syscalc_path : None (user set)
> rivet_path : None
> yoda_path : None
> lhapdf : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
> lhapdf_py2 : None
> lhapdf_py3 : /global/u1/r/rmastand/HEP_installs/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
> timeout : 60
> f2py_compiler : None
> f2py_compiler_py2 : None
> f2py_compiler_py3 : None
> web_browser : None
> eps_viewer : None
> text_editor : None
> fortran_compiler : None
> cpp_compiler : None
> auto_update : 7 (user set)
> cluster_type : condor
> cluster_status_update : (600, 30)
> cluster_nb_retry : 1
> cluster_local_path : None
> cluster_retry_wait : 300
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.5.1 2023-07-11 *
> #* *
> #* 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 low_mem_multicore_nlo_generation False
> set complex_mass_scheme False
> set include_lepton_initiated_processes False
> set gauge unitary
> set loop_optimized_output True
> set loop_color_flows False
> set max_npoint_for_channel 0
> set default_unset_couplings 99
> set max_t_for_channel 99
> set zerowidth_tchannel True
> set nlo_mixed_expansion True
> import model SMEFTsim_U35_MwScheme_UFO-full
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> generate p p > w- h, h > b b~, w- > l- vl~
> output ./mg_processes/signal1
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.180000e+00 # MB
> 6 1.727600e+02 # MT
> 23 9.118760e+01 # MZ
> 25 1.250900e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 15 0.000000e+00 # ta- : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 9000005 9.118760e+01 # z1 : MZ
> 9000006 8.038700e+01 # w1+ : MWsm
> 9000007 1.727600e+02 # t1 : MT
> 9000008 1.250900e+02 # h1 : MH
> 24 8.038700e+01 # w+ : MW
>
> ###################################
> ## INFORMATION FOR SMEFT
> ###################################
> Block smeft
> 1 0.000000e+00 # cG
> 2 0.000000e+00 # cW
> 3 0.000000e+00 # cH
> 4 0.000000e+00 # cHbox
> 5 0.0 # MadMiner
> 6 0.000000e+00 # cHG
> 7 0.0 # MadMiner
> 8 0.000000e+00 # cHB
> 9 0.000000e+00 # cHWB
> 10 0.000000e+00 # ceHRe11
> 11 0.000000e+00 # ceHRe22
> 12 0.000000e+00 # ceHRe33
> 19 0.000000e+00 # cuHRe11
> 20 0.000000e+00 # cuHRe22
> 21 0.000000e+00 # cuHRe33
> 28 0.000000e+00 # cdHRe11
> 29 0.000000e+00 # cdHRe22
> 30 0.000000e+00 # cdHRe33
> 37 0.000000e+00 # ceWRe11
> 38 0.000000e+00 # ceWRe22
> 39 0.000000e+00 # ceWRe33
> 46 0.000000e+00 # ceBRe11
> 47 0.000000e+00 # ceBRe22
> 48 0.000000e+00 # ceBRe33
> 55 0.000000e+00 # cuGRe11
> 56 0.000000e+00 # cuGRe22
> 57 0.000000e+00 # cuGRe33
> 64 0.000000e+00 # cuWRe11
> 65 0.000000e+00 # cuWRe22
> 66 0.000000e+00 # cuWRe33
> 73 0.000000e+00 # cuBRe11
> 74 0.000000e+00 # cuBRe22
> 75 0.000000e+00 # cuBRe33
> 82 0.000000e+00 # cdGRe11
> 83 0.000000e+00 # cdGRe22
> 84 0.000000e+00 # cdGRe33
> 91 0.000000e+00 # cdWRe11
> 92 0.000000e+00 # cdWRe22
> 93 0.000000e+00 # cdWRe33
> 100 0.000000e+00 # cdBRe11
> 101 0.000000e+00 # cdBRe22
> 102 0.000000e+00 # cdBRe33
> 109 0.000000e+00 # cHl1Re11
> 112 0.000000e+00 # cHl1Re22
> 114 0.000000e+00 # cHl1Re33
> 115 0.000000e+00 # cHl3Re11
> 118 0.000000e+00 # cHl3Re22
> 120 0.000000e+00 # cHl3Re33
> 121 0.000000e+00 # cHq1Re11
> 124 0.000000e+00 # cHq1Re22
> 126 0.000000e+00 # cHq1Re33
> 127 0.000000e+00 # cHq3Re11
> 130 0.000000e+00 # cHq3Re22
> 132 0.000000e+00 # cHq3Re33
> 133 0.000000e+00 # cHeRe11
> 136 0.000000e+00 # cHeRe22
> 138 0.000000e+00 # cHeRe33
> 139 0.000000e+00 # cHuRe11
> 142 0.000000e+00 # cHuRe22
> 144 0.000000e+00 # cHuRe33
> 145 0.000000e+00 # cHdRe11
> 148 0.000000e+00 # cHdRe22
> 150 0.000000e+00 # cHdRe33
> 151 0.000000e+00 # cHudRe11
> 152 0.000000e+00 # cHudRe22
> 153 0.000000e+00 # cHudRe33
> 160 0.000000e+00 # cllRe1111
> 161 0.000000e+00 # cllRe1122
> 162 0.000000e+00 # cllRe1221
> 163 0.000000e+00 # cllRe1133
> 164 0.000000e+00 # cllRe1331
> 165 0.000000e+00 # cllRe2222
> 166 0.000000e+00 # cllRe2233
> 167 0.000000e+00 # cllRe2332
> 168 0.000000e+00 # cllRe3333
> 172 0.000000e+00 # cllRe1212
> 179 0.000000e+00 # cllRe1313
> 185 0.000000e+00 # cllRe2323
> 187 0.000000e+00 # cqq1Re1111
> 188 0.000000e+00 # cqq1Re1122
> 189 0.000000e+00 # cqq1Re1221
> 190 0.000000e+00 # cqq1Re1133
> 191 0.000000e+00 # cqq1Re1331
> 192 0.000000e+00 # cqq1Re2222
> 193 0.000000e+00 # cqq1Re2233
> 194 0.000000e+00 # cqq1Re2332
> 195 0.000000e+00 # cqq1Re3333
> 199 0.000000e+00 # cqq1Re1212
> 206 0.000000e+00 # cqq1Re1313
> 212 0.000000e+00 # cqq1Re2323
> 214 0.000000e+00 # cqq3Re1111
> 215 0.000000e+00 # cqq3Re1122
> 216 0.000000e+00 # cqq3Re1221
> 217 0.000000e+00 # cqq3Re1133
> 218 0.000000e+00 # cqq3Re1331
> 219 0.000000e+00 # cqq3Re2222
> 220 0.000000e+00 # cqq3Re2233
> 221 0.000000e+00 # cqq3Re2332
> 222 0.000000e+00 # cqq3Re3333
> 226 0.000000e+00 # cqq3Re1212
> 233 0.000000e+00 # cqq3Re1313
> 239 0.000000e+00 # cqq3Re2323
> 241 0.000000e+00 # clq1Re1111
> 245 0.000000e+00 # clq1Re1122
> 246 0.000000e+00 # clq1Re1133
> 265 0.000000e+00 # clq1Re2211
> 268 0.000000e+00 # clq1Re2222
> 270 0.000000e+00 # clq1Re2233
> 280 0.000000e+00 # clq1Re3311
> 283 0.000000e+00 # clq1Re3322
> 284 0.000000e+00 # clq1Re3333
> 286 0.000000e+00 # clq3Re1111
> 290 0.000000e+00 # clq3Re1122
> 291 0.000000e+00 # clq3Re1133
> 310 0.000000e+00 # clq3Re2211
> 313 0.000000e+00 # clq3Re2222
> 315 0.000000e+00 # clq3Re2233
> 325 0.000000e+00 # clq3Re3311
> 328 0.000000e+00 # clq3Re3322
> 329 0.000000e+00 # clq3Re3333
> 331 0.000000e+00 # ceeRe1111
> 332 0.000000e+00 # ceeRe1122
> 333 0.000000e+00 # ceeRe1133
> 334 0.000000e+00 # ceeRe2222
> 335 0.000000e+00 # ceeRe2233
> 336 0.000000e+00 # ceeRe3333
> 341 0.000000e+00 # ceeRe1212
> 344 0.000000e+00 # ceeRe1313
> 351 0.000000e+00 # ceeRe2323
> 352 0.000000e+00 # cuuRe1111
> 353 0.000000e+00 # cuuRe1122
> 354 0.000000e+00 # cuuRe1221
> 355 0.000000e+00 # cuuRe1133
> 356 0.000000e+00 # cuuRe1331
> 357 0.000000e+00 # cuuRe2222
> 358 0.000000e+00 # cuuRe2233
> 359 0.000000e+00 # cuuRe2332
> 360 0.000000e+00 # cuuRe3333
> 364 0.000000e+00 # cuuRe1212
> 371 0.000000e+00 # cuuRe1313
> 377 0.000000e+00 # cuuRe2323
> 379 0.000000e+00 # cddRe1111
> 380 0.000000e+00 # cddRe1122
> 381 0.000000e+00 # cddRe1221
> 382 0.000000e+00 # cddRe1133
> 383 0.000000e+00 # cddRe1331
> 384 0.000000e+00 # cddRe2222
> 385 0.000000e+00 # cddRe2233
> 386 0.000000e+00 # cddRe2332
> 387 0.000000e+00 # cddRe3333
> 391 0.000000e+00 # cddRe1212
> 398 0.000000e+00 # cddRe1313
> 404 0.000000e+00 # cddRe2323
> 406 0.000000e+00 # ceuRe1111
> 410 0.000000e+00 # ceuRe1122
> 411 0.000000e+00 # ceuRe1133
> 430 0.000000e+00 # ceuRe2211
> 433 0.000000e+00 # ceuRe2222
> 434 0.000000e+00 # ceuRe2223
> 435 0.000000e+00 # ceuRe2233
> 445 0.000000e+00 # ceuRe3311
> 448 0.000000e+00 # ceuRe3322
> 449 0.000000e+00 # ceuRe3333
> 451 0.000000e+00 # cedRe1111
> 455 0.000000e+00 # cedRe1122
> 456 0.000000e+00 # cedRe1133
> 475 0.000000e+00 # cedRe2211
> 478 0.000000e+00 # cedRe2222
> 480 0.000000e+00 # cedRe2233
> 490 0.000000e+00 # cedRe3311
> 493 0.000000e+00 # cedRe3322
> 494 0.000000e+00 # cedRe3333
> 496 0.000000e+00 # cud1Re1111
> 500 0.000000e+00 # cud1Re1122
> 501 0.000000e+00 # cud1Re1133
> 520 0.000000e+00 # cud1Re2211
> 523 0.000000e+00 # cud1Re2222
> 525 0.000000e+00 # cud1Re2233
> 535 0.000000e+00 # cud1Re3311
> 538 0.000000e+00 # cud1Re3322
> 539 0.000000e+00 # cud1Re3333
> 541 0.000000e+00 # cud8Re1111
> 545 0.000000e+00 # cud8Re1122
> 546 0.000000e+00 # cud8Re1133
> 565 0.000000e+00 # cud8Re2211
> 568 0.000000e+00 # cud8Re2222
> 569 0.000000e+00 # cud8Re2223
> 570 0.000000e+00 # cud8Re2233
> 580 0.000000e+00 # cud8Re3311
> 583 0.000000e+00 # cud8Re3322
> 584 0.000000e+00 # cud8Re3333
> 586 0.000000e+00 # cleRe1111
> 590 0.000000e+00 # cleRe1122
> 591 0.000000e+00 # cleRe1133
> 610 0.000000e+00 # cleRe2211
> 613 0.000000e+00 # cleRe2222
> 615 0.000000e+00 # cleRe2233
> 625 0.000000e+00 # cleRe3311
> 628 0.000000e+00 # cleRe3322
> 629 0.000000e+00 # cleRe3333
> 631 0.000000e+00 # cluRe1111
> 635 0.000000e+00 # cluRe1122
> 636 0.000000e+00 # cluRe1133
> 655 0.000000e+00 # cluRe2211
> 658 0.000000e+00 # cluRe2222
> 660 0.000000e+00 # cluRe2233
> 670 0.000000e+00 # cluRe3311
> 673 0.000000e+00 # cluRe3322
> 674 0.000000e+00 # cluRe3333
> 676 0.000000e+00 # cldRe1111
> 680 0.000000e+00 # cldRe1122
> 681 0.000000e+00 # cldRe1133
> 700 0.000000e+00 # cldRe2211
> 703 0.000000e+00 # cldRe2222
> 705 0.000000e+00 # cldRe2233
> 715 0.000000e+00 # cldRe3311
> 718 0.000000e+00 # cldRe3322
> 719 0.000000e+00 # cldRe3333
> 721 0.000000e+00 # cqeRe1111
> 725 0.000000e+00 # cqeRe1122
> 726 0.000000e+00 # cqeRe1133
> 745 0.000000e+00 # cqeRe2211
> 748 0.000000e+00 # cqeRe2222
> 750 0.000000e+00 # cqeRe2233
> 760 0.000000e+00 # cqeRe3311
> 763 0.000000e+00 # cqeRe3322
> 764 0.000000e+00 # cqeRe3333
> 766 0.000000e+00 # cqu1Re1111
> 770 0.000000e+00 # cqu1Re1122
> 771 0.000000e+00 # cqu1Re1133
> 790 0.000000e+00 # cqu1Re2211
> 793 0.000000e+00 # cqu1Re2222
> 795 0.000000e+00 # cqu1Re2233
> 805 0.000000e+00 # cqu1Re3311
> 808 0.000000e+00 # cqu1Re3322
> 809 0.000000e+00 # cqu1Re3333
> 811 0.000000e+00 # cqu8Re1111
> 815 0.000000e+00 # cqu8Re1122
> 816 0.000000e+00 # cqu8Re1133
> 835 0.000000e+00 # cqu8Re2211
> 838 0.000000e+00 # cqu8Re2222
> 840 0.000000e+00 # cqu8Re2233
> 850 0.000000e+00 # cqu8Re3311
> 853 0.000000e+00 # cqu8Re3322
> 854 0.000000e+00 # cqu8Re3333
> 856 0.000000e+00 # cqd1Re1111
> 860 0.000000e+00 # cqd1Re1122
> 861 0.000000e+00 # cqd1Re1133
> 880 0.000000e+00 # cqd1Re2211
> 883 0.000000e+00 # cqd1Re2222
> 885 0.000000e+00 # cqd1Re2233
> 895 0.000000e+00 # cqd1Re3311
> 898 0.000000e+00 # cqd1Re3322
> 899 0.000000e+00 # cqd1Re3333
> 901 0.000000e+00 # cqd8Re1111
> 905 0.000000e+00 # cqd8Re1122
> 906 0.000000e+00 # cqd8Re1133
> 925 0.000000e+00 # cqd8Re2211
> 928 0.000000e+00 # cqd8Re2222
> 930 0.000000e+00 # cqd8Re2233
> 940 0.000000e+00 # cqd8Re3311
> 943 0.000000e+00 # cqd8Re3322
> 944 0.000000e+00 # cqd8Re3333
> 946 0.000000e+00 # cledqRe1111
> 950 0.000000e+00 # cledqRe1122
> 954 0.000000e+00 # cledqRe1133
> 955 0.000000e+00 # cledqRe2211
> 959 0.000000e+00 # cledqRe2222
> 963 0.000000e+00 # cledqRe2233
> 964 0.000000e+00 # cledqRe3311
> 968 0.000000e+00 # cledqRe3322
> 972 0.000000e+00 # cledqRe3333
> 1027 0.000000e+00 # cquqd1Re1111
> 1031 0.000000e+00 # cquqd1Re1122
> 1035 0.000000e+00 # cquqd1Re1133
> 1036 0.000000e+00 # cquqd1Re2211
> 1040 0.000000e+00 # cquqd1Re2222
> 1044 0.000000e+00 # cquqd1Re2233
> 1045 0.000000e+00 # cquqd1Re3311
> 1049 0.000000e+00 # cquqd1Re3322
> 1053 0.000000e+00 # cquqd1Re3333
> 1108 0.000000e+00 # cquqd8Re1111
> 1112 0.000000e+00 # cquqd8Re1122
> 1116 0.000000e+00 # cquqd8Re1133
> 1117 0.000000e+00 # cquqd8Re2211
> 1121 0.000000e+00 # cquqd8Re2222
> 1125 0.000000e+00 # cquqd8Re2233
> 1126 0.000000e+00 # cquqd8Re3311
> 1130 0.000000e+00 # cquqd8Re3322
> 1134 0.000000e+00 # cquqd8Re3333
> 1189 0.000000e+00 # clequ1Re1111
> 1193 0.000000e+00 # clequ1Re1122
> 1197 0.000000e+00 # clequ1Re1133
> 1198 0.000000e+00 # clequ1Re2211
> 1202 0.000000e+00 # clequ1Re2222
> 1206 0.000000e+00 # clequ1Re2233
> 1207 0.000000e+00 # clequ1Re3311
> 1211 0.000000e+00 # clequ1Re3322
> 1215 0.000000e+00 # clequ1Re3333
> 1270 0.000000e+00 # clequ3Re1111
> 1274 0.000000e+00 # clequ3Re1122
> 1278 0.000000e+00 # clequ3Re1133
> 1279 0.000000e+00 # clequ3Re2211
> 1283 0.000000e+00 # clequ3Re2222
> 1287 0.000000e+00 # clequ3Re2233
> 1288 0.000000e+00 # clequ3Re3311
> 1292 0.000000e+00 # clequ3Re3322
> 1296 0.000000e+00 # clequ3Re3333
>
> ###################################
> ## INFORMATION FOR SMEFTCUTOFF
> ###################################
> Block smeftcutoff
> 1 1.000000e+03 # LambdaSMEFT
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 8.038700e+01 # MW
> 2 1.166379e-05 # Gf
> 3 1.179000e-01 # aS
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.180000e+00 # ymb
> 6 1.727600e+02 # ymt
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.330000e+00 # WT
> DECAY 23 2.495200e+00 # WZ
> DECAY 24 2.085000e+00 # WW
> DECAY 25 4.070000e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.000000e+00 # d : 0.0
> DECAY 2 0.000000e+00 # u : 0.0
> DECAY 3 0.000000e+00 # s : 0.0
> DECAY 4 0.000000e+00 # c : 0.0
> DECAY 5 0.000000e+00 # b : 0.0
> DECAY 11 0.000000e+00 # e- : 0.0
> DECAY 12 0.000000e+00 # ve : 0.0
> DECAY 13 0.000000e+00 # mu- : 0.0
> DECAY 14 0.000000e+00 # vm : 0.0
> DECAY 15 0.000000e+00 # ta- : 0.0
> DECAY 16 0.000000e+00 # vt : 0.0
> DECAY 21 0.000000e+00 # g : 0.0
> DECAY 22 0.000000e+00 # a : 0.0
> DECAY 9000005 2.495200e+00 # z1 : WZ
> DECAY 9000006 2.085000e+00 # w1+ : WW
> DECAY 9000007 1.330000e+00 # t1 : WT
> DECAY 9000008 4.070000e-03 # h1 : WH
> #===========================================================
> # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
> #===========================================================
>
> Block QNUMBERS 9000005 # z1
> 1 0 # 3 times electric charge
> 2 3 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # Particle/Antiparticle distinction (0=own anti)
> Block QNUMBERS 9000006 # w1+
> 1 3 # 3 times electric charge
> 2 3 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # Particle/Antiparticle distinction (0=own anti)
> Block QNUMBERS 9000007 # t1
> 1 2 # 3 times electric charge
> 2 2 # number of spin states (2S+1)
> 3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # Particle/Antiparticle distinction (0=own anti)
> Block QNUMBERS 9000008 # h1
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2S+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # Particle/Antiparticle distinction (0=own anti)#*********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat MadEvent *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a '# ' are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # To display more options, you can type the command: *
> # update full_run_card *
> #*********************************************************************
>
>
> #*******************************
> # Running parameters *
> #*******************************
>
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_1 = run_tag ! name of the run
>
> #*********************************************************************
> # Number of events and rnd seed *
> # Warning: Do not generate more than 1M events in a single run *
> # If you want to run Pythia, avoid more than 50k events in a run. *
> #*********************************************************************
> 10000 = nevents ! Number of unweighted events requested
> 0 = iseed ! rnd seed (0=assigned automatically=default))
>
> #*********************************************************************
> # Collider type and energy *
> # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
> # 3=photon from electron *
> #*********************************************************************
> 1 = lpp1 ! beam 1 type
> 1 = lpp2 ! beam 2 type
> 6500.0 = ebeam1 ! beam 1 total energy in GeV
> 6500.0 = ebeam2 ! beam 2 total energy in GeV
>
> #*********************************************************************
> # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> #*********************************************************************
> nn23lo1 = pdlabel ! PDF set
> 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
>
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> 1.0 = scalefact ! scale factor for event-by-event scales
>
> #*********************************************************************
> # Type and output format *
> #*********************************************************************
> False = gridpack ! True = setting up the grid pack
> -1.0 = time_of_flight ! threshold (in mm) below which the invariant lifetime is not written (-1 means not written)
> 3.0 = lhe_version ! Change the way clustering information pass to shower.
> True = clusinfo ! include clustering tag in output
> sum = event_norm ! average/sum. Normalization of the weight in the LHEF
>
> #*********************************************************************
> # Matching parameter (MLM only) *
> #*********************************************************************
> 0 = ickkw ! 0 no matching, 1 MLM
> 1.0 = alpsfact ! scale factor for QCD emission vx
> False = chcluster ! cluster only according to channel diag
> 4 = asrwgtflavor ! highest quark flavor for a_s reweight
> False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 (turn off for VBF and single top processes)
> 0.0 = xqcut ! minimum kt jet measure between partons
>
> #*********************************************************************
> # handling of the helicities: *
> # 0: sum over all helicities *
> # 1: importance sampling over helicities *
> #*********************************************************************
> 0 = nhel ! using helicities importance sampling or not.
>
> #*********************************************************************
> # Generation bias, check the wiki page below for more information: *
> # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
> #*********************************************************************
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
>
>
> #*******************************
> # Parton level cuts definition *
> #*******************************
>
> #*********************************************************************
> # BW cutoff. Define on/off-shell for "$" and decay *
> #*********************************************************************
> 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
>
> #*********************************************************************
> # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #*********************************************************************
> False = cut_decays ! Cut decay products
>
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20.0 = ptj ! minimum pt for the jets
> 0.0 = ptb ! minimum pt for the b
> 10.0 = pta ! minimum pt for the photons
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptbmax ! maximum pt for the b
> -1.0 = ptamax ! maximum pt for the photons
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
>
> #*********************************************************************
> # Minimum and maximum E's (in the center of mass frame) *
> #*********************************************************************
> 0.0 = ej ! minimum E for the jets
> 0.0 = eb ! minimum E for the b
> 0.0 = ea ! minimum E for the photons
> 0.0 = el ! minimum E for the charged leptons
> -1.0 = ejmax ! maximum E for the jets
> -1.0 = ebmax ! maximum E for the b
> -1.0 = eamax ! maximum E for the photons
> -1.0 = elmax ! maximum E for the charged leptons
> {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
>
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 5.0 = etaj ! max rap for the jets
> -1.0 = etab ! max rap for the b
> 2.5 = etaa ! max rap for the photons
> 2.5 = etal ! max rap for the charged leptons
> 0.0 = etajmin ! min rap for the jets
> 0.0 = etabmin ! min rap for the b
> 0.0 = etaamin ! min rap for the photons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
>
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.4 = drjj ! min distance between jets
> 0.0 = drbb ! min distance between b's
> 0.4 = drll ! min distance between leptons
> 0.4 = draa ! min distance between gammas
> 0.0 = drbj ! min distance between b and jet
> 0.4 = draj ! min distance between gamma and jet
> 0.4 = drjl ! min distance between jet and lepton
> 0.0 = drab ! min distance between gamma and b
> 0.0 = drbl ! min distance between b and lepton
> 0.4 = dral ! min distance between gamma and lepton
> -1.0 = drjjmax ! max distance between jets
> -1.0 = drbbmax ! max distance between b's
> -1.0 = drllmax ! max distance between leptons
> -1.0 = draamax ! max distance between gammas
> -1.0 = drbjmax ! max distance between b and jet
> -1.0 = drajmax ! max distance between gamma and jet
> -1.0 = drjlmax ! max distance between jet and lepton
> -1.0 = drabmax ! max distance between gamma and b
> -1.0 = drblmax ! max distance between b and lepton
> -1.0 = dralmax ! max distance 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
> 120.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
> 130.0 = mmaamax ! max invariant mass of gamma gamma pair
> -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
> ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
>
> #*********************************************************************
> # Minimum and maximum invariant mass for all leptons *
> #*********************************************************************
> 0.0 = mmnl ! min invariant mass for all leptons (l+- and vl)
> -1.0 = mmnlmax ! max invariant mass for all leptons (l+- and vl)
>
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons *
> #*********************************************************************
> 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
>
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0.0 = ptheavy ! minimum pt for at least one heavy final state
> 0.0 = xptj ! minimum pt for at least one jet
> 0.0 = xptb ! minimum pt for at least one b
> 0.0 = xpta ! minimum pt for at least one photon
> 0.0 = xptl ! minimum pt for at least one charged lepton
>
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0.0 = ptj1min ! minimum pt for the leading jet in pt
> 0.0 = ptj2min ! minimum pt for the second jet in pt
> 0.0 = ptj3min ! minimum pt for the third jet in pt
> 0.0 = ptj4min ! minimum pt for the fourth jet in pt
> -1.0 = ptj1max ! maximum pt for the leading jet in pt
> -1.0 = ptj2max ! maximum pt for the second jet in pt
> -1.0 = ptj3max ! maximum pt for the third jet in pt
> -1.0 = ptj4max ! maximum pt for the fourth jet in pt
> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
>
> #*********************************************************************
> # Control the pt's of leptons sorted by pt *
> #*********************************************************************
> 0.0 = ptl1min ! minimum pt for the leading lepton in pt
> 0.0 = ptl2min ! minimum pt for the second lepton in pt
> 0.0 = ptl3min ! minimum pt for the third lepton in pt
> 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1.0 = ptl1max ! maximum pt for the leading lepton in pt
> -1.0 = ptl2max ! maximum pt for the second lepton in pt
> -1.0 = ptl3max ! maximum pt for the third lepton in pt
> -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
>
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
> 0.0 = ihtmin ! inclusive Ht for all partons (including b)
> -1.0 = ihtmax ! inclusive Ht for all partons (including b)
> 0.0 = ht2min ! minimum Ht for the two leading jets
> 0.0 = ht3min ! minimum Ht for the three leading jets
> 0.0 = ht4min ! minimum Ht for the four leading jets
> -1.0 = ht2max ! maximum Ht for the two leading jets
> -1.0 = ht3max ! maximum Ht for the three leading jets
> -1.0 = ht4max ! maximum Ht for the four leading jets
>
> #*********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442 *
> # When ptgmin=0, all the other parameters are ignored *
> # When ptgmin>0, pta and draj are not going to be used *
> #*********************************************************************
> 0.0 = ptgmin ! Min photon transverse momentum
> 0.4 = R0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> True = isoEM ! isolate photons from EM energy (photons and leptons)
>
> #*********************************************************************
> # WBF cuts *
> #*********************************************************************
> 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
> 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
>
> #*********************************************************************
> # Turn on either the ktdurham or ptlund cut to activate *
> # CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
> #*********************************************************************
> -1.0 = ktdurham
> 0.4 = dparameter
> -1.0 = ptlund
> 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above
>
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
>
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: Do not use for interference type of computation *
> #*********************************************************************
> # False = use_syst ! Enable systematics studies # Commented out by MadMiner
> # systematics = systematics_program ! none, systematics [python], SysCalc [deprecated, C++] # Commented out by MadMiner
> # ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule # Commented out by MadMiner
>
>
>
> #*********************************************************************
> # MadMiner systematics setup *
> #*********************************************************************
> False = use_syst
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

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

To post a message you must log in.