Problem computing widths for vertices with form factors

Asked by Lara Mason on 2019-08-29

Hi experts,

I am attempting to calculate branching ratios in model where I have introduced some form factors in the vertices. I manage to generate some events with the model with no problem, except when I do
compute_widths __ or set the width of the particle to Auto.

When I do this, I get the error

Command "generate_events run_01" interrupted in sub-command:
"set max_npoint_for_channel 0" with error:
NameError : name 'Fgg' is not defined
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

(Fgg is the name of one of my form factors)

Is there a known issue with calculating decays widths for models including form factors, and is there a way to work around it?

I'm using Madgraph5 VERSION 2.6.6 .

Thanks very much in advance,
Lara

The contents of MG5_debug below:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.6 2018-06-28 *
#* *
#* 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 loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0

Traceback (most recent call last):
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2461, in do_generate_events
    switch_mode = self.ask_run_configuration(mode, args)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 6172, in ask_run_configuration
    self.ask_edit_cards(cards, plot=False, first_cmd=first_cmd)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 966, in ask_edit_cards
    banner=banner)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 1038, in ask_edit_card_static
    cards=cards, mode=mode, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1128, in ask
    fct=question_instance, fct_timeout=fct_timeout)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1728, in timed_input
    result = fct(question)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 2091, in __call__
    return self.cmdloop()
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 2259, in cmdloop
    super(SmartQuestion,self).cmdloop(intro)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 170, in cmdloop
    stop = self.postcmd(stop, line)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 5956, in postcmd
    self.do_update('dependent', timer=20)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 6002, in do_update
    self.do_compute_widths('')
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 6372, in do_compute_widths
    out = self.mother_interface.do_compute_widths(line)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 2288, in do_compute_widths
    out = cmd.exec_cmd(line, model=opts['model'])
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/master_interface.py", line 343, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madgraph_interface.py", line 8165, in do_compute_widths
    ), skip_2body=skip_2body, model=decaymodel)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/master_interface.py", line 331, in do_decay_diagram
    return self.cmd.do_decay_diagram(self, *args, **opts)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madgraph_interface.py", line 8439, in do_decay_diagram
    model.find_all_channels(2,generate_abstract=False)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/mg5decay/decay_objects.py", line 2845, in find_all_channels
    part.find_channels_nextlevel(self, min_br)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/mg5decay/decay_objects.py", line 914, in find_channels_nextlevel
    self['apx_decaywidth'] += temp_channel.get_apx_decaywidth(model)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/mg5decay/decay_objects.py", line 4424, in get_apx_decaywidth
    self.get_apx_psarea(model)/ \
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/mg5decay/decay_objects.py", line 4097, in get_apx_matrixelement_sq
    lor_value = eval(new_structure % q_dict_lor)
  File "<string>", line 1, in <module>
NameError: name 'Fgg' is not defined

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2019-08-29
Last reply:
2019-09-04
Lara Mason (laram) said : #1

Contd:

I found a related question
https://answers.launchpad.net/mg5amcnlo/+question/408658

and tried to implement a similar patch fix at line 4097, where the issue seems to start, but I then get other errors (see following). As I wasn't confident in my patch test, I wanted to confirm whether there is a known fix for this issue.

Thanks very much

--------------------------
Traceback (most recent call last):
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/master_interface.py", line 343, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madgraph_interface.py", line 8181, in do_compute_widths
    self.exec_cmd('output %s -f' % decay_dir,child=False)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/master_interface.py", line 296, in do_output
    self.cmd.do_output(self, line, *args, **opts)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madgraph_interface.py", line 7681, in do_output
    self.finalize(nojpeg, flaglist=flaglist)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/interface/madgraph_interface.py", line 7969, in finalize
    wanted_couplings)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/iolibs/export_v4.py", line 3553, in convert_model
    wanted_lorentz, wanted_couplings)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/madgraph/iolibs/export_v4.py", line 906, in convert_model
    aloha_model.compute_subset(wanted_lorentz)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/aloha/create_aloha.py", line 915, in compute_subset
    self.compute_aloha(builder, routines=routines)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/aloha/create_aloha.py", line 994, in compute_aloha
    wavefunction = builder.compute_routine(outgoing, tag)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/aloha/create_aloha.py", line 162, in compute_routine
    self.expr = self.compute_aloha_high_kernel(mode, factorize)
  File "/Users/laramason/Downloads/MG5_aMC_v2_6_5/aloha/create_aloha.py", line 289, in compute_aloha_high_kernel
    raise ALOHAERROR, 'unknow type in Lorentz Evaluation: %s ' % str(error)
ALOHAERROR: unknow type in Lorentz Evaluation: name 'Error' is not defined

Hi,

I guess I can not help on this without the model and without the command that you use (such that I can reproduce that issue).
Depending on "how" you implemented the form-factors, it can indeed not be compatible with that features.

Can you put the model online somewhere (or send it to me by email).

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Lara Mason for more information if necessary.

To post a message you must log in.