compute_widths fails with Higgs_Effective_Couplings_FormFact model

Asked by Emanuel Gouveia

Hello,

I am using the Higgs_Effective_Couplings_FormFact model (available here: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Models/ggHFullLoop) with version 2.5.4 of MadGraph5_aMC@NLO.

When I try to automatically obtain the width of a particle, either using compute_widths or setting the width to Auto in the param_card.dat file, I get this output:

Please note that the automatic computation of the width is
only valid in narrow-width approximation and at tree-level.
INFO: Get two body decay from FeynRules formula
partial width of particle t lower than QCD scale:0.00246492181726. Set it to zero. ([24, 3])
partial width of particle t lower than QCD scale:0.000108817188062. Set it to zero. ([24, 1])
Results written to /home/emanuel/MG5_aMC_v2_5_4/models/Higgs_Effective_Couplings_FormFact/param_card.dat
INFO: get decay diagram for t
Vertexlist of this model has not been searched.Automatically run the model.find_vertexlist()
Found 1 stable particles
Command "compute_widths t" interrupted with error:
NameError : name 'myGGH' 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.

This is the content of MG5_debug:

compute_widths t
Traceback (most recent call last):
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/extended_cmd.py", line 1415, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/extended_cmd.py", line 1369, in onecmd_orig
    return func(arg, **opt)
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/master_interface.py", line 334, in do_compute_widths
    return self.cmd.do_compute_widths(self, *args, **opts)
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/madgraph_interface.py", line 7925, in do_compute_widths
    ), skip_2body=skip_2body, model=decaymodel)
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/master_interface.py", line 322, in do_decay_diagram
    return self.cmd.do_decay_diagram(self, *args, **opts)
  File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/madgraph_interface.py", line 8185, in do_decay_diagram
    model.find_all_channels(2,generate_abstract=False)
  File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 2842, in find_all_channels
    part.find_channels_nextlevel(self, min_br)
  File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 914, in find_channels_nextlevel
    self['apx_decaywidth'] += temp_channel.get_apx_decaywidth(model)
  File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 4421, in get_apx_decaywidth
    self.get_apx_psarea(model)/ \
  File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 4094, in get_apx_matrixelement_sq
    lor_value = eval(new_structure % q_dict_lor)
  File "<string>", line 1, in <module>
NameError: name 'myGGH' is not defined
                          MadGraph5_aMC@NLO Options
                          ----------------
        complex_mass_scheme : False
                      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 : ./HEPTools/lib
               cpp_compiler : None
               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 : lhapdf-config
          madanalysis5_path : None (user set)
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : None (user set)
                      ninja : ./HEPTools/lib
        output_dependencies : external
                      pjfry : None (user set)
            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

Thank you very much!

Cheers,
Emanuel

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

This model is not a valid UFO model, it is indeed supported by MG5 for the main functionality but it is indeed incompatible with the auto width computation.

Cheers,

Olivier

> On 25 Apr 2017, at 21:17, Emanuel Gouveia <email address hidden> wrote:
>
> New question #628365 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/628365
>
> Hello,
>
> I am using the Higgs_Effective_Couplings_FormFact model (available here: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Models/ggHFullLoop) with version 2.5.4 of MadGraph5_aMC@NLO.
>
> When I try to automatically obtain the width of a particle, either using compute_widths or setting the width to Auto in the param_card.dat file, I get this output:
>
> Please note that the automatic computation of the width is
> only valid in narrow-width approximation and at tree-level.
> INFO: Get two body decay from FeynRules formula
> partial width of particle t lower than QCD scale:0.00246492181726. Set it to zero. ([24, 3])
> partial width of particle t lower than QCD scale:0.000108817188062. Set it to zero. ([24, 1])
> Results written to /home/emanuel/MG5_aMC_v2_5_4/models/Higgs_Effective_Couplings_FormFact/param_card.dat
> INFO: get decay diagram for t
> Vertexlist of this model has not been searched.Automatically run the model.find_vertexlist()
> Found 1 stable particles
> Command "compute_widths t" interrupted with error:
> NameError : name 'myGGH' 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.
>
> This is the content of MG5_debug:
>
> compute_widths t
> Traceback (most recent call last):
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/extended_cmd.py", line 1415, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/extended_cmd.py", line 1369, in onecmd_orig
> return func(arg, **opt)
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/master_interface.py", line 334, in do_compute_widths
> return self.cmd.do_compute_widths(self, *args, **opts)
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/madgraph_interface.py", line 7925, in do_compute_widths
> ), skip_2body=skip_2body, model=decaymodel)
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/master_interface.py", line 322, in do_decay_diagram
> return self.cmd.do_decay_diagram(self, *args, **opts)
> File "/home/emanuel/MG5_aMC_v2_5_4/madgraph/interface/madgraph_interface.py", line 8185, in do_decay_diagram
> model.find_all_channels(2,generate_abstract=False)
> File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 2842, in find_all_channels
> part.find_channels_nextlevel(self, min_br)
> File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 914, in find_channels_nextlevel
> self['apx_decaywidth'] += temp_channel.get_apx_decaywidth(model)
> File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 4421, in get_apx_decaywidth
> self.get_apx_psarea(model)/ \
> File "/home/emanuel/MG5_aMC_v2_5_4/mg5decay/decay_objects.py", line 4094, in get_apx_matrixelement_sq
> lor_value = eval(new_structure % q_dict_lor)
> File "<string>", line 1, in <module>
> NameError: name 'myGGH' is not defined
> MadGraph5_aMC@NLO Options
> ----------------
> complex_mass_scheme : False
> 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 : ./HEPTools/lib
> cpp_compiler : None
> 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 : lhapdf-config
> madanalysis5_path : None (user set)
> madanalysis_path : ./MadAnalysis
> mg5amc_py8_interface_path : None (user set)
> ninja : ./HEPTools/lib
> output_dependencies : external
> pjfry : None (user set)
> 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
>
> Thank you very much!
>
> Cheers,
> Emanuel
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Emanuel Gouveia (eddygouveia91) said :
#2

Thanks Olivier Mattelaer, that solved my question.