Re: a question?

 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