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