EM continuous process not reporting final results at info-level
In standard end-user configuration, ie eith debug level set to info
, the example program em_shower
(the same applies to all other examples) for a 1TeV shower returns:
[corsika:info (Cascade.hpp:83)] This is CORSIKA 8.0.0.0
[corsika:info (Cascade.hpp:85)] Tracking algorithm: LeapFrog-curved (version 1.0.0)
[corsika:info (Cascade.hpp:88)] Stack - with full cascade HISTORY.
[corsika:info (ParticleCut.inl:193)]
******************************
kinetic energy removed by cut of electromagnetic (GeV): 0 (number: 0)
kinetic energy removed by cut of invisible (GeV): 0 (number: 0)
kinetic energy removed by kinetic energy cut (GeV): 195.0162809527017 (number: 65649)
kinetic energy removed by time cut (GeV): 0
******************************
[corsika:info (ObservationPlane.inl:95)] ******************************
ObservationPlane:
energy an ground (GeV) : 4.29224044041467
no. of particles at ground : 132
******************************
total cut energy (GeV): 976.563
relative difference (%): -2.34369
This is very confusing output as most of the energy seems to go missing.. the reason is that the continuous process does not report the total losses. One fix would be to replace CORSIKA_LOG_DEBUG
with CORSIKA_LOG_INFO
in ContinuousProcess::showResults
of the PROPOSAL interface.