**From:** Francesco Cerutti (*Francesco.Cerutti@cern.ch*)

**Date:** Tue Nov 13 2007 - 17:50:44 CET

**Previous message:**Stefan Roesler: "Re: Fluence-to-Effective Dose/Ambient Dose Equivalent Calculations"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Dear Denis,

*> I have the following experiment to simulate:
*

*> Projectile : A gold beam particle at 25 GeV/nucleons
*

*> Fixed target: 1% interaction Au target (corresponds ~ 250 micrometers )
*

*> Iron beam dump at the end of the cave.
*

*> Now i would like to study the fluences of electron, protons, pions and
*

*> mostly neutrons inside the cave.
*

*> For minimum bias event i use Fluka heavy ion definition with DpmJetIII.
*

*> For Au-Au collision i use a external source ( using source.f ) reading
*

*> UqQmd output file and filling the stack.
*

*> Now i want to estimate the fluence of all particles in the cave.
*

*> My question is how to modify source.f in order to have the correct
*

*> normalization for the fluence?
*

*> In source.f there is the line calculating total weight for primaries:
*

*> WEIPRI = WEIPRI + WTFLKA( NPFLKA )
*

*> should i modify it to
*

*> WEIPRI = 1/100 ?
*

if you want to use your own list of secondary particles - so that each

of them represents a primary history - and to assume an interaction

probability of 1% for the incident Au projectiles, you should in

principle randomly select one Au+Au event in your list and hence load all

the respective secondary products one after the other, each as a primary

particle.

At the same time I would keep unchanged the weight definition in source.f

* Wt is the weight of the particle

WTFLK (NPFLKA) = ONEONE

WEIPRI = WEIPRI + WTFLK (NPFLKA)

[by the way, WEIPRI = 1/100 means WEIPRI = 0 since you have a division

between two integers]

and I would rescale fluences for each cycle eventually multiplying them by

(SCAFAC * 1D-2) where 1d-2 is the assumed ratio between the interacting

gold ions and the incident ones, and SCAFAC is the ratio between the

number of transported reaction products (i.e. your number of histories N)

and the number of considered collisions.

Then, merging results from different cycles, you should NOT weigh them

according to their number of histories N_i, but according to N_i/SCAFAC_i.

On the other hand, to properly overcome the low interaction probability

due to your very thin target, you might want to apply biasing: the

LAM-BIAS card allows you to reduce the hadronic inelastic interaction

length (by the factor WHAT(2) with SDUM blank) for a given material

(WHAT(3), gold in your case).

Hope this helps

Francesco

**************************************************

Francesco Cerutti

CERN-AB

CH-1211 Geneva 23

Switzerland

tel. ++41 22 7678962

fax ++41 22 7668854

**Previous message:**Stefan Roesler: "Re: Fluence-to-Effective Dose/Ambient Dose Equivalent Calculations"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Tue Nov 13 2007 - 23:29:32 CET
*