From ea615e57bfd9a9196ac2003fefa12bf72d52cc3e Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Sat, 23 Jan 2021 14:39:43 +0100
Subject: [PATCH] fixed energy deposit profile
---
corsika/detail/modules/energy_loss/BetheBlochPDG.inl | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
index fbde1729c..776531540 100644
--- a/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
+++ b/corsika/detail/modules/energy_loss/BetheBlochPDG.inl
@@ -201,9 +201,9 @@ namespace corsika {
GrammageType const grammageEnd = shower_axis_.getProjectedX(vTrack.getPosition(1));
const auto deltaX = grammageEnd - grammageStart;
- int binStart = grammageStart / dX_;
+ int const binStart = grammageStart / dX_;
if (binStart < 0) return;
- int binEnd = grammageEnd / dX_;
+ int const binEnd = grammageEnd / dX_;
if (binEnd > int(profile_.size() - 1)) return;
if (deltaX < dX_threshold_) return;
@@ -226,7 +226,7 @@ namespace corsika {
} else {
fill(binStart, ((1 + binStart) * dX_ - grammageStart) / deltaX);
fill(binEnd, (grammageEnd - binEnd * dX_) / deltaX);
- for (int bin = binStart + 1; bin < binEnd; ++bin) { fill(bin, 1); }
+ for (int bin = binStart + 1; bin < binEnd; ++bin) { fill(bin, dX_ / deltaX); }
}
CORSIKA_LOG_DEBUG("total energy added to histogram: {} ", energyCount);
--
GitLab