Re: Birks' correction

From: Anna Ferrari <a.ferrari_at_hzdr.de>
Date: Fri, 22 Mar 2013 18:01:53 +0100

2013 18:01:57 +0100 (CET)
X-Virus-Scanned: clamav-milter 0.97.7 at smtp1.mi.infn.it
X-Virus-Status: Clean
Sender: owner-fluka-discuss_at_mi.infn.it

Dear Joe Confort,
thank you, first of all, for pointing out a point, that maybe could be
clarified a bit more in the manual.
The use of the USERDUMP card with SDUM=UDQUENCH enables –as you know-the
quenching correction to the deposited energy in the track reconstruction. You
have the possibility to load different set of parameters (always supposed as
two-order corrections, as in the case of the Birks parameters) until a limit
of three sets. This is what you set in the card. The use of these different
sets in association with different materials is up to you (we will see how
below), but in fact in each material the code computes the quenching
corrections for each set. In the array DTQUEN the quenched values are stored
and you can select the set of parameters that you want to use. DTQUEN is not
used anymore (the authors in case can correct me on this point!) because is
there to give to the user the possibility to retrieve – via mgdraw - the
quenched deposited energy. This can be done at two different levels, along the
track steps in MGDRAW or for a "spot" deposition in ENDRAW:
1) in MGDRAW you have to sum, instead of the original energy depositions
dtrack(k)(with k=1,...mtrack), the corresponding quenched values given by the
dtquen(k,i), where i refers to the set of parameters that you want to use.

   2) in ENDRAW you have to sum, instead of the variable ‘rull’, the value of
dtquen(1,i) (here the first 1 means that we have only one "spot" deposit).

I hope this answers your doubts…
I would like to add two observations:
1) the treatment of the quenching in case of spot depositions is not trivial
to compute, because we have always to do assumptions, depending of the kind of
spot.
2) if you want to apply the quenching corrections to the binned energy, that
is scored via USRBIN or EVENTBIN, you have to use the TCQUENCH card (see the
manual). In this case the correction proceeds through a different, independent
way.

Hope it helps,
Regards,

Anna




Am Tue, 12 Mar 2013 01:25:35 -0700 (MST) schrieb Joseph Comfort
<Joseph.Comfort_at_asu.edu>:
> What is the meaning of the first, second, and third 'material' on the
> USERDUMP card with SDUM=UDQUENCH? I can not find a connection to a
> material, either in the manual or the source code.
>
> It sppears to me that the Birk's parameters are used to fill the array
> DTQUEN in MGDRAW. They can be printed out there, but I do not see where
> that array is used anywhere.
>
> Thank you,
> Joe Comfort
>
>

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Anna Ferrari
Institute of Radiation Physics
Helmholtz-Zentrum Dresden-Rossendorf e.V.
Bautzner Landstraße 400
D - 01328 Dresden (Germany)
Tel. +49 351 260 2872
Received on Fri Mar 22 2013 - 20:59:23 CET

This archive was generated by hypermail 2.3.0 : Fri Mar 22 2013 - 20:59:24 CET