**From:** Sebastien WURTH (*wurth@ipno.in2p3.fr*)

**Date:** Wed Mar 05 2008 - 09:21:30 CET

**Previous message:**¿µçô: "fluka & cygwin"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

040910

X-Accept-Language: fr, en

MIME-Version: 1.0

To: hamide jalali <jalali.hb@gmail.com>

CC: fluka-discuss@fluka.org

Subject: Re: deq99c.f & dose equivalent

References: <200803040821.m248L7As029152@smtp2.mi.infn.it>=09 <47CD0CCB.409=

0509@ipno.in2p3.fr> <64c3dcd00803041240v4bc7ffbcmf89f4bec5898ce3e@mail.gmai=

l.com>

In-Reply-To: <64c3dcd00803041240v4bc7ffbcmf89f4bec5898ce3e@mail.gmail.com>

Content-Type: text/plain; charset=3DISO-8859-1; format=3Dflowed

Content-Transfer-Encoding: 8bit

Sender: owner-fluka-discuss@smtp2.mi.infn.it

Hello Hamide,

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=

=B2/kg.

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.

Best regards.

Sebastien.

hamide jalali wrote :

*> Dear Sebastien
*

*> =20
*

*> 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?
*

*> =20
*

*> Any help would be appreciated in advance.
*

*>
*

*> Kind Regards
*

*> Hamide
*

*> =20
*

*> =20
*

*> On Tue, Mar 4, 2008 at 12:48 AM, Sebastien WURTH <wurth@ipno.in2p3.fr=20
*

*> <mailto:wurth@ipno.in2p3.fr>> 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=
*

s

*> included in the routine is the right thing to do, I used it very
*

*> often.
*

*> 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.
*

*> Regards.
*

*> Sebastien.
*

*>
*

*>
*

*> 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=
*

s

*> >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
*

*> equivalent
*

*> >for each secondary separately, then add them up? How about
*

*> generalised
*

*> >particle 208 (energy)?
*

*> >And the other question, among USRBIN,USRBDX and USRTRACK which one i=
*

s

*> >recommended to calculate dose equivalent?
*

*> >
*

*> >Best Regards
*

*> >Hamide
*

*> >
*

*> >
*

*> >
*

*>
*

*>
*

**Previous message:**¿µçô: "fluka & cygwin"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Wed Mar 05 2008 - 12:08:28 CET
*