MadGraph5_aMC@NLO

Running mg5_amc@NLO with pythia

Asked by Andrew McAdams

Hi! I am trying to simulate the process

generate p p > b b~ [QCD]

When trying to launch it using the hadronization option pythia8, with the default options in run_card, I get empty .hempc file. This does not happen if launching the process p p > b b~.

From the logs that I include below, it follows the following error:

PYTHIA Error in LHAPDF::init: libpythia8lhapdf6.so: cannot open shared object file: No such file or directory

However, I have previously successfully installed lhapdf6 inside madgraph.

Could you please tell me what may be a problem?

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Andrew McAdams
Solved:
Last query:
Last reply:
Revision history for this message
Andrew McAdams (mcadams1994) said : 		#1

Below, there is some additional info:

mcatnlo.log

 USING PYTHIA 8, version >= 8.200
g++ -O -I/home/name/Downloads/mg5_new/HEPTools/pythia8/include \
  -I/home/name/Downloads/mg5_new/bin/b-production/MCatNLO/include Pythia82.cc -o Pythia8.exe \
  -L/home/name/Downloads/mg5_new/HEPTools/pythia8/lib -lpythia8 \
  -I/home/name/Downloads/mg5_new/HEPTools/hepmc/include -L/home/name/Downloads/mg5_new/HEPTools/hepmc/lib -lHepMC -L../lib -L/home/name/Downloads/mg5_new/HEPTools/pythia8//lib -L/home/name/Downloads/mg5_new/HEPTools/zlib/lib -lpythia8 -lz -ldl -lstdc++ \

***** Now ready for showering 10000 events with PYTHIA8

Revision history for this message
Andrew McAdams (mcadams1994) said : 		#2

From terminal:

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 2.9.3 2021-03-25 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* http://amcatnlo.cern.ch *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/name/Downloads/mg5_new/bin/b-production/Cards/amcatnlo_configuration.txt
INFO: load configuration from /home/name/Downloads/mg5_new/input/mg5_configuration.txt
INFO: load configuration from /home/name/Downloads/mg5_new/bin/b-production/Cards/amcatnlo_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
launch auto
The following switches determine which programs are run:
/=============================================================================\
1. Type of perturbative computation order = NLO
2. No MC@[N]LO matching / event generation fixed_order = OFF
3. Shower the generated events shower = PYTHIA6Q
4. Decay onshell particles madspin = OFF
5. Add weights to events for new hypp. reweight = OFF
6. Run MadAnalysis5 on the events generated madanalysis = Not Avail.
\=============================================================================/
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:
/=============================================================================\
1. Type of perturbative computation order = NLO
2. No MC@[N]LO matching / event generation fixed_order = OFF
3. Shower the generated events shower = PYTHIA6PT
4. Decay onshell particles madspin = OFF
5. Add weights to events for new hypp. reweight = OFF
6. Run MadAnalysis5 on the events generated madanalysis = Not Avail.
\=============================================================================/
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:
/=============================================================================\
1. Type of perturbative computation order = NLO
2. No MC@[N]LO matching / event generation fixed_order = OFF
3. Shower the generated events shower = PYTHIA8
4. Decay onshell particles madspin = OFF
5. Add weights to events for new hypp. reweight = OFF
6. Run MadAnalysis5 on the events generated madanalysis = Not Avail.
\=============================================================================/
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.
>0
 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/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, shower, enter path][90s to answer]
>0
 INFO: modify parameter extrapaths of the shower_card.dat to ../lib /home/name/Downloads/mg5_new/HEPTools/pythia8//lib /home/name/Downloads/mg5_new/HEPTools/zlib/lib
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
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 16 cores
INFO: Compiling P0_gg_bbx...
INFO: Compiling P0_uux_bbx...
INFO: Compiling P0_uxu_bbx...
INFO: P0_uxu_bbx done.
INFO: P0_uux_bbx done.
INFO: P0_gg_bbx done.
INFO: Checking test output:
INFO: P0_gg_bbx
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_uux_bbx
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_uxu_bbx
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: 5, Completed: 0 [ current time: 13h46 ]
INFO: Idle: 0, Running: 4, Completed: 1 [ current time: 13h46 ]
INFO: Idle: 0, Running: 3, Completed: 2 [ 0.034s ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 1.1s ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 1.1s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 5s ]
sum of cpu time of last step: 0 second
INFO: Determining the number of unweighted events per channel

      Intermediate results:
      Random seed: 38
      Total cross section: 3.849e+08 +- 1.1e+07 pb
      Total abs(cross section): 8.224e+08 +- 1.2e+07 pb

INFO: Computing upper envelope
INFO: Idle: 0, Running: 5, Completed: 0 [ current time: 13h46 ]
INFO: Idle: 0, Running: 4, Completed: 1 [ 0.66s ]
INFO: Idle: 0, Running: 3, Completed: 2 [ 0.74s ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 2s ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 2.1s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 4.8s ]
sum of cpu time of last step: 0 second
INFO: Updating the number of unweighted events per channel

      Intermediate results:
      Random seed: 38
      Total cross section: 3.741e+08 +- 6.7e+06 pb
      Total abs(cross section): 8.664e+08 +- 7.0e+06 pb

INFO: Generating events
INFO: Idle: 0, Running: 3, Completed: 2 [ current time: 13h46 ]
INFO: Idle: 0, Running: 2, Completed: 3 [ 0.35s ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 7.1s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 9.2s ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 2, Completed: 3 [ current time: 13h47 ]
INFO: Idle: 0, Running: 1, Completed: 4 [ 0.61s ]
INFO: Idle: 0, Running: 0, Completed: 5 [ 0.65s ]
INFO: Collecting events
INFO:
   --------------------------------------------------------------
      Summary:
      Process p p > b b~ [QCD]
      Run at p-p collider (6500.0 + 6500.0 GeV)
      Number of events generated: 10000
      Total cross section: 3.741e+08 +- 6.7e+06 pb
   --------------------------------------------------------------
      Scale variation (computed from LHE events):
          Dynamical_scale_choice -1 (envelope of 9 values):
              3.644e+08 pb +19.4% -45.2%
   --------------------------------------------------------------

INFO: The /home/name/Downloads/mg5_new/bin/b-production/Events/run_06/events.lhe.gz file has been generated.

INFO: Events generated
reweight -from_cards
decay_events -from_cards
INFO: Preparing MCatNLO run
INFO: Using LHAPDF v6.3.0 interface for PDFs
INFO: Compiling MCatNLO for PYTHIA8...
INFO: ... done
INFO: Showering events...
INFO: (Running in /home/name/Downloads/mg5_new/bin/b-production/MCatNLO/RUN_PYTHIA8_4)
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 13h47 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 0.29s ]
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 13h47 ]
INFO: The file /home/name/Downloads/mg5_new/bin/b-production/Events/run_06/events_PYTHIA8_0.hepmc.gz has been generated.
It contains showered and hadronized events in the HEPMC format obtained showering the parton-level event file /home/name/Downloads/mg5_new/bin/b-production/Events/run_06/events.lhe.gz with PYTHIA8
INFO: Run complete
INFO:
quit
INFO:

Revision history for this message
Andrew McAdams (mcadams1994) said : 		#3

File mcatnlo_run.log from folder RUN_PYTHIA8_1

PYTHIA Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off
 PYTHIA Warning in SLHAinterface::initSLHA: ignoring QNUMBERS for id = 82 (already exists)
 PYTHIA Info from SLHAinterface::initSLHA: importing MASS entries for id = {25}
 PYTHIA Warning in SLHAinterface::initSLHA: ignoring MASS entries for id = {1,2,3,4,5,6,11,12,13,14,15,16,21,22,23,24} (SLHA:keepSM. Use id > 1000000 for new particles)
 PYTHIA Warning in SLHAinterface::initSLHA: ignoring DECAY tables for id = {6,23,24,1,2,3,4,5,11,12,13,14,15,16,21,22} (SLHA:keepSM. Use id > 1000000 for new particles)
 PYTHIA Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables for id = {25} (total width provided but no Branching Ratios)
 PYTHIA Error in LHAPDF::init: libpythia8lhapdf6.so: cannot open shared object file: No such file or directory
 PYTHIA Error in Pythia::init: could not set up PDF for beam A
 PYTHIA Abort from Pythia::init: PDF initialization failed
 PYTHIA Abort from Pythia::next: not properly initialized so cannot generate events

 *------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------*
 | |
 | Subprocess Code | Number of events | sigma +- delta |
 | | Tried Selected Accepted | (estimated) (mb) |
 | | | |
 |-----------------------------------------------------------------------------------------------------------------|
 | | | |
 | | | |
 | sum | 0 0 0 | 0.000e+00 0.000e+00 |
 | |
 *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------*

 *------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 1 Abort from Pythia::init: PDF initialization failed |
 | 10 Abort from Pythia::next: not properly initialized so cannot generate events |
 | 1 Error in LHAPDF::init: libpythia8lhapdf6.so: cannot open shared object file: No such file or directory |
 | 1 Error in Pythia::init: could not set up PDF for beam A |
 | 1 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 1 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 1 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring QNUMBERS |
 | 1 Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#4

Hi,

I do not have any clue, this sounds an installation issue with pythia8 and the only suggestion that I can give is to contact Pythia8 authors.

Cheers,

Olivier

Revision history for this message
Andrew McAdams (mcadams1994) said : 		#5

An update: the problem has been solved. I have had to reinstall pythia, as lhapdf6 has been installed after the initial pythia installation.