add properties of first interaction to output
The differences between any two showers are mostly due to the characteristics of the first interaction. For debugging and for air shower phenomenology it would be useful to have the properties of the first interaction stored in the output. This could be either implemented as a separate process, a wrapper to the interaction model (a la InteractionCounter
) or a writer
for the interaction process. Or some other ingenious way...