Error: End code 21.0

Asked by Ruoquan Wang on 2021-05-17

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?

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

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.

Revision history for this message
Ruoquan Wang (isaac95) said :
#2

Thanks Olivier Mattelaer, that solved my question.