systematics_program input type error

Hi,

You can use the following patch to make it work, this is going to be fixed in 2.6.1.

Cheers,

Olivier

=== modified file 'MadSpin/interface_madspin.py'
--- MadSpin/interface_madspin.py 2017-04-12 13:04:38 +0000
+++ MadSpin/interface_madspin.py 2017-09-07 13:43:52 +0000
@@ -1021,7 +1021,7 @@

                    run_card["iseed"] = self.seed
                    run_card['gridpack'] = True
- run_card['systematics_program'] = False
+ run_card['use_syst'] = False
                    run_card.write(pjoin(decay_dir, "Cards", "run_card.dat"))
                    param_card = self.banner['slha']
                    open(pjoin(decay_dir, "Cards", "param_card.dat"),"w").write(param_card)

> On 7 Oct 2017, at 22:27, Sihyun Jeon <email address hidden> wrote:
>
> New question #659116 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/659116
>
> Hi MadGraph team,
>
> I am having trouble making gridpacks with the following error.
> The generation stops when madspin starts to run.
>
> I put in
> =====================================================================
> set ms_dir ./madspingrid
> set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> set max_running_process 1
>
> set spinmode onshell
> decay n1 > mu+ j j
> decay n1 > mu- j j
>
> # running the actual code
> launch
> =======================================================
>
> in the madspin.
> It's strange that this only happens when m(n1) < m(W).
> This does not happen when m(n1) > m(W) (in this case, decay n1 > mu+ w-, w- > j j i sused).
> The error log is below.
>
> Thank you in advance.
> Sincerely, Sihyun Jeon
>
>
>
>
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.2.6.0 20xx-xx-xx *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* and *
> #* http://amcatnlo.cern.ch *
> #* *
> #************************************************************
> #* *
> #* Command File for aMCatNLO *
> #* *
> #* run as ./bin/aMCatNLO.py filename *
> #* *
> #************************************************************
> launch -n pilotrun
> Traceback (most recent call last):
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 1438, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
> return func(arg, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/amcatnlo_run_interface.py", line 1243, in do_launch
> self.exec_cmd('decay_events -from_cards', postcmd=False)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 1465, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
> return func(arg, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/processtmp/bin/internal/common_run_interface.py", line 3463, in do_decay_events
> madspin_cmd.import_command_file(path)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1580, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1465, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
> return func(arg, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/various/misc.py", line 99, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/MadSpin/interface_madspin.py", line 538, in do_launch
> return self.run_onshell(line)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/MadSpin/interface_madspin.py", line 1192, in run_onshell
> evt_decayfile[pdg], pwidth = self.generate_events(pdg, nb_needed, mg5, output_width=True)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/MadSpin/interface_madspin.py", line 1024, in generate_events
> run_card['systematics_program'] = False
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/various/banner.py", line 1118, in __setitem__
> value = self.format_variable(value, targettype, name=name)
> File "/afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/madgraph/various/banner.py", line 1191, in format_variable
> (name, type(value), targettype, value)
> Exception: Wrong input type for systematics_program found <type 'bool'> and expecting <type 'str'> for value False
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_local_path : /cvmfs/cms.cern.ch/slc6_amd64_gcc481/external/lhapdf/6.1.6-ddibom/share/LHAPDF
> cluster_queue : local
> madanalysis_path : None
> lhapdf : /cvmfs/cms.cern.ch/slc6_amd64_gcc481/external/lhapdf/6.1.6-ddibom/share/LHAPDF/../../bin/lhapdf-config
> cluster_size : 100
> cluster_memory : None
> pjfry : None
> cluster_status_update : (600, 30)
> cluster_time : None
> f2py_compiler : None
> ninja : /afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/HEPTools/lib
> hepmc_path : None
> pythia8_path : None
> hwpp_path : None
> automatic_html_opening : False
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> mg5_path : /afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0
> td_path : None
> delphes_path : None
> thepeg_path : None
> cluster_type : condor
> madanalysis5_path : None
> exrootanalysis_path : None
> fortran_compiler : None
> nb_core : 16
> collier : /afs/cern.ch/work/s/shjeon/genproductions_2017_mg26x/bin/MadGraph5_aMCatNLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO/HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO_gridpack/work/MG5_aMC_v2_6_0/HEPTools/lib
> auto_update : 0
> cluster_nb_retry : 1
> eps_viewer : None
> syscalc_path : None
> fastjet : None
> golem : None
> cpp_compiler : None
> notification_center : True
> run_mode : 2
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.0 2017-08-16 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model SM_HeavyN_NLO
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
> define j = p
> define p = u c d s b u~ c~ d~ s~ b~ g
> define j = u c d s b u~ c~ d~ s~ b~ g
> define mu = mu+ mu-
> define e = e+ e-
> define w = w+ w-
> generate p p > n1 mu [QCD]
> output HeavyMajoranaNeutrinoToMuMu_Schannel_M60_NLO
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR LOOP
> ###################################
> BLOCK LOOP #
> 1 9.118800e+01 # mu_r
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 6 1.733000e+02 # mt
> 23 9.118760e+01 # mz
> 25 1.257000e+02 # mh
> 9900012 6.000000e+01 # mn1
> 9900014 5.000000e+02 # mn2
> 9900016 1.000000e+03 # mn3
> 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
> 5 0.000000e+00 # b : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 0.0
> 14 0.000000e+00 # vm : 0.0
> 15 0.000000e+00 # ta- : 0.0
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 7.995123e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
> 9000002 9.118760e+01 # ghz : mz
> 9000003 7.995123e+01 # ghwp : mw
> 9000004 7.995123e+01 # ghwm : mw
> ###################################
> ## INFORMATION FOR NUMIXING
> ###################################
> BLOCK NUMIXING #
> 1 1.000000e-01 # ven1
> 2 0.000000e+00 # ven2
> 3 0.000000e+00 # ven3
> 4 1.000000e-01 # vmun1
> 5 0.000000e+00 # vmun2
> 6 0.000000e+00 # vmun3
> 7 0.000000e+00 # vtan1
> 8 0.000000e+00 # vtan2
> 9 0.000000e+00 # vtan3
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.279400e+02 # aewm1
> 2 1.174560e-05 # gf
> 3 1.184000e-01 # as
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 6 1.733000e+02 # ymt
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> BLOCK QNUMBERS 9900012 # 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 9900014
> ###################################
> BLOCK QNUMBERS 9900014 # 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 9900016
> ###################################
> BLOCK QNUMBERS 9900016 # 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)
> #
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.350000e+00
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> #
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.495200e+00
> #
> # PDG Width
> DECAY 24 2.085000e+00
> #
> # PDG Width
> DECAY 25 4.170000e-03
> #
> # PDG Width
> DECAY 9000002 2.495200e+00
> #
> # PDG Width
> DECAY 9000003 2.085000e+00
> #
> # PDG Width
> DECAY 9000004 2.085000e+00
> #
> # PDG Width
> DECAY 9900012 1.298704e-05
> # BR NDA ID1 ID2 ...
> 6.528047e-02 3 -13 -2 1 # 8.47800075109e-07
> 6.528047e-02 3 -13 -4 3 # 8.47800075109e-07
> 6.469143e-02 3 -11 -2 1 # 8.40150189067e-07
> 6.469143e-02 3 -11 -4 3 # 8.40150189067e-07
> 6.459133e-02 3 -1 2 13 # 8.38850186363e-07
> 6.459133e-02 3 -3 4 13 # 8.38850186363e-07
> 6.436033e-02 3 -1 2 11 # 8.35850180123e-07
> 6.436033e-02 3 -3 4 11 # 8.35850180123e-07
> 2.143791e-02 3 -13 11 14 # 2.78414994686e-07
> 2.143791e-02 3 -15 11 16 # 2.78414994686e-07
> 2.142598e-02 3 -13 -12 11 # 2.78260059299e-07
> 2.142598e-02 3 -16 -13 15 # 2.78260059299e-07
> 2.134975e-02 3 -11 12 13 # 2.7727005724e-07
> 2.134975e-02 3 -15 13 16 # 2.7727005724e-07
> 2.131163e-02 3 -14 -11 13 # 2.76774991275e-07
> 2.131163e-02 3 -16 -11 15 # 2.76774991275e-07
> 1.332790e-02 3 -11 11 12 # 1.73089970416e-07
> 1.329402e-02 3 -14 -13 13 # 1.72649969501e-07
> 1.327400e-02 3 -12 -11 11 # 1.7238996896e-07
> 1.318853e-02 3 -13 13 14 # 1.71279966651e-07
> 1.060391e-02 3 -1 1 14 # 1.37713403326e-07
> 1.060391e-02 3 -3 3 14 # 1.37713403326e-07
> 1.060391e-02 3 -5 5 14 # 1.37713403326e-07
> 1.051818e-02 3 -12 -1 1 # 1.36600024387e-07
> 1.051818e-02 3 -12 -3 3 # 1.36600024387e-07
> 1.051818e-02 3 -12 -5 5 # 1.36600024387e-07
> 1.044349e-02 3 -14 -1 1 # 1.3563002237e-07
> 1.044349e-02 3 -14 -3 3 # 1.3563002237e-07
> 1.044349e-02 3 -14 -5 5 # 1.3563002237e-07
> 1.040011e-02 3 -1 1 12 # 1.35066644574e-07
> 1.040011e-02 3 -3 3 12 # 1.35066644574e-07
> 1.040011e-02 3 -5 5 12 # 1.35066644574e-07
> 8.207029e-03 3 -14 -2 2 # 1.06585013904e-07
> 8.207029e-03 3 -14 -4 4 # 1.06585013904e-07
> 8.178924e-03 3 -12 -2 2 # 1.06220013145e-07
> 8.178924e-03 3 -12 -4 4 # 1.06220013145e-07
> 8.160829e-03 3 -2 2 14 # 1.05985012656e-07
> 8.160829e-03 3 -4 4 14 # 1.05985012656e-07
> 8.121174e-03 3 -2 2 12 # 1.05470011585e-07
> 8.121174e-03 3 -4 4 12 # 1.05470011585e-07
> 4.820961e-03 3 -14 -12 12 # 6.26100133454e-08
> 4.820961e-03 3 -16 -14 16 # 6.26100133454e-08
> 4.779766e-03 3 -12 12 14 # 6.20750122326e-08
> 4.779766e-03 3 -16 14 16 # 6.20750122326e-08
> 4.771681e-03 3 -14 12 14 # 6.19700120142e-08
> 4.771681e-03 3 -16 12 16 # 6.19700120142e-08
> 4.743191e-03 3 -14 -12 14 # 6.16000112446e-08
> 4.743191e-03 3 -16 -12 16 # 6.16000112446e-08
> 2.417641e-03 3 -14 -14 14 # 3.13980003726e-08
> 2.391307e-03 3 -13 12 13 # 3.10559996613e-08
> 2.391307e-03 3 -15 12 15 # 3.10559996613e-08
> 2.388920e-03 3 -14 -11 11 # 3.10249995968e-08
> 2.388920e-03 3 -15 -14 15 # 3.10249995968e-08
> 2.383877e-03 3 -11 11 14 # 3.09595059541e-08
> 2.383877e-03 3 -15 14 15 # 3.09595059541e-08
> 2.369054e-03 3 -14 14 14 # 3.07669990602e-08
> 2.364896e-03 3 -13 -12 13 # 3.07129989478e-08
> 2.364896e-03 3 -15 -12 15 # 3.07129989478e-08
> 2.357273e-03 3 -12 12 12 # 3.06139987419e-08
> 2.351729e-03 3 -12 -12 12 # 3.05419985922e-08
> #
> # PDG Width
> DECAY 9900014 1.500000e+00
> #
> # PDG Width
> DECAY 9900016 1.230000e+01
> #***********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat aMC@NLO *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a hash (#) are info or comments *
> # *
> # mind the format: value = variable ! comment *
> #***********************************************************************
> #
> #*******************
> # Running parameters
> #*******************
> #
> #***********************************************************************
> # Tag name for the run (one word) *
> #***********************************************************************
> tag_1 = run_tag ! name of the run
> #***********************************************************************
> # Number of events (and their normalization) and the required *
> # (relative) accuracy on the Xsec. *
> # These values are ignored for fixed order runs *
> #***********************************************************************
> 1000 = nevents ! Number of unweighted events requested
> 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 100 = nevt_job! Max number of events per job in event generation.
> ! (-1= no split).
> average = event_norm ! Normalize events to sum or average to the X sect.
> #***********************************************************************
> # Number of points per itegration channel (ignored for aMC@NLO runs) *
> #***********************************************************************
> 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
> ! number of points and iter. below)
> # These numbers are ignored except if req_acc_FO is equal to -1
> 5000 = npoints_FO_grid ! number of points to setup grids
> 4 = niters_FO_grid ! number of iter. to setup grids
> 10000 = npoints_FO ! number of points to compute Xsec
> 6 = niters_FO ! number of iter. to compute Xsec
> #***********************************************************************
> # Random number seed *
> #***********************************************************************
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #***********************************************************************
> # Collider type and energy *
> #***********************************************************************
> 1 = lpp1 ! beam 1 type (0 = no PDF)
> 1 = lpp2 ! beam 2 type (0 = no PDF)
> 6500 = ebeam1 ! beam 1 energy in GeV
> 6500 = ebeam2 ! beam 2 energy in GeV
> #***********************************************************************
> # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> #***********************************************************************
> lhapdf = pdlabel ! PDF set
> 90200 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> #***********************************************************************
> # Include the NLO Monte Carlo subtr. terms for the following parton *
> # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> #***********************************************************************
> PYTHIA8 = parton_shower
> 1.0 = shower_scale_factor ! multiply default shower starting
> ! scale by this factor
> #***********************************************************************
> # Renormalization and factorization scales *
> # (Default functional form for the non-fixed scales is the sum of *
> # the transverse masses of all final state particles and partons. This *
> # can be changed in SubProcesses/set_scales.f) *
> #***********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.118 = muR_ref_fixed ! fixed ren reference scale
> 91.118 = muF1_ref_fixed ! fixed fact reference scale for pdf1
> 91.118 = muF2_ref_fixed ! fixed fact reference scale for pdf2
> -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> #***********************************************************************
> # Renormalization and factorization scales (advanced and NLO options) *
> #***********************************************************************
> False = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
> 91.118 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
> 1 = muR_over_ref ! ratio of current muR over reference muR
> 1 = muF1_over_ref ! ratio of current muF1 over reference muF1
> 1 = muF2_over_ref ! ratio of current muF2 over reference muF2
> 1 = QES_over_ref ! ratio of current QES over reference QES
> #***********************************************************************
> # Reweight flags to get scale dependence and PDF uncertainty *
> # For scale dependence: factor rw_scale_up/down around central scale *
> # For PDF uncertainty: use LHAPDF with supported set *
> #***********************************************************************
> True = reweight_scale ! reweight to get scale dependence
> 0.5 = rw_Rscale_down ! lower bound for ren scale variations
> 2.0 = rw_Rscale_up ! upper bound for ren scale variations
> 0.5 = rw_Fscale_down ! lower bound for fact scale variations
> 2.0 = rw_Fscale_up ! upper bound for fact scale variations
> True = reweight_PDF ! reweight to get PDF uncertainty
> 90201 = PDF_set_min ! First of the error PDF sets
> 90303 = PDF_set_max ! Last of the error PDF sets
> #***********************************************************************
> # Merging - WARNING! Applies merging only at the hard-event level. *
> # After showering an MLM-type merging should be applied as well. *
> # See http://amcatnlo.cern.ch/FxFx_merging.htm for more details. *
> #***********************************************************************
> 0 = ickkw ! 0 no merging, 3 FxFx merging
> #***********************************************************************
> #
> #***********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) *
> #***********************************************************************
> 15 = bwcutoff
> #***********************************************************************
> # Cuts on the jets *
> # Jet clustering is performed by FastJet.
> # When matching to a parton shower, these generation cuts should be *
> # considerably softer than the analysis cuts. *
> # (more specific cuts can be specified in SubProcesses/cuts.f) *
> #***********************************************************************
> 1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 0.7 = jetradius ! The radius parameter for the jet algorithm
> 0.0001 = ptj ! Min jet transverse momentum
> -1 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> #***********************************************************************
> # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> # (more specific gen cuts can be specified in SubProcesses/cuts.f) *
> #***********************************************************************
> 0 = ptl ! Min lepton transverse momentum
> -1 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 0 = drll ! Min distance between opposite sign lepton pairs
> 0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 4 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442 *
> # When ptgmin=0, all the other parameters are ignored *
> #***********************************************************************
> 0 = ptgmin ! Min photon transverse momentum
> -1 = etagamma ! Max photon abs(pseudo-rap)
> 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)
> #***********************************************************************
> # Maximal PDG code for quark to be considered a jet when applying cuts.*
> # At least all massless quarks of the model should be included here. *
> #***********************************************************************
> 5 = maxjetflavor
> #***********************************************************************
> # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> #***********************************************************************
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
> #***********************************************************************
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.