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