From 965c8b36a2100039033920807fa4be8b56892de0 Mon Sep 17 00:00:00 2001
From: Maximilian Reininghaus <maximilian.reininghaus@kit.edu>
Date: Wed, 7 Dec 2022 17:01:24 +0100
Subject: [PATCH] a little bit of documentation
---
corsika/modules/thinning/EMThinning.hpp | 12 ++++++++++--
1 file changed, 10 insertions(+), 2 deletions(-)
diff --git a/corsika/modules/thinning/EMThinning.hpp b/corsika/modules/thinning/EMThinning.hpp
index b414be492..d7c580b90 100644
--- a/corsika/modules/thinning/EMThinning.hpp
+++ b/corsika/modules/thinning/EMThinning.hpp
@@ -16,7 +16,9 @@
namespace corsika {
- //! This process implements thinning for EM splitting processes (1 -> 2)
+ /**
+ * This process implements thinning for EM splitting processes (1 -> 2).
+ */
class EMThinning : public SecondariesProcess<EMThinning> {
public:
@@ -32,7 +34,13 @@ namespace corsika {
* Apply thinning to secondaries. Only EM primaries with two EM secondaries are
* considered.
*
- * @tparam TStackView
+ * If the maximum weight is still out of reach, Hillas thinning is applied (i.e.
+ * one of the two secondaries is kept, the other one discarded). If the acceptance
+ * probabilities would lead to a weight factor exceeding the maximum weight, we resort
+ * to statistical thinning (i.e. the secondaries are kept/discared randomly each on is
+ * own). In that case, acceptance probabilities can be assigned without constraints
+ * (sum does not need to be 1) and we increase the acceptance probability such that
+ * the maximum weight is not exceeded.
*/
template <typename TStackView>
void doSecondaries(TStackView&);
--
GitLab