parameter scan fails

Asked by Semyon Obodenko

Dear experts,

I want to make such scan:

import model Pseudoscalar_2HDM
generate p p > h1 Xd Xd~ [noborn=QCD]
output bbMET_mA100
launch bbMET_mA100
set nevents 10000
set param_card sinp 0.35
set param_card mxd 10
set param_card laP1 3
set param_card laP2 3
set param_card lam3 3
set param_card gPXd 1
set param_card sinbma 1
set param_card tanbeta 1
set run_card ebeam1 7000
set run_card ebeam2 7000
set param_card mh1 125
set param_card mh2 100
set param_card mh3 100
set param_card mhc 100
set param_card mh4 scan:[100,200,300,400,500,600,700,800,1000]
launch bbMET_mA100 -i
print_results --path=./cross_section_mA100ma.txt --format=short

but after start it by typing
./mg5_aMC bbMET_mA100.txt
Madgraph gives a number of errors and warnings, the solution of which I tried to find on this site, but unsuccessfully

Madgraph gives:

[username@lxplus909 bin]$ ./mg5_aMC bbMET_mA100.txt
************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 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 *
* and *
* http://amcatnlo.web.cern.ch/amcatnlo/ *
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
************************************************************
load MG5 configuration from ../input/mg5_configuration.txt
fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
 MG5_aMC> set fastjet /PATH/TO/fastjet-config

set ninja to /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/HEPTools/lib
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
Note that you can still compile and run aMC@NLO with the built-in PDFs
 MG5_aMC> set lhapdf /PATH/TO/lhapdf-config

set lhapdf to /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/HEPTools/lhapdf6_py3/bin/lhapdf-config
Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
import /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100.txt
The import format was not given, so we guess it as command
import model Pseudoscalar_2HDM
INFO: Restrict model Pseudoscalar_2HDM with file ../models/Pseudoscalar_2HDM/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Defined multiparticle all = g ghg ghg~ u c d s u~ c~ d~ s~ a gha gha~ ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ ghz ghwp ghwm h+ h1 h2 h3 h4 w- ghz~ ghwp~ ghwm~ h- ta- xd ta+ xd~
generate p p > h1 Xd Xd~ [noborn=QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS QCD+2*QED+DS
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS QCD+2*QED+DS
INFO: Contributing diagrams generated: 0 Born, 60(+4) loops, 0 R2, 0 UV
INFO: Contributing diagrams generated: 0 Born, 2 loops, 2 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, 2 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, 2 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, 2 R2, 0 UV
INFO: Process u~ u > h1 xd xd~ added to mirror process u u~ > h1 xd xd~
INFO: Process c~ c > h1 xd xd~ added to mirror process c c~ > h1 xd xd~
INFO: Process d~ d > h1 xd xd~ added to mirror process d d~ > h1 xd xd~
INFO: Process s~ s > h1 xd xd~ added to mirror process s s~ > h1 xd xd~
5 processes with 68 diagrams generated in 8.021 s
Total: 5 processes with 68 diagrams
output bbMET_mA100
INFO: directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100 already exists.
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n]
found line : launch bbMET_mA100
This answer is not valid for current question. Keep it for next question and use here default: y
INFO: initialize a new directory: bbMET_mA100
INFO: remove old information in bbMET_mA100
WARNING: The collier reduction library could not be found with PATH:/afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/HEPTools/lib specified in mg5_configuration.txt. It will not be available.
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: u u~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: c c~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Helas calls for process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Processing color information for loop process: d d~ > h1 xd xd~ [ noborn = QCD ]
INFO: Creating color matrix loop process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Combined process s s~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ] with process d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_0_1_gg_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_h1xdxdx
INFO: Generating Feynman diagrams for Process: g g > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_h1xdxdx
INFO: Creating files in directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_1_uux_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_2_ccx_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_3_ddx_h1xdxdx
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Creating files in directory P0_qq_h1xdxdx
INFO: Generating Feynman diagrams for Process: u u~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Feynman diagrams for Process: c c~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Generating Feynman diagrams for Process: d d~ > h1 xd xd~ WEIGHTED<=21 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group qq_h1xdxdx
Generated helas calls for 4 subprocesses (66 diagrams) in 0.359 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 35 routines in 1.346 s
save configuration file to /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100 done.
Type "launch" to generate events from this process, or see
/afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/README
Run "open index.html" to see more information about this process.
launch bbMET_mA100
************************************************************
* *
* 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 3.4.2 2023-01-20 *
* *
* 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 /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/Cards/me5_configuration.txt
INFO: load configuration from /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/input/mg5_configuration.txt
INFO: load configuration from /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/Cards/me5_configuration.txt
generate_events run_01
set nevents 10000
INFO: modify parameter nevents of the run_card.dat to 10000
set param_card sinp 0.35
INFO: modify param_card information BLOCK higgs with id (5,) set to 0.35
set param_card mxd 10
INFO: modify param_card information BLOCK mass with id (52,) set to 10.0
set param_card laP1 3
INFO: modify param_card information BLOCK higgs with id (2,) set to 3.0
set param_card laP2 3
INFO: modify param_card information BLOCK higgs with id (3,) set to 3.0
set param_card lam3 3
INFO: modify param_card information BLOCK higgs with id (1,) set to 3.0
set param_card gPXd 1
INFO: modify param_card information BLOCK dminputs with id (1,) set to 1.0
set param_card sinbma 1
INFO: modify param_card information BLOCK frblock with id (3,) set to 1.0
set param_card tanbeta 1
INFO: modify param_card information BLOCK frblock with id (2,) set to 1.0
set run_card ebeam1 7000
INFO: modify parameter ebeam1 of the run_card.dat to 7000.0
set run_card ebeam2 7000
INFO: modify parameter ebeam2 of the run_card.dat to 7000.0
set param_card mh1 125
INFO: modify param_card information BLOCK mass with id (25,) set to 125.0
set param_card mh2 100
INFO: modify param_card information BLOCK mass with id (35,) set to 100.0
set param_card mh3 100
INFO: modify param_card information BLOCK mass with id (36,) set to 100.0
set param_card mhc 100
INFO: modify param_card information BLOCK mass with id (37,) set to 100.0
set param_card mh4 scan:[100,200,300,400,500,600,700,800,1000]
INFO: modify param_card information BLOCK mass with id (55,) set to scan:[100,200,300,400,500,600,700,800,1000]
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]
found line : launch bbMET_mA100 -i
This answer is not valid for current question. Keep it for next question and use here default: 0
INFO: Create the next param_card in the scan definition
INFO: change parameter mass with code (55,) to 100
INFO: Update the dependent parameter of the param_card.dat
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
restorecon: SELinux: Could not get canonical path for /afs/cern.ch/user/s/username/.mozilla/firefox/* restorecon: No such file or directory.
initMadLoop -r -f
libGL error: No matching fbConfigs or visuals found
libGL error: failed to load driver: swrast
Initializing MadLoop loop-induced matrix elements (this can take some time)...
root: Error while executing ./check in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_1_uux_h1xdxdx
root: Failed at running the process in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_1_uux_h1xdxdx.
root: Error while executing ./check in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_2_ccx_h1xdxdx
root: Failed at running the process in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_2_ccx_h1xdxdx.
root: Error while executing ./check in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_3_ddx_h1xdxdx
root: Failed at running the process in /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/SubProcesses/PV0_1_3_ddx_h1xdxdx.
Error detected in "generate_events run_01"
write debug file /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : Failed the initialization of loop-induced matrix element 'PV0_1_1_uux_h1xdxdx' (trying with a maximum of 18 PS points).
INFO:
command not executed: print_results --path=./cross_section_mA100ma.txt --format=short
quit
INFO:
launch bbMET_mA100 -i
************************************************************
* *
* 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 3.4.2 2023-01-20 *
* *
* 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 /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/Cards/me5_configuration.txt
INFO: load configuration from /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/input/mg5_configuration.txt
INFO: load configuration from /afs/cern.ch/user/s/username/private/MG5_aMC_v3_4_1/bin/bbMET_mA100/Cards/me5_configuration.txt
quit
INFO:
quit
[username@lxplus909 bin]$

...
What am I doing wrong?

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

This is an error that I have never fully understood and which is transient.
The issue is that we fail to initialize the computation of the loop-induced process (which is non trivial due to the lack of born to get a scale). My advise is typically to retry multiple times.
If it is consistent, the question is then if your loop is not zero due to some symmetry.

I assign this to Valentin who might want to add something.

Cheers,

Olivier

Can you help with this problem?

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

To post a message you must log in.