diff --git a/Processes/ParticleCut/ParticleCut.cc b/Processes/ParticleCut/ParticleCut.cc index 8a350f259ecd3267ac2af968377118a40584f8e1..28e7013d08d864471de08c742250b3ecb4357760 100644 --- a/Processes/ParticleCut/ParticleCut.cc +++ b/Processes/ParticleCut/ParticleCut.cc @@ -27,13 +27,11 @@ namespace corsika::process { auto const energyLab = vP.GetEnergy(); // nuclei if (vP.GetPID() == particles::Code::Nucleus) { + // calculate energy per nucleon auto const ElabNuc = energyLab / vP.GetNuclearA(); - auto const EcmNN = sqrt(2. * units::constants::nucleonMass * ElabNuc); - return (ElabNuc < fECut || EcmNN < 10_GeV); + return (ElabNuc < fECut); } else { - // TODO: center-of-mass energy hard coded - const HEPEnergyType Ecm = sqrt(2. * units::constants::nucleonMass * energyLab); - return (energyLab < fECut || Ecm < 10_GeV); + return (energyLab < fECut); } }