Dear Felicia,
> I don't think I need to restrict the region where it deposits its energy.
Do you check that energy deposition has happened within TARGET region?
It looks like the user routine scores energy deposition in both TARGET
and BLKBODY regions, in that case all kinetic energy of the particle
will be scored as energy deposition, and the result will be always the same.
For example in Geant4 a sensitive detector is created for a specific
logical volume, like a region in FLUKA.
Best regards,
Mikhail
On 23.08.2018 17:03, Felicia Sutanto wrote:
> 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:26:41 CEST