error in writing subprocesses with MadSTR

Asked by Hesham El Faham on 2020-06-15

Hi,

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]
output
Traceback (most recent call last):
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/master_interface.py", 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/madgraph_interface.py", line 5450, in do_import
    self.import_command_file(args[1])
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", 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/extended_cmd.py", 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/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/madstr_interface.py", 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/madstr_interface.py", line 369, in export
    path,self.options['OLP'])
  File "/home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/PLUGIN/MadSTR/madstr_exporter.py", 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/export_fks.py", 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/madstr_exporter.py", 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/madstr_exporter.py", 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?
Thanks,
Hesham

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
marco zaro
Solved:
2020-06-16
Last query:
2020-06-16
Last reply:
2020-06-16
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,

Marco

Hesham El Faham (helfaham) said : #2

Thanks marco zaro, that solved my question.