parameter scan fails
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,
launch bbMET_mA100 -i
print_results --path=
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:/
* and *
* http://
* *
* 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/
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/
MG5_aMC> set fastjet /PATH/TO/
set ninja to /afs/cern.
lhapdf-config does not seem to correspond to a valid lhapdf-config executable.
Please set the 'lhapdf' variable to the (absolute) /PATH/TO/
Note that you can still compile and run aMC@NLO with the built-in PDFs
MG5_aMC> set lhapdf /PATH/TO/
set lhapdf to /afs/cern.
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
import /afs/cern.
The import format was not given, so we guess it as command
import model Pseudoscalar_2HDM
INFO: Restrict model Pseudoscalar_2HDM with file ../models/
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.
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/
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.
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.
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.
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.
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.
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /afs/cern.
Type "launch" to generate events from this process, or see
/afs/cern.
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:/
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
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,
INFO: modify param_card information BLOCK mass with id (55,) set to scan:[100,
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/
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.
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.
root: Failed at running the process in /afs/cern.
root: Error while executing ./check in /afs/cern.
root: Failed at running the process in /afs/cern.
root: Error while executing ./check in /afs/cern.
root: Failed at running the process in /afs/cern.
Error detected in "generate_events run_01"
write debug file /afs/cern.
If you need help with this issue please contact us on https:/
str : Failed the initialization of loop-induced matrix element 'PV0_1_
INFO:
command not executed: print_results --path=
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:/
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
quit
INFO:
quit
[username@lxplus909 bin]$
...
What am I doing wrong?
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- Valentin Hirschi Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Semyon Obodenko for more information if necessary.