MadGraph5_aMC@NLO

Problem generating loop induced process with loop_sm_MSbar_yb_yc model

Asked by Tiziano

Dear MG experts,

I'm trying to generate a loop induced process using the loop_sm_MSbar_yb_yc model (we partially discussed this topic with some of you in an e-mail thread). I would like to generate some gridpacks to use in the CMS workflow but for now I have tried a MadGraph standalone approach. I've tested this setup with MG5_aMC_v2_6_5 and MG5_aMC_v2_9_18, both times without success.

The process I'm trying to generate is the following:

define lj = g u u~ d d~ s s~

generate g g > h [noborn=QCD]
add process lj lj > h lj [noborn=QCD]

If I try to generate this using the default loop_sm model the process runs through, while as soon as I add the line:

import model loop_sm_MSbar_yb_yc-yc4FS

the generation doesn't run through, stopping with the following output:

MG5_aMC>import model loop_sm_MSbar_yb_yc-yc4FS
INFO: Restrict model loop_sm_MSbar_yb_yc with file models/loop_sm_MSbar_yb_yc/restrict_yc4FS.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
MG5_aMC>define lj = g u u~ d d~ s s~
Defined multiparticle lj = g u d s u~ d~ s~
MG5_aMC>generate g g > h [noborn=QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS 2*QED+QCD
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
1 processes with 2 diagrams generated in 0.113 s
Total: 1 processes with 2 diagrams
MG5_aMC>add process lj lj > h lj [noborn=QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=6: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=7: WEIGTHED IS 2*QED+QCD
INFO: Contributing diagrams generated: 0 Born, 12 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
INFO: Process u g > h u added to mirror process g u > h u
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
INFO: Process d g > h d added to mirror process g d > h d
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
INFO: Process s g > h s added to mirror process g s > h s
WARNING: Some loop diagrams contributing to this process are discarded because they are not pure (QCD)-perturbation.
Make sure you did not want to include them.
INFO: Contributing diagrams generated: 0 Born, 2 loops, 0 R2, 0 UV
INFO: Process u~ g > h u~ added to mirror process g u~ > h u~
INFO: Process u~ u > h g added to mirror process u u~ > h g
INFO: Process d~ g > h d~ added to mirror process g d~ > h d~
INFO: Process d~ d > h g added to mirror process d d~ > h g
INFO: Process s~ g > h s~ added to mirror process g s~ > h s~
INFO: Process s~ s > h g added to mirror process s s~ > h g
10 processes with 30 diagrams generated in 3.930 s
Total: 11 processes with 32 diagrams
MG5_aMC>output yc-loop_2
INFO: initialize a new directory: yc-loop_2
INFO: remove old information in yc-loop_2
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h g [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Generating Helas calls for process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Generating Helas calls for process: g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Processing color information for loop process: g u > h u [ noborn = QCD ]
INFO: Creating color matrix loop process: g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process g d > h d WEIGHTED<=16 [ noborn = QCD ] with process g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process g s > h s WEIGHTED<=16 [ noborn = QCD ] with process g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Generating Helas calls for process: g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Processing color information for loop process: g u~ > h u~ [ noborn = QCD ]
INFO: Creating color matrix loop process: g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process g d~ > h d~ WEIGHTED<=16 [ noborn = QCD ] with process g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process g s~ > h s~ WEIGHTED<=16 [ noborn = QCD ] with process g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Generating Helas calls for process: u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Processing color information for loop process: u u~ > h g [ noborn = QCD ]
INFO: Creating color matrix loop process: u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process d d~ > h g WEIGHTED<=16 [ noborn = QCD ] with process u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Combined process s s~ > h g WEIGHTED<=16 [ noborn = QCD ] with process u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Creating files in directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_0_1_gg_hg
ALOHA: aloha creates FFS1 set of routines with options: L2
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_hg
INFO: Generating Feynman diagrams for Process: g g > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_hg
INFO: Creating files in directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_1_1_gg_h
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_h
INFO: Generating Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_h
INFO: Creating files in directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_1_gu_hu
ALOHA: aloha creates FFV1 set of routines with options: L2
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Creating files in directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_2_gux_hux
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Creating files in directory P0_gq_hq
INFO: Generating Feynman diagrams for Process: g u > h u WEIGHTED<=16 [ noborn = QCD ]
INFO: Generating Feynman diagrams for Process: g u~ > h u~ WEIGHTED<=16 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gq_hq
INFO: Creating files in directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_3_1_uux_hg
ALOHA: aloha creates FFS1 set of routines with options: L2
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Creating files in directory P0_qq_hg
INFO: Generating Feynman diagrams for Process: u u~ > h g WEIGHTED<=16 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group qq_hg
Generated helas calls for 5 subprocesses (20 diagrams) in 0.185 s
save configuration file to /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_18/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
INFO: Generate web pages
Output to directory /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2 done.
Type "launch" to generate events from this process, or see
/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch
************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.6.5 2018-02-03 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/Cards/me5_configuration.txt
INFO: load configuration from /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/input/mg5_configuration.txt
INFO: load configuration from /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
No valid eps viewer found. Please set in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|============= values ==============|======== other options ========\
| 1. Choose the shower/hadronization program | shower = Not Avail. | Please install module |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\=======================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. MadLoopParams : MadLoopParams.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, MadLoopParams, enter path][90s to answer]
>2
open /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/Cards/run_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. MadLoopParams : MadLoopParams.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, MadLoopParams, enter path]
>1
open /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/Cards/param_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. MadLoopParams : MadLoopParams.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, MadLoopParams, enter path]
>
INFO: Update the dependent parameter of the param_card.dat
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
INFO: Using LHAPDF v6.2.1 interface for PDFs
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
root: Error while executing ./check in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_3_1_uux_hg
root: Failed at running the process in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_3_1_uux_hg.
root: Error while executing ./check in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_2_gux_hux
root: Failed at running the process in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_2_gux_hux.
root: Error while executing ./check in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_1_1_gg_h
root: Failed at running the process in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_1_1_gg_h.
root: Error while executing ./check in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_1_gu_hu
root: Failed at running the process in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_2_1_gu_hu.
root: Error while executing ./check in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_0_1_gg_hg
root: Failed at running the process in /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/SubProcesses/PV0_0_1_gg_hg.
Error detected in "generate_events run_01"
write debug file /work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/yc-loop_2/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
MadGraph5Error : Failed the initialization of loop-induced matrix element 'PV0_0_1_gg_hg' (trying with a maximum of 18 PS points).
quit

If I open the log file I can get this further piece information:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/common_run_interface.py", line 6869, in new_fct
    original_fct(obj, *args, **opts)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2508, in run_generate_events
    postcmd=False)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 3277, in do_survey
    self.configure_directory()
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 5601, in configure_directory
    self.do_treatcards('')
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 3254, in do_treatcards
    self.exec_cmd('initMadLoop -r -f')
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 2613, in do_initMadLoop
    subproc_prefix='PV', MG_options=self.options, interface=self)
  File "/work/bevila_t/HpC_Analysis/CMSSW_10_6_18/src/MG_playground/MG5_aMC_v2_6_5/madgraph/interface/madevent_interface.py", line 7061, in init_MadLoop
    %(max_mult*n_PS))
MadGraph5Error: Failed the initialization of loop-induced matrix element 'PV0_0_1_gg_hg' (trying with a maximum of 18 PS points).
                              Run Options
                              -----------
               stdout_level : 20 (user set)

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : None (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 : lhapdf-config
              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

If I try to run with the newer release MG5_aMC_v2_9_18 the outcome is the same, with the additional issue that I have to run with python3.10 due to this warning which otherwise would generate an error:

WARNING: [Errno 2] No such file or directory: '/Users/tiziano-bevilacqua/madgraph/MG5_aMC_v2_9_18/Template/loop_material/StandAlone/SubProcesses/MadLoop5_resources'

so in this case I also have to convert the model using the command:

convert model /Users/tiziano-bevilacqua/madgraph/MG5_aMC_v2_9_18/models/loop_sm_MSbar_yb_yc

Do you have any suggestion on how to solve this problem?

Best,
Tiziano

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Last query:
Last reply:
Revision history for this message
Tiziano (tbevilac) said : 		#1

Actually the generation also goes through if one does:

import model loop_sm_MSbar_yb_yc

instead of:

import model loop_sm_MSbar_yb_yc-yc4FS

so the problem seems to originate from the yc4FS restriction

Revision history for this message
marco zaro (marco-zaro) said : 		#2

Hi Tiziano,
this is strange; the only differences between the two restrictions are that restrict_default has a tau yukawa !=0 (which I can bet it is irrelevant), and has ymb=mb. In restrict_yb4FS ymb and mb differ.
Can you check what happens putting them both to the same value?

Thanks,

Marco

Revision history for this message
Tiziano (tbevilac) said : 		#3

Hi Marco,
Isn't the restricted card for loop_sm_MSbar_yb_yc-yc4FS restrict_yc4FS.dat? here I have that ymb=ytau=0 and that ymc=mc. If I add back ymb=mb together with ymc=mc the generation indeed works fine again. It appears that if I only have ymc!=0 it generates some problem.

Thanks,
Tiziano

Revision history for this message
Tiziano (tbevilac) said : 		#4

Hi Marco,
I haven't manage to fix the issue, do you have any further suggestion?

Thanks,
Tiziano

Revision history for this message
marco zaro (marco-zaro) said : 		#5

Hi,
the diagram is zero unless mc (not only yb ) is set to !=0.

Can you help with this problem?

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

To post a message you must log in.