wrong support of epsilon_ijk color structure

[Olivier Mattelaer]

Hi

I am trying to generate a signal and my jobs are failling with this error:

Command "generate_events RPVStop100tojj_upto2j_13TeV --cluster -f" interrupted with error:
MadEventError : Error detected Stop running: Error: Color combination

Here is the log file:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.2.2 2014-11-06 *
#* *
#* 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 RPVStop100tojj_upto2j_13TeV --cluster -f
Traceback (most recent call last):
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/extended_cmd.py", line 879, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/madevent_interface.py", line 1923, in do_generate_events
    postcmd=False)
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/cms/gomez/MG5_aMC_v2_2_2/madgraph/interface/madevent_interface.py", line 2361, in do_survey
    open(pjoin(self.me_dir,'error')).read()
MadEventError: Error detected Stop running: Error: Color combination 1 2 -> 0 not implemented in addmothers.f

                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
          cluster_temp_path : None
              cluster_queue : None (user set)
             cluster_memory : None (user set)
               cluster_time : None (user set)
                    nb_core : 8 (user set)
                   run_mode : 1 (user set)

                      Configuration Options
                      ---------------------
                text_editor : 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
        exrootanalysis_path : None (user set)
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : lhapdf-config
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 0 (user set)
           cluster_nb_retry : 1
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.2.2 2014-11-06 *
#* *
#* 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 loop_optimized_output True
set group_subprocesses Auto
set ignore_six_quark_processes False
set gauge unitary
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model RPVMSSM_UFO -modelname
define j = g u c d s u~ c~ d~ s~ b b~
generate p p > su1 su1~ /su2 QED=0, su1 > d~ s~, su1~ > d s @0
add process p p > su1 su1~ j /su2 QED=0, su1 > d~ s~, su1~ > d s @1
add process p p > su1 su1~ j j /su2 QED=0, su1 > d~ s~, su1~ > d s @2
output RPVStop_UDD312_ISR_upto2j
######################################################################
## 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 DSQMIX
###################################
Block dsqmix
    1 3 9.862151e-01 # RRd1x3
    1 6 1.654681e-01 # RRd1x6
    2 3 -1.654681e-01 # RRd2x3
    2 6 9.862151e-01 # RRd2x6
    3 5 1.000000e+00 # RRd3x5
    4 4 1.000000e+00 # RRd4x4
    5 2 1.000000e+00 # RRd5x2
    6 1 1.000000e+00 # RRd6x1

###################################
## INFORMATION FOR FRALPHA
###################################
Block fralpha
    1 -1.075727e-01 # alp

###################################
## INFORMATION FOR HMIX
###################################
Block hmix
    1 4.642435e+02 # RMUH
    2 9.749107e+00 # tb
    4 3.251216e+05 # MA2

###################################
## INFORMATION FOR IMVCKM
###################################
Block imvckm
    1 3 -3.349042e-03 # ICKM1x3
    2 1 -1.377152e-04 # ICKM2x1
    2 2 -3.212913e-05 # ICKM2x2
    3 1 -3.258550e-03 # ICKM3x1
    3 2 -7.602238e-04 # ICKM3x2

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.832388e+00 # MB
    6 1.731000e+02 # MT
   15 1.776990e+00 # Mta
   23 9.118760e+01 # MZ
   24 7.982901e+01 # MW
   25 1.139668e+02 # MH01
   35 5.643097e+02 # MH02
   36 5.646897e+02 # MA0
   37 5.704586e+02 # MH
  1000001 5.378740e+02 # Msd1
  1000002 100 # Msu1
  1000003 6.105113e+02 # Msd2
  1000004 6.133992e+02 # Msu2
  1000005 6.155144e+02 # Msd3
  1000006 6.157041e+02 # Msu3
  1000011 3.037530e+02 # Msl1
  1000012 3.336805e+02 # Msn1
  1000013 3.168669e+02 # Msl2
  1000014 3.370975e+02 # Msn2
  1000015 3.169185e+02 # Msl3
  1000016 3.371096e+02 # Msn3
  1000021 6.161578e+02 # Mgo
  1000022 9.948601e+01 # Mneu1
  1000023 1.895047e+02 # Mneu2
  1000024 1.894680e+02 # Mch1
  1000025 -4.695446e+02 # Mneu3
  1000035 4.773438e+02 # Mneu4
  1000037 4.805668e+02 # Mch2
  2000001 6.155223e+02 # Msd4
  2000002 6.157276e+02 # Msu4
  2000003 6.349333e+02 # Msd5
  2000004 6.302198e+02 # Msu5
  2000005 6.349335e+02 # Msd6
  2000006 6.302357e+02 # Msu6
  2000011 3.464667e+02 # Msl4
  2000013 3.464780e+02 # Msl5
  2000015 3.481667e+02 # Msl6
## 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

###################################
## INFORMATION FOR MSD2
###################################
Block msd2
    1 1 3.543979e+05 # RmD21x1
    2 2 3.543929e+05 # RmD22x2
    3 3 3.461351e+05 # RmD23x3

###################################
## INFORMATION FOR MSE2
###################################
Block mse2
    1 1 9.819512e+04 # RmE21x1
    2 2 9.817883e+04 # RmE22x2
    3 3 9.361293e+04 # RmE23x3

###################################
## INFORMATION FOR MSL2
###################################
Block msl2
    1 1 1.174128e+05 # RmL21x1
    2 2 1.174047e+05 # RmL22x2
    3 3 1.151411e+05 # RmL23x3

###################################
## INFORMATION FOR MSOFT
###################################
Block msoft
    1 1.021024e+02 # RMx1
    2 1.926534e+02 # RMx2
    3 5.884243e+02 # RMx3
   21 1.021411e+05 # mHu2
   22 -2.154519e+05 # mHd2

###################################
## INFORMATION FOR MSQ2
###################################
Block msq2
    1 1 3.752466e+05 # RmQ21x1
    2 2 3.752418e+05 # RmQ22x2
    3 3 2.699498e+05 # RmQ23x3

###################################
## INFORMATION FOR MSU2
###################################
Block msu2
    1 1 3.564881e+05 # RmU21x1
    2 2 3.564835e+05 # RmU22x2
    3 3 1.494114e+05 # RmU23x3

###################################
## INFORMATION FOR NMIX
###################################
Block nmix
    1 1 -9.927936e-01 # RNN1x1
    1 2 3.701124e-02 # RNN1x2
    1 3 -1.087277e-01 # RNN1x3
    1 4 3.419427e-02 # RNN1x4
    2 1 -6.145626e-02 # RNN2x1
    2 2 -9.731574e-01 # RNN2x2
    2 3 1.992759e-01 # RNN2x3
    2 4 -9.734984e-02 # RNN2x4
    3 1 4.928328e-02 # RNN3x1
    3 2 -7.521552e-02 # RNN3x2
    3 3 -6.987909e-01 # RNN3x3
    3 4 -7.096514e-01 # RNN3x4
    4 1 9.030553e-02 # RNN4x1
    4 2 -2.143305e-01 # RNN4x2
    4 3 -6.783499e-01 # RNN4x3
    4 4 6.969568e-01 # RNN4x4

###################################
## INFORMATION FOR RVLAMLLE
###################################
Block rvlamlle
    1 2 1 2.000000e-01 # RLLE1x2x1
    1 2 2 2.000000e-01 # RLLE1x2x2
    1 2 3 2.000000e-01 # RLLE1x2x3
    1 3 1 2.000000e-01 # RLLE1x3x1
    1 3 2 2.000000e-01 # RLLE1x3x2
    1 3 3 2.000000e-01 # RLLE1x3x3
    2 1 1 -2.000000e-01 # RLLE2x1x1
    2 1 2 -2.000000e-01 # RLLE2x1x2
    2 1 3 -2.000000e-01 # RLLE2x1x3
    2 3 1 2.000000e-01 # RLLE2x3x1
    2 3 2 2.000000e-01 # RLLE2x3x2
    2 3 3 2.000000e-01 # RLLE2x3x3
    3 1 1 -2.000000e-01 # RLLE3x1x1
    3 1 2 -2.000000e-01 # RLLE3x1x2
    3 1 3 -2.000000e-01 # RLLE3x1x3
    3 2 1 -2.000000e-01 # RLLE3x2x1
    3 2 2 -2.000000e-01 # RLLE3x2x2
    3 2 3 -2.000000e-01 # RLLE3x2x3

###################################
## INFORMATION FOR RVLAMLQD
###################################
Block rvlamlqd
    1 1 1 1.000000e-01 # RLQD1x1x1
    1 1 2 1.000000e-01 # RLQD1x1x2
    1 1 3 1.000000e-01 # RLQD1x1x3
    1 2 1 1.000000e-01 # RLQD1x2x1
    1 2 2 1.000000e-01 # RLQD1x2x2
    1 2 3 1.000000e-01 # RLQD1x2x3
    1 3 1 1.000000e-01 # RLQD1x3x1
    1 3 2 1.000000e-01 # RLQD1x3x2
    1 3 3 1.000000e-01 # RLQD1x3x3
    2 1 1 1.000000e-01 # RLQD2x1x1
    2 1 2 1.000000e-01 # RLQD2x1x2
    2 1 3 1.000000e-01 # RLQD2x1x3
    2 2 1 1.000000e-01 # RLQD2x2x1
    2 2 2 1.000000e-01 # RLQD2x2x2
    2 2 3 1.000000e-01 # RLQD2x2x3
    2 3 1 1.000000e-01 # RLQD2x3x1
    2 3 2 1.000000e-01 # RLQD2x3x2
    2 3 3 1.000000e-01 # RLQD2x3x3
    3 1 1 1.000000e-01 # RLQD3x1x1
    3 1 2 1.000000e-01 # RLQD3x1x2
    3 1 3 1.000000e-01 # RLQD3x1x3
    3 2 1 1.000000e-01 # RLQD3x2x1
    3 2 2 1.000000e-01 # RLQD3x2x2
    3 2 3 1.000000e-01 # RLQD3x2x3
    3 3 1 1.000000e-01 # RLQD3x3x1
    3 3 2 1.000000e-01 # RLQD3x3x2
    3 3 3 1.000000e-01 # RLQD3x3x3

###################################
## INFORMATION FOR RVLAMUDD
###################################
Block rvlamudd
    1 1 2 2.000000e-01 # RLDD1x1x2
    1 1 3 2.000000e-01 # RLDD1x1x3
    1 2 1 -2.000000e-01 # RLDD1x2x1
    1 2 3 2.000000e-01 # RLDD1x2x3
    1 3 1 -2.000000e-01 # RLDD1x3x1
    1 3 2 -2.000000e-01 # RLDD1x3x2
    2 1 2 2.000000e-01 # RLDD2x1x2
    2 1 3 2.000000e-01 # RLDD2x1x3
    2 2 1 -2.000000e-01 # RLDD2x2x1
    2 2 3 2.000000e-01 # RLDD2x2x3
    2 3 1 -2.000000e-01 # RLDD2x3x1
    2 3 2 -2.000000e-01 # RLDD2x3x2
    3 1 2 2.000000e-01 # RLDD3x1x2
    3 1 3 2.000000e-01 # RLDD3x1x3
    3 2 1 -2.000000e-01 # RLDD3x2x1
    3 2 3 2.000000e-01 # RLDD3x2x3
    3 3 1 -2.000000e-01 # RLDD3x3x1
    3 3 2 -2.000000e-01 # RLDD3x3x2

###################################
## INFORMATION FOR RVTLLE
###################################
Block rvtlle
    1 2 1 2.000000e-01 # RTLE1x2x1
    1 2 2 2.000000e-01 # RTLE1x2x2
    1 2 3 2.000000e-01 # RTLE1x2x3
    1 3 1 2.000000e-01 # RTLE1x3x1
    1 3 2 2.000000e-01 # RTLE1x3x2
    1 3 3 2.000000e-01 # RTLE1x3x3
    2 1 1 -2.000000e-01 # RTLE2x1x1
    2 1 2 -2.000000e-01 # RTLE2x1x2
    2 1 3 -2.000000e-01 # RTLE2x1x3
    2 3 1 2.000000e-01 # RTLE2x3x1
    2 3 2 2.000000e-01 # RTLE2x3x2
    2 3 3 2.000000e-01 # RTLE2x3x3
    3 1 1 -2.000000e-01 # RTLE3x1x1
    3 1 2 -2.000000e-01 # RTLE3x1x2
    3 1 3 -2.000000e-01 # RTLE3x1x3
    3 2 1 -2.000000e-01 # RTLE3x2x1
    3 2 2 -2.000000e-01 # RTLE3x2x2
    3 2 3 -2.000000e-01 # RTLE3x2x3

###################################
## INFORMATION FOR RVTLQD
###################################
Block rvtlqd
    1 1 1 1.000000e-01 # RTQD1x1x1
    1 1 2 1.000000e-01 # RTQD1x1x2
    1 1 3 1.000000e-01 # RTQD1x1x3
    1 2 1 1.000000e-01 # RTQD1x2x1
    1 2 2 1.000000e-01 # RTQD1x2x2
    1 2 3 1.000000e-01 # RTQD1x2x3
    1 3 1 1.000000e-01 # RTQD1x3x1
    1 3 2 1.000000e-01 # RTQD1x3x2
    1 3 3 1.000000e-01 # RTQD1x3x3
    2 1 1 1.000000e-01 # RTQD2x1x1
    2 1 2 1.000000e-01 # RTQD2x1x2
    2 1 3 1.000000e-01 # RTQD2x1x3
    2 2 1 1.000000e-01 # RTQD2x2x1
    2 2 2 1.000000e-01 # RTQD2x2x2
    2 2 3 1.000000e-01 # RTQD2x2x3
    2 3 1 1.000000e-01 # RTQD2x3x1
    2 3 2 1.000000e-01 # RTQD2x3x2
    2 3 3 1.000000e-01 # RTQD2x3x3
    3 1 1 1.000000e-01 # RTQD3x1x1
    3 1 2 1.000000e-01 # RTQD3x1x2
    3 1 3 1.000000e-01 # RTQD3x1x3
    3 2 1 1.000000e-01 # RTQD3x2x1
    3 2 2 1.000000e-01 # RTQD3x2x2
    3 2 3 1.000000e-01 # RTQD3x2x3
    3 3 1 1.000000e-01 # RTQD3x3x1
    3 3 2 1.000000e-01 # RTQD3x3x2
    3 3 3 1.000000e-01 # RTQD3x3x3

###################################
## INFORMATION FOR RVTUDD
###################################
Block rvtudd
    1 1 2 2.000000e-01 # RTDD1x1x2
    1 1 3 2.000000e-01 # RTDD1x1x3
    1 2 1 -2.000000e-01 # RTDD1x2x1
    1 2 3 2.000000e-01 # RTDD1x2x3
    1 3 1 -2.000000e-01 # RTDD1x3x1
    1 3 2 -2.000000e-01 # RTDD1x3x2
    2 1 2 2.000000e-01 # RTDD2x1x2
    2 1 3 2.000000e-01 # RTDD2x1x3
    2 2 1 -2.000000e-01 # RTDD2x2x1
    2 2 3 2.000000e-01 # RTDD2x2x3
    2 3 1 -2.000000e-01 # RTDD2x3x1
    2 3 2 -2.000000e-01 # RTDD2x3x2
    3 1 2 2.000000e-01 # RTDD3x1x2
    3 1 3 2.000000e-01 # RTDD3x1x3
    3 2 1 -2.000000e-01 # RTDD3x2x1
    3 2 3 2.000000e-01 # RTDD3x2x3
    3 3 1 -2.000000e-01 # RTDD3x3x1
    3 3 2 -2.000000e-01 # RTDD3x3x2

###################################
## INFORMATION FOR SELMIX
###################################
Block selmix
    1 3 3.474325e-01 # RRl1x3
    1 6 9.377050e-01 # RRl1x6
    2 5 1.000000e+00 # RRl2x5
    3 4 1.000000e+00 # RRl3x4
    4 1 1.000000e+00 # RRl4x1
    5 2 1.000000e+00 # RRl5x2
    6 3 -9.377050e-01 # RRl6x3
    6 6 3.474325e-01 # RRl6x6

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.279340e+02 # aEWM1
    3 1.176000e-01 # aS

###################################
## INFORMATION FOR SNUMIX
###################################
Block snumix
    1 3 1.000000e+00 # RRn1x3
    2 2 1.000000e+00 # RRn2x2
    3 1 1.000000e+00 # RRn3x1

###################################
## INFORMATION FOR TD
###################################
Block td
    1 1 -2.551888e-01 # Rtd1x1
    2 2 -4.374642e+00 # Rtd2x2
    3 3 -1.677180e+02 # Rtd3x3

###################################
## INFORMATION FOR TE
###################################
Block te
    1 1 -2.155743e-02 # Rte1x1
    2 2 -4.457265e+00 # Rte2x2
    3 3 -7.436274e+01 # Rte3x3

###################################
## INFORMATION FOR TU
###################################
Block tu
    1 1 -8.900702e-03 # Rtu1x1
    2 2 -3.560256e+00 # Rtu2x2
    3 3 -5.887617e+02 # Rtu3x3

###################################
## INFORMATION FOR UMIX
###################################
Block umix
    1 1 -9.576948e-01 # RUU1x1
    1 2 2.877859e-01 # RUU1x2
    2 1 2.877859e-01 # RUU2x1
    2 2 9.576948e-01 # RUU2x2

###################################
## INFORMATION FOR UPMNS
###################################
Block upmns
    1 1 1.000000e+00 # RMNS1x1
    2 2 1.000000e+00 # RMNS2x2
    3 3 1.000000e+00 # RMNS3x3

###################################
## INFORMATION FOR USQMIX
###################################
Block usqmix
    1 3 5.111614e-01 # RRu1x3
    1 6 8.594848e-01 # RRu1x6
    2 3 -8.594848e-01 # RRu2x3
    2 6 5.111614e-01 # RRu2x6
    3 5 1.000000e+00 # RRu3x5
    4 4 1.000000e+00 # RRu4x4
    5 1 1.000000e+00 # RRu5x1
    6 2 1.000000e+00 # RRu6x2

###################################
## INFORMATION FOR VCKM
###################################
Block vckm
    1 1 9.738404e-01 # RCKM1x1
    1 2 2.271982e-01 # RCKM1x2
    1 3 2.173989e-03 # RCKM1x3
    2 1 -2.270868e-01 # RCKM2x1
    2 2 9.729587e-01 # RCKM2x2
    2 3 4.222469e-02 # RCKM2x3
    3 1 7.478279e-03 # RCKM3x1
    3 2 -4.161426e-02 # RCKM3x2
    3 3 9.991002e-01 # RCKM3x3

###################################
## INFORMATION FOR VMIX
###################################
Block vmix
    1 1 -9.899065e-01 # RVV1x1
    1 2 1.417221e-01 # RVV1x2
    2 1 1.417221e-01 # RVV2x1
    2 2 9.899065e-01 # RVV2x2

###################################
## INFORMATION FOR YD
###################################
Block yd
    1 1 1.884671e-04 # Ryd1x1
    2 2 3.230863e-03 # Ryd2x2
    3 3 1.360783e-01 # Ryd3x3

###################################
## INFORMATION FOR YE
###################################
Block ye
    1 1 2.882278e-05 # Rye1x1
    2 2 5.959638e-03 # Rye2x2
    3 3 1.002358e-01 # Rye3x3

###################################
## INFORMATION FOR YU
###################################
Block yu
    1 1 8.789391e-06 # Ryu1x1
    2 2 3.515756e-03 # Ryu2x2
    3 3 8.829349e-01 # Ryu3x3

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.463305e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.141000e+00 # WW
DECAY 25 3.518961e-03 # WH01
DECAY 35 1.389974e+00 # WH02
DECAY 36 1.742413e+00 # WA0
DECAY 37 1.176094e+00 # WH
DECAY 1000001 7.060640e+00 # Wsd1
DECAY 1000002 0.1 # Wsu1
DECAY 1000003 6.043630e-01 # Wsd2
DECAY 1000004 8.637400e+00 # Wsu2
DECAY 1000005 3.208630e-01 # Wsd3
DECAY 1000006 1.283740e+00 # Wsu3
DECAY 1000011 1.333692e+00 # Wsl1
DECAY 1000012 2.205383e+00 # Wsn1
DECAY 1000013 1.315543e+00 # Wsl2
DECAY 1000014 2.263504e+00 # Wsn2
DECAY 1000015 1.315543e+00 # Wsl3
DECAY 1000016 2.263504e+00 # Wsn3
DECAY 1000021 2.524968e+00 # Wgo
DECAY 1000023 4.107365e-05 # Wneu2
DECAY 1000024 1.056560e-03 # Wch1
DECAY 1000025 3.091347e+00 # Wneu3
DECAY 1000035 3.193022e+00 # Wneu4
DECAY 1000037 5.260644e+00 # Wch2
DECAY 2000001 3.208630e-01 # Wsd4
DECAY 2000002 1.283740e+00 # Wsu4
DECAY 2000003 6.304480e+00 # Wsd5
DECAY 2000004 6.376850e+00 # Wsu5
DECAY 2000005 6.304480e+00 # Wsd6
DECAY 2000006 6.376850e+00 # Wsu6
DECAY 2000011 2.332189e+00 # Wsl4
DECAY 2000013 2.332189e+00 # Wsl5
DECAY 2000015 2.187202e+00 # Wsl6
## 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 # tau- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
DECAY 1000022 0.000000 # n1 : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 1000022 # n1
        1 0 # 3 times electric charge
        2 2 # 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 1000023 # n2
        1 0 # 3 times electric charge
        2 2 # 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 1000025 # n3
        1 0 # 3 times electric charge
        2 2 # 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 1000035 # n4
        1 0 # 3 times electric charge
        2 2 # 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 1000024 # x1+
        1 3 # 3 times electric charge
        2 2 # 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 1000037 # x2+
        1 3 # 3 times electric charge
        2 2 # 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 1000021 # go
        1 0 # 3 times electric charge
        2 2 # number of spin states (2S+1)
        3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 0 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 35 # h02
        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)
Block QNUMBERS 36 # a0
        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)
Block QNUMBERS 37 # h+
        1 3 # 3 times electric charge
        2 1 # 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 1000012 # sv1
        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 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 1000014 # sv2
        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 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 1000016 # sv3
        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 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 1000011 # sl1-
        1 -3 # 3 times electric charge
        2 1 # 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 1000013 # sl2-
        1 -3 # 3 times electric charge
        2 1 # 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 1000015 # sl3-
        1 -3 # 3 times electric charge
        2 1 # 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 2000011 # sl4-
        1 -3 # 3 times electric charge
        2 1 # 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 2000013 # sl5-
        1 -3 # 3 times electric charge
        2 1 # 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 2000015 # sl6-
        1 -3 # 3 times electric charge
        2 1 # 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 1000002 # su1
        1 2 # 3 times electric charge
        2 1 # 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 1000004 # su2
        1 2 # 3 times electric charge
        2 1 # 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 1000006 # su3
        1 2 # 3 times electric charge
        2 1 # 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 2000002 # su4
        1 2 # 3 times electric charge
        2 1 # 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 2000004 # su5
        1 2 # 3 times electric charge
        2 1 # 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 2000006 # su6
        1 2 # 3 times electric charge
        2 1 # 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 1000001 # sd1
        1 -1 # 3 times electric charge
        2 1 # 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 1000003 # sd2
        1 -1 # 3 times electric charge
        2 1 # 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 1000005 # sd3
        1 -1 # 3 times electric charge
        2 1 # 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 2000001 # sd4
        1 -1 # 3 times electric charge
        2 1 # 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 2000003 # sd5
        1 -1 # 3 times electric charge
        2 1 # 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 2000005 # sd6
        1 -1 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # 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 *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Run to generate the grid pack *
#*********************************************************************
  .false. = gridpack !True = setting up the grid pack
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  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 = ebeam1 ! beam 1 total energy in GeV
     6500 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
        0 = polbeam1 ! beam polarization for beam 1
        0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
 'nn23lo1' = pdlabel ! PDF set
 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 F = fixed_ren_scale ! if .true. use fixed ren scale
 F = fixed_fac_scale ! if .true. use fixed fac scale
 91.1880 = scale ! fixed ren scale
 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 1 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Matching - Warning! ickkw > 1 is still beta
#*********************************************************************
 1 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
 1 = highestmult ! for ickkw=2, highest mult group
 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
 1 = alpsfact ! scale factor for QCD emission vx
 F = chcluster ! cluster only according to channel diag
 T = pdfwgt ! for ickkw=1, perform pdf reweighting
 5 = asrwgtflavor ! highest quark flavor for a_s reweight
 T = clusinfo ! include clustering tag in output
 3.0 = lhe_version ! Change the way clustering information pass to shower.
#*********************************************************************
#**********************************************************
#
#**********************************************************
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#**********************************************************
   T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#**********************************************************
#
#**********************************
# BW cutoff (M+/-bwcutoff*Gamma)
#**********************************
  15 = bwcutoff ! (M+/-bwcutoff*Gamma)
#**********************************************************
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#**********************************************************
   T = cut_decays ! Cut decay products
#*************************************************************
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#*************************************************************
   0 = nhel ! Number of helicities used per event
#*******************
# Standard Cuts
#*******************
#
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 20 = ptj ! minimum pt for the jets
  0 = ptb ! minimum pt for the b
 10 = pta ! minimum pt for the photons
 10 = ptl ! minimum pt for the charged leptons
  0 = misset ! minimum missing Et (sum of neutrino's momenta)
  0 = ptheavy ! minimum pt for one heavy final state
 1.0 = ptonium ! minimum pt for the quarkonium states
 -1 = ptjmax ! maximum pt for the jets
 -1 = ptbmax ! maximum pt for the b
 -1 = ptamax ! maximum pt for the photons
 -1 = ptlmax ! maximum pt for the charged leptons
 -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0 = ej ! minimum E for the jets
  0 = eb ! minimum E for the b
  0 = ea ! minimum E for the photons
  0 = el ! minimum E for the charged leptons
 -1 = ejmax ! maximum E for the jets
 -1 = ebmax ! maximum E for the b
 -1 = eamax ! maximum E for the photons
 -1 = elmax ! maximum E for the charged leptons
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
   5 = etaj ! max rap for the jets
  -1 = etab ! max rap for the b
 2.5 = etaa ! max rap for the photons
 2.5 = etal ! max rap for the charged leptons
 0.6 = etaonium ! max rap for the quarkonium states
   0 = etajmin ! min rap for the jets
   0 = etabmin ! min rap for the b
   0 = etaamin ! min rap for the photons
   0 = etalmin ! main rap for the charged leptons
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.4 = drjj ! min distance between jets
 0 = drbb ! min distance between b's
 0.4 = drll ! min distance between leptons
 0.4 = draa ! min distance between gammas
 0 = drbj ! min distance between b and jet
 0.4 = draj ! min distance between gamma and jet
 0.4 = drjl ! min distance between jet and lepton
 0 = drab ! min distance between gamma and b
 0 = drbl ! min distance between b and lepton
 0.4 = dral ! min distance between gamma and lepton
 -1 = drjjmax ! max distance between jets
 -1 = drbbmax ! max distance between b's
 -1 = drllmax ! max distance between leptons
 -1 = draamax ! max distance between gammas
 -1 = drbjmax ! max distance between b and jet
 -1 = drajmax ! max distance between gamma and jet
 -1 = drjlmax ! max distance between jet and lepton
 -1 = drabmax ! max distance between gamma and b
 -1 = drblmax ! max distance between b and lepton
 -1 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
 0 = mmjj ! min invariant mass of a jet pair
 0 = mmbb ! min invariant mass of a b pair
 0 = mmaa ! min invariant mass of gamma gamma pair
 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1 = mmjjmax ! max invariant mass of a jet pair
 -1 = mmbbmax ! max invariant mass of a b pair
 -1 = mmaamax ! max invariant mass of gamma gamma pair
 -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0 = xptj ! minimum pt for at least one jet
 0 = xptb ! minimum pt for at least one b
 0 = xpta ! minimum pt for at least one photon
 0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0 = ptj1min ! minimum pt for the leading jet in pt
 0 = ptj2min ! minimum pt for the second jet in pt
 0 = ptj3min ! minimum pt for the third jet in pt
 0 = ptj4min ! minimum pt for the fourth jet in pt
 -1 = ptj1max ! maximum pt for the leading jet in pt
 -1 = ptj2max ! maximum pt for the second jet in pt
 -1 = ptj3max ! maximum pt for the third jet in pt
 -1 = ptj4max ! maximum pt for the fourth jet in pt
 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
 0 = ptl1min ! minimum pt for the leading lepton in pt
 0 = ptl2min ! minimum pt for the second lepton in pt
 0 = ptl3min ! minimum pt for the third lepton in pt
 0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1 = ptl1max ! maximum pt for the leading lepton in pt
 -1 = ptl2max ! maximum pt for the second lepton in pt
 -1 = ptl3max ! maximum pt for the third lepton in pt
 -1 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1 = htjmax ! maximum jet HT=Sum(jet pt)
 0 = ihtmin !inclusive Ht for all partons (including b)
 -1 = ihtmax !inclusive Ht for all partons (including b)
 0 = ht2min ! minimum Ht for the two leading jets
 0 = ht3min ! minimum Ht for the three leading jets
 0 = ht4min ! minimum Ht for the four leading jets
 -1 = ht2max ! maximum Ht for the two leading jets
 -1 = ht3max ! maximum Ht for the three leading jets
 -1 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
   0 = ptgmin ! Min photon transverse momentum
 0.4 = R0gamma ! Radius of isolation code
 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
 .true. = isoEM ! isolate photons from EM energy (photons and leptons)
#*********************************************************************
# WBF cuts *
#*********************************************************************
 0 = xetamin ! minimum rapidity for two jets in the WBF case
 0 = deltaeta ! minimum rapidity for two jets in the WBF case
#*********************************************************************
# KT DURHAM CUT *
#*********************************************************************
 -1 = ktdurham
 0.4 = dparameter
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
# Jet measure cuts *
#*********************************************************************
 0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: If use_syst is T, matched Pythia output is *
# meaningful ONLY if plotted taking matchscale *
# reweighting into account! *
#*********************************************************************
   F = use_syst ! Enable systematics studies
#
#**************************************
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#**************************************
#
0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
30 50 = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
CT10nlo.LHgrid = sys_pdf # matching scales
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf

any recommendation?

thank you so much.
Alejandro

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Here is a patch that seems to improve the situation.
I would strongly suggest that you check the color flow associate to each events by running the following code

./bin/madevent check_events PATH_TO_LHE_FILE

I will ask some other author of MG5_aMC to check my patch. But if you can try it and run the above check and tell us the result it will be very usefull.

Cheers,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

Just fixing the patch, which was wrongly created

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

looks like a second error flow error occurs in top of the one fixed.
So this patch is not enough

Cheers,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Here is a new patch,

I did not had the time to fully test it, but since I will not have the time before next week and that the current test are sucesfull, I think that you can be interested.

Cheers,

Olivier

[alefisico (alefisico)]

Hi

I tested it and it works, it creates the gridpack. However, when I am try to generate the events with the run.sh script, I am getting this error:
Traceback (most recent call last):
  File "./bin/gridrun", line 84, in <module>
    cmd_line = cmd_interface.GridPackCmd(me_dir=root_path, nb_event=args[0], seed=args[1])
  File "/afs/cern.ch/work/a/algomez/Generation/tmp/CMSSW_7_1_14/src/process/madevent/bin/internal/madevent_interface.py", line 4553, in __init__
    self.launch(nb_event, seed)
  File "/afs/cern.ch/work/a/algomez/Generation/tmp/CMSSW_7_1_14/src/process/madevent/bin/internal/madevent_interface.py", line 4575, in launch
    self.refine4grid(nb_event)
  File "/afs/cern.ch/work/a/algomez/Generation/tmp/CMSSW_7_1_14/src/process/madevent/bin/internal/madevent_interface.py", line 4619, in refine4grid
    cwd=Pdir)
  File "/afs/cern.ch/work/a/algomez/Generation/tmp/CMSSW_7_1_14/src/process/madevent/bin/internal/misc.py", line 551, in deco_f
    % arg
Exception: ['./../../bin/internal/gen_ximprove'] fails with no such file or directory
Error: event file not found !

I checked here https://bugs.launchpad.net/mg5amcnlo/+bug/1449125
that you had a similar bug before, but It does not work for me. Do you have any idea?

thanks for all your help and time!

cheers,

Alejandro

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

I have found the problem (my previous fix was not fine in all situation and one more pathological case was present)

I will not provide a patch, but directly the new file. You have to put it in madgraph/iolibs/template_files (with name addmothers.f)

I have check it and succeeds to generate events without any error for the cards that you send to me.
I know that you have a more complex process, can you check it as well?

Cheers,

Olivier

[alefisico (alefisico)]

Hi Olivier

Thanks for all your help, but I still cannot fix the issue. Can you please give me the steps that you are doing? In case I am doing something wrong.

In any case I found two issues:

- when I try to run the creation of the gridpack in a cluster (lsf but I also try it in condor) I got this error:
generate_events pilotrun
Traceback (most recent call last):
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 879, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/madevent_interface.py", line 1983, in do_generate_events
    postcmd=False)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 919, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/madevent_interface.py", line 2654, in do_survey
    self.monitor(run_type='All jobs submitted for survey', html=True)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/madevent_interface.py", line 3763, in monitor
    self.cluster.wait(self.me_dir, update_status, update_first=update_first)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/cluster.py", line 57, in deco_f_interupt
    return f(self, *args, **opt)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/cluster.py", line 304, in wait
    idle, run, finish, fail = self.control(me_dir)
  File "/afs/cern.ch/work/a/algomez/Generation/genproductions/bin/MadGraph5_aMCatNLO/RPVSt100_UDD312_ISR2j_LO/RPVSt100_UDD312_ISR2j_LO_gridpack/work/processtmp/bin/internal/misc.py", line 244, in deco_f_retry
    raise error.__class__, '[Fail %i times] \n %s ' % (i+1, error)
ClusterManagmentError: [Fail 5 times]
 Job exhausted all retries.

- In the case that I try to run locally, it seems like it is running fine. I got a gridpack at the end, but I found these messages in t...