IAP GITLAB

Skip to content
Snippets Groups Projects
Commit 43dacbfd authored by Felix Riehn's avatar Felix Riehn
Browse files

removed CoM energy from energy cut

parent 11419201
No related branches found
No related tags found
1 merge request!115Resolve "create ParticleCut process"
...@@ -27,13 +27,11 @@ namespace corsika::process { ...@@ -27,13 +27,11 @@ namespace corsika::process {
auto const energyLab = vP.GetEnergy(); auto const energyLab = vP.GetEnergy();
// nuclei // nuclei
if (vP.GetPID() == particles::Code::Nucleus) { if (vP.GetPID() == particles::Code::Nucleus) {
// calculate energy per nucleon
auto const ElabNuc = energyLab / vP.GetNuclearA(); auto const ElabNuc = energyLab / vP.GetNuclearA();
auto const EcmNN = sqrt(2. * units::constants::nucleonMass * ElabNuc); return (ElabNuc < fECut);
return (ElabNuc < fECut || EcmNN < 10_GeV);
} else { } else {
// TODO: center-of-mass energy hard coded return (energyLab < fECut);
const HEPEnergyType Ecm = sqrt(2. * units::constants::nucleonMass * energyLab);
return (energyLab < fECut || Ecm < 10_GeV);
} }
} }
......
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