<div dir="ltr"><div><div><div><div>Hi,<br></div>I have implemented this. I have not changed the method because there is a deterministic part. I see the problem as is: you have a measure which varies stochastically in each pixel (because there is a stochastic part in the computation) and I estimate the uncertainy using the history-by-history method.<br></div>So N is the number of particles during the run (mNumberOfEventsInRun) and X is the measure (Compton scatter, Rayleigh scatter, total scatter etc.) in one pixel. This is computed in <a href="https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/include/GateFixedForcedDetectionFunctors.hh#L1016">this function</a>. Note that all the squared and sums are computed with ITK filters in the code, look for <span class="gmail-pl-smi">mIsSecondarySquaredImageEnabled in <a href="https://github.com/OpenGATE/Gate/blob/develop/source/digits_hits/src/GateFixedForcedDetectionActor.cc">GateFixedForcedDetectionActor.cc</a>. The code is not very easy to read but that's just mimicking the algorithm in the reference you sent for complete projections (all pixels).<br></span></div>I would suggest to use the batch approach to validate it if you have any doubt.<br></div>Simon<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 24, 2017 at 11:59 AM, Triltsch, Nicolas <span dir="ltr"><<a href="mailto:nicolas.triltsch@tum.de" target="_blank">nicolas.triltsch@tum.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Gate fellows,<br>
<br>
My question aims at the history-by-history method which is used to calculate statistical uncertainties for the scatter images (Compton and Rayleigh) in the fixed forced detection actor. I have troubles to understand how the method  is applied within this algorithm which combines statistical calculation of scatter points within the object using classical Monte Carlo (MC) simulations (stage 1) with deterministic calculations of the probability to reach a certain detector pixel (stage 2).<br>
<br>
What are the parameters N and X (see the original paper I have appended) in our case which are required to calculate the variance for the new method?<br>
<br>
Uncertainty should only play a role in the first stage of the simulation where classical MC is used, however, in the second stage purely deterministic calculations are used. How is this considered within the uncertainty method?<br>
<br>
I hope that someone can help me to understand this in more detail!<br>
<br>
Thanks<br>
<br>
Nicolas Triltsch<br>
<br>
-- <br>
B.Sc. Nicolas Triltsch<br>
Masterand<br>
<br>
Technische Universität München<br>
Physik-Department<br>
Lehrstuhl für Biomedizinische Physik E17<br>
<br>
James-Franck-Straße 1<br>
85748 Garching b. München<br>
<br>
Tel: <a href="tel:%2B49%2089%20289%2012591" value="+498928912591" target="_blank">+49 89 289 12591</a><br>
<br>
<a href="mailto:nicolas.triltsch@tum.de" target="_blank">nicolas.triltsch@tum.de</a><br>
<a href="http://www.e17.ph.tum.de" rel="noreferrer" target="_blank">www.e17.ph.tum.de</a><br>
<br>
<br>______________________________<wbr>_________________<br>
Gate-users mailing list<br>
<a href="mailto:Gate-users@lists.opengatecollaboration.org">Gate-users@lists.<wbr>opengatecollaboration.org</a><br>
<a href="http://lists.opengatecollaboration.org/mailman/listinfo/gate-users" rel="noreferrer" target="_blank">http://lists.<wbr>opengatecollaboration.org/<wbr>mailman/listinfo/gate-users</a><br></blockquote></div><br></div>