From e7e34470abd006031fcd29f9168c40d3b98b3a54 Mon Sep 17 00:00:00 2001
From: ralfulrich <ralf.ulrich@kit.edu>
Date: Thu, 9 Dec 2021 20:46:26 +0100
Subject: [PATCH] there is no momentum update in bethe-bloch, just kinetic
energy.
---
corsika/detail/modules/energy_loss/BetheBlochPDG.inl | 10 +---------
1 file changed, 1 insertion(+), 9 deletions(-)
diff --git a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
index c840de798..066997364 100644
--- a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
+++ b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
@@ -163,8 +163,7 @@ namespace corsika {
auto Enew = E + dE;
CORSIKA_LOG_TRACE("EnergyLoss dE={} MeV, E={} GeV, Ekin={} GeV, Enew={} GeV",
dE / 1_MeV, E / 1_GeV, Ekin / 1_GeV, Enew / 1_GeV);
- p.setEnergy(Enew);
- updateMomentum(p, Enew);
+ p.setEnergy(Enew); // kinetic energy on stack
fillProfile(t, dE);
return ProcessReturn::Ok;
}
@@ -194,13 +193,6 @@ namespace corsika {
vTrack, maxGrammage);
}
- template <typename TParticle>
- inline void BetheBlochPDG::updateMomentum(TParticle& vP, HEPEnergyType Enew) {
- HEPMomentumType Pnew = elab2plab(Enew, vP.getMass());
- auto pnew = vP.getMomentum();
- vP.setMomentum(pnew * Pnew / pnew.getNorm());
- }
-
template <typename TTrajectory>
inline void BetheBlochPDG::fillProfile(TTrajectory const& vTrack,
const HEPEnergyType dE) {
--
GitLab