diff --git a/corsika/detail/modules/ParticleCut.inl b/corsika/detail/modules/ParticleCut.inl index a73d0c1c10c585e775ede3f36114ddfd88047908..551234fc9b9579e9b312770623f78fd05a59c06d 100644 --- a/corsika/detail/modules/ParticleCut.inl +++ b/corsika/detail/modules/ParticleCut.inl @@ -104,18 +104,15 @@ namespace corsika { } } - bool ParticleCut::isInvisible(Code const &vCode) const { - return is_neutrino(vCode); - } - + bool ParticleCut::isInvisible(Code const& vCode) const { return is_neutrino(vCode); } + template <typename TParticle> bool ParticleCut::checkCutParticle(const TParticle& particle) { const Code pid = particle.getPID(); HEPEnergyType energy = particle.getEnergy(); - CORSIKA_LOG_DEBUG("ParticleCut: checking {}, E= {} GeV, EcutTot={} GeV", - pid, energy / 1_GeV, - (em_energy_ + inv_energy_ + energy_) / 1_GeV); + CORSIKA_LOG_DEBUG("ParticleCut: checking {}, E= {} GeV, EcutTot={} GeV", pid, + energy / 1_GeV, (em_energy_ + inv_energy_ + energy_) / 1_GeV); if (doCutEm_ && is_em(pid)) { CORSIKA_LOG_DEBUG("removing em. particle..."); em_energy_ += energy; diff --git a/corsika/modules/ParticleCut.hpp b/corsika/modules/ParticleCut.hpp index 61f832898616e4bb8fa98dd8d972b3827d20d1bf..c9771c9367a505db62c8c57c23a240a0e7edcdfa 100644 --- a/corsika/modules/ParticleCut.hpp +++ b/corsika/modules/ParticleCut.hpp @@ -79,8 +79,8 @@ namespace corsika { bool isBelowEnergyCut(TParticle const&) const; //! defines which particles are invisible, by default only neutrinos - bool isInvisible(Code const &) const; - + bool isInvisible(Code const&) const; + private: bool doCutEm_; bool doCutInv_;