Skip to content

Commit 5858cfd

Browse files
rewording after re-reading
1 parent 1372eb2 commit 5858cfd

File tree

1 file changed

+14
-7
lines changed

1 file changed

+14
-7
lines changed

src/physics/stats/electrons-on-target.md

Lines changed: 14 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -49,15 +49,22 @@ where
4949
- \\(\sum w\\) is the sum of the event weights of those same events
5050
- \\(N_\text{attempt}\\) is the number of events that were attempted when simulating
5151

52+
Another way to phrase this equation (one that is more appropriate for the statistical analysis below) is
53+
\\[
54+
N_\text{EoT}^\text{equiv} = \frac{N_\mathrm{attempt}}{\langle{w}\rangle}
55+
\\]
56+
where \\(\langle{w}\rangle\\) is the mean event weight over the sampled events.
57+
Dividing by this mean event weight is like multiplying by an "effective biasing factor" since
58+
the event weights being with a value of \\(1/B\\).
59+
5260
~~~admonish note title="Finding Number of Attempted Events"
5361
Currently, the number of attempted events is stored as the `num_tries_` member of
5462
the `RunHeader`
5563
56-
Samples created with ldmx-sw verions newer than v3.3.4 (>= v3.3.5) have an update
57-
to the processing framework to store this information more directly
58-
(in the `numTries_` member of the RunHeader or <= v4.4.7 and the `num_tries_` member
59-
for newer).
60-
Samples created with ldmx-sw versions older than v3.1.12 (<= v3.1.11) have access
64+
Samples created with ldmx-sw verions >= v3.3.5 and <= v4.4.7 store this number
65+
in the `numTries_` member of the RunHeader.
66+
67+
Samples created with ldmx-sw versions <= v3.1.11 have access
6168
to the "Events Began" field of the `intParameters_` member of the RunHeader.
6269
6370
The easiest way to know for certain the number of tries is to just set the maximum
@@ -110,7 +117,7 @@ In the thin-target regime (where nothing happens to a biased particle besides th
110117
biased process), this equation reduces to the simpler \\(B N_\mathrm{attempt}\\) used in other
111118
analyses since biased tracks in Geant4 begin with a weight of \\(1/B\\).
112119
\\[
113-
N_\text{EoT}^\text{equiv} = \frac{N_\mathrm{sampled}}{\sum w}N_\mathrm{attempt}
120+
N_\text{EoT}^\text{equiv} = \frac{N_\mathrm{sampled}}{\sum w}N_\mathrm{attempt} = \frac{N_\mathrm{attempt}}{\langle{w}\rangle}
114121
\\]
115122

116123
## Event Yield Estimation
@@ -145,7 +152,7 @@ how the probability in the biased sample differs from one in an inclusive sample
145152
Using "i" to stand for "inclusive" and "b" to stand for "biased".
146153
There are two options that we have used in LDMX.
147154
1. \\(p_\text{b} = B p_\text{i}\\) where \\(B\\) is the biasing factor.
148-
2. \\(p_\text{b} = W p_\text{i}\\) where \\(W\\) is the ratio of the average event weights between the two samples. Since the inclusive sample has all event weights equal to one, \\(W = \sum_\text{b} w / N\\) so it represents the EoT estimate described above.
155+
2. \\(p_\text{b} = W p_\text{i}\\) where \\(W\\) is the ratio of the average event weights between the two samples. Since the inclusive sample has all event weights equal to one, \\(W = \langle{w}\rangle_b\\) so it represents the EoT estimate described above.
149156

150157
~~~admonish note title="Binomial Basics"
151158
- Binomials are valid for distributions corresponding to some number of binary yes/no questions.

0 commit comments

Comments
 (0)