From: <Anna.Ferrari@lnf.infn.it>

Date: Mon Jul 28 2008 - 09:17:39 CEST

Date: Mon Jul 28 2008 - 09:17:39 CEST

Hi Elham,

to obtain energy distribution from secondary ions produced by a proton

beam on a target (like your C12H18O7 compound) you have to follow several

steps, to:

a) activate ion production and transport;

b) obtain energy distribution of different heavy fragments.

a) - build your FLUKA executable by linking with DPMJET and RQMD codes.

See for that the script $FLUPRO/flutil/ldp3qmd

- activate in your input card heavy fragment evaporation with the card

PHYSICS with sdum=EVAPORAT and what(1) = 3.0

- activate the transport of all heavy recoils and ions with energy loss

and multiple scattering via the card EVENTYPE with what(3)=2.0.

By setting sdum=DPMJET you activate heavy ion nuclear interactions,

starting from 100 MeV/u up to cosmic ray energies.

Therefore.. be careful: you must look at your energy range if you want

to take into account heavy ion nuclear interactions.

b) To obtain energy distributions of produced heavy fragments, the

simplest way is to use the existing scoring cards like USRTRACK,

USRBDX, etc, together with the user routine fluscw.f, that allows the

heavy ion scoring to be appropriately set (see the answer of Vasilis on

the fluka-list of july, 21: it's exactly your case).

Suppose for example that you want to calculate the energy distribution

of C-12,N-15 and O-16 at the exit of your target. The fluence boundary

crossing estimator is appropriate in this case: you will use three

subsequent USRBDX cards (three different detectors), each one

corresponding to a different ion scoring. The right association is made

in the FLUSCW routine, where you set for the three detectors (JSCRNG =

1,2,3), that are associated to a boundary crossing estimator (ISCRNG =

1), the right baryon number IBARCH(IJ) and the right charge

ICHRGE(IJ) (from the (PAPROP) include file).

Your fluscw.f will look as:

......

INCLUDE '(PAPROP)'

.....

* Check for the relevant boundary crossing estimator:

IF (ISCRNG .EQ. 1 .AND. JSCRNG .EQ. 1) THEN

LSCZER = .NOT. (IBARCH(IJ) .EQ. 12 .AND. ICHRGE(IJ) .EQ. 6)

ELSEIF (ISCRNG .EQ. 1 .AND. JSCRNG .EQ. 2) THEN

LSCZER = .NOT. (IBARCH(IJ) .EQ. 15 .AND. ICHRGE(IJ) .EQ. 7)

ELSEIF (ISCRNG .EQ. 1 .AND. JSCRNG .EQ. 3) THEN

LSCZER = .NOT. (IBARCH(IJ) .EQ. 16 .AND. ICHRGE(IJ) .EQ. 8)

END IF

...........

Last, remember that you activate the FLUSCW routine via the USERWEIGht

card, with what(3)> 0.0

Regards,

Anna

*> Dear Fluka-users:
*

*>
*

*> How can we obtain to energy distribution of secondary particles like
*

*> O-16,N-15,Be-8,B-11,C-12, Li-5(their fluence or flux) or at least their
*

*> energy range on (C12H18O7) target by 30MeV proton beam?
*

*>
*

*>
*

*> Best regards,
*

*> E.Salmasi
*

*>
*

*>
*

Received on Mon Jul 28 09:25:31 2008

*
This archive was generated by hypermail 2.1.8
: Mon Jul 28 2008 - 09:25:32 CEST
*