IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 89deba43 authored by Felix Riehn's avatar Felix Riehn Committed by ralfulrich
Browse files

removed fmt and reintroduce isInvisible

parent 77e941d6
No related branches found
No related tags found
1 merge request!303Resolve "advanced ParticleCut process"
......@@ -104,14 +104,18 @@ namespace corsika {
}
}
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(fmt::format("ParticleCut: checking {}, E= {} GeV, EcutTot={} GeV",
CORSIKA_LOG_DEBUG("ParticleCut: checking {}, E= {} GeV, EcutTot={} GeV",
pid, energy / 1_GeV,
(em_energy_ + inv_energy_ + energy_) / 1_GeV));
(em_energy_ + inv_energy_ + energy_) / 1_GeV);
if (doCutEm_ && is_em(pid)) {
CORSIKA_LOG_DEBUG("removing em. particle...");
em_energy_ += energy;
......
......@@ -78,6 +78,9 @@ namespace corsika {
template <typename TParticle>
bool isBelowEnergyCut(TParticle const&) const;
//! defines which particles are invisible, by default only neutrinos
bool isInvisible(Code const &) const;
private:
bool doCutEm_;
bool doCutInv_;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment