MadGraph sometimes crashes during automatic width computation

Asked by Leonard Schwarze

Dear MadGraph experts,

I am studying an extension of the SM introducing an additional scalar sector. The automatic width computation sometimes causes MadGraph to crash when running compute_widths. The crash only occurs when calculating three-body decays, although there are certain particle masses for which the three-body decays are calculated as expected.
The error message is the following:

Command "generate_events run_01" interrupted in sub-command:
"set nlo_mixed_expansion True" with error:
ValueError : -24 is not in list

I am appending the MG5_debug file at the end of this post.
I have checked that this error occurs in MadGraph versions 3.3.1, 3.3.0 and 2.9.7.

I would be thankful for any help.
Leonard

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.3.0 2021-11-12 *
#* *
#* 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
Traceback (most recent call last):
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 5838, in generate_matrix_elements
    me_index = amplitude_tags.index(amplitude_tag)
ValueError: [False, 0, False, 2, 0, [array('H', [2, 0, 1, 0, 1, 3]), array('H', [2, 3, 0, 1, 0, 1])], [(((((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 2),), (((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 4),)), (((1, 1, False, 'mdl_MS12', 'mdl_WS12', 6100003), (((0, 0), 'GC_36'), ((0, 1), 'GC_37'), '1 ', 'VVS1', 'VVS2', 'QED')),); 1, (((1, 0, 1, None, False, False, 'mdl_MS11', 'mdl_WS11', 1), 1),), (((0, 0, 3, None, False, False, 'mdl_MW', 'mdl_WW', 1), 3),)), (((((0, 0), 'GC_45'), ((0, 1), 'GC_48'), '1 ', 'VSS1', 'VSS2', 'QED'),),); 3, (((((1, 0, 1, None, False, False, 'mdl_MS11', 'mdl_WS11', 1), 1),), (((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 2),)), (((3, 1, True, 'mdl_MZ', 'mdl_WZ', 23), (((0, 0), 'GC_28'), '1 ', 'VVS1', 'QED')),); 1, (((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 4),), (((0, 0, 3, None, False, False, 'mdl_MW', 'mdl_WW', 1), 3),)), (((((0, 6), 'GC_60'), ((0, 7), 'GC_59'), '1 ', 'VVV1', 'VVV2', 'VVV3', 'VVV4', 'VVV5', 'VVV6', 'VVV7', 'VVV8', 'QED'),),); 3, (((((1, 0, 1, None, False, False, 'mdl_MS11', 'mdl_WS11', 1), 1),), (((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 3),)), (((3, 1, True, 'mdl_MZ', 'mdl_WZ', 23), (((0, 0), 'GC_28'), '1 ', 'VVS1', 'QED')),); 1, (((0, 0, 3, None, True, False, 'mdl_MW', 'mdl_WW', 1), 2),), (((0, 0, 3, None, False, False, 'mdl_MW', 'mdl_WW', 1), 4),)), (((((0, 6), 'GC_60'), ((0, 7), 'GC_59'), '1 ', 'VVV1', 'VVV2', 'VVV3', 'VVV4', 'VVV5', 'VVV6', 'VVV7', 'VVV8', 'QED'),),); 3]] is not in list

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madevent_interface.py", line 2354, in do_generate_events
    switch_mode = self.ask_run_configuration(mode, args)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madevent_interface.py", line 6164, in ask_run_configuration
    self.ask_edit_cards(cards, plot=False, first_cmd=first_cmd)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 1037, in ask_edit_cards
    self.ask_edit_card_static(cards, mode, plot, self.options['timeout'],
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 1109, in ask_edit_card_static
    out = ask(question, '0', possible_answer, timeout=int(1.5*timeout),
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1135, in ask
    question_instance.postcmd(answer, line)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 6294, in postcmd
    self.do_update('dependent', timer=20)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 6340, in do_update
    self.do_compute_widths('')
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 6711, in do_compute_widths
    out = self.mother_interface.do_compute_widths(line)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/common_run_interface.py", line 2380, in do_compute_widths
    out = cmd.exec_cmd(line, model=opts['model'])
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/master_interface.py", line 346, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madgraph_interface.py", line 8700, in do_compute_widths
    self.exec_cmd('output madevent %s -f' % decay_dir,child=False)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/master_interface.py", line 299, in do_output
    self.cmd.do_output(self, line, *args, **opts)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madgraph_interface.py", line 8188, in do_output
    self.export(nojpeg, main_file_name, group_processes, args)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/master_interface.py", line 313, in export
    return self.cmd.export(self, *args, **opts)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madgraph_interface.py", line 8318, in export
    ndiags, cpu_time = generate_matrix_elements(self,group_processes)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/interface/madgraph_interface.py", line 8298, in generate_matrix_elements
    self._curr_matrix_elements = HelasMultiProcessClass(
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 5606, in __init__
    self.set('matrix_elements', self.generate_matrix_elements(argument,
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 5847, in generate_matrix_elements
    matrix_element_list = [cls.matrix_element_class(amplitude,
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 3447, in __init__
    self.generate_helas_diagrams(amplitude, optimization, decay_ids)
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 3638, in generate_helas_diagrams
    wf, wf_number = wf.check_and_fix_fermion_flow(\
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 1192, in check_and_fix_fermion_flow
    self.set('mothers', self.get('mothers').sort_by_pdg_codes(\
  File "/home/leonard/Uni/programs/MG5_aMC_v3_3_0/madgraph/core/helas_objects.py", line 2395, in sort_by_pdg_codes
    index = mother_codes.index(code)
ValueError: -24 is not in list
                          MadGraph5_aMC@NLO Options
                          ----------------
        complex_mass_scheme : False
    default_unset_couplings : 99 (user set)
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
  include_lepton_initiated_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0 (user set)
          max_t_for_channel : 99 (user set)
        nlo_mixed_expansion : False (user set)
               stdout_level : 20 (user set)
         zerowidth_tchannel : True

                         MadEvent Options
                          ----------------
     automatic_html_opening : False (user set)
                    nb_core : 8 (user set)
        notification_center : True
                   run_mode : 2

                      Configuration Options
                      ---------------------
                        OLP : MadLoop
   acknowledged_v3.1_syntax : False
         auto_convert_model : True (user set)
                auto_update : 7 (user set)
         cluster_local_path : None
           cluster_nb_retry : 1
              cluster_queue : None (user set)
         cluster_retry_wait : 300
               cluster_size : 100
      cluster_status_update : (600, 30)
          cluster_temp_path : None
               cluster_type : condor
                    collier : ./HEPTools/lib
               cpp_compiler : None
             crash_on_error : False
               delphes_path : None (user set)
                 eps_viewer : None
        exrootanalysis_path : None (user set)
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                    fastjet : None (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : /home/leonard/Uni/programs/MG5_aMC_v3_3_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /home/leonard/Uni/programs/MG5_aMC_v3_3_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
          madanalysis5_path : None (user set)
           madanalysis_path : None (user set)
  mg5amc_py8_interface_path : /home/leonard/Uni/programs/MG5_aMC_v3_3_0/HEPTools/MG5aMC_PY8_interface (user set)
                      ninja : ./HEPTools/lib
        output_dependencies : external
                   pineappl : pineappl
            pythia-pgs_path : None (user set)
               pythia8_path : /home/leonard/Uni/programs/MG5_aMC_v3_3_0/HEPTools/pythia8 (user set)
                    samurai : None
               syscalc_path : None (user set)
                    td_path : None (user set)
                text_editor : None
                thepeg_path : None (user set)
                    timeout : 60
                web_browser : None

Question information

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

see bug report for the patch

Can you help with this problem?

Provide an answer of your own, or ask Leonard Schwarze for more information if necessary.

To post a message you must log in.