Dear Riccardo,
You may explore the RULL in comscw.f instead. Not recommended though, as
you won't get apportioning. Are you absolutely sure that USRBIN cannot
answer your needs? I don't know what 'other quantities' you need. USRBIN
does score by region (just note that the normalisation would be
different from other choices of WHAT(1)s).
Energy deposition is about the last thing to intercept with user
routines. Not a good idea.
:) mary
On 2017-02-24 17:15, dal_bello_at_stud.uni-heidelberg.de wrote:
> Dear FLUKA users,
>
> I have a question regarding the use of the entry ENDRAW in the mgdraw
> routine.
>
> My personalized mgdraw routine scores, among other quantities, the
> energy deposition in different regions of my detector (primary beam:
> gamma between 0. and 10. MeV).
> I tired to score the energy deposition accumulating the value of RULL
> in ENDRAW, however, I do not get correct value (lower values).
>
> I worked around the problem accumulating for each event the energy of
> gamma, e- and e+ entering and exiting the volumes (done in BXDRAW),
> but of course it is not a clean way to solve the problem because it
> slows the simulation and involves considerations about the rest mass
> of e+/e-, which may lead to errors in subtractions.
>
> I accumulate the values of RULL simply with:
>
> ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
> CALL GEOR2N ( MREG, NRGNAM, IERR2 )
>
> IF (NRGNAM(1:6).EQ.'SCINT0') THEN
> ERULL0 = ERULL0 + RULL
> ENDIF
>
> RETURN
>
> Where ERULL0 is a DOUBLE PRECISION variable initialized to 0 and
> re-set to 0 at the end of each event (EEDRAW).
>
> Being the values of energy deposition with RULL smaller than the
> actual ones, I suspect that RULL includes just some types of energy
> deposition.
> Is there a way to score in ENDRAW the correct energy deposition? I
> looked in previous documentation but unsuccessfully.
>
>
> Thank you,
> Best regards,
> Riccardo Dal Bello
>
> -------------------
>
>
> Attached: few energy deposition events in SCINT0 from a 4.5 MeV
> primary gamma beam
>
> 1. column: ID event
> 2. column: true energy deposition (done in BXDRAW)
> 3. column: ERULL0 value called in EEDRAW
>
> 4 0.39547E-02 0.33276E-03
> 6 0.45000E-02 0.49065E-03
> 9 0.37780E-02 0.24186E-03
> 11 0.28722E-02 0.33501E-03
> 24 0.39805E-02 0.96697E-04
> 27 0.13814E-02 0.19710E-03
> 30 0.41481E-02 0.68106E-03
> 33 0.26909E-02 0.28380E-03
> 47 0.34780E-02 0.29714E-03
> 50 0.39890E-02 0.61295E-03
> 55 0.67689E-03 0.22618E-03
> 71 0.43555E-03 0.22666E-03
> 73 0.39474E-02 0.67036E-03
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id¬c_info
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Feb 24 2017 - 13:08:36 CET