Dear Micha,
for the energy balance you should use the current (USRBDX, WHAT(1) =1) and not the fluence (USRBDX, WHAT(1) =101).
The latter is weighted by 1/cos(theta).
Hope this could help.
Best luigi
On 2 Jun 2019, at 07:10, micha weiss <michaweiss89_at_gmail.com<mailto:michaweiss89_at_gmail.com>> wrote:
Thanks for the replay Luigi,
I run this minimal inp file with different gamma energies: 0.1, 0.5, 5, 10 MeV.
Thanks again
Micha Weiss
On Fri, May 31, 2019 at 12:55 PM Luigi Salvatore Esposito <luigi.salvatore.esposito_at_cern.ch<mailto:luigi.salvatore.esposito_at_cern.ch>> wrote:
Dear Micha,
could you possibly send your input?
Thanks, luigi
On 31 May 2019, at 09:11, micha weiss <michaweiss89_at_gmail.com<mailto:michaweiss89_at_gmail.com>> wrote:
Hey FLUKA experts..
I am trying to compare dose-EM in region binning (path length estimator) to the difference of energy of outgoing to incoming gamma and E+E- estimated by USRBDX (surface estimator).
I have a 5^3 cm^3 LIF crystal and I use a 0.1,0.5,5,10 MeV gamma source.
I take the current (DI/DE) passing the surface multiply by the width of the bin and multyply by the energy per bin (DI/DE)*DE*E the I sum the bins. My understanding that this is the energy passing the surface. I then take the difference of the outgoing to the incoming gamma and E+-E+ energy to the LIF crystal and divide it by density of LIF. This is the same as dose-EM.
For dos-EM per region in USRBIN I just divide by the volume to get correct units.
When I divide what I get for dose-EM from USRBIN by USRBDX they are almost identical but I get that there is a deviation that depends on the gamma source energy.
(dose-EM)_USRBDX/(dose-EM)_USRBIN = 0.90 for 0.1 MeV gamaa
(dose-EM)_USRBDX/(dose-EM)_USRBIN = 0.99 for 0.5 MeV gamaa
(dose-EM)_USRBDX/(dose-EM)_USRBIN = 1.05 for 5 MeV gamaa
(dose-EM)_USRBDX/(dose-EM)_USRBIN = 1.15 for 10 MeV gamaa
What am I doing incorrectly?
Thank you for the help.
[
https://ssl.gstatic.com/ui/v1/icons/mail/images/cleardot.gif]
<LiFDoseCheck.inp>
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Mon Jun 03 2019 - 11:03:48 CEST
