Not able to run NLO process in models in Madgraph 3.4.2

Bug #2011588 reported by Esteban Chalbaud
6
This bug affects 1 person
Affects Status Importance Assigned to Milestone
MadGraph5_aMC@NLO
New
Undecided
Kentarou Mawatari

Bug Description

Dear all I am having problems to run the following instructions in MadGraph:
import model DMsimp_s_spin1
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > xd xd~ j [QCD]

After I get:

INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS 2*DMV+QCD+2*QED
Setting the born squared orders automatically in the process definition to QED^2<=4 QCD^2<=2 DMV=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 > xd xd~ j [QCD]" 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.

The report says:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.2 2023-01-20 *
#* *
#* 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 DMsimp_s_spin1
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > xd xd~ j [QCD]
Traceback (most recent call last):
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/fks/fks_base.py", line 173, in __init__
    super(FKSMultiProcess, self).__init__(procdef, **options)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/core/diagram_generation.py", line 1611, in __init__
    self.get('amplitudes')
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/core/diagram_generation.py", line 1655, in get
    diagram_filter=self['diagram_filter']))
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/core/diagram_generation.py", line 1904, in generate_multi_amplitudes
    process_definition.nice_string())
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~ > xd xd~ g/u/c/d/s/b/u~/c~/d~/s~/b~ DMV=2 [ all = QED DMV QCD ] QED^2=4 QCD^2=4 DMV^2=0. Please enter a valid process

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/master_interface.py", line 229, in do_add
    return self.cmd.do_add(self, line, *args, **opts)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/amcatnlo_interface.py", line 667, in do_add
    fksproc =fks_base.FKSMultiProcess(myprocdef,fks_options)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/fks/fks_base.py", line 182, in __init__
    " For this, use the 'virt=' mode, without multiparticle labels.")
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~ > xd xd~ g/u/c/d/s/b/u~/c~/d~/s~/b~ DMV=2 [ all = QED DMV QCD ] QED^2=4 QCD^2=4 DMV^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 "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/master_interface.py", line 278, in do_generate
    return self.cmd.do_generate(self, line, *args, **opts)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/madgraph_interface.py", line 4712, in do_generate
    self.do_add(" ".join(args))
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/master_interface.py", line 236, in do_add
    return self.cmd.create_loop_induced(self, line, *args, **opts)
  File "/lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/madgraph/interface/madgraph_interface.py", line 5290, 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 : True (user set)
         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
                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 : /lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/HEPTools/fastjet/bin/fastjet-config (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : lhapdf-config
                 lhapdf_py2 : None
                 lhapdf_py3 : /lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/HEPTools/lhapdf6_py3/bin/lhapdf-config (user set)
          madanalysis5_path : None (user set)
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : /lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/HEPTools/MG5aMC_PY8_interface (user set)
                      ninja : ./HEPTools/lib
        output_dependencies : external
                   pineappl : pineappl
            pythia-pgs_path : ./pythia-pgs
               pythia8_path : /lstore/titan/atlas/emogollon/MG5_aMC_v3_4_2/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

Any guesses why this is happening?
Best,
Esteban Chalbaud

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) wrote : Re: [Bug 2011588] [NEW] Not able to run NLO process in models in Madgraph 3.4.2
Download full text (20.5 KiB)

So the issue is the default definition of the order that MG5aMC is taking as default.
This force the loop induced mode (that you likely do not want)
and then MG5aMC detects an inconsistency.

So the solution is likely to specify what you want.
My guess is that you want to do is:
> generate p p > xd xd~ j DMV<=2 [QCD]
(but obviously I can be wrong)

Cheers,

Olivier

> On 14 Mar 2023, at 15:17, Esteban Chalbaud <email address hidden> wrote:
>
> Public bug reported:
>
> Dear all I am having problems to run the following instructions in MadGraph:
> import model DMsimp_s_spin1
> define p = g u c d s b u~ c~ d~ s~ b~
> define j = g u c d s b u~ c~ d~ s~ b~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> generate p p > xd xd~ j [QCD]
>
> After I get:
>
> INFO: Checking for minimal orders which gives processes.
> INFO: Please specify coupling orders to bypass this step.
> INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS 2*DMV+QCD+2*QED
> Setting the born squared orders automatically in the process definition to QED^2<=4 QCD^2<=2 DMV=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 > xd xd~ j [QCD]" 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.
>
>
> The report says:
>
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.4.2 2023-01-20 *
> #* *
> #* 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...

Revision history for this message
Esteban Chalbaud (echalbau) wrote :

I think I cannot use this notation DMV<=2 [QCD] in this version of madgraph I get this error:

Command "generate p p > t t~ y1 DMV<=2 [QCD]" interrupted with error:
Exception : Potentially ambigious syntax detected. Note that the syntax of paper 1804.10017 (used in 3.0.x) is not used anymore (since version 3.1.0).
 If you want to follow the syntax of that paper, you can just replace "QED" by "aEW" and "QCD" by "aS".
 More information here: http://amcatnlo.cern.ch/co.htm
 If you know the current meaning of the syntax you can bypass this crash by running (once per machine) this command:
  set acknowledged_v3.1_syntax True --global
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) wrote : Re: [Bug 2011588] Not able to run NLO process in models in Madgraph 3.4.2

You have to run once the following command:
set acknowledged_v3.1_syntax True --global

and then it should work

Cheers,

Olivier

> On 22 Mar 2023, at 15:15, Esteban Chalbaud <email address hidden> wrote:
>
> set acknowledged_v3.1_syntax True --global

Revision history for this message
Esteban Chalbaud (echalbau) wrote :
Download full text (44.3 KiB)

Now it is working but i get a very long error after launch:

INFO: Starting run
Not able to open file /lstore/titan/atlas/emogollon/MG5_aMC_v3_1_1/bin/direct_xd/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 8 cores
INFO: Compiling P0_gu_xdxdxu...
INFO: Compiling P0_gc_xdxdxc...
INFO: Compiling P0_gd_xdxdxd...
INFO: Compiling P0_gd_xdxdxb...
INFO: Compiling P0_gs_xdxdxs...
INFO: Compiling P0_gux_xdxdxux...
INFO: Compiling P0_gcx_xdxdxcx...
INFO: Compiling P0_gdx_xdxdxdx...
STOP 1
INFO: P0_gd_xdxdxb done.
INFO: Compiling P0_gdx_xdxdxbx...
INFO: P0_gu_xdxdxu done.
INFO: Compiling P0_gsx_xdxdxsx...
INFO: P0_gcx_xdxdxcx done.
INFO: Compiling P0_gb_xdxdxd...
INFO: P0_gc_xdxdxc done.
INFO: Compiling P0_gb_xdxdxb...
INFO: P0_gux_xdxdxux done.
INFO: Compiling P0_gbx_xdxdxdx...
INFO: P0_gd_xdxdxd done.
INFO: Compiling P0_gbx_xdxdxbx...
INFO: P0_gs_xdxdxs done.
INFO: Compiling P0_ug_xdxdxu...
INFO: P0_gdx_xdxdxdx done.
INFO: Compiling P0_uux_xdxdxg...
STOP 1
STOP 1
STOP 1
INFO: P0_uux_xdxdxg done.
INFO: Compiling P0_cg_xdxdxc...
INFO: P0_gdx_xdxdxbx done.
INFO: Compiling P0_ccx_xdxdxg...
INFO: P0_ug_xdxdxu done.
INFO: Compiling P0_dg_xdxdxd...
INFO: P0_gb_xdxdxd done.
INFO: Compiling P0_dg_xdxdxb...
INFO: P0_gb_xdxdxb done.
INFO: Compiling P0_ddx_xdxdxg...
INFO: P0_gbx_xdxdxbx done.
INFO: Compiling P0_dbx_xdxdxg...
INFO: P0_gbx_xdxdxdx done.
INFO: Compiling P0_sg_xdxdxs...
INFO: P0_gsx_xdxdxsx done.
INFO: Compiling P0_ssx_xdxdxg...
STOP 1
STOP 1
INFO: P0_ccx_xdxdxg done.
INFO: Compiling P0_uxg_xdxdxux...
INFO: P0_dbx_xdxdxg done.
INFO: Compiling P0_uxu_xdxdxg...
INFO: P0_cg_xdxdxc done.
INFO: Compiling P0_cxg_xdxdxcx...
INFO: P0_dg_xdxdxb done.
INFO: Compiling P0_cxc_xdxdxg...
INFO: P0_ddx_xdxdxg done.
INFO: Compiling P0_dxg_xdxdxdx...
INFO: P0_dg_xdxdxd done.
INFO: Compiling P0_dxg_xdxdxbx...
INFO: P0_sg_xdxdxs done.
INFO: Compiling P0_dxd_xdxdxg...
INFO: P0_ssx_xdxdxg done.
INFO: Compiling P0_dxb_xdxdxg...
STOP 1
INFO: P0_uxu_xdxdxg done.
INFO: P0_cxc_xdxdxg done.
INFO: Compiling P0_sxg_xdxdxsx...
INFO: Compiling P0_sxs_xdxdxg...
INFO: P0_dxg_xdxdxbx done.
INFO: Compiling P0_bg_xdxdxd...
INFO: P0_dxg_xdxdxdx done.
INFO: Compiling P0_bg_xdxdxb...
INFO: P0_cxg_xdxdxcx done.
INFO: Compiling P0_bdx_xdxdxg...
STOP 1
INFO: P0_uxg_xdxdxux done.
INFO: Compiling P0_bbx_xdxdxg...
INFO: P0_dxb_xdxdxg done.
INFO: Compiling P0_bxg_xdxdxdx...
INFO: P0_dxd_xdxdxg done.
INFO: Compiling P0_bxg_xdxdxbx...
STOP 1
STOP 1
INFO: P0_sxs_xdxdxg done.
INFO: Compiling P0_bxd_xdxdxg...
INFO: P0_sxg_xdxdxsx done.
INFO: Compiling P0_bxb_xdxdxg...
INFO: P0_bdx_xdxdxg done.
INFO: P0_bg_xdxdxd done.
INFO: P0_bg_xdxdxb done.
STOP 1
INFO: P0_bxg_xdxdxdx done.
INFO: P0_bbx_xdxdxg done.
INFO: P0_bxg_xdxdxbx done.
STOP 1
INFO: P0_bxd_xdxdxg done.
INFO...

Changed in mg5amcnlo:
assignee: nobody → Kentarou Mawatari (kentarou-mawatari)
To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.