max() iterable argument

Hi,

I guess that the issue is that the code that determines the channel of integration is failing to link to a library or something similar.

Can you go to do SubProcesses directory
inside you should have a directory with a name starting with P
something like
P0_vmar_vmara
go inside that directory as well and run
./gensym
(if you do not have such executable --unlikely--) then do first "make gensym"

I guess that this will reveal an hardware related issue (or maybe a model specific one?)

Cheers,

Olivier

PS: Do you have other model/process working?

> On 30 Sep 2024, at 02:25, Sam Carey <email address hidden> wrote:
>
> New question #818864 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/818864
>
> Hi,
> I'm trying to generate events for an (NP) interaction on Argon nucleus: vm ar > vm ar a (with mass 37.52 GeV and spin 0). I created a vertex for p.ar p.ar p.a like that of a scalar-photon interaction. But I get this error when trying to run the events for ebeam1 = 1 GeV, ebeam2 = 37.52 GeV.: ValueError : max() iterable argument is empty.
>
>
> generate_events run_01
> Traceback (most recent call last):
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1543, in onecmd
> return self.onecmd_orig(line, **opt)
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
> return func(arg, **opt)
> ^^^^^^^^^^^^^^^^
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2403, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/common_run_interface.py", line 7656, in new_fct
> original_fct(obj, *args, **opts)
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 2597, in run_generate_events
> self.exec_cmd('survey %s %s' % (self.run_name,' '.join(args)),
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1572, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/extended_cmd.py", line 1492, in onecmd_orig
> return func(arg, **opt)
> ^^^^^^^^^^^^^^^^
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/interface/madevent_interface.py", line 3456, in do_survey
> jobs, P_zero_result = ajobcreator.get_helicity()
> ^^^^^^^^^^^^^^^^^^^^^^^^^^
> File "/Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/madgraph/madevent/gen_ximprove.py", line 159, in get_helicity
> nb_channel = max([math.floor(float(d)) for d in stdout.split()])
> ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
> ValueError: max() iterable argument is empty
> 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 : 8 (user set)
> cluster_temp_path : None
>
> Configuration Options
> ---------------------
> pythia8_path : /Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/HEPTools/pythia8 (user set)
> hwpp_path : None (user set)
> thepeg_path : None (user set)
> hepmc_path : None (user set)
> madanalysis_path : None (user set)
> madanalysis5_path : /Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/HEPTools/madanalysis5/madanalysis5 (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> exrootanalysis_path : None (user set)
> syscalc_path : None (user set)
> rivet_path : None
> yoda_path : None
> lhapdf : /Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
> lhapdf_py2 : None
> lhapdf_py3 : /Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
> 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.5.1 2023-07-11 *
> #* *
> #* 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 /Users/samcarey/Downloads/HEPSoftware/MG5_aMC_v3_5_1/mode\
> ls/Raylnulightscalar_Form_Factor_UFO
> generate vm ar > vm ar a
> output aaaa
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR FRBLOCK
> ###################################
> BLOCK FRBLOCK #
> 1 1.000000e+00 # c71
> 2 1.000000e+00 # c1
> 3 1.000000e+00 # c2
> 4 5.000000e-02 # msigma
> 5 1.000000e+00 # mn
> 7 3.752000e+01 # mar
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 1 5.040000e-03 # md
> 2 2.550000e-03 # mu
> 3 1.010000e-01 # ms
> 4 1.270000e+00 # mc
> 5 4.700000e+00 # mb
> 6 1.720000e+02 # mt
> 11 5.110000e-04 # me
> 13 1.056600e-01 # mmu
> 15 1.777000e+00 # mta
> 23 9.118760e+01 # mz
> 25 1.250000e+02 # mh
> 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.041823e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
> 9000005 5.000000e-02 # sigma : msigma
> 9000006 1.000000e+00 # n : mn
> 9000008 3.752000e+01 # ar : mar
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.325070e+02 # aewm1
> 2 1.166370e-05 # gf
> 3 1.300000e-01 # as (note that parameter not used if you use a pdf set)
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 1 5.040000e-03 # ymdo
> 2 2.550000e-03 # ymup
> 3 1.010000e-01 # yms
> 4 1.270000e+00 # ymc
> 5 4.700000e+00 # ymb
> 6 1.720000e+02 # ymt
> 11 5.110000e-04 # yme
> 13 1.056600e-01 # ymm
> 15 1.777000e+00 # ymtau
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.508336e+00 # wt
> DECAY 23 2.495200e+00 # wz
> DECAY 24 2.085000e+00 # ww
> DECAY 25 4.070000e-03 # wh
> DECAY 9000005 1.000000e+00 # wsigma
> DECAY 9000006 1.000000e+00 # wn
> DECAY 9000008 1.000000e+00 # war
> DECAY 1 0.000000e+00 # d : 0.0
> DECAY 2 0.000000e+00 # u : 0.0
> DECAY 3 0.000000e+00 # s : 0.0
> DECAY 4 0.000000e+00 # c : 0.0
> DECAY 5 0.000000e+00 # b : 0.0
> DECAY 11 0.000000e+00 # e- : 0.0
> DECAY 12 0.000000e+00 # ve : 0.0
> DECAY 13 0.000000e+00 # mu- : 0.0
> DECAY 14 0.000000e+00 # vm : 0.0
> DECAY 15 0.000000e+00 # ta- : 0.0
> DECAY 16 0.000000e+00 # vt : 0.0
> DECAY 21 0.000000e+00 # g : 0.0
> DECAY 22 0.000000e+00 # a : 0.0
> ###################################
> ## INFORMATION FOR QNUMBERS 9000005
> ###################################
> BLOCK QNUMBERS 9000005 # sigma
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 9000006
> ###################################
> BLOCK QNUMBERS 9000006 # n
> 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 1 # particle/antiparticle distinction (0=own anti)
> ###################################
> ## INFORMATION FOR QNUMBERS 9000008
> ###################################
> BLOCK QNUMBERS 9000008 # ar
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 0 # particle/antiparticle distinction (0=own anti)
>
> #*********************************************************************
> # 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 to_full *
> #*********************************************************************
> #
> #*********************************************************************
> # 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=elastic photon of proton/ion beam *
> # +/-3=PDF of electron/positron beam *
> # +/-4=PDF of muon/antimuon beam *
> #*********************************************************************
> 0 = lpp1 ! beam 1 type
> 1 = lpp2 ! beam 2 type
> 5.0 = ebeam1 ! beam 1 total energy in GeV
> 37.52 = 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 alpha_s and its evol. *
> # pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
> # iww=Improved Weizsaecker-Williams Approx.[hep-ph/9310350] *
> # eva=Effective W/Z/A Approx. [2111.02442] *
> # edff=EDFF in gamma-UPC [eq.(11) in 2207.03012] *
> # chff=ChFF in gamma-UPC [eq.(13) in 2207.03012] *
> # none=No PDF, same as lhapdf with lppx=0 *
> #*********************************************************************
> none = pdlabel1 ! PDF type for beam #1
> nn23lo1 = pdlabel2 ! PDF type for beam #2
> 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> # To see heavy ion options: type "update ion_pdf"
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale1 ! if .true. use fixed fac scale for beam 1
> False = fixed_fac_scale2 ! if .true. use fixed fac scale for beam 2
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> 1.0 = scalefact ! scale factor for event-by-event scales
>
> 0 = ievo_eva ! scale evolution for EW pdfs (eva):
> ! 0 for evo by q^2; 1 for evo by pT^2
>
> #*********************************************************************
> # 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.00773)
> ! 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"
> #*********************************************************************
> # Customization (custom cuts/scale/bias/...) *
> # list of files containing fortran function that overwrite default *
> #*********************************************************************
> = custom_fcts ! List of files containing user hook function
> #*******************************
> # Parton level cuts definition *
> #*******************************
> 0.0 = dsqrt_shat ! minimal shat for full process
> #
> #
> #*********************************************************************
> # 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.0017 = pta ! minimum pt for the photons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptamax ! maximum pt for the photons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = 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 = etaa ! max rap for the photons
> 0.0 = etajmin ! min rap for the jets
> 0.0 = etaamin ! min rap for the photons
> {} = 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 *
> #*********************************************************************
> #*********************************************************************
> # 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 *
> #*********************************************************************
> 0.0 = ptheavy ! minimum pt for at least one heavy final state
> #*********************************************************************
> # 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 *
> #*********************************************************************
> True = 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
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,
I did find a file named P1_vmar_vmara and upon running the gensym, I got the following error message:
At line 125 of file symmetry.f (unit = 25)
Fortran runtime error: Cannot open file 'symfact_orig.dat': No such file or directory

Error termination. Backtrace:
#0 0x100d1ad2f
#1 0x100d1b8d7
#2 0x100d1c1a7
#3 0x100e0b0a3
#4 0x100e0b5bf
#5 0x10023c72f
#6 0x10023ed97

P.S. Yes, the model file does work for a vm e- scattering. But the scattering upon Ar seems to be having this issue. Also, I have only specified a vertex for the ar-ar-a and a Lorentz structure SSV for the interaction (including a form factor).

So the issue is (as reported) that the file symfact_orig.dat is not present.
This file is normally created when you create the code (so at the "output" command stage).
Can you re-generate the code to check if such file is indeed generated?

This is one of the last step of the code creation, but sounds weird that the code will not write it without complaining...

Cheers,

Olivier

Hi,
So, I tried running the code again for scattering on Ar again (Also for a process that works vm scattering on electron). In the first case the file symfact_orig.dat wasn't created, while in the second it was created. But when I tried running the gensym from SubProcesses, I got the same error message again (Cannot open file 'symfact_orig.dat': No such file or directory) for both the events.

Ok can you send me the model to this email adress <email address hidden> (with a link to this thread please such that I know what the issue is). I will try to duplicate the issue and to understand why such file is not created.

Cheers,

Olivier

So I have tested your process and it crash but not with the same error.

I have tested 3.5.1 (with python3.10) and the (recently release 3.5.6, with python3.12)
and in both case, the file symfact_orig.dat is generated. So I'm sorry to not reproduce your error and therefore not really able to help you (which version of python do you use?)

I have only try to run the code with 3.5.1 and it crash with the following error:

INFO: P1_vlar_vlara
Command "generate_events run_01" interrupted with error:
Exception : /Users/omattelaer/Documents/git_workspace/version3/PROC_Raylnulightscalar_Form_Factor_UFO_1/error
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/Users/omattelaer/Documents/git_workspace/version3/PROC_Raylnulightscalar_Form_Factor_UFO_1/run_01_tag_1_debug.log'.
Please attach this file to your report.
INFO:
quit
INFO:
MG5_aMC>^D

[version3]$ cat /Users/omattelaer/Documents/git_workspace/version3/PROC_Raylnulightscalar_Form_Factor_UFO_1/error
 Error: PDF not supported for pdg 9000008

Which makes sense, I guess you need to change lpp2 to 0 and not to 1.

Cheers,

Olivier

Hi,
I have python 3.12 installed in the system. Running the model on v 3.5.6 now simply prints out zero cross section with the following error message:
Survey return zero cross section.
   Typical reasons are the following:
   1) A massive s-channel particle has a width set to zero.
   2) The pdf are zero for at least one of the initial state particles
      or you are using maxjetflavor=4 for initial state b:s.
   3) The cuts are too strong.
   Please check/correct your param_card and/or your run_card.

I am also setting lpp2 to 0.

If I compute the amplitude for one PS point:

 -----------------------------------------------------------------------------
 n E px py pz m
 1 0.4992961E+03 0.0000000E+00 0.0000000E+00 0.4992961E+03 0.0000000E+00
 2 0.5007039E+03 0.0000000E+00 0.0000000E+00 -0.4992961E+03 0.3752000E+02
 3 0.4576928E+03 0.1691258E+03 0.3789202E+03 -0.1931286E+03 0.6031566E-05
 4 0.3652952E+03 -0.1829446E+02 -0.3470325E+03 0.1061441E+03 0.3752000E+02
 5 0.1770119E+03 -0.1508314E+03 -0.3188764E+02 0.8698449E+02 0.1651812E-05
 -----------------------------------------------------------------------------
 Matrix element = 0.0000000000000000 GeV^ -2
 -----------------------------------------------------------------------------

The problem seems that your matrix-element is numerically zero since (here in python but you will have the same issue in many language)

>>> math.exp(-5631*5)
0.0

So I guess that you need to factor out a factor to avoid those numerical factor.

Cheers,

Olivier

Hi,
I understand the issue with the code now.
Thank You very much for the assistance.
Appreciate the help!