Create heavy ion objects in HepMC output
The following runcard simulating UPC with lead beams
EVENTS: 1k
EVENT_GENERATION_MODE: W
MI_HANDLER: None
ALPHAS(MZ): 0.118
ORDER_ALPHAS: 2
ALPHAS: {USE_PDF: 0}
SCALES: METS
SCALE_VARIATIONS: None
BEAMS: [1000822070, 1000822070]
BEAM_ENERGIES: 1500
BEAM_SPECTRA: EPA
SELECTORS:
- NJetFinder:
N: 2
ETMin: 2.5
R: 1.0
Exp: 1
EtaMax: 2.5
MassMax: 5.0
- NJetFinder:
N: 1
ETMin: 4.5
R: 1.0
Exp: 1
EtaMax: 2.5
MassMax: 5.0
ANALYSIS: Rivet
RIVET:
--ignore-beams: 1
--analyses: ATLAS_PBPB_CENTRALITY
USE_HEPMC_SHORT: 1
PDF_LIBRARY: SASGSherpa
MPI_PDF_LIBRARY: SASGSherpa
PDF_SET: SAS1M
MPI_PDF_SET: SAS1M
PROCESSES:
- 93 93 -> 93 93:
Order: {QCD: 2, EW: 0}
does not produce any results in Rivet due to the following warning:
Rivet.Projection.HepMCHeavyIon: WARN Could not find the HepMC HeavyIon object
i.e. the heavy ions don't seem to properly be saved in HepMC.