GEANT321/fluka/sitsao.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:05  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.45  by  S.Giani
  11 *-- Author :
  12 *$ CREATE SITSAO.FOR
  13 *COPY SITSAO
  14 *
  15 *=== sitsao ===========================================================*
  16 *
  17       FUNCTION SITSAO ( E, IZ, A )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *                                                                      *
  25 *     Original from A. Fasso`, CERN-TIS                                *
  26 *                                                                      *
  27 *     Created on 20 september 1991 by    Alfredo Ferrari & Paola Sala  *
  28 *                                                   Infn - Milan       *
  29 *                                                                      *
  30 *     Last change on 20-sep-91     by    Alfredo Ferrari               *
  31 *                                                                      *
  32 *                                                                      *
  33 *----------------------------------------------------------------------*
  34 *
  35 *----------------------------------------------------------------------*
  36 *   Calculates proton-nucleus inelastic cross-section by the formula   *
  37 *   of Letaw, Silberberg and Tsao, Astrophys. J. Suppl. 51, 271 (1983) *
  38 *----------------------------------------------------------------------*
  39 *   E = Energy in MeV, IZ = Atomic number, A = Mass number (At. weight)*
  40 *----------------------------------------------------------------------*
  41 *                                                formula (2)
  42       SITSAO = 45.D+00 * A**0.7D+00 * ( 1.D+00 + 0.016D+00 * SIN
  43      &       ( 5.3D+00 - 2.63D+00 * LOG(A)) )
  44 *                                                formula (5)
  45       IF(E .LT. 2000.D+00) SITSAO = SITSAO * ( 1.D+00 - 0.62D+00
  46      &                            / EXP ( E / 200.D+00 )
  47      &                            * SIN ( 10.9D+00 / E**0.28D+00 ) )
  48       IF (IZ .EQ. 2) THEN
  49          SITSAO=SITSAO*0.8D+00
  50       ELSE IF (IZ .EQ. 4) THEN
  51          SITSAO=SITSAO*(1.D+00 + 0.75D+00/EXP(E/75.D+00))
  52       END IF
  53       RETURN
  54       END