Dear Felicia,
There is a subtle point you have to consider: when a particle reaches the boundary to blackhole, then the code makes a final call to ENDRAW and RULL contains the remaining energy of the particle. So this is not an energy deposition event in your target and the particle is actually considered to be in the blackhole region.
This explains why you end up always with the same energy (i.e. the initial energy).
To make sure that the call to ENDRAW is indeed an energy deposition event in your target, you have to filter events according to the region number (available as MREG in ENDRAW).
Cheers, Anton
________________________________
From: Felicia Sutanto [fsutanto_at_umich.edu]
Sent: 23 August 2018 16:03
To: Anton Lechner; pmk_at_ihep.ru
Cc: fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: NaI detector response, mgdraw, and endraw
Dear Mikhail and Anton,
Thank you for checking my simulation!
The detect card indeed shows the same result as geant, so It must be my way of scoring in mgdraw that's wrong.
But I cannot understand why my way of scoring in mgdraw is wrong...
I'd really appreciate it if you guys can tell me what's wrong with the following:
I commented out all of the write commands in mgdraw.f and
I changed the unformatted collision tape file to a formatted one (in MGDRAW, ENDRAW, and SODRAW):
OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM = 'FORMATTED' )
Then in MGDRAW, I added:
mgIn = 0
WRITE (IODRAW,*) mgIn, ( SNGL (DTRACK (I)), I = 1, MTRACK )
and in ENDRAW, I added:
enIn = 1
WRITE (IODRAW,*) enIn, SNGL (RULL)
and in SODRAW, I added:
soIn = 2;
WRITE (IODRAW, *) soIn, TKEFLK(NPFLKA)
This produces a trakfile with a list of the energy deposited (index: 0 and 1) and the energy of the source particle (index: 2):
2 6.6169999999999998E-004
0
1 6.61699974E-04
2 6.6169999999999998E-004
0
1 1.45529995E-07
0 1.07328451E-05
0 6.06724370E-06
0 1.79876715E-05
0 7.13503368E-06
0 4.00992985E-06
0 6.33815262E-06
0 5.01241266E-06
0 6.40999087E-06
0 6.91474997E-06
0 4.94164897E-06
0 3.67576922E-06
0 4.34047570E-06
0 3.34160836E-06
0 6.53870757E-06
0 1.80250277E-06
1 9.95420487E-05
0
1 4.66763682E-04
In my post-processing, I added the energies with indexes 0 and 1, and tally them every time I go to index 2.
In my simulation, I only have a spherical NaI inside a spherical blackbody. I don't think I need to restrict the region where it deposits its energy.
Sorry this may be a very trivial question, but I am really not sure what's wrong with this way of scoring...
Thank you and best regards,
FS
On Thu, Aug 23, 2018 at 8:52 AM Anton Lechner <Anton.Lechner_at_cern.ch<mailto:Anton.Lechner_at_cern.ch>> wrote:
Dear Felicia,
Indeed, as pointed out by Mikhail, it must be your way of scoring in MGDRAW or the post-processing, not the physics. Attached you can find the results I calculated with your input file using the DETECT card. The results are pretty much the same as you get with Geant4 (for 1E6 primaries, like you used in Geant4).
Cheers, Anton
________________________________________
From: owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it> [owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it>] on behalf of Mikhail Polkovnikov <Михаил Полковников> [pmk_at_ihep.ru<mailto:pmk_at_ihep.ru>]
Sent: 23 August 2018 12:06
To: Felicia Sutanto; fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: NaI detector response, mgdraw, and endraw
On 21.08.2018 22:58, Felicia Sutanto wrote:
> Dear FLUKA users,
>
> I am a newbie and I have a question about MGDRAW and ENDRAW.
>
> I am trying to simulate a very simple case (input file is attached):
> 661.7 keV photon at the center of a spherical NaI detector (radius = 5cm)
>
> The purpose of this exercise is only to familiarize myself with the
> collision tape.
> I try to tally the energy deposited by summing over the DTRACK in
> MGDRAW and RULL in the ENDRAW. The result is attached
> (fluka_nai_result.pdf).
>
> I also try to do the same simulation in Geant4 and the result is also
> attached (geant4_nai_result.pdf). It seems that the Compton continuum
> is missing in FLUKA...
>
> Am I supposed to maybe include certain physics in FLUKA?
>
> Best regards,
> FS
>
>
>
>
>
>
Dear Felicia,
You can check an energy deposition process by adding a DETECT card.
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....
DETECT 0.0 1.E-05 1.E-03 1.E-06 1. TARGETEDEP
How do you process the MGDRAW, ENDRAW and EEDRAW entries within the user
routine, and how do you process the collision tape itself? Both DETECT
and user dump work fine for me.
Best regards,
Mikhail
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Aug 23 2018 - 19:42:42 CEST