lhapdf6 interface

Asked by Guglielmo Coloretti

Dear all,

I am using latest version of ubuntu and latest version of MG5. When I am trying to install lhapdf within MG5, it does not create the bin directory with the executable.

I downloaded and installed lhapdf manually and set its path in the configuration file for MG5 (though I had to set it twice, as you can see from the banner below). I was wondering if pythia8 is working properly and is correctly linked to lhapdf.

I have already set
export PYTHONPATH=/home/workstation/Desktop/MG5_aMC_v3_3_1/HEPTools/lhapdf6_py3/lib/python3.10/site-packages:$PYTHONPATH
export LD_LIBRARY_PATH=/home/workstation/Desktop/MG5_aMC_v3_4_1/HEPTools/lhapdf6_py3/lib:$LD_LIBRARY_PATH

Generations work if I use non-partonic initial states. As soon as I am trying to generate for instance p p > a with the SM, I obtain the following error:

***********************************************************
* *
* 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.1 2022-09-01 *
* *
* 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 /home/workstation/Desktop/MG5_aMC_v3_4_1/scratch/Cards/me5_configuration.txt
INFO: load configuration from /home/workstation/Desktop/MG5_aMC_v3_4_1/input/mg5_configuration.txt
INFO: load configuration from /home/workstation/Desktop/MG5_aMC_v3_4_1/scratch/Cards/me5_configuration.txt
Using default text editor "vi". Set another one 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 = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = ExRoot | OFF |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=======================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' 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]
>1
The following switches determine which programs are run:
/=================== Description ===================|============= values ==============|======== other options ========\
| 1. Choose the shower/hadronization program | shower = Pythia8 | OFF |
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = ExRoot | OFF |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=======================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=OFF' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>2
The following switches determine which programs are run:
/=================== Description ===================|============= values ==============|======== other options ========\
| 1. Choose the shower/hadronization program | shower = Pythia8 | OFF |
| 2. Choose the detector simulation program | detector = Delphes | OFF |
| 3. Choose an analysis package (plot/convert) | analysis = ExRoot | OFF |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=======================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=OFF' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
| 4. delphes : delphes_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, pythia8, 4, enter path, ... ][90s to answer]
>set pdlabel lhapdf
INFO: modify parameter pdlabel of the run_card.dat to lhapdf
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
| 4. delphes : delphes_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, pythia8, 4, enter path, ... ]
>
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.1190025
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
INFO: Using LHAPDF v6.5.3 interface for PDFs
Error: no DISPLAY environment variable specified
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: Compiling for process 1/1.
INFO: P1_qq_a
INFO: P1_qq_a
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h44 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 0.037s ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 0.037s ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
sum of cpu time of last step: 0 second
INFO: Effective Luminosity nan pb^-1
INFO: need to improve 0 channels
- Current estimate of cross-section: nan +- nan
    P1_qq_a
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 12h44 ]
INFO: Combining runs
sum of cpu time of last step: 0 second
INFO: finish refine
refine 10000 --treshold=0.9
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
sum of cpu time of last step: 0 second
INFO: Effective Luminosity nan pb^-1
INFO: need to improve 0 channels
- Current estimate of cross-section: nan +- nan
    P1_qq_a
INFO: Idle: 0, Running: 0, Completed: 0 [ current time: 12h44 ]
INFO: Combining runs
sum of cpu time of last step: 0 second
INFO: finish refine
INFO: Combining Events
  === Results Summary for run: run_01 tag: tag_1 ===

     Cross-section : nan +- nan pb
     Nb of events : 0

INFO: Running Systematics computation
Command "generate_events run_01" interrupted with error:
InvalidCmd : Invalid run name. Please retry
quit
INFO: storing files of previous run
INFO: Done
INFO:
more information in /home/workstation/Desktop/MG5_aMC_v3_4_1/scratch/index.html
MG5_aMC>exit

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Guglielmo Coloretti
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

The process p p > a indeed does not make sense and the crash is expected.
The issue is that you ask a single photon onshell out of two onshell quark.
Since the mass of a photon is zero, it means that you ask the two quark to also have a zero energy ...

Cheers,

Olivier

Revision history for this message
Guglielmo Coloretti (guglielmocoloretti) said :
#2

Dear Olivier,

Apologies for having posted the answer, it was quite straightforward.... but thank you very much for your kind reply!!

Best regards,
Guglielmo Coloretti