MadGraph/Pythia8 hepmc and Rivet
Hi,
I'm trying to generate a hepmc file with MadGraph and showering by Pythia8 that I can afterwards use with Rivet. However, there are two problems. As I'm not sure whether they are unrelated or not, I'll post them here together.
1)
When I'm running MG5_aMC_v2_7_2 with the following script:
generate p p > t t~
output mininmal_ex
launch
shower=Pythia8
set nevents 1000
Pythia8 is not "visibly" switched on in the options, as was the case for MG5_aMC_v2_6_7:
MG5_aMC_v2_7_2 logfile snapshot:
[...]
*******
INFO: load configuration from /lustre/
INFO: load configuration from /lustre/
INFO: load configuration from /lustre/
Using default text editor "vi". Set another one in ./input/
generate_events run_01
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program | ^[[1m detector^[[0m = ^[[01mNot Avail.^[[0m | Please install module |
| 3. Choose an analysis package (plot/convert) | ^[[1m analysis^[[0m = ^[[01mNot Avail.^[[0m | Please install module |
| 4. Decay onshell particles | ^[[1m madspin^[[0m = ^[[31mOFF^[[0m | ON|onshell |
| 5. Add weights to events for new hypp. | ^[[1m reweight^[[0m = ^[[31mOFF^[[0m | ON |
\======
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.
set nevents 1000
^[[1mINFO: modify parameter nevents of the run_card.dat to 1000 ^[[0m
Do you want to edit a card (press enter to bypass editing)?
/------
| ^[[31m1^[[0m. param : ^[[32mparam_
| ^[[31m2^[[0m. run : ^[[32mrun_
| ^[[31m3^[[0m. pythia8 : ^[[32mpythia8_
\------
you can also
[...]
MG5_aMC_v2_6_7 logfile snapshot:
[...]
*******
INFO: load configuration from /lustre/
INFO: load configuration from /lustre/
INFO: load configuration from /lustre/
Using default text editor "vi". Set another one in ./input/
generate_events run_01
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program | ^[[1m detector^[[0m = ^[[01mNot Avail.^[[0m | Please install module |
| 3. Choose an analysis package (plot/convert) | ^[[1m analysis^[[0m = ^[[01mNot Avail.^[[0m | Please install module |
| 4. Decay onshell particles | ^[[1m madspin^[[0m = ^[[31mOFF^[[0m | ON|onshell |
| 5. Add weights to events for new hypp. | ^[[1m reweight^[[0m = ^[[31mOFF^[[0m | ON |
\======
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.
The following switches determine which programs are run:
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program | ^[[1m detector^[[0m = ^[[31mOFF^[[0m | Please install module |
| 3. Choose an analysis package (plot/convert) | ^[[1m analysis^[[0m = ^[[31mOFF^[[0m | Please install module |
| 4. Decay onshell particles | ^[[1m madspin^[[0m = ^[[31mOFF^[[0m | ON|onshell |
| 5. Add weights to events for new hypp. | ^[[1m reweight^[[0m = ^[[31mOFF^[[0m | ON |
\======
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.
set nevents 1000
^[[1mINFO: modify parameter nevents of the run_card.dat to 1000 ^[[0m
Do you want to edit a card (press enter to bypass editing)?
/------
| ^[[31m1^[[0m. param : ^[[32mparam_
| ^[[31m2^[[0m. run : ^[[32mrun_
| ^[[31m3^[[0m. pythia8 : ^[[32mpythia8_
\------
you can also
[...]
Is that expected and something that was changed from MG5 2.6.7 to MG5 2.7.2? As far as I can tell, the actual showering with Pythia8 *does* happen in both cases.
2)
Depending on the MadGraph and Rivet version, interplay between those two does or does not work. What I found out so far is that
MG5 2.6.7 + Rivet 3.1.0: does not work (a)
MG5 2.6.7 + Rivet 3.1.1: works
MG5 2.7.2 + Rivet 3.1.0: does not work (b)
MG5 2.7.2 + Rivet 3.1.1: does not work (c)
Hereby, "works" means that after running the short MG5 script from above events can be processed as expected while "does not work" means that I end up with one of the following error messages:
case (a)
» rivet -a ATLAS_2016_I1469071 minimal_
Rivet 3.1.0 running on machine wgs32.zeuthen.
Could not identify nominal weight. Please check your HEPMC file.
»
case (b)
» rivet -a ATLAS_2016_I1469071 minimal_
Rivet 3.1.0 running on machine wgs32.zeuthen.
streaming input: end key not found setting badbit.
Could not read from 'minimal_
»
case (c)
» rivet -a ATLAS_2016_I1469071 minimal_
Rivet 3.1.1 running on machine wgs32.zeuthen.
streaming input: end key not found setting badbit.
Could not read from 'minimal_
»
Do you any idea why that is?
Thanks,
Martin
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply: