Error: Color combination 1 2 -> 0 not implemented in addmothers.f

Would you mind to check with the 2.3 version of the code?

https://launchpad.net/mg5amcnlo/2.0/2.2.0/+download/MG5_aMC_v2.3.0.beta.tar.gz

I’m not sure but i think that I solve one similar problem since the release of 2.2.2

Cheers,

Olivier

On 19 Apr 2015, at 18:01, alefisico <email address hidden> wrote:

> New question #265473 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/265473
>
> 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
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi

well I managed to produce this lhe file. I think I was including some process in the calculation that I shouldn't. I run it in 2.2 and it works fine.

However, just for the sake of this question, I just run the same process in 2.3 and I am getting other error:

enerate_events RPVStop100tojj_upto2j_13TeV --cluster -f
Traceback (most recent call last):
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/extended_cmd.py", line 879, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/madevent_interface.py", line 1996, in do_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/interface/madevent_interface.py", line 2740, in do_refine
    x_improve.launch() # create the ajob for the refinment.
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/madevent/gen_ximprove.py", line 763, in launch
    self.results = sum_html.collect_result(self.cmd, None)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/madevent/sum_html.py", line 672, in collect_result
    P_comb.add_results(name, pjoin(P_path,name,'results.dat'), mfactor)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/madevent/sum_html.py", line 398, in add_results
    oneresult.read_results(filepath)
  File "/cms/gomez/MG5_aMC_v2_3_0_beta/madgraph/madevent/sum_html.py", line 267, in read_results
    finput = open(filepath)
IOError: [Errno 2] No such file or directory: '/cms/gomez/MG5_aMC_v2_3_0_beta/RPVStop_UDD312_ISR_upto2j/SubProcesses/P2_qq_su1su1xgg_su1_qq_su1x_qq/G1/results.dat'

I dont know if it helps.

cheers,
Alejandro