Re: [fluka-discuss]: mgdraw user routine, entry ENDRAW, ICODE 22, RULL issue

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Tue, 7 Feb 2017 21:47:25 +0100

Dear Andrea,

you are overlooking the statistical weight of your photon that in the
printed line appears to be 1.3E-5 and transforms the residual kinetic
energy of 232 keV (you wrote MeV instead of GeV) into the indicated RULL
value.
The mentioned statistical weight seems to indicate that you adopted a
quite strong biasing, which may have not innocent consequences and anyway
prevents a clean event-by-event analysis. The interest in the latter can
justify the use of mgdraw, otherwise local energy deposition can be safely
obtained by USRBIN (also per region).
Moreover, your photon transport threshold (>232 keV) looks pretty high.
[iarg=2 means that the photon energy is lower than both the transport AND
production threshold, which is obvious if the two - as expected - are the
same].

Best regards

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Tue, 7 Feb 2017, Andrea Celentano wrote:

> Dear FLUKA users,
> I am experiencing the following issue when using the MGDRAW user routine, in
> particular the entry ENDRAW, that I want to use to score local energy
> deposition.
> In my code, I simply dump to a file the data from the ENDRAW call, each time
> the region name is equal to those of my interest:
>
>    CALL GEOR2N(MREG,REGNAME,IREGERR)
>     
>       if ((REGNAME.eq.'CRS1').OR.(REGNAME.eq.'CRS2').
>      &   or.(REGNAME.eq.'CRS3')) then
>          write (IODRAW,*) NCASE,REGNAME,
>      &         JTRACK,WTRACK,WSCRNG,
>      &         ETRACK,ICODE,'1',
>      &         RULL
>       endif
>
> After the simulation, If I open the file, I see lines as the one that
> follows:
>
>  65 CRS1               7   1.3337934658997385E-005  
> 1.3337934658997385E-005   2.3205639141605591E-004          22 1  
> 3.0951529859100757E-009
>
>
>
> Event n. 65, region is CRS1, RULL (local energy deposition) is 3.09515E-9
> GeV, ICODE is 22, i.e.  below threshold, iarg=2 (what is iarg, btw?).
>
>
> However, I do not understand why this RULL value:
>
> * If I can use JTRACK,WTRACK,ETRACK from (TRACK), then the energy of this
> photon is 2.32E-004 MeV: why isn't RULL equal to this value? I expect that,
> if the particle is killed being below threshold, all the kinetic energy is
> deposited locally.
>
> * If, instead, I can't use JTRACK,WTRACK,ETRACK from (TRACK), how can I
> access to the information of the particle that is depositing energy?
>
>
> Actually, all the lines where ICODE==22 have this issue.
>
>
> Thanks in advance
>
> Andrea Celentano
>
>

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Feb 07 2017 - 23:31:05 CET

This archive was generated by hypermail 2.3.0 : Tue Feb 07 2017 - 23:31:06 CET