TopEFT Gridpack Creation

Asked by Douglas Burns on 2018-04-16

Dear Madgraph experts,

I am trying to create some gridpacks for ttbar production using an eft model (dim6top: http://feynrules.irmp.ucl.ac.be/wiki/dim6top).

I am running the process :

p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~ > w- b~ , w- > wdecay wdecay)
where wdecay = g u c d s u~ c~ d~ s~ e+ mu+ tau+ e- mu- tau- ve vm vt ve~ vm~ vt~

I am however running into problems. When trying to create the gridpacks i end up with:

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

     Cross-section : 0 +- 0 pb
     Nb of events : 0

before failing with the error:

Command "generate_events " interrupted with error:
IndexError : list index out of range

and in the log file (run_01_tag_1_debug.log):

  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1467, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 2487, in do_generate_events
    self.exec_cmd('combine_events', postcmd=False)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1494, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 3593, in do_combine_events
    proc_charac=self.proc_characteristic)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 1114, in unweight
    return super(MultiEventFile, self).unweight(outputpath, get_wgt_multi, **opts)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 398, in unweight
    max_wgt = all_wgt[-1]
IndexError: list index out of range

Looking into some of the SubProcesses logs, I can see that none of the events passed the phase space cuts and it says to loosen cuts or increase max_events. I do not know where to increase the max_events from and the phase space cuts i think are the default ones which work fine with SM ttbar production.

Cheers,
Doug

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2018-04-17
Last query:
2018-04-17
Last reply:
2018-04-16
Douglas Burns (doburns) said : #1

Here is the full error log:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events
Traceback (most recent call last):
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1467, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 2487, in do_generate_events
    self.exec_cmd('combine_events', postcmd=False)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1494, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 3593, in do_combine_events
    proc_charac=self.proc_characteristic)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 1114, in unweight
    return super(MultiEventFile, self).unweight(outputpath, get_wgt_multi, **opts)
  File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 398, in unweight
    max_wgt = all_wgt[-1]
IndexError: list index out of range
                              Run Options
                              -----------
               stdout_level : None

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : lhapdf-config
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 0 (user set)
        exrootanalysis_path : None (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.1 2017-12-12 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model dim6top_LO_UFO
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~
define l = l+ l- ta+ ta-
define nu = vl vl~
define wdecay = j l nu
generate p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~\
 > w- b~ , w- > wdecay wdecay)
output process -nojpeg
######################################################################
## 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 DIM6
###################################
Block dim6
    1 1.000000e+03 # Lambda
    2 0.000000e+00 # ctp
    3 0.000000e+00 # ctpI
    4 0.000000e+00 # cpQM
    5 0.000000e+00 # cpQ3
    6 0.000000e+00 # cpt
    7 0.000000e+00 # cpb
    8 0.000000e+00 # cptb
    9 0.000000e+00 # cptbI
   10 0.000000e+00 # ctW
   11 0.000000e+00 # ctZ
   12 0.000000e+00 # ctWI
   13 0.000000e+00 # ctZI
   14 0.000000e+00 # cbW
   15 0.000000e+00 # cbWI
   16 0.000000e+00 # ctG
   17 0.000000e+00 # ctGI
   18 0.000000e+00 # cQlM1
   19 0.000000e+00 # cQlM2
   20 0.000000e+00 # cQlM3
   21 0.000000e+00 # cQl31
   22 0.000000e+00 # cQl32
   23 0.000000e+00 # cQl33
   24 0.000000e+00 # cQe1
   25 0.000000e+00 # cQe2
   26 0.000000e+00 # cQe3
   27 0.000000e+00 # ctl1
   28 0.000000e+00 # ctl2
   29 0.000000e+00 # ctl3
   30 0.000000e+00 # cte1
   31 0.000000e+00 # cte2
   32 0.000000e+00 # cte3
   33 0.000000e+00 # ctlS1
   34 0.000000e+00 # ctlSI1
   35 0.000000e+00 # ctlS2
   36 0.000000e+00 # ctlSI2
   37 0.000000e+00 # ctlS3
   38 0.000000e+00 # ctlSI3
   39 0.000000e+00 # ctlT1
   40 0.000000e+00 # ctlTI1
   41 0.000000e+00 # ctlT2
   42 0.000000e+00 # ctlTI2
   43 0.000000e+00 # ctlT3
   44 0.000000e+00 # ctlTI3
   45 0.000000e+00 # cblS1
   46 0.000000e+00 # cblSI1
   47 0.000000e+00 # cblS2
   48 0.000000e+00 # cblSI2
   49 0.000000e+00 # cblS3
   50 0.000000e+00 # cblSI3
   51 0.000000e+00 # cQq83
   52 0.000000e+00 # cQq81
   53 0.000000e+00 # cQu8
   54 0.000000e+00 # cQd8
   55 0.000000e+00 # ctq8
   56 0.000000e+00 # ctu8
   57 0.000000e+00 # ctd8
   58 0.000000e+00 # cQq13
   59 0.000000e+00 # cQq11
   60 0.000000e+00 # cQu1
   61 0.000000e+00 # cQd1
   62 0.000000e+00 # ctq1
   63 0.000000e+00 # ctu1
   64 0.000000e+00 # ctd1
   65 0.000000e+00 # cQQ1
   66 0.000000e+00 # cQQ8
   67 0.000000e+00 # cQt1
   68 0.000000e+00 # cQb1
   69 0.000000e+00 # ctt1
   70 0.000000e+00 # ctb1
   71 0.000000e+00 # cQt8
   72 0.000000e+00 # cQb8
   73 0.000000e+00 # ctb8
   74 0.000000e+00 # ctQqu1
   75 0.000000e+00 # ctQqu1I
   76 0.000000e+00 # ctQqu8
   77 0.000000e+00 # ctQqu8I
   78 0.000000e+00 # cbQqd1
   79 0.000000e+00 # cbQqd1I
   80 0.000000e+00 # cbQqd8
   81 0.000000e+00 # cbQqd8I
   82 0.000000e+00 # cQtqd1
   83 0.000000e+00 # cQtqd1I
   84 0.000000e+00 # cQtqd8
   85 0.000000e+00 # cQtqd8I
   86 0.000000e+00 # cQbqu1
   87 0.000000e+00 # cQbqu1I
   88 0.000000e+00 # cQbqu8
   89 0.000000e+00 # cQbqu8I
   90 0.000000e+00 # cbtud1
   91 0.000000e+00 # cbtud1I
   92 0.000000e+00 # cbtud8
   93 0.000000e+00 # cbtud8I
   94 0.000000e+00 # cQtQb1
   95 0.000000e+00 # cQtQb8
   96 0.000000e+00 # cQtQb1I
   97 0.000000e+00 # cQtQb8I

###################################
## INFORMATION FOR FCNC
###################################
Block fcnc
    1 0.000000e+00 # ctpx13
    2 0.000000e+00 # ctpx23
    3 0.000000e+00 # ctpx31
    4 0.000000e+00 # ctpx32
    5 0.000000e+00 # ctpIx13
    6 0.000000e+00 # ctpIx23
    7 0.000000e+00 # ctpIx31
    8 0.000000e+00 # ctpIx32
    9 0.000000e+00 # cpQMx31
   10 0.000000e+00 # cpQMx32
   11 0.000000e+00 # cpQMIx31
   12 0.000000e+00 # cpQMIx32
   13 0.000000e+00 # cpQ3x31
   14 0.000000e+00 # cpQ3x32
   15 0.000000e+00 # cpQ3Ix31
   16 0.000000e+00 # cpQ3Ix32
   17 0.000000e+00 # cptx31
   18 0.000000e+00 # cptx32
   19 0.000000e+00 # cptIx31
   20 0.000000e+00 # cptIx32
   21 0.000000e+00 # cptbx13
   22 0.000000e+00 # cptbx23
   23 0.000000e+00 # cptbx31
   24 0.000000e+00 # cptbx32
   25 0.000000e+00 # cptbIx13
   26 0.000000e+00 # cptbIx23
   27 0.000000e+00 # cptbIx31
   28 0.000000e+00 # cptbIx32
   29 0.000000e+00 # ctWx13
   30 0.000000e+00 # ctWx23
   31 0.000000e+00 # ctWx31
   32 0.000000e+00 # ctWx32
   33 0.000000e+00 # ctWIx13
   34 0.000000e+00 # ctWIx23
   35 0.000000e+00 # ctWIx31
   36 0.000000e+00 # ctWIx32
   37 0.000000e+00 # ctZx13
   38 0.000000e+00 # ctZx23
   39 0.000000e+00 # ctZx31
   40 0.000000e+00 # ctZx32
   41 0.000000e+00 # ctZIx13
   42 0.000000e+00 # ctZIx23
   43 0.000000e+00 # ctZIx31
   44 0.000000e+00 # ctZIx32
   45 0.000000e+00 # cbWx13
   46 0.000000e+00 # cbWx23
   47 0.000000e+00 # cbWx31
   48 0.000000e+00 # cbWx32
   49 0.000000e+00 # cbWIx13
   50 0.000000e+00 # cbWIx23
   51 0.000000e+00 # cbWIx31
   52 0.000000e+00 # cbWIx32
   53 0.000000e+00 # ctGx13
   54 0.000000e+00 # ctGx23
   55 0.000000e+00 # ctGx31
   56 0.000000e+00 # ctGx32
   57 0.000000e+00 # ctGIx13
   58 0.000000e+00 # ctGIx23
   59 0.000000e+00 # ctGIx31
   60 0.000000e+00 # ctGIx32
   61 0.000000e+00 # cQl3x1x31
   62 0.000000e+00 # cQl3x1x32
   63 0.000000e+00 # cQl3x2x31
   64 0.000000e+00 # cQl3x2x32
   65 0.000000e+00 # cQl3x3x31
   66 0.000000e+00 # cQl3x3x32
   67 0.000000e+00 # cQl3Ix1x31
   68 0.000000e+00 # cQl3Ix1x32
   69 0.000000e+00 # cQl3Ix2x31
   70 0.000000e+00 # cQl3Ix2x32
   71 0.000000e+00 # cQl3Ix3x31
   72 0.000000e+00 # cQl3Ix3x32
   73 0.000000e+00 # cQlMx1x31
   74 0.000000e+00 # cQlMx1x32
   75 0.000000e+00 # cQlMx2x31
   76 0.000000e+00 # cQlMx2x32
   77 0.000000e+00 # cQlMx3x31
   78 0.000000e+00 # cQlMx3x32
   79 0.000000e+00 # cQlMIx1x31
   80 0.000000e+00 # cQlMIx1x32
   81 0.000000e+00 # cQlMIx2x31
   82 0.000000e+00 # cQlMIx2x32
   83 0.000000e+00 # cQlMIx3x31
   84 0.000000e+00 # cQlMIx3x32
   85 0.000000e+00 # cQex1x31
   86 0.000000e+00 # cQex1x32
   87 0.000000e+00 # cQex2x31
   88 0.000000e+00 # cQex2x32
   89 0.000000e+00 # cQex3x31
   90 0.000000e+00 # cQex3x32
   91 0.000000e+00 # cQeIx1x31
   92 0.000000e+00 # cQeIx1x32
   93 0.000000e+00 # cQeIx2x31
   94 0.000000e+00 # cQeIx2x32
   95 0.000000e+00 # cQeIx3x31
   96 0.000000e+00 # cQeIx3x32
   97 0.000000e+00 # ctlx1x31
   98 0.000000e+00 # ctlx1x32
   99 0.000000e+00 # ctlx2x31
  100 0.000000e+00 # ctlx2x32
  101 0.000000e+00 # ctlx3x31
  102 0.000000e+00 # ctlx3x32
  103 0.000000e+00 # ctlIx1x31
  104 0.000000e+00 # ctlIx1x32
  105 0.000000e+00 # ctlIx2x31
  106 0.000000e+00 # ctlIx2x32
  107 0.000000e+00 # ctlIx3x31
  108 0.000000e+00 # ctlIx3x32
  109 0.000000e+00 # ctex1x31
  110 0.000000e+00 # ctex1x32
  111 0.000000e+00 # ctex2x31
  112 0.000000e+00 # ctex2x32
  113 0.000000e+00 # ctex3x31
  114 0.000000e+00 # ctex3x32
  115 0.000000e+00 # cteIx1x31
  116 0.000000e+00 # cteIx1x32
  117 0.000000e+00 # cteIx2x31
  118 0.000000e+00 # cteIx2x32
  119 0.000000e+00 # cteIx3x31
  120 0.000000e+00 # cteIx3x32
  121 0.000000e+00 # ctlSx1x13
  122 0.000000e+00 # ctlSx1x23
  123 0.000000e+00 # ctlSx1x31
  124 0.000000e+00 # ctlSx1x32
  125 0.000000e+00 # ctlSx2x13
  126 0.000000e+00 # ctlSx2x23
  127 0.000000e+00 # ctlSx2x31
  128 0.000000e+00 # ctlSx2x32
  129 0.000000e+00 # ctlSx3x13
  130 0.000000e+00 # ctlSx3x23
  131 0.000000e+00 # ctlSx3x31
  132 0.000000e+00 # ctlSx3x32
  133 0.000000e+00 # ctlSIx1x13
  134 0.000000e+00 # ctlSIx1x23
  135 0.000000e+00 # ctlSIx1x31
  136 0.000000e+00 # ctlSIx1x32
  137 0.000000e+00 # ctlSIx2x13
  138 0.000000e+00 # ctlSIx2x23
  139 0.000000e+00 # ctlSIx2x31
  140 0.000000e+00 # ctlSIx2x32
  141 0.000000e+00 # ctlSIx3x13
  142 0.000000e+00 # ctlSIx3x23
  143 0.000000e+00 # ctlSIx3x31
  144 0.000000e+00 # ctlSIx3x32
  145 0.000000e+00 # ctlTx1x13
  146 0.000000e+00 # ctlTx1x23
  147 0.000000e+00 # ctlTx1x31
  148 0.000000e+00 # ctlTx1x32
  149 0.000000e+00 # ctlTx2x13
  150 0.000000e+00 # ctlTx2x23
  151 0.000000e+00 # ctlTx2x31
  152 0.000000e+00 # ctlTx2x32
  153 0.000000e+00 # ctlTx3x13
  154 0.000000e+00 # ctlTx3x23
  155 0.000000e+00 # ctlTx3x31
  156 0.000000e+00 # ctlTx3x32
  157 0.000000e+00 # ctlTIx1x13
  158 0.000000e+00 # ctlTIx1x23
  159 0.000000e+00 # ctlTIx1x31
  160 0.000000e+00 # ctlTIx1x32
  161 0.000000e+00 # ctlTIx2x13
  162 0.000000e+00 # ctlTIx2x23
  163 0.000000e+00 # ctlTIx2x31
  164 0.000000e+00 # ctlTIx2x32
  165 0.000000e+00 # ctlTIx3x13
  166 0.000000e+00 # ctlTIx3x23
  167 0.000000e+00 # ctlTIx3x31
  168 0.000000e+00 # ctlTIx3x32
  169 0.000000e+00 # cqq11x3331
  170 0.000000e+00 # cqq11x3332
  171 0.000000e+00 # cqq11Ix3331
  172 0.000000e+00 # cqq11Ix3332
  173 0.000000e+00 # cqq13x3331
  174 0.000000e+00 # cqq13x3332
  175 0.000000e+00 # cqq13Ix3331
  176 0.000000e+00 # cqq13Ix3332
  177 0.000000e+00 # cuu1x3331
  178 0.000000e+00 # cuu1x3332
  179 0.000000e+00 # cuu1Ix3331
  180 0.000000e+00 # cuu1Ix3332
  181 0.000000e+00 # cqu1x3331
  182 0.000000e+00 # cqu1x3332
  183 0.000000e+00 # cqu1Ix3331
  184 0.000000e+00 # cqu1Ix3332
  185 0.000000e+00 # cqu1x3133
  186 0.000000e+00 # cqu1x3233
  187 0.000000e+00 # cqu1Ix3133
  188 0.000000e+00 # cqu1Ix3233
  189 0.000000e+00 # cqu8x3331
  190 0.000000e+00 # cqu8x3332
  191 0.000000e+00 # cqu8Ix3331
  192 0.000000e+00 # cqu8Ix3332
  193 0.000000e+00 # cqu8x3233
  194 0.000000e+00 # cqu8x3133
  195 0.000000e+00 # cqu8Ix3133
  196 0.000000e+00 # cqu8Ix3233
  197 0.000000e+00 # cqd1x3331
  198 0.000000e+00 # cqd1x3332
  199 0.000000e+00 # cqd1Ix3331
  200 0.000000e+00 # cqd1Ix3332
  201 0.000000e+00 # cqd1x3133
  202 0.000000e+00 # cqd1x3233
  203 0.000000e+00 # cqd1Ix3133
  204 0.000000e+00 # cqd1Ix3233
  205 0.000000e+00 # cqd8x3331
  206 0.000000e+00 # cqd8x3332
  207 0.000000e+00 # cqd8Ix3331
  208 0.000000e+00 # cqd8Ix3332
  209 0.000000e+00 # cqd8x3133
  210 0.000000e+00 # cqd8x3233
  211 0.000000e+00 # cqd8Ix3133
  212 0.000000e+00 # cqd8Ix3233
  213 0.000000e+00 # cud1x3331
  214 0.000000e+00 # cud1x3332
  215 0.000000e+00 # cud1Ix3331
  216 0.000000e+00 # cud1Ix3332
  217 0.000000e+00 # cud1x3133
  218 0.000000e+00 # cud1x3233
  219 0.000000e+00 # cud1Ix3133
  220 0.000000e+00 # cud1Ix3233
  221 0.000000e+00 # cud8x3331
  222 0.000000e+00 # cud8x3332
  223 0.000000e+00 # cud8Ix3331
  224 0.000000e+00 # cud8Ix3332
  225 0.000000e+00 # cud8x3133
  226 0.000000e+00 # cud8x3233
  227 0.000000e+00 # cud8Ix3133
  228 0.000000e+00 # cud8Ix3233
  229 0.000000e+00 # cquqd1x3331
  230 0.000000e+00 # cquqd1x3332
  231 0.000000e+00 # cquqd1x3313
  232 0.000000e+00 # cquqd1x3323
  233 0.000000e+00 # cquqd1x3133
  234 0.000000e+00 # cquqd1x3233
  235 0.000000e+00 # cquqd1x1333
  236 0.000000e+00 # cquqd1x2333
  237 0.000000e+00 # cquqd1Ix3331
  238 0.000000e+00 # cquqd1Ix3332
  239 0.000000e+00 # cquqd1Ix3313
  240 0.000000e+00 # cquqd1Ix3323
  241 0.000000e+00 # cquqd1Ix3133
  242 0.000000e+00 # cquqd1Ix3233
  243 0.000000e+00 # cquqd1Ix1333
  244 0.000000e+00 # cquqd1Ix2333
  245 0.000000e+00 # cquqd8x3331
  246 0.000000e+00 # cquqd8x3332
  247 0.000000e+00 # cquqd8x3313
  248 0.000000e+00 # cquqd8x3323
  249 0.000000e+00 # cquqd8x3133
  250 0.000000e+00 # cquqd8x3233
  251 0.000000e+00 # cquqd8x1333
  252 0.000000e+00 # cquqd8x2333
  253 0.000000e+00 # cquqd8Ix3331
  254 0.000000e+00 # cquqd8Ix3332
  255 0.000000e+00 # cquqd8Ix3313
  256 0.000000e+00 # cquqd8Ix3323
  257 0.000000e+00 # cquqd8Ix3133
  258 0.000000e+00 # cquqd8Ix3233
  259 0.000000e+00 # cquqd8Ix1333
  260 0.000000e+00 # cquqd8Ix2333
  261 0.000000e+00 # cqq11x31ii
  262 0.000000e+00 # cqq11x32ii
  263 0.000000e+00 # cqq11Ix31ii
  264 0.000000e+00 # cqq11Ix32ii
  265 0.000000e+00 # cqq13x31ii
  266 0.000000e+00 # cqq13x32ii
  267 0.000000e+00 # cqq13Ix31ii
  268 0.000000e+00 # cqq13Ix32ii
  269 0.000000e+00 # cqq81x31ii
  270 0.000000e+00 # cqq81x32ii
  271 0.000000e+00 # cqq81Ix31ii
  272 0.000000e+00 # cqq81Ix32ii
  273 0.000000e+00 # cqq83x31ii
  274 0.000000e+00 # cqq83x32ii
  275 0.000000e+00 # cqq83Ix31ii
  276 0.000000e+00 # cqq83Ix32ii
  277 0.000000e+00 # cuu1x31ii
  278 0.000000e+00 # cuu1x32ii
  279 0.000000e+00 # cuu1Ix31ii
  280 0.000000e+00 # cuu1Ix32ii
  281 0.000000e+00 # cuu8x31ii
  282 0.000000e+00 # cuu8x32ii
  283 0.000000e+00 # cuu8Ix31ii
  284 0.000000e+00 # cuu8Ix32ii
  285 0.000000e+00 # cud1x31ii
  286 0.000000e+00 # cud1x32ii
  287 0.000000e+00 # cud1Ix31ii
  288 0.000000e+00 # cud1Ix32ii
  289 0.000000e+00 # cud8x31ii
  290 0.000000e+00 # cud8x32ii
  291 0.000000e+00 # cud8Ix31ii
  292 0.000000e+00 # cud8Ix32ii
  293 0.000000e+00 # cqu1x31ii
  294 0.000000e+00 # cqu1x32ii
  295 0.000000e+00 # cqu1Ix31ii
  296 0.000000e+00 # cqu1Ix32ii
  297 0.000000e+00 # cqu8x31ii
  298 0.000000e+00 # cqu8x32ii
  299 0.000000e+00 # cqu8Ix31ii
  300 0.000000e+00 # cqu8Ix32ii
  301 0.000000e+00 # cqu1xii31
  302 0.000000e+00 # cqu1xii32
  303 0.000000e+00 # cqu1Ixii31
  304 0.000000e+00 # cqu1Ixii32
  305 0.000000e+00 # cqu8xii31
  306 0.000000e+00 # cqu8xii32
  307 0.000000e+00 # cqu8Ixii31
  308 0.000000e+00 # cqu8Ixii32
  309 0.000000e+00 # cqd1x31ii
  310 0.000000e+00 # cqd1x32ii
  311 0.000000e+00 # cqd1Ix31ii
  312 0.000000e+00 # cqd1Ix32ii
  313 0.000000e+00 # cqd8Ix31ii
  314 0.000000e+00 # cqd8Ix32ii
  315 0.000000e+00 # cqd8x31ii
  316 0.000000e+00 # cqd8x32ii
  317 0.000000e+00 # cpbIx31
  318 0.000000e+00 # cpbIx32
  319 0.000000e+00 # cpbx31
  320 0.000000e+00 # cpbx32
  321 0.000000e+00 # cblSx1x13
  322 0.000000e+00 # cblSx1x23
  323 0.000000e+00 # cblSx1x31
  324 0.000000e+00 # cblSx1x32
  325 0.000000e+00 # cblSx2x13
  326 0.000000e+00 # cblSx2x23
  327 0.000000e+00 # cblSx2x31
  328 0.000000e+00 # cblSx2x32
  329 0.000000e+00 # cblSx3x13
  330 0.000000e+00 # cblSx3x23
  331 0.000000e+00 # cblSx3x31
  332 0.000000e+00 # cblSx3x32
  333 0.000000e+00 # cblSIx1x13
  334 0.000000e+00 # cblSIx1x23
  335 0.000000e+00 # cblSIx1x31
  336 0.000000e+00 # cblSIx1x32
  337 0.000000e+00 # cblSIx2x13
  338 0.000000e+00 # cblSIx2x23
  339 0.000000e+00 # cblSIx2x31
  340 0.000000e+00 # cblSIx2x32
  341 0.000000e+00 # cblSIx3x13
  342 0.000000e+00 # cblSIx3x23
  343 0.000000e+00 # cblSIx3x31
  344 0.000000e+00 # cblSIx3x32

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.720000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118760e+01 # MZ
   25 1.250000e+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.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 79.824360 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.279000e+02 # aEWM1
    2 1.166370e-05 # Gf
    3 1.184000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.700000e+00 # ymb
    6 1.720000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.085000e+00 # WW
DECAY 25 4.070000e-03 # WH
## 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.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
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
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
     0.0 = polbeam1 ! beam polarization for beam 1
     0.0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.188 = scale ! fixed ren scale
 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
.true. = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  3.0 = lhe_version ! Change the way clustering information pass to shower.
  True = clusinfo ! include clustering tag in output
  average = event_norm ! average/sum. Normalization of the weight in the LHEF

#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
1 = 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)
30.0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
   False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20.0 = ptj ! minimum pt for the jets
 0.0 = ptb ! minimum pt for the b
 10.0 = pta ! minimum pt for the photons
 10.0 = ptl ! minimum pt for the charged leptons
 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
 -1.0 = ptjmax ! maximum pt for the jets
 -1.0 = ptbmax ! maximum pt for the b
 -1.0 = ptamax ! maximum pt for the photons
 -1.0 = ptlmax ! maximum pt for the charged leptons
 -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
 {} = 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 = 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 = 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
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
0 = drjlmax ! max distance between jet and lepton
 -1.0 = drabmax ! max distance between gamma and b
 -1.0 = drblmax ! max distance between b and lepton
 -1.0 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
 0.0 = mmjj ! min invariant mass of a jet pair
 0.0 = mmbb ! min invariant mass of a b pair
 0.0 = mmaa ! min invariant mass of gamma gamma pair
 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1.0 = mmjjmax ! max invariant mass of a jet pair
 -1.0 = mmbbmax ! max invariant mass of a b pair
 -1.0 = mmaamax ! max invariant mass of gamma gamma pair
 -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
 {} = 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 letpons *
#*********************************************************************
 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = ptheavy ! minimum pt for at least one heavy final state
 0.0 = xptj ! minimum pt for at least one jet
 0.0 = xptb ! minimum pt for at least one b
 0.0 = xpta ! minimum pt for at least one photon
 0.0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0.0 = ptj1min ! minimum pt for the leading jet in pt
 0.0 = ptj2min ! minimum pt for the second jet in pt
 0.0 = ptj3min ! minimum pt for the third jet in pt
 0.0 = ptj4min ! minimum pt for the fourth jet in pt
 -1.0 = ptj1max ! maximum pt for the leading jet in pt
 -1.0 = ptj2max ! maximum pt for the second jet in pt
 -1.0 = ptj3max ! maximum pt for the third jet in pt
 -1.0 = ptj4max ! maximum pt for the fourth jet in pt
 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
 0.0 = ptl1min ! minimum pt for the leading lepton in pt
 0.0 = ptl2min ! minimum pt for the second lepton in pt
 0.0 = ptl3min ! minimum pt for the third lepton in pt
 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1.0 = ptl1max ! maximum pt for the leading lepton in pt
 -1.0 = ptl2max ! maximum pt for the second lepton in pt
 -1.0 = ptl3max ! maximum pt for the third lepton in pt
 -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
 0.0 = ihtmin !inclusive Ht for all partons (including b)
 -1.0 = ihtmax !inclusive Ht for all partons (including b)
 0.0 = ht2min ! minimum Ht for the two leading jets
 0.0 = ht3min ! minimum Ht for the three leading jets
 0.0 = ht4min ! minimum Ht for the four leading jets
 -1.0 = ht2max ! maximum Ht for the two leading jets
 -1.0 = ht3max ! maximum Ht for the three leading jets
 -1.0 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
 0.0 = ptgmin ! Min photon transverse momentum
 0.4 = R0gamma ! Radius of isolation code
 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
 True = isoEM ! isolate photons from EM energy (photons and leptons)
#*********************************************************************
# WBF cuts *
#*********************************************************************
 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#***********************************************************************
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#***********************************************************************
 -1.0 = ktdurham
 0.4 = dparameter
 -1.0 = ptlund
 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
#**************************************
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#**************************************
#
0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_lo_as_0130_qed = sys_pdf # list of pdf sets
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf
#

Hi,

Looks like that the model request you to do something like:
 generate p p > t t~ FCNC=0 DIM6=1, (t > w+ b DIM6=0, w+ > wdecay wdecay DIM6=0), (t~ > w- b~ DIM6=0, w- > wdecay wdecay DIM6=0)

Cheers,

Olivier

> On 16 Apr 2018, at 14:13, Douglas Burns <email address hidden> wrote:
>
> New question #667989 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/667989
>
> Dear Madgraph experts,
>
> I am trying to create some gridpacks for ttbar production using an eft model (dim6top: http://feynrules.irmp.ucl.ac.be/wiki/dim6top).
>
> I am running the process :
>
> p p > t t~ FCNC=0 DIM6=1, (t > w+ b, w+ > wdecay wdecay), (t~ > w- b~ , w- > wdecay wdecay)
> where wdecay = g u c d s u~ c~ d~ s~ e+ mu+ tau+ e- mu- tau- ve vm vt ve~ vm~ vt~
>
> I am however running into problems. When trying to create the gridpacks i end up with:
>
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 0 +- 0 pb
> Nb of events : 0
>
> before failing with the error:
>
> Command "generate_events " interrupted with error:
> IndexError : list index out of range
>
> and in the log file (run_01_tag_1_debug.log):
>
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1467, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
> return func(arg, **opt)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 2487, in do_generate_events
> self.exec_cmd('combine_events', postcmd=False)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1494, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
> return func(arg, **opt)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/madevent_interface.py", line 3593, in do_combine_events
> proc_charac=self.proc_characteristic)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 1114, in unweight
> return super(MultiEventFile, self).unweight(outputpath, get_wgt_multi, **opts)
> File "/scratch/db0268/EFT/TTJetsEFT/process/bin/internal/lhe_parser.py", line 398, in unweight
> max_wgt = all_wgt[-1]
> IndexError: list index out of range
>
> Looking into some of the SubProcesses logs, I can see that none of the events passed the phase space cuts and it says to loosen cuts or increase max_events. I do not know where to increase the max_events from and the phase space cuts i think are the default ones which work fine with SM ttbar production.
>
> Cheers,
> Doug
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Douglas Burns (doburns) said : #3

Thanks Olivier Mattelaer, that solved my question.