diff --git a/corsika/detail/modules/thinning/EMThinning.inl b/corsika/detail/modules/thinning/EMThinning.inl index 1a9982f7974ea4196f0f0d19126714a321791209..daef9d1c9b7f7b2d693d1ba15c3bbe270e462133 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);