From: Sebastien WURTH (firstname.lastname@example.org)
Date: Wed Mar 05 2008 - 09:21:30 CET
X-Accept-Language: fr, en
To: hamide jalali <email@example.com>
Subject: Re: deq99c.f & dose equivalent
References: <200803040821.m248L7As029152@smtp2.mi.infn.it>=09 <47CD0CCB.409=
Content-Type: text/plain; charset=3DISO-8859-1; format=3Dflowed
I reply to the all discussion group, I hope you don't mind, some people=20
could just have better answers than mine.
I'm not sure I really understood your question. In shielding problems=20
you want to calculate dose equivalent for alpha particles ?
The dose should be huge yes, but on very short distances !.
Then, if you want to calculate dose from secondaries induced by=20
reactions of these alpha, you have the case we discussed yesterday.
For heavy ions, I am no expert.
But, in general, you can always use FLUSCW routine without deq99c, and=20
convert your fluence into dose with your own conversion factors,
I did it before I started to use Stefan Roesler's routine.
For instance for photons (monoenergetic example here), you can find this=20
relation : D' =3D F'*E*=B5en/rh=F4, where :
D' is the absorbed dose (=3D equivalent dose for photons) rate in Gy per=20
second, E the energy in J and =B5en/rh=F4 the absorption coefficient in m=
F' is the fluence rate in particles per m=B2 per second : F' =3D A* Nu /=20
4*pi*d=B2 where A is the activity in Bq, Nu the branching ratio of the=20
energy considered, d the distance considered in meters. This is given by=20
FLUKA, in fact, FLUKA gives F'/A in particules/cm=B2/primary.
Of course, some work on the units must be done to use MeV or GeV for=20
energy, cm=B2/g for absorption coefficient, etc.
You calculate your factor to be applied to fluence during run time, you=20
modify FLUSCW routine as explained in the manual (user routines=20
section), you compile, link with your executable and run it. You could=20
use such relations for alpha or heavy ions...
Hope it helps.
hamide jalali wrote :
> Dear Sebastien
> I am very thankful to you for your reply. Now I can be sure that my=20
> calculation is correct.
> I use deq99c.f for the same purpose as you, I mean shielding problem.=20
> In a part of my problem I should calculate dose equivalent when the=20
> source is alpha particles and also heavy ions. In this case FLUKA=20
> don't support my calculations because it hasn't any routine for this=20
> kind of source so my calculation will remain half-done.
> Have you ever faced this problem? What is your suggestion when source=20
> is heavy nuclei?
> Any help would be appreciated in advance.
> Kind Regards
> On Tue, Mar 4, 2008 at 12:48 AM, Sebastien WURTH <firstname.lastname@example.org=20
> <mailto:email@example.com>> wrote:
> Hello Hamide,
> The code to include in your detector's SDUM is explained at the
> beginning of the deq routine (maybe AMB74 here).
> Using particle code 201 to score the dose equivalent for all particle=
> included in the routine is the right thing to do, I used it very
> You can calculate separately and add up, if everything goes OK, it
> should be the same.
> Generalised particle is not an option in my opinion. We are talking
> about conversion from fluence to dose equivalent here (using FLUSCW
> routine), not from energy deposited to dose (for this one, COMSCW
> routine is the suitable one).
> You can use the three kinds of detectors to calculate dose, I usually
> use USRBIN to plot spatial distributions for shielding purposes, it
> depends of what you need.
> Hope it helps.
> hamide jalali wrote :
> >Dear FLUKA users
> >I'd like to get ambient dose equivalent in ICRU sphere by using
> >deq99c.f while the source is Proton with a spectrum of energy
> which is
> >defined by source.f.
> >I want to know how can get total dose equivalent from all secondarie=
> >which are included deq99c.f ? Is it enough to set the particle type
> >to be scored=3D 201 (all particle) or I should calculate dose
> >for each secondary separately, then add them up? How about
> >particle 208 (energy)?
> >And the other question, among USRBIN,USRBDX and USRTRACK which one i=
> >recommended to calculate dose equivalent?
> >Best Regards
This archive was generated by hypermail 2.1.6 : Wed Mar 05 2008 - 12:08:28 CET