Invalid configuration used in corsika example
The lines currently used in the corsika example initialize the modules describing the electromagnetic component in this way:
HEPEnergyType heHadronModelThreshold = 63.1_GeV;
corsika::proposal::Interaction emCascade(env, sibyll.getHadronInteractionModel(),
heHadronModelThreshold);
// NOT available for PROPOSAL due to interface trouble:
// InteractionCounter emCascadeCounted(emCascade);
// corsika::proposal::ContinuousProcess<SubWriter<decltype(dEdX)>> emContinuous(env);
BetheBlochPDG<SubWriter<decltype(dEdX)>> emContinuous{dEdX};
This means that for continuous processes, BetheBloch is used, while for stochastic processes, PROPOSAL is used. This makes no sense, as BetheBloch is (from my understanding) supposed to take all energy losses into account on its own. When you use proposal::Interaction in addition, this module while do its own energy losses, which means that energy losses will be taken into account twice. Results from this will be completely unphysical.
For the corsika example, one should either use proposal::Interaction in combination with proposal::ContinuousProcess, or only use BetheBloch.