error in writing subprocesses with MadSTR

Asked by Hesham El Faham on 2020-06-15


I am trying to generate (p p > t z w b) at NLO in QCD in the 5F scheme. I use MadSTR to account for interference with ttZ. Everything works fine in the computation part but then the program fails in writing the subprocesses giving the following error (attached along the rest of MG5_debug):
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm-no_b_mass
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
import model loop_sm-no_b_mass
define p = g u c d s u~ c~ d~ s~ b b~
define j = g u c d s u~ c~ d~ s~ b b~
define w = w+ w-
define tt = t t~
define bb = b b~
generate p p > tt w z bb [QCD]
Traceback (most recent call last):
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 278, in do_import
    self.cmd.do_import(self, *args, **opts)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 5450, in do_import
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 1661, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 1544, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/", line 239, in do_output
    self.export(nojpeg, main_file_name, group_processes=group_processes)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/", line 369, in export
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/", line 94, in generate_directories_fks
    calls = super(MadSTRExporter, self).generate_directories_fks(matrix_elements, *args)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/iolibs/", line 489, in generate_directories_fks
    self.write_real_matrix_elements(matrix_element, fortran_model)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/", line 203, in write_real_matrix_elements
    fksreal.matrix_element, os_me, suffix, fortran_model)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/", line 223, in write_os_wrapper
    permutation.append(ninitial + real_ids.index(os_id))
ValueError: -4 is not in list
                          MadGraph5_aMC@NLO Options
        complex_mass_scheme : False
    default_unset_couplings : 99
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0
               stdout_level : 20 (user set)

                         MadEvent Options
     automatic_html_opening : True
                    nb_core : None
        notification_center : True
                   run_mode : 2
                     Configuration Options
                        OLP : MadLoop
                    amcfast : amcfast-config
                   applgrid : applgrid-config
                auto_update : 7
         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 : None (user set)
               cpp_compiler : None
             crash_on_error : False
               delphes_path : ./Delphes
                 eps_viewer : None
        exrootanalysis_path : ./ExRootAnalysis
              f2py_compiler : None
                    fastjet : None (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/lhapdf6/bin/lhapdf-config (user set)
          madanalysis5_path : ./HEPTools/madanalysis5/madanalysis5
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : None (user set)
                      ninja : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/lib (user set)
        output_dependencies : external
            pythia-pgs_path : ./pythia-pgs
               pythia8_path : None (user set)
                    samurai : None
               syscalc_path : ./SysCalc
                    td_path : ./td
                text_editor : None
                thepeg_path : None (user set)
                    timeout : 60
                web_browser : None
May you please help with that?

Question information

English Edit question
MadGraph5_aMC@NLO Edit question
marco zaro Edit question
Solved by:
marco zaro
Last query:
Last reply:
Best marco zaro (marco-zaro) said : #1

Dear Hesham,
the process you are trying
generate p p > tt w z bb [QCD]
has resonant contributions which appear already at the Born. This is not supported by MadSTR.
If instead you wanted to generate the t w z final state, the correct syntax is

generate p p > tt w z [QCD]

Let me know

Best wishes,


Hesham El Faham (helfaham) said : #2

Thanks marco zaro, that solved my question.