From fec782a4dfc1b0ce2c6ec228227a9feec497d4da Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@tu-dortmund.de>
Date: Fri, 24 Mar 2023 21:15:53 +0100
Subject: [PATCH] update weight correctly
---
corsika/detail/modules/thinning/EMThinning.inl | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/corsika/detail/modules/thinning/EMThinning.inl b/corsika/detail/modules/thinning/EMThinning.inl
index 1a9982f79..daef9d1c9 100644
--- a/corsika/detail/modules/thinning/EMThinning.inl
+++ b/corsika/detail/modules/thinning/EMThinning.inl
@@ -57,10 +57,10 @@ namespace corsika {
if (w1 <= maxWeightFactor && w2 <= maxWeightFactor) { // apply Hillas thinning
if (uniform_(rng_) <= p1) { // keep 1st with probability p1
particle2.setWeight(0);
- particle1.setWeight(w1 * parentWeight);
+ particle1.setWeight(w1);
} else { // keep 2nd
particle1.setWeight(0);
- particle2.setWeight(w2 * parentWeight);
+ particle2.setWeight(w2);
}
} else { // weight limitation kicks in, do statistical thinning
double const w1prime = std::min(w1, maxWeightFactor);
--
GitLab