Dear FLUKA Experts,
Sorry for resurrecting this topic (Uncertainty: Overall Value) but after going through the wiki page, I found that I could implement the general formula for uncertainty propagation. That is,
to answer your first question, when scoring ENERGY with a Region mesh, the USRBIN result is reported in GeV/primary and the error, as written in the .lis file, is a percentage error.
Looking at your DETECT file, the first bin has the highest counts comparing it to the others: excluding it from the calcualtion you are performing causes the deviation from the average energy deposited computed with USRBIN .
Simply including the first bin counts in the average, I get 2.48E-07 GeV/primary, closer to the USRBIN result.
In order to obtain an even similar result, you should set more appropriately the DETECT bins since, as it is now, the majority of the bins give 0.
For example, with the following DETECT card, I get, using your formulas, an average energy deposited of 1.38E-07 GeV/primary with uncertainty 6.9E-01 %
I will let you set it accordingly to you needs but in general the USRBIN gives you already the average energy deposited in the region that you are interested with its own error without the need of any post processing.
Finally to point out that, if I understand correctly the structure of your DETECT file, it looks like that you are not considering Ei = LBi + 0.5*(UPi-LBi), rather Ei = UBi.
I hope this helps you in solving your doubts!
Giuseppe
Dear Giuseppe,
Thank you very much for your reply. Is there any reason you chose 22.78 keV instead of 500 keV as your E_max for the DETECT card? Or is this based on the readings from my “…tab.lis” file? By checking, I see that the counts reduce to zero beyond this value.
as you correctly understood, I have chosen the Emax for DETECT simply looking at which was the last bin with significant counts in the .lis file, no physic reason behind. From this come also my statement to set this value depending on your needs.
yes, USRBIN (with selected mesh Region) gives the average energy deposition (in GeV) per primary.
It is done summing, for each primary of a cycle, the energy deposited in the region and finally dividing it by the number of primaries per cycle (ni). If a primary gives 0 as energy deposited, it is taken into account in the sum.