Problem with Higgs characterization model, gluon fusion

Asked by Laura Buonincontri

Dear experts,

I am writing since I am trying to run the Higgs characterization model to calculate the p p > h, p p > h j at NLO.
I followed the README list of recommendations and looked at the launchpad, and everybody seems to have no issues in running the process by doing:

 >import model HC_NLO_X0_UFO-heft
>generate p p > x0 / t [QCD] @0
>add process p p > x0 j / t [QCD] @1

However, I am not able to do it. When I type the first generate:

MG5_aMC>generate p p > x0 / t [QCD] @0
Setting the born squared orders automatically in the process definition to QED^2<=2 QCD=0 QNP=0.
If this is not what you need, please regenerate with the correct orders.
---------------------------------------------------------------------------
 No Born diagrams found. Now switching to the loop-induced mode.
 Please also cite ref. 'arXiv:1507.00020' when using results from this mode.
---------------------------------------------------------------------------
Command "generate p p > x0 / t [QCD] @0" interrupted with error:
Exception :
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

Here I attached the MG5_debug report here:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.0 2022-05-06 *
#* *
#* 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 group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model sm
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~
import model HC_NLO_X0_UFO-heft
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > x0 / t [QCD] @0
Traceback (most recent call last):
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/fks/fks_base.py", line 173, in __init__
    super(FKSMultiProcess, self).__init__(procdef, **options)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/core/diagram_generation.py", line 1611, in __init__
    self.get('amplitudes')
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/core/diagram_generation.py", line 1650, in get
    self.generate_multi_amplitudes(process_def,
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/core/diagram_generation.py", line 1903, in generate_multi_amplitudes
    raise NoDiagramException("No amplitudes generated from process %s. Please enter a valid process" % \
madgraph.core.diagram_generation.NoDiagramException: No amplitudes generated from process Process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > x0 / t QNP=1 [ all = QED QNP QCD ] QED^2=2 QCD^2=2 QNP^2=0. Please enter a valid process

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/master_interface.py", line 229, in do_add
    return self.cmd.do_add(self, line, *args, **opts)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/amcatnlo_interface.py", line 667, in do_add
    fksproc =fks_base.FKSMultiProcess(myprocdef,fks_options)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/fks/fks_base.py", line 177, in __init__
    raise NoBornException("Born diagrams could not be generated for the "+\
madgraph.fks.fks_base.NoBornException: Born diagrams could not be generated for the process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > x0 / t QNP=1 [ all = QED QNP QCD ] QED^2=2 QCD^2=2 QNP^2=0. Notice that aMC@NLO does not handle loop-induced processes yet, but you can still use MadLoop if you want to only generate them. For this, use the 'virt=' mode, without multiparticle labels.

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/master_interface.py", line 278, in do_generate
    return self.cmd.do_generate(self, line, *args, **opts)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/madgraph_interface.py", line 4704, in do_generate
    self.do_add(" ".join(args))
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/master_interface.py", line 236, in do_add
    return self.cmd.create_loop_induced(self, line, *args, **opts)
  File "/home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/madgraph/interface/madgraph_interface.py", line 5282, in create_loop_induced
    raise Exception
Exception
                          MadGraph5_aMC@NLO Options
                          ----------------
                auto_update : 7
        complex_mass_scheme : False
    default_unset_couplings : 99
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
  include_lepton_initiated_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0
          max_t_for_channel : 99
        nlo_mixed_expansion : True
               stdout_level : 20 (user set)
         zerowidth_tchannel : True

                         MadEvent Options
                          ----------------
     automatic_html_opening : True
                    nb_core : None
        notification_center : True
                   run_mode : 2

                      Configuration Options
                      ---------------------
                        OLP : MadLoop
   acknowledged_v3.1_syntax : False
         auto_convert_model : False
         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 : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/lib (user set)
                contur_path : ./HEPTools/contur
               cpp_compiler : None
             crash_on_error : False
               delphes_path : ./Delphes
                 eps_viewer : None
        exrootanalysis_path : ./ExRootAnalysis
              f2py_compiler : None
          f2py_compiler_py2 : None
          f2py_compiler_py3 : None
                    fastjet : None (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
                 lhapdf_py2 : None
                 lhapdf_py3 : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
          madanalysis5_path : None (user set)
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/MG5aMC_PY8_interface (user set)
                      ninja : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/lib (user set)
        output_dependencies : external
                   pineappl : pineappl
            pythia-pgs_path : ./pythia-pgs
               pythia8_path : /home/user/madgraph/MG5_aMC_v3.4.0/MG5_aMC_v3_4_0/HEPTools/pythia8 (user set)
                 rivet_path : ./HEPTools/rivet
                    samurai : None
               syscalc_path : ./SysCalc
                    td_path : ./td
                text_editor : None
                thepeg_path : None (user set)
                    timeout : 60
                web_browser : None
                  yoda_path : ./HEPTools/yoda

Can you see if I am missing something?

Thank you,

Laura

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
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

I guess that those instructions where for version2 of the code.
Therefore I would suggest to use our Long Term Stable which is still on version 2.

For version 3 to work you likely need to specify exactly which coupling order you want for this tricky case:
I do not know that model, but I would guess that you would need:
generate p p > x0 / t QCD<=2 QNP<=1 [QCD] @0

Now, be carefull that some of those models do have issue with scale variation in version 3 since some of them are wrong on the power of alphas in the effective coupling. (in version 2 we where fitting the power of alphas which was limiting us in many way but was avoiding that particular issue)

Cheers,

Olivier

Revision history for this message
Laura Buonincontri (laurabuonincontri) said :
#2

Hi Olivier,

yes that is the problem for the generation of the process, thank you very much!
I downloaded the LTS version 2 and now it works.
I typed the commands written in the README file.

 import model HC_NLO_X0-heft
 generate p p > x0 / t [QCD]

I had a problem after the launch command.

INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
Error detected in "launch auto "
write debug file /home/user/madgraph/MG5aMC_LTS_v2.9.16/MG5_aMC_v2_9_16/prova_ggf/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : Poles do not cancel, run cannot continue

But I solved it by putting to -1 #IRPoleCheckThreshold and #PrecisionVirtualAtRunTime, like suggested in the FAQ (https://answers.launchpad.net/mg5amcnlo/+faq/2720).

I still have one issue when I run madspin on the x0. I would like to make it decay into b b~.
When Madspin is run it shows these available decays:

INFO: decay channels for x0 : ( width = 0.0004512381 GeV )
INFO: BR d1 d2
INFO: 5.734223e-01 ta+ ta-
INFO: 4.058390e-01 g g
INFO: 1.397089e-02 a a
INFO: 6.767824e-03 a z
INFO:

So obviously, it crashes by saying:

Command "launch auto " interrupted with error:
NoDiagramException : No amplitudes generated from process Process: x0 > b b~ QCD=99 QNP=1. Please enter a valid process

How can I make x0 decaying into b b~? Do I have to specify something in the generation, like QNP=1? I saw that in the coupling.py and parameters.py the yukhawa coupling is coded.

Thank you,

Laura Buonincontri

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#3

I guess that you might need to change the model for the decay within madspin.
since it looks like your model (or likely the heft restriction of the model) does not have a coupling from x0 to b b~ (or it is set to zero).

So you either need to change the model, the restriction or just the param_card to let madspin use a benchmark where x0 can decay to b b~

Cheers,

Olivier

Revision history for this message
Laura Buonincontri (laurabuonincontri) said :
#4

Thanks Olivier Mattelaer, that solved my question.