AW: [fluka-discuss]: track averaged LET in a volume

From: Anna Ferrari <a.ferrari_at_hzdr.de>
Date: Fri, 12 Dec 2014 20:15:32 +0100

Dear Brandon,
Sender: owner-fluka-discuss_at_smtp2.mi.infn.it



to obtain the LET in a material you can use the FLUKA function GETLET.

The calling arguments are:



     DOUBLE PRECISION FUNCTION GETLET ( IJ, EKIN, PLA, TDELTA, MATLET )



where:

IJ =3D particle index

EKIN =3D particle kinetic energy (GeV)

PLA =3D particle momentum (GeV/c)

TDELTA =3D maximum secondary electron energy (GeV) (unrestricted if =3D< 0)=


MATLET =3D material index for which LET is requested



You can call this function from a user routine. Since you would like to cal=
culate the track averaged LET in a volume,

I would suggest that you define a USRTRACK for that volume, and use the ro=
utine FLUSCW to multiply the score by GETLET



Hope it helps,

Regards,

Anna



-----------------------------------------------------

Anna Ferrari

Institute of Radiation Physics

Helmholtz-Zentrum Dresden-Rossendorf

PF 510119 - 01314 Dresden (Germany)



Von: owner-fluka-discuss_at_mi.infn.it [mailto:owner-fluka-discuss_at_mi.infn.it]=
  Im Auftrag von REDDELL, BRANDON D. (JSC-EV511)
Gesendet: Mittwoch, 3. Dezember 2014 16:48
An: Vittorio Boccone
Cc: fluka-discuss_at_fluka.org
Betreff: RE: [fluka-discuss]: track averaged LET in a volume



Ok..



One last question: can I get access to the method/routine that FLUKA uses t=
o compute DOSEQLET?



Thanks again for your help,



Brandon



From: Vittorio Boccone [mailto:dr.vittorio.boccone_at_ieee.org]
Sent: Wednesday, December 03, 2014 8:25 AM
To: REDDELL, BRANDON D. (JSC-EV511)
Cc: fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: track averaged LET in a volume



I'm afraid yes.

V.



On Wed, Dec 3, 2014 at 3:02 PM, REDDELL, BRANDON D. (JSC-EV511) <brandon.d.=
reddell_at_nasa.gov> wrote:

Yes, that is what I am trying to get. After writing my last response, I was=
  thinking this may be very complicated if I use something like mgdraw.



Is this the user routine you are referring to or is there a more efficient =
one?



Thanks for your help!


On Dec 3, 2014, at 6:30 AM, Vittorio Boccone <dr.vittorio.boccone_at_ieee.org>=
  wrote:

Hi Brandon,

  them I misunderstood what you want to obtain. You want actually to obtain =
the LET averaged on all the primary and secondary tracks for all the partic=
les of the shower generated by the primary event. Am I correct?



Using DOSEQLET you will get back the Q(LET) relation which we use (ICRP60) =
where LET=3D LET(inf) in water.

What you request can be done of course with an user routine.



Vittorio







On Wed, Dec 3, 2014 at 4:01 AM, REDDELL, BRANDON D. (JSC-EV511) <brandon.d.=
reddell_at_nasa.gov> wrote:

Thanks for the reply Vittorio. Yes, I have obtained the stopping power tabl=
es as you describe, but I don't think that helps me for a volume where I ha=
ve a shower of charged particles passing through a volume, each particle ha=
ving different initial energies and path lengths through the volume.



Shouldn't this calculation be available somewhere since FLUKA can provide D=
OSEQLET in a volume? If so, how can I get this LET? Could I score both DOSE=
  and DOSEQLET in the same volume, and use these values to get average IRCP6=
0 quality factor, then determine track-averaged LET from the Q(LET) equatio=
ns?



Brandon



Sent from my iPhone


On Dec 2, 2014, at 3:11 PM, Vittorio Boccone <dr.vittorio.boccone_at_ieee.org>=
  wrote:

Dear Brandon,

  the track averaged LET is in principle the stopping power.

You can get the stopping power calculated by FLUKA for each
material using the DELTARAY. (See Manual 7.1.3.7 of the pdf,
or 7.12 of the ASCII + references therein)

Best,

Vittorio



On Tue, Dec 2, 2014 at 4:30 PM, REDDELL, BRANDON D. (JSC-EV511) <brandon.d.=
reddell_at_nasa.gov> wrote:

Hello,



I know how to use USRYIELD to get LET across a surface, is there a way to g=
et track averaged LET over a volume (i.e. Delta-E/pathlength) in FLUKA?



Brandon
Received on Sat Dec 13 2014 - 00:26:47 CET