How to use a Pythia afterburner for O2

Dear O2 experts,

For my studies I need to add a pythia afterburner to the simulation.
This afterburner will calculate the amount of Internal bremsstrahlung created when a collision takes places. The amount of photons produced (and in particular soft photons) is predicted via the Gellman and Low’s theorem.

These photons will be added to the event (since they are soft, they will have a negligible impact on the production processes, so an afterburner seems to be the best option).

How do I integrate such an afterburner?
Do I make it like a user hook and place it in AliceO2/Generators/share/egconfig at dev · AliceO2Group/AliceO2 · GitHub ? Or is there a better way I can tackle this?

Kind regards,

Pinging @rpreghen as expert for generation.

hi @Cas ,
the user hooks approach is IMO not the right one.
I would rather implement a custom event generator that inherits from Pythia.

an starting point can be this one

an approach can be that you let Pythia generate the event like in the Bool_t generateEvent() override { function. You analyse the event when generated and add particles to the Pythia stack as per your needs.


1 Like