FCNC key error when using dim6op and MadSpin

Asked by Hesham El Faham on 2021-02-18

Hi,
I am trying to generate the following simple EFT process using 3.0.3-betternlops (I also tried using MG5_aMC_v2_7_2):
-->
import model dim6top_LO_UFO_each_coupling_order
generate g g > t t~ h FCNC=0 DIM6=1
<--
Everything works fine until I turn MadSpin on, that's when I get the following error (log file pasted here):
-->
generate_events
Traceback (most recent call last):
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1518, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1467, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/madevent_interface.py", line 2468, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/common_run_interface.py", line 7015, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/madevent_interface.py", line 2567, in run_generate_events
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1547, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1467, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/common_run_interface.py", line 3708, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1664, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/interface_madspin.py", line 634, in do_launch
    self.options)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2046, in __init__
    self.generate_all_matrix_element()
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 1934, in deco_f_set
    out = f(*args, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2663, in generate_all_matrix_element
    self.get_branching_ratio()
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2921, in get_branching_ratio
    width.extract_br(self.list_branches, self.mgcmd)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1351, in extract_br
    self.launch_width_evaluation(to_decay, mgcmd)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1618, in launch_width_evaluation
    self.compute_widths(model, argument)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1656, in compute_widths
    cmd.exec_cmd(line)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 345, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8419, in do_compute_widths
    self.exec_cmd('output %s -f' % decay_dir,child=False)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
    return func(arg, **opt)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 298, in do_output
    self.cmd.do_output(self, line, *args, **opts)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 7916, in do_output
    self.export(nojpeg, main_file_name, group_processes, args)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 312, in export
    return self.cmd.export(self, *args, **opts)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8042, in export
    ndiags, cpu_time = generate_matrix_elements(self,group_processes)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8024, in generate_matrix_elements
    matrix_element_opts=mode)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5626, in __init__
    compute_loop_nc = compute_loop_nc))
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5902, in generate_matrix_elements
    compute_loop_nc=compute_loop_nc)
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5696, in process_color
    new_amp = matrix_element.get_base_amplitude()
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 4634, in get_base_amplitude
    diag.calculate_orders(self.get('processes')[0].get('model'))
  File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/base_objects.py", line 2573, in calculate_orders
    (c,n) in couplings.items()])
KeyError: 'FCNC'
                             Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : None
             cluster_memory : None
               cluster_size : 100
              cluster_queue : None
                    nb_core : 8 (user set)
               cluster_time : None
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
               web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/lhapdf6/bin/lhapdf-config (user set)
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : None (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.0.3 2020-06-21 *
#* BZR 3.0.3-betternlops 979 *
#* *
#* 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 default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set include_lepton_initiated_processes False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model dim6top_LO_UFO_each_coupling_order
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~
generate g g > t t~ h FCNC=0 DIM6=1
output test_ttH_EFT
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DIM6
###################################
BLOCK DIM6 #
      1 1.000000e+03 # lambda
      2 0.000000e+00 # ctp
      3 0.000000e+00 # ctpi
      4 0.000000e+00 # cpqm
      5 0.000000e+00 # cpq3
      6 1.000000e+00 # cpt
      7 0.000000e+00 # cpb
      8 0.000000e+00 # cptb
      9 0.000000e+00 # cptbi
      10 0.000000e+00 # ctw
      11 0.000000e+00 # ctz
      12 0.000000e+00 # ctwi
      13 0.000000e+00 # ctzi
      14 0.000000e+00 # cbw
      15 0.000000e+00 # cbwi
      16 0.000000e+00 # ctg
      17 0.000000e+00 # ctgi
      18 0.000000e+00 # cqlm1
      19 0.000000e+00 # cqlm2
      20 0.000000e+00 # cqlm3
      21 0.000000e+00 # cql31
      22 0.000000e+00 # cql32
      23 0.000000e+00 # cql33
      24 0.000000e+00 # cqe1
      25 0.000000e+00 # cqe2
      26 0.000000e+00 # cqe3
      27 0.000000e+00 # ctl1
      28 0.000000e+00 # ctl2
      29 0.000000e+00 # ctl3
      30 0.000000e+00 # cte1
      31 0.000000e+00 # cte2
      32 0.000000e+00 # cte3
      33 0.000000e+00 # ctls1
      34 0.000000e+00 # ctlsi1
      35 0.000000e+00 # ctls2
      36 0.000000e+00 # ctlsi2
      37 0.000000e+00 # ctls3
      38 0.000000e+00 # ctlsi3
      39 0.000000e+00 # ctlt1
      40 0.000000e+00 # ctlti1
      41 0.000000e+00 # ctlt2
      42 0.000000e+00 # ctlti2
      43 0.000000e+00 # ctlt3
      44 0.000000e+00 # ctlti3
      45 0.000000e+00 # cbls1
      46 0.000000e+00 # cblsi1
      47 0.000000e+00 # cbls2
      48 0.000000e+00 # cblsi2
..
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.466900e+00
# BR NDA ID1 ID2 ...
   1.000000e+00 2 5 24 # 1.4669
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495000e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 6.526578e-03
# BR NDA ID1 ID2 ...
   8.259611e-01 2 -5 5 # 0.00539069954412
   3.964405e-02 2 -15 15 # 0.000258739984561
   2.110294e-02 3 -4 3 24 # 0.000137729983939
   2.108915e-02 3 -24 -1 2 # 0.000137639982429
   2.092674e-02 3 -2 1 24 # 0.000136580000896
   2.079497e-02 3 -24 -3 4 # 0.000135719993713
   7.045269e-03 3 -12 11 24 # 4.59814976595e-05
   7.045269e-03 3 -14 13 24 # 4.59814976595e-05
   7.020447e-03 3 -24 -11 12 # 4.58194949404e-05
   7.020447e-03 3 -24 -13 14 # 4.58194949404e-05
   6.982066e-03 3 -24 -15 16 # 4.55689983501e-05
   6.971800e-03 3 -16 15 24 # 4.55019965004e-05
   1.925281e-03 3 -1 1 23 # 1.25654966184e-05
   1.925281e-03 3 -3 3 23 # 1.25654966184e-05
   1.545067e-03 3 -5 5 23 # 1.00840002907e-05
   1.516798e-03 3 -4 4 23 # 9.89950045724e-06
   1.483304e-03 3 -2 2 23 # 9.68089925371e-06
<--
I tried running with the default numbers for the decay widths and with the option 'auto' but both lead to the same error.
May you please help with this?
Best,
Hesham

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2021-02-19
Last query:
2021-02-19
Last reply:
2021-02-19

Hi,

This is a bug in the model definition, they do define coupling to be of type FCNC.
But they do not define the property of what FCNC is.

The model authors need to modify their UFO model to follow the UFO convention.
(I have put them in copy of this message)

Cheers,

Olivier

>
> New question #695632 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/695632 <https://answers.launchpad.net/mg5amcnlo/+question/695632>
>
> Hi,
> I am trying to generate the following simple EFT process using 3.0.3-betternlops (I also tried using MG5_aMC_v2_7_2):
> -->
> import model dim6top_LO_UFO_each_coupling_order
> generate g g > t t~ h FCNC=0 DIM6=1
> <--
> Everything works fine until I turn MadSpin on, that's when I get the following error (log file pasted here):
> -->
> generate_events
> Traceback (most recent call last):
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1518, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1467, in onecmd_orig
> return func(arg, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/madevent_interface.py", line 2468, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/common_run_interface.py", line 7015, in new_fct
> original_fct(obj, *args, **opts)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/madevent_interface.py", line 2567, in run_generate_events
> self.exec_cmd('decay_events -from_cards', postcmd=False)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1547, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/extended_cmd.py", line 1467, in onecmd_orig
> return func(arg, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/ttH_EFT_reparam/test_ttH_EFT/bin/internal/common_run_interface.py", line 3708, in do_decay_events
> madspin_cmd.import_command_file(path)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1664, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
> return func(arg, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 100, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/interface_madspin.py", line 634, in do_launch
> self.options)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2046, in __init__
> self.generate_all_matrix_element()
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 100, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/various/misc.py", line 1934, in deco_f_set
> out = f(*args, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2663, in generate_all_matrix_element
> self.get_branching_ratio()
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 2921, in get_branching_ratio
> width.extract_br(self.list_branches, self.mgcmd)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1351, in extract_br
> self.launch_width_evaluation(to_decay, mgcmd)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1618, in launch_width_evaluation
> self.compute_widths(model, argument)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/MadSpin/decay.py", line 1656, in compute_widths
> cmd.exec_cmd(line)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
> return func(arg, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 345, in do_compute_widths
> return self.cmd.do_compute_widths(self, *args, **opts)
> return self.cmd.do_compute_widths(self, *args, **opts)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8419, in do_compute_widths
> self.exec_cmd('output %s -f' % decay_dir,child=False)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1547, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/extended_cmd.py", line 1467, in onecmd_orig
> return func(arg, **opt)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 298, in do_output
> self.cmd.do_output(self, line, *args, **opts)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 7916, in do_output
> self.export(nojpeg, main_file_name, group_processes, args)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/master_interface.py", line 312, in export
> return self.cmd.export(self, *args, **opts)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8042, in export
> ndiags, cpu_time = generate_matrix_elements(self,group_processes)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/interface/madgraph_interface.py", line 8024, in generate_matrix_elements
> matrix_element_opts=mode)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5626, in __init__
> compute_loop_nc = compute_loop_nc))
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5902, in generate_matrix_elements
> compute_loop_nc=compute_loop_nc)
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 5696, in process_color
> new_amp = matrix_element.get_base_amplitude()
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/helas_objects.py", line 4634, in get_base_amplitude
> diag.calculate_orders(self.get('processes')[0].get('model'))
> File "/home/elfaham/Downloads/tzw/working_dir/3.0.3-betternlops/madgraph/core/base_objects.py", line 2573, in calculate_orders
> (c,n) in couplings.items()])
> KeyError: 'FCNC'
> Run Options
> -----------
> stdout_level : None
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None
> cluster_size : 100
> cluster_queue : None
> nb_core : 8 (user set)
> cluster_time : None
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> pythia8_path : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/pythia8 (user set)
> hwpp_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> thepeg_path : None (user set)
> cluster_type : condor
> madanalysis5_path : None (user set)
> cluster_nb_retry : 1
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : /home/elfaham/Downloads/tzw/working_dir/MG5_aMC_v2_7_2/HEPTools/lhapdf6/bin/lhapdf-config (user set)
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : None (user set)
> timeout : 60
> cpp_compiler : None
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.0.3 2020-06-21 *
> #* BZR 3.0.3-betternlops 979 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph <https://server06.fynu.ucl.ac.be/projects/madgraph> *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set default_unset_couplings 99
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set include_lepton_initiated_processes False
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model dim6top_LO_UFO_each_coupling_order
> 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~
> generate g g > t t~ h FCNC=0 DIM6=1
> output test_ttH_EFT
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR DIM6
> ###################################
> BLOCK DIM6 #
> 1 1.000000e+03 # lambda
> 2 0.000000e+00 # ctp
> 3 0.000000e+00 # ctpi
> 4 0.000000e+00 # cpqm
> 5 0.000000e+00 # cpq3
> 6 1.000000e+00 # cpt
> 7 0.000000e+00 # cpb
> 8 0.000000e+00 # cptb
> 9 0.000000e+00 # cptbi
> 10 0.000000e+00 # ctw
> 11 0.000000e+00 # ctz
> 12 0.000000e+00 # ctwi
> 13 0.000000e+00 # ctzi
> 14 0.000000e+00 # cbw
> 15 0.000000e+00 # cbwi
> 16 0.000000e+00 # ctg
> 17 0.000000e+00 # ctgi
> 18 0.000000e+00 # cqlm1
> 19 0.000000e+00 # cqlm2
> 20 0.000000e+00 # cqlm3
> 21 0.000000e+00 # cql31
> 22 0.000000e+00 # cql32
> 23 0.000000e+00 # cql33
> 24 0.000000e+00 # cqe1
> 25 0.000000e+00 # cqe2
> 26 0.000000e+00 # cqe3
> 27 0.000000e+00 # ctl1
> 28 0.000000e+00 # ctl2
> 29 0.000000e+00 # ctl3
> 30 0.000000e+00 # cte1
> 31 0.000000e+00 # cte2
> 32 0.000000e+00 # cte3
> 33 0.000000e+00 # ctls1
> 34 0.000000e+00 # ctlsi1
> 35 0.000000e+00 # ctls2
> 36 0.000000e+00 # ctlsi2
> 37 0.000000e+00 # ctls3
> 38 0.000000e+00 # ctlsi3
> 39 0.000000e+00 # ctlt1
> 40 0.000000e+00 # ctlti1
> 41 0.000000e+00 # ctlt2
> 42 0.000000e+00 # ctlti2
> 43 0.000000e+00 # ctlt3
> 44 0.000000e+00 # ctlti3
> 45 0.000000e+00 # cbls1
> 46 0.000000e+00 # cblsi1
> 47 0.000000e+00 # cbls2
> 48 0.000000e+00 # cblsi2
> ..
> #*************************
> # Decay widths *
> #*************************
> #
> # PDG Width
> DECAY 1 0.000000e+00
> #
> # PDG Width
> DECAY 2 0.000000e+00
> #
> # PDG Width
> DECAY 3 0.000000e+00
> #
> # PDG Width
> DECAY 4 0.000000e+00
> #
> # PDG Width
> DECAY 5 0.000000e+00
> #
> # PDG Width
> DECAY 6 1.466900e+00
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 5 24 # 1.4669
> #
> # PDG Width
> DECAY 11 0.000000e+00
> #
> # PDG Width
> DECAY 12 0.000000e+00
> #
> # PDG Width
> DECAY 13 0.000000e+00
> # PDG Width
> DECAY 14 0.000000e+00
> #
> # PDG Width
> DECAY 15 0.000000e+00
> #
> # PDG Width
> DECAY 16 0.000000e+00
> #
> # PDG Width
> DECAY 21 0.000000e+00
> #
> # PDG Width
> DECAY 22 0.000000e+00
> #
> # PDG Width
> DECAY 23 2.495000e+00
> #
> # PDG Width
> DECAY 24 2.085000e+00
> #
> # PDG Width
> DECAY 25 6.526578e-03
> # BR NDA ID1 ID2 ...
> 8.259611e-01 2 -5 5 # 0.00539069954412
> 3.964405e-02 2 -15 15 # 0.000258739984561
> 2.110294e-02 3 -4 3 24 # 0.000137729983939
> 2.108915e-02 3 -24 -1 2 # 0.000137639982429
> 2.092674e-02 3 -2 1 24 # 0.000136580000896
> 2.079497e-02 3 -24 -3 4 # 0.000135719993713
> 7.045269e-03 3 -12 11 24 # 4.59814976595e-05
> 7.045269e-03 3 -14 13 24 # 4.59814976595e-05
> 7.020447e-03 3 -24 -11 12 # 4.58194949404e-05
> 7.020447e-03 3 -24 -13 14 # 4.58194949404e-05
> 6.982066e-03 3 -24 -15 16 # 4.55689983501e-05
> 6.971800e-03 3 -16 15 24 # 4.55019965004e-05
> 1.925281e-03 3 -1 1 23 # 1.25654966184e-05
> 1.925281e-03 3 -3 3 23 # 1.25654966184e-05
> 1.545067e-03 3 -5 5 23 # 1.00840002907e-05
> 1.516798e-03 3 -4 4 23 # 9.89950045724e-06
> 1.483304e-03 3 -2 2 23 # 9.68089925371e-06
> <--
> I tried running with the default numbers for the decay widths and with the option 'auto' but both lead to the same error.
> May you please help with this?
> Best,
> Hesham
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hesham El Faham (helfaham) said : #2

Thanks, Olivier. I will wait for their response then.

Best,
Hesham

Hi,

Loooks like their answer did not went trough... here it is:

Thanks Olivier for the notice!

Hesham, please add:
---
FCNC = CouplingOrder(name = 'FCNC',
                     expansion_order = 99,
                     hierarchy = 1)
---
in coupling_orders.py of the model directory.

I've just made the addition in the UFO available from the model webpage
(and made it python3 compatible using MG's conversion tool).

Cheers,
Gauthier

Hesham El Faham (helfaham) said : #4

Thanks Olivier Mattelaer, that solved my question.