NLO generation
Dear all,
I would like to simulate NLO process using MadGraph5 version 3.4.1. I would like to shower using Pythia8 and use Delphes as a detector.
First thing, I cannot see any switch for the detector simulation. I was wondering if there is a way to use Delphes within the MadGraph5 interface or I should run it independently.
As for the error I am getting, I am just following the tutorial using SM and generating
generate p p > e+ ve [QCD]
but getting a "Configuration error" that I cannot understand. I point out that all LO generations work fine and with any model (also UFO written by myself).
generate
I paste here the output of the terminal and the log files.
*******
Terminal output:
*******
gcolor@
*******
* *
* 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/mg5_
set fastjet to /home/gcolor/
set ninja to /home/gcolor/
set collier to /home/gcolor/
set lhapdf to /home/gcolor/
set lhapdf to /home/gcolor/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
Loading default model: sm
INFO: Restrict model sm 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 u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
MG5_aMC>generate p p > e+ ve [QCD]
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
import model loop_sm
INFO: Restrict model loop_sm with file models/
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / 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+
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
Setting the born squared orders automatically in the process definition to QED^2<=4 QCD=0.
If this is not what you need, please regenerate with the correct orders.
WARNING: Use of multiparticles is non-trivial for NLO process generation and depends on the orders included, the process considered, as well as the PDF set chosen. See appendix D of arXiv:1804.10017 [hep-ph] for some guidance.
INFO: Generating FKS-subtracted matrix elements for born process: u d~ > e+ ve QCD^2=2 QED^2=4 (1 / 4)
INFO: Generating FKS-subtracted matrix elements for born process: c s~ > e+ ve QCD^2=2 QED^2=4 (2 / 4)
INFO: Generating FKS-subtracted matrix elements for born process: d~ u > e+ ve QCD^2=2 QED^2=4 (3 / 4)
INFO: Generating FKS-subtracted matrix elements for born process: s~ c > e+ ve QCD^2=2 QED^2=4 (4 / 4)
INFO: Generating virtual matrix element with MadLoop for process: u d~ > ve e+ QCD^2=2 QED^2=4 (1 / 4)
INFO: Generating virtual matrix element with MadLoop for process: c s~ > ve e+ QCD^2=2 QED^2=4 (2 / 4)
INFO: Generating virtual matrix element with MadLoop for process: d~ u > ve e+ QCD^2=2 QED^2=4 (3 / 4)
INFO: Generating virtual matrix element with MadLoop for process: s~ c > ve e+ QCD^2=2 QED^2=4 (4 / 4)
INFO: Generated 4 subprocesses with 24 real emission diagrams, 4 born diagrams and 4 virtual diagrams
MG5_aMC>output MY_FIRST_
INFO: Writing out the aMC@NLO code, using optimized Loops
INFO: directory /home/gcolor/
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n][60s to answer]
>
INFO: initialize a new directory: MY_FIRST_
INFO: remove old information in MY_FIRST_
INFO: Generating real emission matrix-elements...
INFO: Generating Helas calls for FKS process: u d~ > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4 (1 / 4)
INFO: Processing color information for process: u d~ > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
INFO: Generating Helas calls for FKS process: c s~ > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4 (2 / 4)
INFO: Generating Helas calls for FKS process: d~ u > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4 (3 / 4)
INFO: Processing color information for process: d~ u > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
INFO: Generating Helas calls for FKS process: s~ c > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4 (4 / 4)
INFO: ... Done
INFO: Writing files in P0_udx_veep (1 / 2)
INFO: Creating files in directory V0_udx_veep
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: u d~ > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
INFO: Generating born Feynman diagrams for Process: u d~ > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
INFO: Writing files in P0_dxu_veep (2 / 2)
INFO: Creating files in directory V0_dxu_veep
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: d~ u > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
INFO: Generating born Feynman diagrams for Process: d~ u > ve e+ [ all = QCD QED ] QCD^2=2 QED^2=4
History written to /home/gcolor/
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 8 routines in 0.130 s
save configuration file to /home/gcolor/
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
INFO: Generating MadAnalysis5 default cards tailored to this process
INFO: Cards created in 1.80s
Type "launch" to generate events from this process, or see
/home/gcolor/
Run "open index.html" to see more information about this process.
MG5_aMC>launch
INFO: *******
* *
* 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 *
* http://
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /home/gcolor/
INFO: load configuration from /home/gcolor/
INFO: load configuration from /home/gcolor/
Using default text editor "vi". Set another one in ./input/
launch auto
The following switches determine which programs are run:
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = HERWIG6 | OFF|PYTHIA6Q|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' 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]
>3
The following switches determine which programs are run:
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = OFF | PYTHIA6Q|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA6Q | PYTHIA6PT|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA6PT | PYTHIA8|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>3
The following switches determine which programs are run:
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA8 | HERWIG6|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Either type the switch number (1 to 6) to change its setting,
Set any switch explicitly (e.g. type 'order=LO' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>
INFO: will run in mode: aMC@NLO
INFO: modify parameter parton_shower of the run_card.dat to PYTHIA8
Do you want to edit a card (press enter to bypass editing)?
/------
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. shower : shower_card.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, shower, enter path][90s to answer]
>
INFO: modify parameter extralibs of the shower_card.dat to pythia8 z dl stdc++
INFO: modify parameter extrapaths of the shower_card.dat to ../lib /home/gcolor/
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.119
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
Gtk-Message: 11:10:24.480: Not loading module "atk-bridge": The functionality is provided by GTK natively. Please try to not load it.
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_udx_veep...
INFO: Compiling P0_dxu_veep...
INFO: P0_udx_veep done.
INFO: P0_dxu_veep done.
INFO: Checking test output:
INFO: P0_udx_veep
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxu_veep
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grids
INFO: Idle: 0, Running: 0, Completed: 2 [ current time: 11h10 ]
sum of cpu time of last step: 0 second
INFO: Determining the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 1.083e+04 +- 7.3e+01 pb
Total abs(cross section): 1.263e+04 +- 7.4e+01 pb
INFO: Computing upper envelope
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h10 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.29s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 0.29s ]
sum of cpu time of last step: 0 second
INFO: Updating the number of unweighted events per channel
Intermediate results:
Random seed: 33
Total cross section: 1.088e+04 +- 4.5e+01 pb
Total abs(cross section): 1.272e+04 +- 4.4e+01 pb
INFO: Generating events
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h10 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 1.4s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 1.6s ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 11h10 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.2s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 0.26s ]
INFO: Collecting events
INFO:
----
Summary:
Process p p > e+ ve [QCD]
Run at p-p collider (6500.0 + 6500.0 GeV)
Number of events generated: 10000
Total cross section: 1.088e+04 +- 4.5e+01 pb
----
Scale variation (computed from LHE events):
----
INFO: The /home/gcolor/
INFO: Events generated
reweight -from_cards
decay_events -from_cards
INFO: Preparing MCatNLO run
INFO: Using LHAPDF v6.5.3 interface for PDFs
INFO: Compiling MCatNLO for PYTHIA8...
INFO: ... done
INFO: Showering events...
INFO: (Running in /home/gcolor/
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 11h11 ]
mv: cannot stat 'Pythia8.hep': No such file or directory
gzip: events.hepmc: No such file or directory
INFO: Idle: 0, Running: 0, Completed: 1 [ 0.41s ]
Error detected in "launch auto "
write debug file /home/gcolor/
If you need help with this issue, please, contact us on https:/
str : No file has been generated, an error occurred. More information in /home/gcolor/
quit
INFO:
MG5_aMC>
*******
ME5_debug:
*******
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.4.2 2023-01-20 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* and *
#* http://
#* *
#******
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#******
launch auto
Traceback (most recent call last):
File "/usr/lib/
os.rename(src, real_dst)
FileNotFoundError: [Errno 2] No such file or directory: '/home/
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/home/
files.
File "/home/
shutil.
File "/usr/lib/
copy_
File "/usr/lib/
copyfile(src, dst, follow_
File "/usr/lib/
with open(src, 'rb') as fsrc:
FileNotFoundError: [Errno 2] No such file or directory: '/home/
During handling of the above exception, another exception occurred:
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
raise aMCatNLOError('No file has been generated, an error occurred.'+\
madgraph.
Value of current Options:
exrootana
cluster_
mg5amc_
output_
acknowledged_
ignore_
low_mem_
complex_
include_
loop_
max_
default_
nlo_
automatic_
notificat
*******
mcatnlo_run.log:
*******
*-----
| |
| *------
| | | |
| | | |
| | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | |
| | P P Y Y T H H I A A This is PYTHIA version 8.306 | |
| | PPP Y T HHHHH I AAAAA Last date of change: 28 Jun 2021 | |
| | P Y T H H I A A | |
| | P Y T H H III A A Now is 09 Jun 2023 at 11:11:06 | |
| | | |
| | Program documentation and an archive of historic versions is found on: | |
| | | |
| | https:/
| | | |
| | PYTHIA is authored by a collaboration consisting of: | |
| | | |
| | Christian Bierlich, Nishita Desai, Leif Gellersen, Ilkka Helenius, Philip | |
| | Ilten, Leif Lönnblad, Stephen Mrenna, Stefan Prestel, Christian Preuss, | |
| | Torbjörn Sjöstrand, Peter Skands, Marius Utheim and Rob Verheyen. | |
| | | |
| | The complete list of authors, including contact information and | |
| | affiliations, can be found on https:/
| | Problems or bugs should be reported on email at <email address hidden>. | |
| | | |
| | The main program reference is 'An Introduction to PYTHIA 8.2', | |
| | T. Sjöstrand et al, Comput. Phys. Commun. 191 (2015) 159 | |
| | [arXiv:1410.3012 [hep-ph]] | |
| | | |
| | The main physics reference is the 'PYTHIA 6.4 Physics and Manual', | |
| | T. Sjöstrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175] | |
| | | |
| | PYTHIA is released under the GNU General Public Licence version 2 or later.| |
| | Please respect the MCnet Guidelines for Event Generator Authors and Users. | |
| | | |
| | Disclaimer: this program comes without any guarantees. | |
| | Beware of errors and use common sense when interpreting results. | |
| | | |
| | Copyright (C) 2021 Torbjörn Sjöstrand | |
| | | |
| | | |
| *------
| |
*-----
*-----
| Last Change 12 Apr 2017 - P. Skands
| Parsing: events.lhe
| (SLHA::readFile) line 86 - storing non-SLHA(2) block: yukawa
| (SLHA::readFile) line 112 - storing QNUMBERS for id = 82 gh
*-----
PYTHIA Info from SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Warning in SLHAinterface:
PYTHIA Error in Plugin::Plugin: libpythia8lhapd
PYTHIA Error in Pythia::init: could not set up PDF for beam A
PYTHIA Abort from Pythia::init: PDF initialization failed
Error: could not initialise Pythia
I thank you very much for your help!!
Best regards,
Guglielmo Coloretti
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Guglielmo Coloretti
- Solved:
- Last query:
- Last reply: