Error: End code 21.0

Hi,

You set the dynamical scale to 0 "user custom function" but you did not define such function
(or not correctly) since the code crash with
" User-defined scale not set
 STOP 21
"

Cheers,

Olivier

> On 17 May 2021, at 21:30, Ruoquan Wang <email address hidden> wrote:
>
> New question #697107 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/697107
>
> A fixed target experiment simulation, the script for mg5 is here:
> -----------------------------------------------------------------------------------
> import model sm-lepton_masses
> define p g u c d s u~ c~ d~ s~
> define j g u c d s u~ c~ d~ s~
> generate vm p > mu- j
> output ~/Work/tryvm
>
> launch
> set lpp1 0
> set polbeam1 = -100
> set ebeam2 0.
> set no_parton_cut
> set nevents 10000.
> set dynamical_scale_choice 0
> set use_syst False
> set nb_proton2 1 #!
> set nb_neutron2 0
> set mass_ion2 0.93827231 #!
> set ebeam1 1e6
> -----------------------------------------------------------------------------------
>
> My mg5 version is 3.1.0, on MacOS 10.15.
> I get:
>
> Error when reading /Users/isaac/Work/tryvm/SubProcesses/P1_vlq_mumq/G1/results.dat
> Command "generate_events run_01" interrupted with error:
> Exception : Reported error: End code 21.0
> Full associated log:
> Process in group number 1
> Setting beam polarization -100.00000000000000 0.0000000000000000
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> ****************************************
>
> NNPDFDriver version 1.0.3
> Grid: NNPDF23_lo_as_0130_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23lo1: 0.13000000000000000
> Define smin to 1.1164035600000000E-002
> Define smin to 1.1164035600000000E-002
> Warning: pt or E min of a jet should in general be >0
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_yme = 5.1099999999999995E-004
> mdl_ymm = 0.10566000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MH = 125.00000000000000
> mdl_Me = 5.1099999999999995E-004
> mdl_MM = 0.10566000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WT = 1.4915000000000000
> mdl_WH = 6.3823389999999999E-003
> mdl_WTau = 2.2699999999999998E-012
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM1x1 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_complexi = (0.0000000000000000,1.0000000000000000)
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_ye = 2.9350485578697581E-006
> mdl_ym = 6.0688303449025181E-004
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_I1x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_I2x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I3x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I4x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm-lepton_masses
> ---------------------------------
>
> GC_100 = 0.00000E+00 0.46191E+00
>
> Collider parameters:
> --------------------
>
> Running at e P machine @ 1937.2891472364161 GeV
> PDF set = nn23lo1
> alpha_s(Mz)= 0.1300 running at 2 loops.
> alpha_s(Mz)= 0.1300 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> getting user params
> Enter number of events and max and min iterations:
> Number of events and iterations 1000 5 3
> Enter desired fractional accuracy:
> Desired fractional accuracy: 0.10000000000000001
> Enter 0 for fixed, 2 for adjustable grid:
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Explicitly summing over helicities
> Enter Configuration Number:
> Running Configuration Number: 1
> Not subdividing B.W.
> Attempting mappinvarients 1 4
> Determine nb_t
> T-channel found: 1
> Completed mapping 4
> about to integrate 3 1000 5 3 3 1
> Using non-zero grid deformation.
> 3 dimensions 1000 events 3 invarients 5 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 16
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 1369.87
> Mapping Graph 1 to config 1
> Determine nb_t
> T-channel found: 1
> Warning: No cutoff for shat integral found
> Minimum set to 1.0000000000000000E-004
> Setting grid 1 0.10000E-03 1
> Warning: No minimum found for integration
> Setting minimum to 5.3289408640017579E-007
> Transforming s_hat 1/s 3 5.3289408640017579E-007 1.1164035600000000E-002 1876545.5003549277
> 1 1 2 3
> Masses: 0.000E+00 0.000E+00 0.106E+00 0.000E+00
> Using random seed offsets 1 : 1
> with seed 24
> Ranmar initialization seeds 31059 9398
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Particle 3 4
> Et > 0.0 0.0
> E > 0.0 0.0
> Eta < -1.0 -1.0
> xqcut: 0.0 0.0
> d R # 3 > -0.0 0.0
> s min # 3> 0.0 0.0
> xqcutij # 3> 0.0 0.0
> User-defined scale not set
> STOP 21
>
> ls status:
> input_app.txt
> run1_app.log
>
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/Users/isaac/Work/tryvm/run_01_tag_1_debug.log'.
> Please attach this file to your report.
>
> The log file is here:
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.1.0 2021-03-30 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadEvent *
> #* *
> #* run as ./bin/madevent.py filename *
> #* *
> #************************************************************
> generate_events run_01
> Traceback (most recent call last):
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/madevent/sum_html.py", line 307, in read_results
> self.axsec, self.xerru, self.xerrc, self.nevents, self.nw,\
> ValueError: not enough values to unpack (expected 10, got 5)
>
> During handling of the above exception, another exception occurred:
>
> Traceback (most recent call last):
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1544, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/madevent_interface.py", line 2474, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/common_run_interface.py", line 7114, in new_fct
> original_fct(obj, *args, **opts)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/madevent_interface.py", line 2517, in run_generate_events
> self.exec_cmd('refine %s' % nb_event, postcmd=False)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/interface/madevent_interface.py", line 3512, in do_refine
> x_improve.launch() # create the ajob for the refinment.
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/madevent/gen_ximprove.py", line 1048, in launch
> self.results = sum_html.collect_result(self.cmd,
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/madevent/sum_html.py", line 754, in collect_result
> P_comb.add_results(os.path.basename(G), path, mfactors[G])
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/madevent/sum_html.py", line 433, in add_results
> oneresult.read_results(filepath)
> File "/Users/isaac/Work/MG5_aMC_v3_1_0/madgraph/madevent/sum_html.py", line 316, in read_results
> raise Exception("Reported error: End code %s \n Full associated log: \n%s"\
> Exception: Reported error: End code 21.0
> Full associated log:
> Process in group number 1
> Setting beam polarization -100.00000000000000 0.0000000000000000
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> ****************************************
>
> NNPDFDriver version 1.0.3
> Grid: NNPDF23_lo_as_0130_qed_mem0.grid
> ****************************************
> New value of alpha_s from PDF nn23lo1: 0.13000000000000000
> Define smin to 1.1164035600000000E-002
> Define smin to 1.1164035600000000E-002
> Warning: pt or E min of a jet should in general be >0
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_yme = 5.1099999999999995E-004
> mdl_ymm = 0.10566000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MT = 173.00000000000000
> mdl_MB = 4.7000000000000002
> mdl_MH = 125.00000000000000
> mdl_Me = 5.1099999999999995E-004
> mdl_MM = 0.10566000000000000
> mdl_MTA = 1.7769999999999999
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WT = 1.4915000000000000
> mdl_WH = 6.3823389999999999E-003
> mdl_WTau = 2.2699999999999998E-012
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM1x1 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_complexi = (0.0000000000000000,1.0000000000000000)
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_ye = 2.9350485578697581E-006
> mdl_ym = 6.0688303449025181E-004
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_I1x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_I2x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I3x33 = (0.99366614581500623,0.0000000000000000)
> mdl_I4x33 = (2.69955542504654901E-002,0.0000000000000000)
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm-lepton_masses
> ---------------------------------
>
> GC_100 = 0.00000E+00 0.46191E+00
>
> Collider parameters:
> --------------------
>
> Running at e P machine @ 1937.2891472364161 GeV
> PDF set = nn23lo1
> alpha_s(Mz)= 0.1300 running at 2 loops.
> alpha_s(Mz)= 0.1300 running at 2 loops.
> Renormalization scale set on event-by-event basis
> Factorization scale set on event-by-event basis
>
>
> getting user params
> Enter number of events and max and min iterations:
> Number of events and iterations 1000 5 3
> Enter desired fractional accuracy:
> Desired fractional accuracy: 0.10000000000000001
> Enter 0 for fixed, 2 for adjustable grid:
> Suppress amplitude (0 no, 1 yes)?
> Using suppressed amplitude.
> Exact helicity sum (0 yes, n = number/event)?
> Explicitly summing over helicities
> Enter Configuration Number:
> Running Configuration Number: 1
> Not subdividing B.W.
> Attempting mappinvarients 1 4
> Determine nb_t
> T-channel found: 1
> Completed mapping 4
> about to integrate 3 1000 5 3 3 1
> Using non-zero grid deformation.
> 3 dimensions 1000 events 3 invarients 5 iterations 1 config(s), (0.99)
> Using h-tuple random number sequence.
> Error opening grid
> Using Uniform Grid! 16
> Using uniform alpha 1.0000000000000000
> Grid defined OK
> Set CM energy to 1369.87
> Mapping Graph 1 to config 1
> Determine nb_t
> T-channel found: 1
> Warning: No cutoff for shat integral found
> Minimum set to 1.0000000000000000E-004
> Setting grid 1 0.10000E-03 1
> Warning: No minimum found for integration
> Setting minimum to 5.3289408640017579E-007
> Transforming s_hat 1/s 3 5.3289408640017579E-007 1.1164035600000000E-002 1876545.5003549277
> 1 1 2 3
> Masses: 0.000E+00 0.000E+00 0.106E+00 0.000E+00
> Using random seed offsets 1 : 1
> with seed 24
> Ranmar initialization seeds 31059 9398
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Particle 3 4
> Et > 0.0 0.0
> E > 0.0 0.0
> Eta < -1.0 -1.0
> xqcut: 0.0 0.0
> d R # 3 > -0.0 0.0
> s min # 3> 0.0 0.0
> xqcutij # 3> 0.0 0.0
> User-defined scale not set
> STOP 21
>
> ls status:
> input_app.txt
> run1_app.log
>
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> run_mode : 2
> cluster_queue : None (user set)
> cluster_time : None (user set)
> cluster_size : 100
> cluster_memory : 100 (user set)
> nb_core : 12 (user set)
> cluster_temp_path : None
>
> Configuration Options
> ---------------------
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> thepeg_path : None (user set)
> hepmc_path : None (user set)
> madanalysis_path : None (user set)
> madanalysis5_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> exrootanalysis_path : /Users/isaac/Work/MG5_aMC_v3_1_0/ExRootAnalysis (user set)
> syscalc_path : None (user set)
> lhapdf : /Users/isaac/Work/lhaPDF/bin/lhapdf-config (user set)
> lhapdf_py2 : None
> lhapdf_py3 : None
> timeout : 60
> f2py_compiler : None
> f2py_compiler_py2 : None
> f2py_compiler_py3 : None
> web_browser : None
> eps_viewer : None
> text_editor : None
> fortran_compiler : None
> cpp_compiler : None
> auto_update : 7 (user set)
> cluster_type : condor
> cluster_status_update : (600, 30)
> cluster_nb_retry : 1
> cluster_local_path : None
> cluster_retry_wait : 300
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.1.0 2021-03-30 *
> #* *
> #* 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 low_mem_multicore_nlo_generation False
> set complex_mass_scheme False
> set include_lepton_initiated_processes False
> set gauge unitary
> set loop_optimized_output True
> set loop_color_flows False
> set max_npoint_for_channel 0
> set default_unset_couplings 99
> set max_t_for_channel 99
> set zerowidth_tchannel True
> set nlo_mixed_expansion True
> 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 sm-lepton_masses
> generate vm p > mu- j
> set lhapdf /Users/isaac/Work/lhaPDF/bin/lhapdf-config
> output ~/Work/tryvm
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 11 5.110000e-04 # Me
> 13 1.056600e-01 # MM
> 15 1.777000e+00 # MTA
> 23 9.118800e+01 # MZ
> 25 1.250000e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.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
> 12 0.000000e+00 # ve : 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
> 24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 1.325070e+02 # aEWM1
> 2 1.166390e-05 # Gf
> 3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set)
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 11 5.110000e-04 # yme
> 13 1.056600e-01 # ymm
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 15 2.270000e-12 # WTau
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.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 16 0.000000e+00 # vt : 0.0
> DECAY 21 0.000000e+00 # g : 0.0
> DECAY 22 0.000000e+00 # a : 0.0
> #*********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat MadEvent *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a '# ' are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # To display more options, you can type the command: *
> # update full_run_card *
> #*********************************************************************
> #
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_1 = run_tag ! name of the run
> #*********************************************************************
> # Number of events and rnd seed *
> # Warning: Do not generate more than 1M events in a single run *
> #*********************************************************************
> 10000 = nevents ! Number of unweighted events requested
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #*********************************************************************
> # Collider type and energy *
> # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
> # 3=photon from electron, 4=photon from muon *
> #*********************************************************************
> 0 = lpp1 ! beam 1 type
> 1 = lpp2 ! beam 2 type
> 1000000.0 = ebeam1 ! beam 1 total energy in GeV
> 0.0 = ebeam2 ! beam 2 total energy in GeV
> #*********************************************************************
> # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> #*********************************************************************
> -100.0 = polbeam1 ! beam polarization for beam 1
> 0.0 = polbeam2 ! beam polarization for beam 2
>
>
> #*********************************************************************
> # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> #*********************************************************************
> nn23lo1 = pdlabel ! PDF set
> 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
>
> #*********************************************************************
> # Heavy ion PDF / rescaling of PDF *
> #*********************************************************************
> 1 = nb_proton1 # number of proton for the first beam
> 0 = nb_neutron1 # number of neutron for the first beam
> -1.0 = mass_ion1 # mass of the heavy ion (first beam)
> # Note that seting differently the two beams only work if you use
> # group_subprocess=False when generating your matrix-element
> 1 = nb_proton2 # number of proton for the second beam
> 0 = nb_neutron2 # number of neutron for the second beam
> 0.93827231 = mass_ion2 # mass of the heavy ion (second beam)
>
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> 0 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> 1.0 = scalefact ! scale factor for event-by-event scales
> #*********************************************************************
> # Type and output format
> #*********************************************************************
> False = gridpack !True = setting up the grid pack
> -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
> average = event_norm ! average/sum. Normalization of the weight in the LHEF
> # To see MLM/CKKW merging options: type "update MLM" or "update CKKW"
>
> #*********************************************************************
> #
> #*********************************************************************
> # Phase-Space Optimization strategy (basic options)
> #*********************************************************************
> 0 = nhel ! using helicities importance sampling or not.
> ! 0: sum over helicity, 1: importance sampling
> 2 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx)
> ! 1 is old strategy (using amp square)
> ! 2 is new strategy (using only the denominator)
> # To see advanced option for Phase-Space optimization: type "update psoptim"
> #*********************************************************************
> # Generation bias, check the wiki page below for more information: *
> # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
> #*********************************************************************
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
> #
> #*******************************
> # Parton level cuts definition *
> #*******************************
> #
> #
> #*********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
> #*********************************************************************
> 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 0.0 = ptj ! minimum pt for the jets
> 0.0 = ptl ! minimum pt for the charged leptons
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptlmax ! maximum pt for the charged leptons
> {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
> #
> # For display option for energy cut in the partonic center of mass frame type 'update ecut'
> #
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> -1.0 = etaj ! max rap for the jets
> -1.0 = etal ! max rap for the charged leptons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.0 = drjl ! min distance between jet and lepton
> -1.0 = drjlmax ! max distance between jet and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
> ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: Do not use for interference type of computation *
> #*********************************************************************
> False = use_syst ! Enable systematics studies
> #
> systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
> ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
> # Syscalc is deprecated but to see the associate options type'update syscalc'
>
>
> Any hints about this?
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.