madspin for heavy resonance decay

Hi,

Could you check if you still have this bug with 2.3.3, Your version seems quite old.

Cheers,

Olivier
> On Jan 21, 2016, at 06:08, toodles <email address hidden> wrote:
>
> New question #281261 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/281261
>
> I am now using MG5 and mssm-full model to generate h2 in 5FS, then h2 decay to ttbar and then jj+lvl.
>
> The pattern I am using in madspin is : decay h2 > t t~, (t > w+ b, w+ > j j), (t~ > w- b~, w- > l- vl~)
>
> However I got this error. I copied the run_01_tag_1_debug.log content in the following:
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.2.1.2 *
> #* *
> #* 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 *
> #* *
> #************************************************************
> launch
> Traceback (most recent call last):
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 879, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
> return func(arg, **opt)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2139, in do_launch
> self.do_generate_events(line, *args, **opt)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2128, in do_generate_events
> self.exec_cmd('decay_events -from_cards', postcmd=False)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 919, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
> return func(arg, **opt)
> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/common_run_interface.py", line 1554, in do_decay_events
> madspin_cmd.import_command_file(path)
> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 1035, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
> return func(arg, **opt)
> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/interface_madspin.py", line 461, in do_launch
> generate_all.run()
> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2015, in run
> decay_mapping = self.get_identical_decay()
> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2360, in get_identical_decay
> True, self.pid2width, self.pid2mass, BW_cut,self.Ecollider)
> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 460, in generate_momenta
> tree = self["tree"][res]
> KeyError: -4
> Run Options
> -----------
> stdout_level : None
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> cluster_temp_path : None
> cluster_queue : madgraph
> cluster_memory : None
> cluster_time : None
> nb_core : 8 (user set)
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_status_update : (600, 30)
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> pythia-pgs_path : /home/toodles/madness/MG5_aMC_v2_1_2/pythia-pgs (user set)
> td_path : /home/toodles/madness/MG5_aMC_v2_1_2/td (user set)
> delphes_path : /home/toodles/madness/MG5_aMC_v2_1_2/Delphes (user set)
> thepeg_path : None (user set)
> cluster_type : condor
> exrootanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/ExRootAnalysis (user set)
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/MadAnalysis (user set)
> lhapdf : lhapdf-config
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> cluster_nb_retry : 1
> timeout : 60
> cpp_compiler : None
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.1.2 2014-07-03 *
> #* *
> #* 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 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 mssm-full
> define p = p b b~
> generate p p > h2
> output test
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR HMIX
> ###################################
> BLOCK HMIX Q= 1.000000e+00 #
> 1 3.576810e+02 # rmuh
> 2 9.748624e+00 # tb
> 4 1.664391e+05 # ma2
> 3 2.449353e+02 # higgs vev(Q) MSSM DRb
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 5 4.889917e+00 # mb
> 6 1.750000e+02 # mt
> 15 1.777000e+00 # mta
> 23 9.118760e+01 # mz
> 24 7.982901e+01 # mw
> 25 1.108991e+02 # mh01
> 35 3.999601e+02 # mh02
> 36 3.995839e+02 # ma0
> 37 4.078790e+02 # mh
> 1000001 5.684411e+02 # msd1
> 1000002 5.611190e+02 # msu1
> 1000003 5.684411e+02 # msd2
> 1000004 5.611190e+02 # msu2
> 1000005 5.130652e+02 # msd3
> 1000006 3.996685e+02 # msu3
> 1000011 2.029157e+02 # msl1
> 1000012 1.852583e+02 # msn1
> 1000013 2.029157e+02 # msl2
> 1000014 1.852583e+02 # msn2
> 1000015 1.344909e+02 # msl3
> 1000016 1.847085e+02 # msn3
> 1000021 6.077137e+02 # mgo
> 1000022 9.668807e+01 # mneu1
> 1000023 1.810882e+02 # mneu2
> 1000024 1.816965e+02 # mch1
> 1000025 -3.637560e+02 # mneu3
> 1000035 3.817294e+02 # mneu4
> 1000037 3.799393e+02 # mch2
> 2000001 5.452285e+02 # msd4
> 2000002 5.492593e+02 # msu4
> 2000003 5.452285e+02 # msd5
> 2000004 5.492593e+02 # msu5
> 2000005 5.437267e+02 # msd6
> 2000006 5.857858e+02 # msu6
> 2000011 1.441028e+02 # msl4
> 2000013 1.441028e+02 # msl5
> 2000015 2.068678e+02 # msl6
> 1 0.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> ###################################
> ## INFORMATION FOR MSOFT
> ###################################
> BLOCK MSOFT Q= 1.000000e+00 #
> 1 1.013965e+02 # rmx1
> 2 1.915042e+02 # rmx2
> 3 5.882630e+02 # rmx3
> 21 3.233749e+04 # mhd2
> 22 -1.288001e+05 # mhu2
> 31 1.953348e+02 # rml21x1
> 32 1.953348e+02 # rml22x2
> 33 1.944960e+02 # rml23x3
> 34 1.364941e+02 # rme21x1
> 35 1.364941e+02 # rme22x2
> 36 1.340434e+02 # rme23x3
> 41 5.475735e+02 # rmq21x1
> 42 5.475735e+02 # rmq22x2
> 43 4.987639e+02 # rmq23x3
> 44 5.295112e+02 # rmu21x1
> 45 5.295112e+02 # rmu22x2
> 46 4.232459e+02 # rmu23x3
> 47 5.231488e+02 # rmd21x1
> 48 5.231488e+02 # rmd22x2
> 49 5.198673e+02 # rmd23x3
> ###################################
> ## INFORMATION FOR NMIX
> ###################################
> BLOCK NMIX #
> 1 1 9.863644e-01 # rnn1x1
> 1 2 -5.311036e-02 # rnn1x2
> 1 3 1.464340e-01 # rnn1x3
> 1 4 -5.311861e-02 # rnn1x4
> 2 1 9.935054e-02 # rnn2x1
> 2 2 9.449493e-01 # rnn2x2
> 2 3 -2.698467e-01 # rnn2x3
> 2 4 1.561507e-01 # rnn2x4
> 3 1 -6.033880e-02 # rnn3x1
> 3 2 8.770049e-02 # rnn3x2
> 3 3 6.958775e-01 # rnn3x3
> 3 4 7.102270e-01 # rnn3x4
> 4 1 -1.165071e-01 # rnn4x1
> 4 2 3.107390e-01 # rnn4x2
> 4 3 6.492260e-01 # rnn4x3
> 4 4 -6.843778e-01 # rnn4x4
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.279340e+02 # aewm1
> 3 1.180000e-01 # as
> 6 1.750000e+02 # mt
> 7 1.777000e+00 # mta
> 4 9.118760e+01 # mz
> 2 1.166370e-05 # G_F [GeV^-2]
> ###################################
> ## INFORMATION FOR UMIX
> ###################################
> BLOCK UMIX #
> 1 1 9.168349e-01 # ruu1x1
> 1 2 -3.992666e-01 # ruu1x2
> 2 1 3.992666e-01 # ruu2x1
> 2 2 9.168349e-01 # ruu2x2
> ###################################
> ## INFORMATION FOR VMIX
> ###################################
> BLOCK VMIX #
> 1 1 9.725578e-01 # rvv1x1
> 1 2 -2.326612e-01 # rvv1x2
> 2 1 2.326612e-01 # rvv2x1
> 2 2 9.725578e-01 # rvv2x2
> ###################################
> ## INFORMATION FOR YD
> ###################################
> BLOCK YD Q= 1.000000e+00 #
> 3 3 1.388402e-01 # ryd3x3
> ###################################
> ## INFORMATION FOR YE
> ###################################
> BLOCK YE Q= 1.000000e+00 #
> 3 3 1.008908e-01 # rye3x3
> ###################################
> ## INFORMATION FOR YU
> ###################################
> BLOCK YU Q= 1.000000e+00 #
> 3 3 8.928445e-01 # ryu3x3
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.561950e+00 # wt
> DECAY 23 2.411433e+00 # wz
> DECAY 24 2.002822e+00 # ww
> DECAY 25 1.986108e-03 # wh01
> DECAY 35 5.748014e-01 # wh02
> DECAY 36 6.321785e-01 # wa0
> DECAY 37 5.469628e-01 # wh
> DECAY 1000001 5.312788e+00 # wsd1
> DECAY 1000002 5.477195e+00 # wsu1
> DECAY 1000003 5.312788e+00 # wsd2
> DECAY 1000004 5.477195e+00 # wsu2
> DECAY 1000005 3.736276e+00 # wsd3
> DECAY 1000006 2.021596e+00 # wsu3
> DECAY 1000011 2.136822e-01 # wsl1
> DECAY 1000012 1.498816e-01 # wsn1
> DECAY 1000013 2.136822e-01 # wsl2
> DECAY 1000014 1.498816e-01 # wsn2
> DECAY 1000015 1.483273e-01 # wsl3
> DECAY 1000016 1.475190e-01 # wsn3
> DECAY 1000021 5.506754e+00 # wgo
> DECAY 1000023 2.077700e-02 # wneu2
> DECAY 1000024 1.704145e-02 # wch1
> DECAY 1000025 1.915985e+00 # wneu3
> DECAY 1000035 2.585851e+00 # wneu4
> DECAY 1000037 2.486895e+00 # wch2
> DECAY 2000001 2.858123e-01 # wsd4
> DECAY 2000002 1.152973e+00 # wsu4
> DECAY 2000003 2.858123e-01 # wsd5
> DECAY 2000004 1.152973e+00 # wsu5
> DECAY 2000005 8.015663e-01 # wsd6
> DECAY 2000006 7.373133e+00 # wsu6
> DECAY 2000011 2.161216e-01 # wsl4
> DECAY 2000013 2.161216e-01 # wsl5
> DECAY 2000015 2.699061e-01 # wsl6
> DECAY 1 0.000000e+00 # d : 0.0
> DECAY 2 0.000000e+00 # u : 0.0
> DECAY 3 0.000000e+00 # s : 0.0
> DECAY 4 0.000000e+00 # c : 0.0
> DECAY 5 0.000000e+00 # b : 0.0
> DECAY 11 0.000000e+00 # e- : 0.0
> DECAY 12 0.000000e+00 # ve : 0.0
> DECAY 13 0.000000e+00 # mu- : 0.0
> DECAY 14 0.000000e+00 # vm : 0.0
> DECAY 15 0.000000e+00 # ta- : 0.0
> DECAY 16 0.000000e+00 # vt : 0.0
> DECAY 21 0.000000e+00 # g : 0.0
> DECAY 22 0.000000e+00 # a : 0.0
> DECAY 1000022 0.000000e+00 # n1 : 0.0
> ###################################
> ## INFORMATION FOR QNUMBERS 1000022
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000023
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000025
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000035
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000024
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000037
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000021
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 35
> ###################################
> BLOCK QNUMBERS 35 # h2
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 36
> ###################################
> BLOCK QNUMBERS 36 # h3
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 37
> ###################################
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000012
> ###################################
> BLOCK QNUMBERS 1000012 # sve
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000014
> ###################################
> BLOCK QNUMBERS 1000014 # svm
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000016
> ###################################
> BLOCK QNUMBERS 1000016 # svt
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000011
> ###################################
> BLOCK QNUMBERS 1000011 # el-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000013
> ###################################
> BLOCK QNUMBERS 1000013 # mul-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000015
> ###################################
> BLOCK QNUMBERS 1000015 # ta1-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000011
> ###################################
> BLOCK QNUMBERS 2000011 # er-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000013
> ###################################
> BLOCK QNUMBERS 2000013 # mur-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000015
> ###################################
> BLOCK QNUMBERS 2000015 # ta2-
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000002
> ###################################
> BLOCK QNUMBERS 1000002 # ul
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000004
> ###################################
> BLOCK QNUMBERS 1000004 # cl
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000006
> ###################################
> BLOCK QNUMBERS 1000006 # t1
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000002
> ###################################
> BLOCK QNUMBERS 2000002 # ur
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000004
> ###################################
> BLOCK QNUMBERS 2000004 # cr
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000006
> ###################################
> BLOCK QNUMBERS 2000006 # t2
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000001
> ###################################
> BLOCK QNUMBERS 1000001 # dl
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000003
> ###################################
> BLOCK QNUMBERS 1000003 # sl
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 1000005
> ###################################
> BLOCK QNUMBERS 1000005 # b1
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000001
> ###################################
> BLOCK QNUMBERS 2000001 # dr
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000003
> ###################################
> BLOCK QNUMBERS 2000003 # sr
> 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)
> ###################################
> ## INFORMATION FOR QNUMBERS 2000005
> ###################################
> BLOCK QNUMBERS 2000005 # b2
> 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)
> ###################################
> ## INFORMATION FOR MODSEL
> ###################################
> BLOCK MODSEL #
> 1 1 #
> ###################################
> ## INFORMATION FOR STOPMIX
> ###################################
> BLOCK STOPMIX #
> 1 1 5.536450e-01 # rru3x3
> 1 2 8.327528e-01 # rru3x6
> 2 1 8.327528e-01 # rru6x3
> 2 2 -5.536450e-01 # rru6x6
> ###################################
> ## INFORMATION FOR SBOTMIX
> ###################################
> BLOCK SBOTMIX #
> 1 1 9.387379e-01 # rrd3x3
> 1 2 3.446319e-01 # rrd3x6
> 2 1 -3.446319e-01 # rrd6x3
> 2 2 9.387379e-01 # rrd6x6
> ###################################
> ## INFORMATION FOR STAUMIX
> ###################################
> BLOCK STAUMIX #
> 1 1 2.824872e-01 # rrl3x3
> 1 2 9.592711e-01 # rrl3x6
> 2 1 9.592711e-01 # rrl6x3
> 2 2 -2.824872e-01 # rrl6x6
> ###################################
> ## INFORMATION FOR ALPHA
> ###################################
> BLOCK ALPHA #
> -1.138252e-01 # alp
> ###################################
> ## INFORMATION FOR AE
> ###################################
> BLOCK AE Q= 1.000000e+00 #
> 3 3 -2.517769e+02 # A_tau(Q) DRbar
> 1 1 0.000000e+00 # A_e(Q) DRbar
> 2 2 0.000000e+00 # A_mu(Q) DRbar
> ###################################
> ## INFORMATION FOR AU
> ###################################
> BLOCK AU Q= 1.000000e+00 #
> 3 3 -4.981299e+02 # A_t(Q) DRbar
> 1 1 0.000000e+00 # A_u(Q) DRbar
> 2 2 0.000000e+00 # A_c(Q) DRbar
> ###################################
> ## INFORMATION FOR AD
> ###################################
> BLOCK AD Q= 1.000000e+00 #
> 3 3 -7.972741e+02 # A_b(Q) DRbar
> 1 1 0.000000e+00 # A_d(Q) DRbar
> 2 2 0.000000e+00 # A_s(Q) DRbar
> #*********************************************************************
> # 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. *
> #*********************************************************************
> 1000 = 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
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> F = fixed_ren_scale ! if .true. use fixed ren scale
> F = fixed_fac_scale ! if .true. use fixed fac scale
> 91.1880 = scale ! fixed ren scale
> 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> 1 = scalefact ! scale factor for event-by-event scales
> #*********************************************************************
> # Matching - Warning! ickkw > 1 is still beta
> #*********************************************************************
> 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
> 1 = highestmult ! for ickkw=2, highest mult group
> 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
> 1 = alpsfact ! scale factor for QCD emission vx
> F = chcluster ! cluster only according to channel diag
> T = pdfwgt ! for ickkw=1, perform pdf reweighting
> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
> T = clusinfo ! include clustering tag in output
> 3.0 = lhe_version ! Change the way clustering information pass to shower.
> #*********************************************************************
> #**********************************************************
> #
> #**********************************************************
> # Automatic ptj and mjj cuts if xqcut > 0
> # (turn off for VBF and single top processes)
> #**********************************************************
> T = auto_ptj_mjj ! Automatic setting of ptj and mjj
> #**********************************************************
> #
> #**********************************
> # BW cutoff (M+/-bwcutoff*Gamma)
> #**********************************
> 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #**********************************************************
> # Apply pt/E/eta/dr/mij cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #**********************************************************
> T = cut_decays ! Cut decay products
> #*************************************************************
> # Number of helicities to sum per event (0 = all helicities)
> # 0 gives more stable result, but longer run time (needed for
> # long decay chains e.g.).
> # Use >=2 if most helicities contribute, e.g. pure QCD.
> #*************************************************************
> 0 = nhel ! Number of helicities used per event
> #*******************
> # Standard Cuts
> #*******************
> #
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20 = ptj ! minimum pt for the jets
> 0 = ptb ! minimum pt for the b
> 10 = pta ! minimum pt for the photons
> 10 = ptl ! minimum pt for the charged leptons
> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
> 0 = ptheavy ! minimum pt for one heavy final state
> 1.0 = ptonium ! minimum pt for the quarkonium states
> -1 = ptjmax ! maximum pt for the jets
> -1 = ptbmax ! maximum pt for the b
> -1 = ptamax ! maximum pt for the photons
> -1 = ptlmax ! maximum pt for the charged leptons
> -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> #*********************************************************************
> # Minimum and maximum E's (in the center of mass frame) *
> #*********************************************************************
> 0 = ej ! minimum E for the jets
> 0 = eb ! minimum E for the b
> 0 = ea ! minimum E for the photons
> 0 = el ! minimum E for the charged leptons
> -1 = ejmax ! maximum E for the jets
> -1 = ebmax ! maximum E for the b
> -1 = eamax ! maximum E for the photons
> -1 = elmax ! maximum E for the charged leptons
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 5 = etaj ! max rap for the jets
> -1 = etab ! max rap for the b
> 2.5 = etaa ! max rap for the photons
> 2.5 = etal ! max rap for the charged leptons
> 0.6 = etaonium ! max rap for the quarkonium states
> 0 = etajmin ! min rap for the jets
> 0 = etabmin ! min rap for the b
> 0 = etaamin ! min rap for the photons
> 0 = etalmin ! main rap for the charged leptons
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.4 = drjj ! min distance between jets
> 0 = drbb ! min distance between b's
> 0.4 = drll ! min distance between leptons
> 0.4 = draa ! min distance between gammas
> 0 = drbj ! min distance between b and jet
> 0.4 = draj ! min distance between gamma and jet
> 0.4 = drjl ! min distance between jet and lepton
> 0 = drab ! min distance between gamma and b
> 0 = drbl ! min distance between b and lepton
> 0.4 = dral ! min distance between gamma and lepton
> -1 = drjjmax ! max distance between jets
> -1 = drbbmax ! max distance between b's
> -1 = drllmax ! max distance between leptons
> -1 = draamax ! max distance between gammas
> -1 = drbjmax ! max distance between b and jet
> -1 = drajmax ! max distance between gamma and jet
> -1 = drjlmax ! max distance between jet and lepton
> -1 = drabmax ! max distance between gamma and b
> -1 = drblmax ! max distance between b and lepton
> -1 = dralmax ! maxdistance between gamma and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> # WARNING: for four lepton final state mmll cut require to have *
> # different lepton masses for each flavor! *
> #*********************************************************************
> 0 = mmjj ! min invariant mass of a jet pair
> 0 = mmbb ! min invariant mass of a b pair
> 0 = mmaa ! min invariant mass of gamma gamma pair
> 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1 = mmjjmax ! max invariant mass of a jet pair
> -1 = mmbbmax ! max invariant mass of a b pair
> -1 = mmaamax ! max invariant mass of gamma gamma pair
> -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons *
> #*********************************************************************
> 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0 = xptj ! minimum pt for at least one jet
> 0 = xptb ! minimum pt for at least one b
> 0 = xpta ! minimum pt for at least one photon
> 0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0 = ptj1min ! minimum pt for the leading jet in pt
> 0 = ptj2min ! minimum pt for the second jet in pt
> 0 = ptj3min ! minimum pt for the third jet in pt
> 0 = ptj4min ! minimum pt for the fourth jet in pt
> -1 = ptj1max ! maximum pt for the leading jet in pt
> -1 = ptj2max ! maximum pt for the second jet in pt
> -1 = ptj3max ! maximum pt for the third jet in pt
> -1 = ptj4max ! maximum pt for the fourth jet in pt
> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
> #*********************************************************************
> # Control the pt's of leptons sorted by pt *
> #*********************************************************************
> 0 = ptl1min ! minimum pt for the leading lepton in pt
> 0 = ptl2min ! minimum pt for the second lepton in pt
> 0 = ptl3min ! minimum pt for the third lepton in pt
> 0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1 = ptl1max ! maximum pt for the leading lepton in pt
> -1 = ptl2max ! maximum pt for the second lepton in pt
> -1 = ptl3max ! maximum pt for the third lepton in pt
> -1 = ptl4max ! maximum pt for the fourth lepton in pt
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1 = htjmax ! maximum jet HT=Sum(jet pt)
> 0 = ihtmin !inclusive Ht for all partons (including b)
> -1 = ihtmax !inclusive Ht for all partons (including b)
> 0 = ht2min ! minimum Ht for the two leading jets
> 0 = ht3min ! minimum Ht for the three leading jets
> 0 = ht4min ! minimum Ht for the four leading jets
> -1 = ht2max ! maximum Ht for the two leading jets
> -1 = ht3max ! maximum Ht for the three leading jets
> -1 = ht4max ! maximum Ht for the four leading jets
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442 *
> # When ptgmin=0, all the other parameters are ignored *
> # When ptgmin>0, pta and draj are not going to be used *
> #***********************************************************************
> 0 = ptgmin ! Min photon transverse momentum
> 0.4 = R0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> .true. = isoEM ! isolate photons from EM energy (photons and leptons)
> #*********************************************************************
> # WBF cuts *
> #*********************************************************************
> 0 = xetamin ! minimum rapidity for two jets in the WBF case
> 0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #*********************************************************************
> # KT DURHAM CUT *
> #*********************************************************************
> -1 = ktdurham
> 0.4 = dparameter
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> # Jet measure cuts *
> #*********************************************************************
> 0 = xqcut ! minimum kt jet measure between partons
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: If use_syst is T, matched Pythia output is *
> # meaningful ONLY if plotted taking matchscale *
> # reweighting into account! *
> #*********************************************************************
> F = use_syst ! Enable systematics studies
> #
> #**************************************
> # Parameter of the systematics study
> # will be used by SysCalc (if installed)
> #**************************************
> #
> 0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
> 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
> 30 50 = sys_matchscale # variation of merging scale
> # PDF sets and number of members (0 or none for all members).
> CT10nlo.LHgrid = sys_pdf # matching scales
> # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
>
> BTW, if I change the decay pattern to be
> decay h2 > t t~, t > w+ b, w+ > j j, t~ > w- b~, w- > l- vl~
> there is no such error anymore. But the decay get stuck at "Decaying the events.....".
>
> I am really appreciate for your help
>
> Regards,
> Ying-Ying
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi, Olivier

Actually when I use 2.3.3 and use mssm-full, I cannot do output.

Here is the error massage:

Command "output test" interrupted with error:

IOError : [Errno 2] No such file or directory: '/Users/toodles/madness/MG5_aMC_v2_3_3/test/Source/MODEL/param_card_rule.dat'

here is the MG5_debug file:

import model mssm-full
define p = p b b~
generate p p > h2
output test
Traceback (most recent call last):
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 908, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 897, in onecmd_orig
    return func(arg, **opt)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/master_interface.py", line 284, in do_output
    self.cmd.do_output(self, line, *args, **opts)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/madgraph_interface.py", line 6687, in do_output
    self.finalize(nojpeg, flaglist=flaglist)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/madgraph_interface.py", line 6984, in finalize
    wanted_couplings)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 3188, in convert_model_to_mg4
    wanted_lorentz, wanted_couplings)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 753, in convert_model_to_mg4
    model_builder.build(wanted_couplings)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 5149, in build
    self.write_all()
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 5198, in write_all
    self.create_param_card()
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 6248, in create_param_card
    mssm_convert=True)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 6239, in create_param_card_static
    translator.make_valid_param_card(output_path, rule_card_path)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/models/check_param_card.py", line 1461, in make_valid_param_card
    cardrule.load_rule(restrictpath)
  File "/Users/toodles/madness/MG5_aMC_v2_3_3/models/check_param_card.py", line 882, in load_rule
    tree = ET.parse(inputpath)
  File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/xml/etree/ElementTree.py", line 1182, in parse
    tree.parse(source, parser)
  File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/xml/etree/ElementTree.py", line 647, in parse
    source = open(source, "rb")
IOError: [Errno 2] No such file or directory: '/Users/toodles/madness/MG5_aMC_v2_3_3/test/Source/MODEL/param_card_rule.dat'

So I cannot try mad spin in mssm-full model with 2.3.3

Thanks for your reply.

Regards,

Ying-Ying

At 2016-01-21 16:13:17, "Olivier Mattelaer" <email address hidden> wrote:
>Your question #281261 on MadGraph5_aMC@NLO changed:
>https://answers.launchpad.net/mg5amcnlo/+question/281261
>
> Status: Open => Answered
>
>Olivier Mattelaer proposed the following answer:
>Hi,
>
>Could you check if you still have this bug with 2.3.3, Your version
>seems quite old.
>
>Cheers,
>
>Olivier
>> On Jan 21, 2016, at 06:08, toodles <email address hidden> wrote:
>>
>> New question #281261 on MadGraph5_aMC@NLO:
>> https://answers.launchpad.net/mg5amcnlo/+question/281261
>>
>> I am now using MG5 and mssm-full model to generate h2 in 5FS, then h2 decay to ttbar and then jj+lvl.
>>
>> The pattern I am using in madspin is : decay h2 > t t~, (t > w+ b, w+ > j j), (t~ > w- b~, w- > l- vl~)
>>
>> However I got this error. I copied the run_01_tag_1_debug.log content in the following:
>>
>> #************************************************************
>> #* MadGraph5_aMC@NLO/MadEvent *
>> #* *
>> #* * * *
>> #* * * * * *
>> #* * * * * 5 * * * * *
>> #* * * * * *
>> #* * * *
>> #* *
>> #* *
>> #* VERSION 5.2.1.2 *
>> #* *
>> #* 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 *
>> #* *
>> #************************************************************
>> launch
>> Traceback (most recent call last):
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 879, in onecmd
>> return self.onecmd_orig(line, **opt)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
>> return func(arg, **opt)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2139, in do_launch
>> self.do_generate_events(line, *args, **opt)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2128, in do_generate_events
>> self.exec_cmd('decay_events -from_cards', postcmd=False)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 919, in exec_cmd
>> stop = Cmd.onecmd_orig(current_interface, line, **opt)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
>> return func(arg, **opt)
>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/common_run_interface.py", line 1554, in do_decay_events
>> madspin_cmd.import_command_file(path)
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 1035, in import_command_file
>> self.exec_cmd(line, precmd=True)
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
>> stop = Cmd.onecmd_orig(current_interface, line, **opt)
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
>> return func(arg, **opt)
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/interface_madspin.py", line 461, in do_launch
>> generate_all.run()
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2015, in run
>> decay_mapping = self.get_identical_decay()
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2360, in get_identical_decay
>> True, self.pid2width, self.pid2mass, BW_cut,self.Ecollider)
>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 460, in generate_momenta
>> tree = self["tree"][res]
>> KeyError: -4
>> Run Options
>> -----------
>> stdout_level : None
>>
>> MadEvent Options
>> ----------------
>> automatic_html_opening : False (user set)
>> cluster_temp_path : None
>> cluster_queue : madgraph
>> cluster_memory : None
>> cluster_time : None
>> nb_core : 8 (user set)
>> run_mode : 2
>>
>> Configuration Options
>> ---------------------
>> text_editor : None
>> cluster_status_update : (600, 30)
>> pythia8_path : None (user set)
>> hwpp_path : None (user set)
>> pythia-pgs_path : /home/toodles/madness/MG5_aMC_v2_1_2/pythia-pgs (user set)
>> td_path : /home/toodles/madness/MG5_aMC_v2_1_2/td (user set)
>> delphes_path : /home/toodles/madness/MG5_aMC_v2_1_2/Delphes (user set)
>> thepeg_path : None (user set)
>> cluster_type : condor
>> exrootanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/ExRootAnalysis (user set)
>> eps_viewer : None
>> web_browser : None
>> syscalc_path : None (user set)
>> madanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/MadAnalysis (user set)
>> lhapdf : lhapdf-config
>> hepmc_path : None (user set)
>> cluster_retry_wait : 300
>> fortran_compiler : None
>> auto_update : 7 (user set)
>> cluster_nb_retry : 1
>> timeout : 60
>> cpp_compiler : None
>> #************************************************************
>> #* MadGraph5_aMC@NLO *
>> #* *
>> #* * * *
>> #* * * * * *
>> #* * * * * 5 * * * * *
>> #* * * * * *
>> #* * * *
>> #* *
>> #* *
>> #* VERSION 2.1.2 2014-07-03 *
>> #* *
>> #* 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 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 mssm-full
>> define p = p b b~
>> generate p p > h2
>> output test
>> ######################################################################
>> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
>> ######################################################################
>> ###################################
>> ## INFORMATION FOR HMIX
>> ###################################
>> BLOCK HMIX Q= 1.000000e+00 #
>> 1 3.576810e+02 # rmuh
>> 2 9.748624e+00 # tb
>> 4 1.664391e+05 # ma2
>> 3 2.449353e+02 # higgs vev(Q) MSSM DRb
>> ###################################
>> ## INFORMATION FOR MASS
>> ###################################
>> BLOCK MASS #
>> 5 4.889917e+00 # mb
>> 6 1.750000e+02 # mt
>> 15 1.777000e+00 # mta
>> 23 9.118760e+01 # mz
>> 24 7.982901e+01 # mw
>> 25 1.108991e+02 # mh01
>> 35 3.999601e+02 # mh02
>> 36 3.995839e+02 # ma0
>> 37 4.078790e+02 # mh
>> 1000001 5.684411e+02 # msd1
>> 1000002 5.611190e+02 # msu1
>> 1000003 5.684411e+02 # msd2
>> 1000004 5.611190e+02 # msu2
>> 1000005 5.130652e+02 # msd3
>> 1000006 3.996685e+02 # msu3
>> 1000011 2.029157e+02 # msl1
>> 1000012 1.852583e+02 # msn1
>> 1000013 2.029157e+02 # msl2
>> 1000014 1.852583e+02 # msn2
>> 1000015 1.344909e+02 # msl3
>> 1000016 1.847085e+02 # msn3
>> 1000021 6.077137e+02 # mgo
>> 1000022 9.668807e+01 # mneu1
>> 1000023 1.810882e+02 # mneu2
>> 1000024 1.816965e+02 # mch1
>> 1000025 -3.637560e+02 # mneu3
>> 1000035 3.817294e+02 # mneu4
>> 1000037 3.799393e+02 # mch2
>> 2000001 5.452285e+02 # msd4
>> 2000002 5.492593e+02 # msu4
>> 2000003 5.452285e+02 # msd5
>> 2000004 5.492593e+02 # msu5
>> 2000005 5.437267e+02 # msd6
>> 2000006 5.857858e+02 # msu6
>> 2000011 1.441028e+02 # msl4
>> 2000013 1.441028e+02 # msl5
>> 2000015 2.068678e+02 # msl6
>> 1 0.000000e+00 # d : 0.0
>> 2 0.000000e+00 # u : 0.0
>> 3 0.000000e+00 # s : 0.0
>> 4 0.000000e+00 # c : 0.0
>> 11 0.000000e+00 # e- : 0.0
>> 12 0.000000e+00 # ve : 0.0
>> 13 0.000000e+00 # mu- : 0.0
>> 14 0.000000e+00 # vm : 0.0
>> 16 0.000000e+00 # vt : 0.0
>> 21 0.000000e+00 # g : 0.0
>> 22 0.000000e+00 # a : 0.0
>> ###################################
>> ## INFORMATION FOR MSOFT
>> ###################################
>> BLOCK MSOFT Q= 1.000000e+00 #
>> 1 1.013965e+02 # rmx1
>> 2 1.915042e+02 # rmx2
>> 3 5.882630e+02 # rmx3
>> 21 3.233749e+04 # mhd2
>> 22 -1.288001e+05 # mhu2
>> 31 1.953348e+02 # rml21x1
>> 32 1.953348e+02 # rml22x2
>> 33 1.944960e+02 # rml23x3
>> 34 1.364941e+02 # rme21x1
>> 35 1.364941e+02 # rme22x2
>> 36 1.340434e+02 # rme23x3
>> 41 5.475735e+02 # rmq21x1
>> 42 5.475735e+02 # rmq22x2
>> 43 4.987639e+02 # rmq23x3
>> 44 5.295112e+02 # rmu21x1
>> 45 5.295112e+02 # rmu22x2
>> 46 4.232459e+02 # rmu23x3
>> 47 5.231488e+02 # rmd21x1
>> 48 5.231488e+02 # rmd22x2
>> 49 5.198673e+02 # rmd23x3
>> ###################################
>> ## INFORMATION FOR NMIX
>> ###################################
>> BLOCK NMIX #
>> 1 1 9.863644e-01 # rnn1x1
>> 1 2 -5.311036e-02 # rnn1x2
>> 1 3 1.464340e-01 # rnn1x3
>> 1 4 -5.311861e-02 # rnn1x4
>> 2 1 9.935054e-02 # rnn2x1
>> 2 2 9.449493e-01 # rnn2x2
>> 2 3 -2.698467e-01 # rnn2x3
>> 2 4 1.561507e-01 # rnn2x4
>> 3 1 -6.033880e-02 # rnn3x1
>> 3 2 8.770049e-02 # rnn3x2
>> 3 3 6.958775e-01 # rnn3x3
>> 3 4 7.102270e-01 # rnn3x4
>> 4 1 -1.165071e-01 # rnn4x1
>> 4 2 3.107390e-01 # rnn4x2
>> 4 3 6.492260e-01 # rnn4x3
>> 4 4 -6.843778e-01 # rnn4x4
>> ###################################
>> ## INFORMATION FOR SMINPUTS
>> ###################################
>> BLOCK SMINPUTS #
>> 1 1.279340e+02 # aewm1
>> 3 1.180000e-01 # as
>> 6 1.750000e+02 # mt
>> 7 1.777000e+00 # mta
>> 4 9.118760e+01 # mz
>> 2 1.166370e-05 # G_F [GeV^-2]
>> ###################################
>> ## INFORMATION FOR UMIX
>> ###################################
>> BLOCK UMIX #
>> 1 1 9.168349e-01 # ruu1x1
>> 1 2 -3.992666e-01 # ruu1x2
>> 2 1 3.992666e-01 # ruu2x1
>> 2 2 9.168349e-01 # ruu2x2
>> ###################################
>> ## INFORMATION FOR VMIX
>> ###################################
>> BLOCK VMIX #
>> 1 1 9.725578e-01 # rvv1x1
>> 1 2 -2.326612e-01 # rvv1x2
>> 2 1 2.326612e-01 # rvv2x1
>> 2 2 9.725578e-01 # rvv2x2
>> ###################################
>> ## INFORMATION FOR YD
>> ###################################
>> BLOCK YD Q= 1.000000e+00 #
>> 3 3 1.388402e-01 # ryd3x3
>> ###################################
>> ## INFORMATION FOR YE
>> ###################################
>> BLOCK YE Q= 1.000000e+00 #
>> 3 3 1.008908e-01 # rye3x3
>> ###################################
>> ## INFORMATION FOR YU
>> ###################################
>> BLOCK YU Q= 1.000000e+00 #
>> 3 3 8.928445e-01 # ryu3x3
>> ###################################
>> ## INFORMATION FOR DECAY
>> ###################################
>> DECAY 6 1.561950e+00 # wt
>> DECAY 23 2.411433e+00 # wz
>> DECAY 24 2.002822e+00 # ww
>> DECAY 25 1.986108e-03 # wh01
>> DECAY 35 5.748014e-01 # wh02
>> DECAY 36 6.321785e-01 # wa0
>> DECAY 37 5.469628e-01 # wh
>> DECAY 1000001 5.312788e+00 # wsd1
>> DECAY 1000002 5.477195e+00 # wsu1
>> DECAY 1000003 5.312788e+00 # wsd2
>> DECAY 1000004 5.477195e+00 # wsu2
>> DECAY 1000005 3.736276e+00 # wsd3
>> DECAY 1000006 2.021596e+00 # wsu3
>> DECAY 1000011 2.136822e-01 # wsl1
>> DECAY 1000012 1.498816e-01 # wsn1
>> DECAY 1000013 2.136822e-01 # wsl2
>> DECAY 1000014 1.498816e-01 # wsn2
>> DECAY 1000015 1.483273e-01 # wsl3
>> DECAY 1000016 1.475190e-01 # wsn3
>> DECAY 1000021 5.506754e+00 # wgo
>> DECAY 1000023 2.077700e-02 # wneu2
>> DECAY 1000024 1.704145e-02 # wch1
>> DECAY 1000025 1.915985e+00 # wneu3
>> DECAY 1000035 2.585851e+00 # wneu4
>> DECAY 1000037 2.486895e+00 # wch2
>> DECAY 2000001 2.858123e-01 # wsd4
>> DECAY 2000002 1.152973e+00 # wsu4
>> DECAY 2000003 2.858123e-01 # wsd5
>> DECAY 2000004 1.152973e+00 # wsu5
>> DECAY 2000005 8.015663e-01 # wsd6
>> DECAY 2000006 7.373133e+00 # wsu6
>> DECAY 2000011 2.161216e-01 # wsl4
>> DECAY 2000013 2.161216e-01 # wsl5
>> DECAY 2000015 2.699061e-01 # wsl6
>> DECAY 1 0.000000e+00 # d : 0.0
>> DECAY 2 0.000000e+00 # u : 0.0
>> DECAY 3 0.000000e+00 # s : 0.0
>> DECAY 4 0.000000e+00 # c : 0.0
>> DECAY 5 0.000000e+00 # b : 0.0
>> DECAY 11 0.000000e+00 # e- : 0.0
>> DECAY 12 0.000000e+00 # ve : 0.0
>> DECAY 13 0.000000e+00 # mu- : 0.0
>> DECAY 14 0.000000e+00 # vm : 0.0
>> DECAY 15 0.000000e+00 # ta- : 0.0
>> DECAY 16 0.000000e+00 # vt : 0.0
>> DECAY 21 0.000000e+00 # g : 0.0
>> DECAY 22 0.000000e+00 # a : 0.0
>> DECAY 1000022 0.000000e+00 # n1 : 0.0
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000022
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000023
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000025
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000035
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000024
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000037
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000021
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 35
>> ###################################
>> BLOCK QNUMBERS 35 # h2
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 36
>> ###################################
>> BLOCK QNUMBERS 36 # h3
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 37
>> ###################################
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000012
>> ###################################
>> BLOCK QNUMBERS 1000012 # sve
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000014
>> ###################################
>> BLOCK QNUMBERS 1000014 # svm
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000016
>> ###################################
>> BLOCK QNUMBERS 1000016 # svt
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000011
>> ###################################
>> BLOCK QNUMBERS 1000011 # el-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000013
>> ###################################
>> BLOCK QNUMBERS 1000013 # mul-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000015
>> ###################################
>> BLOCK QNUMBERS 1000015 # ta1-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000011
>> ###################################
>> BLOCK QNUMBERS 2000011 # er-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000013
>> ###################################
>> BLOCK QNUMBERS 2000013 # mur-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000015
>> ###################################
>> BLOCK QNUMBERS 2000015 # ta2-
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000002
>> ###################################
>> BLOCK QNUMBERS 1000002 # ul
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000004
>> ###################################
>> BLOCK QNUMBERS 1000004 # cl
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000006
>> ###################################
>> BLOCK QNUMBERS 1000006 # t1
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000002
>> ###################################
>> BLOCK QNUMBERS 2000002 # ur
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000004
>> ###################################
>> BLOCK QNUMBERS 2000004 # cr
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000006
>> ###################################
>> BLOCK QNUMBERS 2000006 # t2
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000001
>> ###################################
>> BLOCK QNUMBERS 1000001 # dl
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000003
>> ###################################
>> BLOCK QNUMBERS 1000003 # sl
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 1000005
>> ###################################
>> BLOCK QNUMBERS 1000005 # b1
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000001
>> ###################################
>> BLOCK QNUMBERS 2000001 # dr
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000003
>> ###################################
>> BLOCK QNUMBERS 2000003 # sr
>> 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)
>> ###################################
>> ## INFORMATION FOR QNUMBERS 2000005
>> ###################################
>> BLOCK QNUMBERS 2000005 # b2
>> 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)
>> ###################################
>> ## INFORMATION FOR MODSEL
>> ###################################
>> BLOCK MODSEL #
>> 1 1 #
>> ###################################
>> ## INFORMATION FOR STOPMIX
>> ###################################
>> BLOCK STOPMIX #
>> 1 1 5.536450e-01 # rru3x3
>> 1 2 8.327528e-01 # rru3x6
>> 2 1 8.327528e-01 # rru6x3
>> 2 2 -5.536450e-01 # rru6x6
>> ###################################
>> ## INFORMATION FOR SBOTMIX
>> ###################################
>> BLOCK SBOTMIX #
>> 1 1 9.387379e-01 # rrd3x3
>> 1 2 3.446319e-01 # rrd3x6
>> 2 1 -3.446319e-01 # rrd6x3
>> 2 2 9.387379e-01 # rrd6x6
>> ###################################
>> ## INFORMATION FOR STAUMIX
>> ###################################
>> BLOCK STAUMIX #
>> 1 1 2.824872e-01 # rrl3x3
>> 1 2 9.592711e-01 # rrl3x6
>> 2 1 9.592711e-01 # rrl6x3
>> 2 2 -2.824872e-01 # rrl6x6
>> ###################################
>> ## INFORMATION FOR ALPHA
>> ###################################
>> BLOCK ALPHA #
>> -1.138252e-01 # alp
>> ###################################
>> ## INFORMATION FOR AE
>> ###################################
>> BLOCK AE Q= 1.000000e+00 #
>> 3 3 -2.517769e+02 # A_tau(Q) DRbar
>> 1 1 0.000000e+00 # A_e(Q) DRbar
>> 2 2 0.000000e+00 # A_mu(Q) DRbar
>> ###################################
>> ## INFORMATION FOR AU
>> ###################################
>> BLOCK AU Q= 1.000000e+00 #
>> 3 3 -4.981299e+02 # A_t(Q) DRbar
>> 1 1 0.000000e+00 # A_u(Q) DRbar
>> 2 2 0.000000e+00 # A_c(Q) DRbar
>> ###################################
>> ## INFORMATION FOR AD
>> ###################################
>> BLOCK AD Q= 1.000000e+00 #
>> 3 3 -7.972741e+02 # A_b(Q) DRbar
>> 1 1 0.000000e+00 # A_d(Q) DRbar
>> 2 2 0.000000e+00 # A_s(Q) DRbar
>> #*********************************************************************
>> # 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. *
>> #*********************************************************************
>> 1000 = 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
>> #*********************************************************************
>> # Renormalization and factorization scales *
>> #*********************************************************************
>> F = fixed_ren_scale ! if .true. use fixed ren scale
>> F = fixed_fac_scale ! if .true. use fixed fac scale
>> 91.1880 = scale ! fixed ren scale
>> 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
>> 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
>> 1 = scalefact ! scale factor for event-by-event scales
>> #*********************************************************************
>> # Matching - Warning! ickkw > 1 is still beta
>> #*********************************************************************
>> 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
>> 1 = highestmult ! for ickkw=2, highest mult group
>> 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
>> 1 = alpsfact ! scale factor for QCD emission vx
>> F = chcluster ! cluster only according to channel diag
>> T = pdfwgt ! for ickkw=1, perform pdf reweighting
>> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
>> T = clusinfo ! include clustering tag in output
>> 3.0 = lhe_version ! Change the way clustering information pass to shower.
>> #*********************************************************************
>> #**********************************************************
>> #
>> #**********************************************************
>> # Automatic ptj and mjj cuts if xqcut > 0
>> # (turn off for VBF and single top processes)
>> #**********************************************************
>> T = auto_ptj_mjj ! Automatic setting of ptj and mjj
>> #**********************************************************
>> #
>> #**********************************
>> # BW cutoff (M+/-bwcutoff*Gamma)
>> #**********************************
>> 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
>> #**********************************************************
>> # Apply pt/E/eta/dr/mij cuts on decay products or not
>> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
>> #**********************************************************
>> T = cut_decays ! Cut decay products
>> #*************************************************************
>> # Number of helicities to sum per event (0 = all helicities)
>> # 0 gives more stable result, but longer run time (needed for
>> # long decay chains e.g.).
>> # Use >=2 if most helicities contribute, e.g. pure QCD.
>> #*************************************************************
>> 0 = nhel ! Number of helicities used per event
>> #*******************
>> # Standard Cuts
>> #*******************
>> #
>> #*********************************************************************
>> # Minimum and maximum pt's (for max, -1 means no cut) *
>> #*********************************************************************
>> 20 = ptj ! minimum pt for the jets
>> 0 = ptb ! minimum pt for the b
>> 10 = pta ! minimum pt for the photons
>> 10 = ptl ! minimum pt for the charged leptons
>> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
>> 0 = ptheavy ! minimum pt for one heavy final state
>> 1.0 = ptonium ! minimum pt for the quarkonium states
>> -1 = ptjmax ! maximum pt for the jets
>> -1 = ptbmax ! maximum pt for the b
>> -1 = ptamax ! maximum pt for the photons
>> -1 = ptlmax ! maximum pt for the charged leptons
>> -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
>> #*********************************************************************
>> # Minimum and maximum E's (in the center of mass frame) *
>> #*********************************************************************
>> 0 = ej ! minimum E for the jets
>> 0 = eb ! minimum E for the b
>> 0 = ea ! minimum E for the photons
>> 0 = el ! minimum E for the charged leptons
>> -1 = ejmax ! maximum E for the jets
>> -1 = ebmax ! maximum E for the b
>> -1 = eamax ! maximum E for the photons
>> -1 = elmax ! maximum E for the charged leptons
>> #*********************************************************************
>> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
>> #*********************************************************************
>> 5 = etaj ! max rap for the jets
>> -1 = etab ! max rap for the b
>> 2.5 = etaa ! max rap for the photons
>> 2.5 = etal ! max rap for the charged leptons
>> 0.6 = etaonium ! max rap for the quarkonium states
>> 0 = etajmin ! min rap for the jets
>> 0 = etabmin ! min rap for the b
>> 0 = etaamin ! min rap for the photons
>> 0 = etalmin ! main rap for the charged leptons
>> #*********************************************************************
>> # Minimum and maximum DeltaR distance *
>> #*********************************************************************
>> 0.4 = drjj ! min distance between jets
>> 0 = drbb ! min distance between b's
>> 0.4 = drll ! min distance between leptons
>> 0.4 = draa ! min distance between gammas
>> 0 = drbj ! min distance between b and jet
>> 0.4 = draj ! min distance between gamma and jet
>> 0.4 = drjl ! min distance between jet and lepton
>> 0 = drab ! min distance between gamma and b
>> 0 = drbl ! min distance between b and lepton
>> 0.4 = dral ! min distance between gamma and lepton
>> -1 = drjjmax ! max distance between jets
>> -1 = drbbmax ! max distance between b's
>> -1 = drllmax ! max distance between leptons
>> -1 = draamax ! max distance between gammas
>> -1 = drbjmax ! max distance between b and jet
>> -1 = drajmax ! max distance between gamma and jet
>> -1 = drjlmax ! max distance between jet and lepton
>> -1 = drabmax ! max distance between gamma and b
>> -1 = drblmax ! max distance between b and lepton
>> -1 = dralmax ! maxdistance between gamma and lepton
>> #*********************************************************************
>> # Minimum and maximum invariant mass for pairs *
>> # WARNING: for four lepton final state mmll cut require to have *
>> # different lepton masses for each flavor! *
>> #*********************************************************************
>> 0 = mmjj ! min invariant mass of a jet pair
>> 0 = mmbb ! min invariant mass of a b pair
>> 0 = mmaa ! min invariant mass of gamma gamma pair
>> 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
>> -1 = mmjjmax ! max invariant mass of a jet pair
>> -1 = mmbbmax ! max invariant mass of a b pair
>> -1 = mmaamax ! max invariant mass of gamma gamma pair
>> -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
>> #*********************************************************************
>> # Minimum and maximum invariant mass for all letpons *
>> #*********************************************************************
>> 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
>> -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
>> #*********************************************************************
>> # Minimum and maximum pt for 4-momenta sum of leptons *
>> #*********************************************************************
>> 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
>> -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
>> #*********************************************************************
>> # Inclusive cuts *
>> #*********************************************************************
>> 0 = xptj ! minimum pt for at least one jet
>> 0 = xptb ! minimum pt for at least one b
>> 0 = xpta ! minimum pt for at least one photon
>> 0 = xptl ! minimum pt for at least one charged lepton
>> #*********************************************************************
>> # Control the pt's of the jets sorted by pt *
>> #*********************************************************************
>> 0 = ptj1min ! minimum pt for the leading jet in pt
>> 0 = ptj2min ! minimum pt for the second jet in pt
>> 0 = ptj3min ! minimum pt for the third jet in pt
>> 0 = ptj4min ! minimum pt for the fourth jet in pt
>> -1 = ptj1max ! maximum pt for the leading jet in pt
>> -1 = ptj2max ! maximum pt for the second jet in pt
>> -1 = ptj3max ! maximum pt for the third jet in pt
>> -1 = ptj4max ! maximum pt for the fourth jet in pt
>> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
>> #*********************************************************************
>> # Control the pt's of leptons sorted by pt *
>> #*********************************************************************
>> 0 = ptl1min ! minimum pt for the leading lepton in pt
>> 0 = ptl2min ! minimum pt for the second lepton in pt
>> 0 = ptl3min ! minimum pt for the third lepton in pt
>> 0 = ptl4min ! minimum pt for the fourth lepton in pt
>> -1 = ptl1max ! maximum pt for the leading lepton in pt
>> -1 = ptl2max ! maximum pt for the second lepton in pt
>> -1 = ptl3max ! maximum pt for the third lepton in pt
>> -1 = ptl4max ! maximum pt for the fourth lepton in pt
>> #*********************************************************************
>> # Control the Ht(k)=Sum of k leading jets *
>> #*********************************************************************
>> 0 = htjmin ! minimum jet HT=Sum(jet pt)
>> -1 = htjmax ! maximum jet HT=Sum(jet pt)
>> 0 = ihtmin !inclusive Ht for all partons (including b)
>> -1 = ihtmax !inclusive Ht for all partons (including b)
>> 0 = ht2min ! minimum Ht for the two leading jets
>> 0 = ht3min ! minimum Ht for the three leading jets
>> 0 = ht4min ! minimum Ht for the four leading jets
>> -1 = ht2max ! maximum Ht for the two leading jets
>> -1 = ht3max ! maximum Ht for the three leading jets
>> -1 = ht4max ! maximum Ht for the four leading jets
>> #***********************************************************************
>> # Photon-isolation cuts, according to hep-ph/9801442 *
>> # When ptgmin=0, all the other parameters are ignored *
>> # When ptgmin>0, pta and draj are not going to be used *
>> #***********************************************************************
>> 0 = ptgmin ! Min photon transverse momentum
>> 0.4 = R0gamma ! Radius of isolation code
>> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
>> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
>> .true. = isoEM ! isolate photons from EM energy (photons and leptons)
>> #*********************************************************************
>> # WBF cuts *
>> #*********************************************************************
>> 0 = xetamin ! minimum rapidity for two jets in the WBF case
>> 0 = deltaeta ! minimum rapidity for two jets in the WBF case
>> #*********************************************************************
>> # KT DURHAM CUT *
>> #*********************************************************************
>> -1 = ktdurham
>> 0.4 = dparameter
>> #*********************************************************************
>> # maximal pdg code for quark to be considered as a light jet *
>> # (otherwise b cuts are applied) *
>> #*********************************************************************
>> 4 = maxjetflavor ! Maximum jet pdg code
>> #*********************************************************************
>> # Jet measure cuts *
>> #*********************************************************************
>> 0 = xqcut ! minimum kt jet measure between partons
>> #*********************************************************************
>> #
>> #*********************************************************************
>> # Store info for systematics studies *
>> # WARNING: If use_syst is T, matched Pythia output is *
>> # meaningful ONLY if plotted taking matchscale *
>> # reweighting into account! *
>> #*********************************************************************
>> F = use_syst ! Enable systematics studies
>> #
>> #**************************************
>> # Parameter of the systematics study
>> # will be used by SysCalc (if installed)
>> #**************************************
>> #
>> 0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
>> 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
>> 30 50 = sys_matchscale # variation of merging scale
>> # PDF sets and number of members (0 or none for all members).
>> CT10nlo.LHgrid = sys_pdf # matching scales
>> # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
>>
>> BTW, if I change the decay pattern to be
>> decay h2 > t t~, t > w+ b, w+ > j j, t~ > w- b~, w- > l- vl~
>> there is no such error anymore. But the decay get stuck at "Decaying the events.....".
>>
>> I am really appreciate for your help
>>
>> Regards,
>> Ying-Ying
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
>--
>If this answers your question, please go to the following page to let us
>know that it is solved:
>https://answers.launchpad.net/mg5amcnlo/+question/281261/+confirm?answer_id=0
>
>If you still need help, you can reply to this email or go to the
>following page to enter your feedback:
>https://answers.launchpad.net/mg5amcnlo/+question/281261
>
>You received this question notification because you asked the question.

Indeed, this bug is fixed in the development version.

So the easiest is that you use the development version:
bzr branch lp:~maddevelopers/mg5amcnlo/2.3.4

Alternatively, you can generate your events with your version, and decay the events with the madspin version of 2.3.3.

Cheers,

Olivier

> On Jan 21, 2016, at 13:42, toodles <email address hidden> wrote:
>
> Question #281261 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/281261
>
> Status: Answered => Open
>
> toodles is still having a problem:
> Hi, Olivier
>
>
> Actually when I use 2.3.3 and use mssm-full, I cannot do output.
>
>
> Here is the error massage:
>
> Command "output test" interrupted with error:
>
> IOError : [Errno 2] No such file or directory:
> '/Users/toodles/madness/MG5_aMC_v2_3_3/test/Source/MODEL/param_card_rule.dat'
>
>
> here is the MG5_debug file:
>
>
> import model mssm-full
> define p = p b b~
> generate p p > h2
> output test
> Traceback (most recent call last):
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 908, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/extended_cmd.py", line 897, in onecmd_orig
> return func(arg, **opt)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/master_interface.py", line 284, in do_output
> self.cmd.do_output(self, line, *args, **opts)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/madgraph_interface.py", line 6687, in do_output
> self.finalize(nojpeg, flaglist=flaglist)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/interface/madgraph_interface.py", line 6984, in finalize
> wanted_couplings)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 3188, in convert_model_to_mg4
> wanted_lorentz, wanted_couplings)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 753, in convert_model_to_mg4
> model_builder.build(wanted_couplings)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 5149, in build
> self.write_all()
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 5198, in write_all
> self.create_param_card()
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 6248, in create_param_card
> mssm_convert=True)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/madgraph/iolibs/export_v4.py", line 6239, in create_param_card_static
> translator.make_valid_param_card(output_path, rule_card_path)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/models/check_param_card.py", line 1461, in make_valid_param_card
> cardrule.load_rule(restrictpath)
> File "/Users/toodles/madness/MG5_aMC_v2_3_3/models/check_param_card.py", line 882, in load_rule
> tree = ET.parse(inputpath)
> File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/xml/etree/ElementTree.py", line 1182, in parse
> tree.parse(source, parser)
> File "/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/xml/etree/ElementTree.py", line 647, in parse
> source = open(source, "rb")
> IOError: [Errno 2] No such file or directory: '/Users/toodles/madness/MG5_aMC_v2_3_3/test/Source/MODEL/param_card_rule.dat'
>
>
> So I cannot try mad spin in mssm-full model with 2.3.3
>
>
> Thanks for your reply.
>
>
> Regards,
>
>
> Ying-Ying
>
>
>
>
>
> At 2016-01-21 16:13:17, "Olivier Mattelaer" <email address hidden> wrote:
>> Your question #281261 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/281261
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi,
>>
>> Could you check if you still have this bug with 2.3.3, Your version
>> seems quite old.
>>
>> Cheers,
>>
>> Olivier
>>> On Jan 21, 2016, at 06:08, toodles <email address hidden> wrote:
>>>
>>> New question #281261 on MadGraph5_aMC@NLO:
>>> https://answers.launchpad.net/mg5amcnlo/+question/281261
>>>
>>> I am now using MG5 and mssm-full model to generate h2 in 5FS, then h2 decay to ttbar and then jj+lvl.
>>>
>>> The pattern I am using in madspin is : decay h2 > t t~, (t > w+ b, w+ > j j), (t~ > w- b~, w- > l- vl~)
>>>
>>> However I got this error. I copied the run_01_tag_1_debug.log content in the following:
>>>
>>> #************************************************************
>>> #* MadGraph5_aMC@NLO/MadEvent *
>>> #* *
>>> #* * * *
>>> #* * * * * *
>>> #* * * * * 5 * * * * *
>>> #* * * * * *
>>> #* * * *
>>> #* *
>>> #* *
>>> #* VERSION 5.2.1.2 *
>>> #* *
>>> #* 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 *
>>> #* *
>>> #************************************************************
>>> launch
>>> Traceback (most recent call last):
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 879, in onecmd
>>> return self.onecmd_orig(line, **opt)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
>>> return func(arg, **opt)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2139, in do_launch
>>> self.do_generate_events(line, *args, **opt)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/madevent_interface.py", line 2128, in do_generate_events
>>> self.exec_cmd('decay_events -from_cards', postcmd=False)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 919, in exec_cmd
>>> stop = Cmd.onecmd_orig(current_interface, line, **opt)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/extended_cmd.py", line 872, in onecmd_orig
>>> return func(arg, **opt)
>>> File "/home/toodles/madness/Project_Heavyhiggs/bbHA/test/bin/internal/common_run_interface.py", line 1554, in do_decay_events
>>> madspin_cmd.import_command_file(path)
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 1035, in import_command_file
>>> self.exec_cmd(line, precmd=True)
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
>>> stop = Cmd.onecmd_orig(current_interface, line, **opt)
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
>>> return func(arg, **opt)
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/interface_madspin.py", line 461, in do_launch
>>> generate_all.run()
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2015, in run
>>> decay_mapping = self.get_identical_decay()
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 2360, in get_identical_decay
>>> True, self.pid2width, self.pid2mass, BW_cut,self.Ecollider)
>>> File "/home/toodles/madness/MG5_aMC_v2_1_2/MadSpin/decay.py", line 460, in generate_momenta
>>> tree = self["tree"][res]
>>> KeyError: -4
>>> Run Options
>>> -----------
>>> stdout_level : None
>>>
>>> MadEvent Options
>>> ----------------
>>> automatic_html_opening : False (user set)
>>> cluster_temp_path : None
>>> cluster_queue : madgraph
>>> cluster_memory : None
>>> cluster_time : None
>>> nb_core : 8 (user set)
>>> run_mode : 2
>>>
>>> Configuration Options
>>> ---------------------
>>> text_editor : None
>>> cluster_status_update : (600, 30)
>>> pythia8_path : None (user set)
>>> hwpp_path : None (user set)
>>> pythia-pgs_path : /home/toodles/madness/MG5_aMC_v2_1_2/pythia-pgs (user set)
>>> td_path : /home/toodles/madness/MG5_aMC_v2_1_2/td (user set)
>>> delphes_path : /home/toodles/madness/MG5_aMC_v2_1_2/Delphes (user set)
>>> thepeg_path : None (user set)
>>> cluster_type : condor
>>> exrootanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/ExRootAnalysis (user set)
>>> eps_viewer : None
>>> web_browser : None
>>> syscalc_path : None (user set)
>>> madanalysis_path : /home/toodles/madness/MG5_aMC_v2_1_2/MadAnalysis (user set)
>>> lhapdf : lhapdf-config
>>> hepmc_path : None (user set)
>>> cluster_retry_wait : 300
>>> fortran_compiler : None
>>> auto_update : 7 (user set)
>>> cluster_nb_retry : 1
>>> timeout : 60
>>> cpp_compiler : None
>>> #************************************************************
>>> #* MadGraph5_aMC@NLO *
>>> #* *
>>> #* * * *
>>> #* * * * * *
>>> #* * * * * 5 * * * * *
>>> #* * * * * *
>>> #* * * *
>>> #* *
>>> #* *
>>> #* VERSION 2.1.2 2014-07-03 *
>>> #* *
>>> #* 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 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 mssm-full
>>> define p = p b b~
>>> generate p p > h2
>>> output test
>>> ######################################################################
>>> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
>>> ######################################################################
>>> ###################################
>>> ## INFORMATION FOR HMIX
>>> ###################################
>>> BLOCK HMIX Q= 1.000000e+00 #
>>> 1 3.576810e+02 # rmuh
>>> 2 9.748624e+00 # tb
>>> 4 1.664391e+05 # ma2
>>> 3 2.449353e+02 # higgs vev(Q) MSSM DRb
>>> ###################################
>>> ## INFORMATION FOR MASS
>>> ###################################
>>> BLOCK MASS #
>>> 5 4.889917e+00 # mb
>>> 6 1.750000e+02 # mt
>>> 15 1.777000e+00 # mta
>>> 23 9.118760e+01 # mz
>>> 24 7.982901e+01 # mw
>>> 25 1.108991e+02 # mh01
>>> 35 3.999601e+02 # mh02
>>> 36 3.995839e+02 # ma0
>>> 37 4.078790e+02 # mh
>>> 1000001 5.684411e+02 # msd1
>>> 1000002 5.611190e+02 # msu1
>>> 1000003 5.684411e+02 # msd2
>>> 1000004 5.611190e+02 # msu2
>>> 1000005 5.130652e+02 # msd3
>>> 1000006 3.996685e+02 # msu3
>>> 1000011 2.029157e+02 # msl1
>>> 1000012 1.852583e+02 # msn1
>>> 1000013 2.029157e+02 # msl2
>>> 1000014 1.852583e+02 # msn2
>>> 1000015 1.344909e+02 # msl3
>>> 1000016 1.847085e+02 # msn3
>>> 1000021 6.077137e+02 # mgo
>>> 1000022 9.668807e+01 # mneu1
>>> 1000023 1.810882e+02 # mneu2
>>> 1000024 1.816965e+02 # mch1
>>> 1000025 -3.637560e+02 # mneu3
>>> 1000035 3.817294e+02 # mneu4
>>> 1000037 3.799393e+02 # mch2
>>> 2000001 5.452285e+02 # msd4
>>> 2000002 5.492593e+02 # msu4
>>> 2000003 5.452285e+02 # msd5
>>> 2000004 5.492593e+02 # msu5
>>> 2000005 5.437267e+02 # msd6
>>> 2000006 5.857858e+02 # msu6
>>> 2000011 1.441028e+02 # msl4
>>> 2000013 1.441028e+02 # msl5
>>> 2000015 2.068678e+02 # msl6
>>> 1 0.000000e+00 # d : 0.0
>>> 2 0.000000e+00 # u : 0.0
>>> 3 0.000000e+00 # s : 0.0
>>> 4 0.000000e+00 # c : 0.0
>>> 11 0.000000e+00 # e- : 0.0
>>> 12 0.000000e+00 # ve : 0.0
>>> 13 0.000000e+00 # mu- : 0.0
>>> 14 0.000000e+00 # vm : 0.0
>>> 16 0.000000e+00 # vt : 0.0
>>> 21 0.000000e+00 # g : 0.0
>>> 22 0.000000e+00 # a : 0.0
>>> ###################################
>>> ## INFORMATION FOR MSOFT
>>> ###################################
>>> BLOCK MSOFT Q= 1.000000e+00 #
>>> 1 1.013965e+02 # rmx1
>>> 2 1.915042e+02 # rmx2
>>> 3 5.882630e+02 # rmx3
>>> 21 3.233749e+04 # mhd2
>>> 22 -1.288001e+05 # mhu2
>>> 31 1.953348e+02 # rml21x1
>>> 32 1.953348e+02 # rml22x2
>>> 33 1.944960e+02 # rml23x3
>>> 34 1.364941e+02 # rme21x1
>>> 35 1.364941e+02 # rme22x2
>>> 36 1.340434e+02 # rme23x3
>>> 41 5.475735e+02 # rmq21x1
>>> 42 5.475735e+02 # rmq22x2
>>> 43 4.987639e+02 # rmq23x3
>>> 44 5.295112e+02 # rmu21x1
>>> 45 5.295112e+02 # rmu22x2
>>> 46 4.232459e+02 # rmu23x3
>>> 47 5.231488e+02 # rmd21x1
>>> 48 5.231488e+02 # rmd22x2
>>> 49 5.198673e+02 # rmd23x3
>>> ###################################
>>> ## INFORMATION FOR NMIX
>>> ###################################
>>> BLOCK NMIX #
>>> 1 1 9.863644e-01 # rnn1x1
>>> 1 2 -5.311036e-02 # rnn1x2
>>> 1 3 1.464340e-01 # rnn1x3
>>> 1 4 -5.311861e-02 # rnn1x4
>>> 2 1 9.935054e-02 # rnn2x1
>>> 2 2 9.449493e-01 # rnn2x2
>>> 2 3 -2.698467e-01 # rnn2x3
>>> 2 4 1.561507e-01 # rnn2x4
>>> 3 1 -6.033880e-02 # rnn3x1
>>> 3 2 8.770049e-02 # rnn3x2
>>> 3 3 6.958775e-01 # rnn3x3
>>> 3 4 7.102270e-01 # rnn3x4
>>> 4 1 -1.165071e-01 # rnn4x1
>>> 4 2 3.107390e-01 # rnn4x2
>>> 4 3 6.492260e-01 # rnn4x3
>>> 4 4 -6.843778e-01 # rnn4x4
>>> ###################################
>>> ## INFORMATION FOR SMINPUTS
>>> ###################################
>>> BLOCK SMINPUTS #
>>> 1 1.279340e+02 # aewm1
>>> 3 1.180000e-01 # as
>>> 6 1.750000e+02 # mt
>>> 7 1.777000e+00 # mta
>>> 4 9.118760e+01 # mz
>>> 2 1.166370e-05 # G_F [GeV^-2]
>>> ###################################
>>> ## INFORMATION FOR UMIX
>>> ###################################
>>> BLOCK UMIX #
>>> 1 1 9.168349e-01 # ruu1x1
>>> 1 2 -3.992666e-01 # ruu1x2
>>> 2 1 3.992666e-01 # ruu2x1
>>> 2 2 9.168349e-01 # ruu2x2
>>> ###################################
>>> ## INFORMATION FOR VMIX
>>> ###################################
>>> BLOCK VMIX #
>>> 1 1 9.725578e-01 # rvv1x1
>>> 1 2 -2.326612e-01 # rvv1x2
>>> 2 1 2.326612e-01 # rvv2x1
>>> 2 2 9.725578e-01 # rvv2x2
>>> ###################################
>>> ## INFORMATION FOR YD
>>> ###################################
>>> BLOCK YD Q= 1.000000e+00 #
>>> 3 3 1.388402e-01 # ryd3x3
>>> ###################################
>>> ## INFORMATION FOR YE
>>> ###################################
>>> BLOCK YE Q= 1.000000e+00 #
>>> 3 3 1.008908e-01 # rye3x3
>>> ###################################
>>> ## INFORMATION FOR YU
>>> ###################################
>>> BLOCK YU Q= 1.000000e+00 #
>>> 3 3 8.928445e-01 # ryu3x3
>>> ###################################
>>> ## INFORMATION FOR DECAY
>>> ###################################
>>> DECAY 6 1.561950e+00 # wt
>>> DECAY 23 2.411433e+00 # wz
>>> DECAY 24 2.002822e+00 # ww
>>> DECAY 25 1.986108e-03 # wh01
>>> DECAY 35 5.748014e-01 # wh02
>>> DECAY 36 6.321785e-01 # wa0
>>> DECAY 37 5.469628e-01 # wh
>>> DECAY 1000001 5.312788e+00 # wsd1
>>> DECAY 1000002 5.477195e+00 # wsu1
>>> DECAY 1000003 5.312788e+00 # wsd2
>>> DECAY 1000004 5.477195e+00 # wsu2
>>> DECAY 1000005 3.736276e+00 # wsd3
>>> DECAY 1000006 2.021596e+00 # wsu3
>>> DECAY 1000011 2.136822e-01 # wsl1
>>> DECAY 1000012 1.498816e-01 # wsn1
>>> DECAY 1000013 2.136822e-01 # wsl2
>>> DECAY 1000014 1.498816e-01 # wsn2
>>> DECAY 1000015 1.483273e-01 # wsl3
>>> DECAY 1000016 1.475190e-01 # wsn3
>>> DECAY 1000021 5.506754e+00 # wgo
>>> DECAY 1000023 2.077700e-02 # wneu2
>>> DECAY 1000024 1.704145e-02 # wch1
>>> DECAY 1000025 1.915985e+00 # wneu3
>>> DECAY 1000035 2.585851e+00 # wneu4
>>> DECAY 1000037 2.486895e+00 # wch2
>>> DECAY 2000001 2.858123e-01 # wsd4
>>> DECAY 2000002 1.152973e+00 # wsu4
>>> DECAY 2000003 2.858123e-01 # wsd5
>>> DECAY 2000004 1.152973e+00 # wsu5
>>> DECAY 2000005 8.015663e-01 # wsd6
>>> DECAY 2000006 7.373133e+00 # wsu6
>>> DECAY 2000011 2.161216e-01 # wsl4
>>> DECAY 2000013 2.161216e-01 # wsl5
>>> DECAY 2000015 2.699061e-01 # wsl6
>>> DECAY 1 0.000000e+00 # d : 0.0
>>> DECAY 2 0.000000e+00 # u : 0.0
>>> DECAY 3 0.000000e+00 # s : 0.0
>>> DECAY 4 0.000000e+00 # c : 0.0
>>> DECAY 5 0.000000e+00 # b : 0.0
>>> DECAY 11 0.000000e+00 # e- : 0.0
>>> DECAY 12 0.000000e+00 # ve : 0.0
>>> DECAY 13 0.000000e+00 # mu- : 0.0
>>> DECAY 14 0.000000e+00 # vm : 0.0
>>> DECAY 15 0.000000e+00 # ta- : 0.0
>>> DECAY 16 0.000000e+00 # vt : 0.0
>>> DECAY 21 0.000000e+00 # g : 0.0
>>> DECAY 22 0.000000e+00 # a : 0.0
>>> DECAY 1000022 0.000000e+00 # n1 : 0.0
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000022
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000023
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000025
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000035
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000024
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000037
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000021
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 35
>>> ###################################
>>> BLOCK QNUMBERS 35 # h2
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 36
>>> ###################################
>>> BLOCK QNUMBERS 36 # h3
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 37
>>> ###################################
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000012
>>> ###################################
>>> BLOCK QNUMBERS 1000012 # sve
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000014
>>> ###################################
>>> BLOCK QNUMBERS 1000014 # svm
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000016
>>> ###################################
>>> BLOCK QNUMBERS 1000016 # svt
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000011
>>> ###################################
>>> BLOCK QNUMBERS 1000011 # el-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000013
>>> ###################################
>>> BLOCK QNUMBERS 1000013 # mul-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000015
>>> ###################################
>>> BLOCK QNUMBERS 1000015 # ta1-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000011
>>> ###################################
>>> BLOCK QNUMBERS 2000011 # er-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000013
>>> ###################################
>>> BLOCK QNUMBERS 2000013 # mur-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000015
>>> ###################################
>>> BLOCK QNUMBERS 2000015 # ta2-
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000002
>>> ###################################
>>> BLOCK QNUMBERS 1000002 # ul
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000004
>>> ###################################
>>> BLOCK QNUMBERS 1000004 # cl
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000006
>>> ###################################
>>> BLOCK QNUMBERS 1000006 # t1
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000002
>>> ###################################
>>> BLOCK QNUMBERS 2000002 # ur
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000004
>>> ###################################
>>> BLOCK QNUMBERS 2000004 # cr
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000006
>>> ###################################
>>> BLOCK QNUMBERS 2000006 # t2
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000001
>>> ###################################
>>> BLOCK QNUMBERS 1000001 # dl
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000003
>>> ###################################
>>> BLOCK QNUMBERS 1000003 # sl
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 1000005
>>> ###################################
>>> BLOCK QNUMBERS 1000005 # b1
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000001
>>> ###################################
>>> BLOCK QNUMBERS 2000001 # dr
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000003
>>> ###################################
>>> BLOCK QNUMBERS 2000003 # sr
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR QNUMBERS 2000005
>>> ###################################
>>> BLOCK QNUMBERS 2000005 # b2
>>> 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)
>>> ###################################
>>> ## INFORMATION FOR MODSEL
>>> ###################################
>>> BLOCK MODSEL #
>>> 1 1 #
>>> ###################################
>>> ## INFORMATION FOR STOPMIX
>>> ###################################
>>> BLOCK STOPMIX #
>>> 1 1 5.536450e-01 # rru3x3
>>> 1 2 8.327528e-01 # rru3x6
>>> 2 1 8.327528e-01 # rru6x3
>>> 2 2 -5.536450e-01 # rru6x6
>>> ###################################
>>> ## INFORMATION FOR SBOTMIX
>>> ###################################
>>> BLOCK SBOTMIX #
>>> 1 1 9.387379e-01 # rrd3x3
>>> 1 2 3.446319e-01 # rrd3x6
>>> 2 1 -3.446319e-01 # rrd6x3
>>> 2 2 9.387379e-01 # rrd6x6
>>> ###################################
>>> ## INFORMATION FOR STAUMIX
>>> ###################################
>>> BLOCK STAUMIX #
>>> 1 1 2.824872e-01 # rrl3x3
>>> 1 2 9.592711e-01 # rrl3x6
>>> 2 1 9.592711e-01 # rrl6x3
>>> 2 2 -2.824872e-01 # rrl6x6
>>> ###################################
>>> ## INFORMATION FOR ALPHA
>>> ###################################
>>> BLOCK ALPHA #
>>> -1.138252e-01 # alp
>>> ###################################
>>> ## INFORMATION FOR AE
>>> ###################################
>>> BLOCK AE Q= 1.000000e+00 #
>>> 3 3 -2.517769e+02 # A_tau(Q) DRbar
>>> 1 1 0.000000e+00 # A_e(Q) DRbar
>>> 2 2 0.000000e+00 # A_mu(Q) DRbar
>>> ###################################
>>> ## INFORMATION FOR AU
>>> ###################################
>>> BLOCK AU Q= 1.000000e+00 #
>>> 3 3 -4.981299e+02 # A_t(Q) DRbar
>>> 1 1 0.000000e+00 # A_u(Q) DRbar
>>> 2 2 0.000000e+00 # A_c(Q) DRbar
>>> ###################################
>>> ## INFORMATION FOR AD
>>> ###################################
>>> BLOCK AD Q= 1.000000e+00 #
>>> 3 3 -7.972741e+02 # A_b(Q) DRbar
>>> 1 1 0.000000e+00 # A_d(Q) DRbar
>>> 2 2 0.000000e+00 # A_s(Q) DRbar
>>> #*********************************************************************
>>> # 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. *
>>> #*********************************************************************
>>> 1000 = 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
>>> #*********************************************************************
>>> # Renormalization and factorization scales *
>>> #*********************************************************************
>>> F = fixed_ren_scale ! if .true. use fixed ren scale
>>> F = fixed_fac_scale ! if .true. use fixed fac scale
>>> 91.1880 = scale ! fixed ren scale
>>> 91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
>>> 91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
>>> 1 = scalefact ! scale factor for event-by-event scales
>>> #*********************************************************************
>>> # Matching - Warning! ickkw > 1 is still beta
>>> #*********************************************************************
>>> 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
>>> 1 = highestmult ! for ickkw=2, highest mult group
>>> 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
>>> 1 = alpsfact ! scale factor for QCD emission vx
>>> F = chcluster ! cluster only according to channel diag
>>> T = pdfwgt ! for ickkw=1, perform pdf reweighting
>>> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
>>> T = clusinfo ! include clustering tag in output
>>> 3.0 = lhe_version ! Change the way clustering information pass to shower.
>>> #*********************************************************************
>>> #**********************************************************
>>> #
>>> #**********************************************************
>>> # Automatic ptj and mjj cuts if xqcut > 0
>>> # (turn off for VBF and single top processes)
>>> #**********************************************************
>>> T = auto_ptj_mjj ! Automatic setting of ptj and mjj
>>> #**********************************************************
>>> #
>>> #**********************************
>>> # BW cutoff (M+/-bwcutoff*Gamma)
>>> #**********************************
>>> 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
>>> #**********************************************************
>>> # Apply pt/E/eta/dr/mij cuts on decay products or not
>>> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
>>> #**********************************************************
>>> T = cut_decays ! Cut decay products
>>> #*************************************************************
>>> # Number of helicities to sum per event (0 = all helicities)
>>> # 0 gives more stable result, but longer run time (needed for
>>> # long decay chains e.g.).
>>> # Use >=2 if most helicities contribute, e.g. pure QCD.
>>> #*************************************************************
>>> 0 = nhel ! Number of helicities used per event
>>> #*******************
>>> # Standard Cuts
>>> #*******************
>>> #
>>> #*********************************************************************
>>> # Minimum and maximum pt's (for max, -1 means no cut) *
>>> #*********************************************************************
>>> 20 = ptj ! minimum pt for the jets
>>> 0 = ptb ! minimum pt for the b
>>> 10 = pta ! minimum pt for the photons
>>> 10 = ptl ! minimum pt for the charged leptons
>>> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
>>> 0 = ptheavy ! minimum pt for one heavy final state
>>> 1.0 = ptonium ! minimum pt for the quarkonium states
>>> -1 = ptjmax ! maximum pt for the jets
>>> -1 = ptbmax ! maximum pt for the b
>>> -1 = ptamax ! maximum pt for the photons
>>> -1 = ptlmax ! maximum pt for the charged leptons
>>> -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
>>> #*********************************************************************
>>> # Minimum and maximum E's (in the center of mass frame) *
>>> #*********************************************************************
>>> 0 = ej ! minimum E for the jets
>>> 0 = eb ! minimum E for the b
>>> 0 = ea ! minimum E for the photons
>>> 0 = el ! minimum E for the charged leptons
>>> -1 = ejmax ! maximum E for the jets
>>> -1 = ebmax ! maximum E for the b
>>> -1 = eamax ! maximum E for the photons
>>> -1 = elmax ! maximum E for the charged leptons
>>> #*********************************************************************
>>> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
>>> #*********************************************************************
>>> 5 = etaj ! max rap for the jets
>>> -1 = etab ! max rap for the b
>>> 2.5 = etaa ! max rap for the photons
>>> 2.5 = etal ! max rap for the charged leptons
>>> 0.6 = etaonium ! max rap for the quarkonium states
>>> 0 = etajmin ! min rap for the jets
>>> 0 = etabmin ! min rap for the b
>>> 0 = etaamin ! min rap for the photons
>>> 0 = etalmin ! main rap for the charged leptons
>>> #*********************************************************************
>>> # Minimum and maximum DeltaR distance *
>>> #*********************************************************************
>>> 0.4 = drjj ! min distance between jets
>>> 0 = drbb ! min distance between b's
>>> 0.4 = drll ! min distance between leptons
>>> 0.4 = draa ! min distance between gammas
>>> 0 = drbj ! min distance between b and jet
>>> 0.4 = draj ! min distance between gamma and jet
>>> 0.4 = drjl ! min distance between jet and lepton
>>> 0 = drab ! min distance between gamma and b
>>> 0 = drbl ! min distance between b and lepton
>>> 0.4 = dral ! min distance between gamma and lepton
>>> -1 = drjjmax ! max distance between jets
>>> -1 = drbbmax ! max distance between b's
>>> -1 = drllmax ! max distance between leptons
>>> -1 = draamax ! max distance between gammas
>>> -1 = drbjmax ! max distance between b and jet
>>> -1 = drajmax ! max distance between gamma and jet
>>> -1 = drjlmax ! max distance between jet and lepton
>>> -1 = drabmax ! max distance between gamma and b
>>> -1 = drblmax ! max distance between b and lepton
>>> -1 = dralmax ! maxdistance between gamma and lepton
>>> #*********************************************************************
>>> # Minimum and maximum invariant mass for pairs *
>>> # WARNING: for four lepton final state mmll cut require to have *
>>> # different lepton masses for each flavor! *
>>> #*********************************************************************
>>> 0 = mmjj ! min invariant mass of a jet pair
>>> 0 = mmbb ! min invariant mass of a b pair
>>> 0 = mmaa ! min invariant mass of gamma gamma pair
>>> 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
>>> -1 = mmjjmax ! max invariant mass of a jet pair
>>> -1 = mmbbmax ! max invariant mass of a b pair
>>> -1 = mmaamax ! max invariant mass of gamma gamma pair
>>> -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
>>> #*********************************************************************
>>> # Minimum and maximum invariant mass for all letpons *
>>> #*********************************************************************
>>> 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
>>> -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
>>> #*********************************************************************
>>> # Minimum and maximum pt for 4-momenta sum of leptons *
>>> #*********************************************************************
>>> 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
>>> -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
>>> #*********************************************************************
>>> # Inclusive cuts *
>>> #*********************************************************************
>>> 0 = xptj ! minimum pt for at least one jet
>>> 0 = xptb ! minimum pt for at least one b
>>> 0 = xpta ! minimum pt for at least one photon
>>> 0 = xptl ! minimum pt for at least one charged lepton
>>> #*********************************************************************
>>> # Control the pt's of the jets sorted by pt *
>>> #*********************************************************************
>>> 0 = ptj1min ! minimum pt for the leading jet in pt
>>> 0 = ptj2min ! minimum pt for the second jet in pt
>>> 0 = ptj3min ! minimum pt for the third jet in pt
>>> 0 = ptj4min ! minimum pt for the fourth jet in pt
>>> -1 = ptj1max ! maximum pt for the leading jet in pt
>>> -1 = ptj2max ! maximum pt for the second jet in pt
>>> -1 = ptj3max ! maximum pt for the third jet in pt
>>> -1 = ptj4max ! maximum pt for the fourth jet in pt
>>> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
>>> #*********************************************************************
>>> # Control the pt's of leptons sorted by pt *
>>> #*********************************************************************
>>> 0 = ptl1min ! minimum pt for the leading lepton in pt
>>> 0 = ptl2min ! minimum pt for the second lepton in pt
>>> 0 = ptl3min ! minimum pt for the third lepton in pt
>>> 0 = ptl4min ! minimum pt for the fourth lepton in pt
>>> -1 = ptl1max ! maximum pt for the leading lepton in pt
>>> -1 = ptl2max ! maximum pt for the second lepton in pt
>>> -1 = ptl3max ! maximum pt for the third lepton in pt
>>> -1 = ptl4max ! maximum pt for the fourth lepton in pt
>>> #*********************************************************************
>>> # Control the Ht(k)=Sum of k leading jets *
>>> #*********************************************************************
>>> 0 = htjmin ! minimum jet HT=Sum(jet pt)
>>> -1 = htjmax ! maximum jet HT=Sum(jet pt)
>>> 0 = ihtmin !inclusive Ht for all partons (including b)
>>> -1 = ihtmax !inclusive Ht for all partons (including b)
>>> 0 = ht2min ! minimum Ht for the two leading jets
>>> 0 = ht3min ! minimum Ht for the three leading jets
>>> 0 = ht4min ! minimum Ht for the four leading jets
>>> -1 = ht2max ! maximum Ht for the two leading jets
>>> -1 = ht3max ! maximum Ht for the three leading jets
>>> -1 = ht4max ! maximum Ht for the four leading jets
>>> #***********************************************************************
>>> # Photon-isolation cuts, according to hep-ph/9801442 *
>>> # When ptgmin=0, all the other parameters are ignored *
>>> # When ptgmin>0, pta and draj are not going to be used *
>>> #***********************************************************************
>>> 0 = ptgmin ! Min photon transverse momentum
>>> 0.4 = R0gamma ! Radius of isolation code
>>> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
>>> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
>>> .true. = isoEM ! isolate photons from EM energy (photons and leptons)
>>> #*********************************************************************
>>> # WBF cuts *
>>> #*********************************************************************
>>> 0 = xetamin ! minimum rapidity for two jets in the WBF case
>>> 0 = deltaeta ! minimum rapidity for two jets in the WBF case
>>> #*********************************************************************
>>> # KT DURHAM CUT *
>>> #*********************************************************************
>>> -1 = ktdurham
>>> 0.4 = dparameter
>>> #*********************************************************************
>>> # maximal pdg code for quark to be considered as a light jet *
>>> # (otherwise b cuts are applied) *
>>> #*********************************************************************
>>> 4 = maxjetflavor ! Maximum jet pdg code
>>> #*********************************************************************
>>> # Jet measure cuts *
>>> #*********************************************************************
>>> 0 = xqcut ! minimum kt jet measure between partons
>>> #*********************************************************************
>>> #
>>> #*********************************************************************
>>> # Store info for systematics studies *
>>> # WARNING: If use_syst is T, matched Pythia output is *
>>> # meaningful ONLY if plotted taking matchscale *
>>> # reweighting into account! *
>>> #*********************************************************************
>>> F = use_syst ! Enable systematics studies
>>> #
>>> #**************************************
>>> # Parameter of the systematics study
>>> # will be used by SysCalc (if installed)
>>> #**************************************
>>> #
>>> 0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
>>> 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
>>> 30 50 = sys_matchscale # variation of merging scale
>>> # PDF sets and number of members (0 or none for all members).
>>> CT10nlo.LHgrid = sys_pdf # matching scales
>>> # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
>>>
>>> BTW, if I change the decay pattern to be
>>> decay h2 > t t~, t > w+ b, w+ > j j, t~ > w- b~, w- > l- vl~
>>> there is no such error anymore. But the decay get stuck at "Decaying the events.....".
>>>
>>> I am really appreciate for your help
>>>
>>> Regards,
>>> Ying-Ying
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https://answers.launchpad.net/mg5amcnlo/+question/281261/+confirm?answer_id=0
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/mg5amcnlo/+question/281261
>>
>> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.