GEANT321/fluka/eexi.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:20:03  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.44  by  S.Giani
  11 *-- Author :
  12 *$ CREATE EEXI.FOR
  13 *COPY EEXI
  14 *                                                                      *
  15 *=== eexi =============================================================*
  16 *                                                                      *
  17       FUNCTION EEXI(IJ,E,A)
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 C********************************************************************
  23 C     VERSION BY                     J. RANFT
  24 C                                    LEIPZIG
  25 C
  26 C     NUCLEAR EXCITATION ENERGY INCLUDING ALL LOW ENERGY SECONDARIES
  27 C     IN COLLISIONS OF PARTICLE IJ OF ENERGY E GEV ON NUCLEI A
  28 C
  29 C     NOTE!!!!!!! NON-STANDARD PARTICLE NUMBERING
  30 C
  31 C********************************************************************
  32 C
  33       IF (E.LE.0.125D0) GO TO 20
  34       B=SQRT(A)/9.D0-.2D0
  35       IF (B.GT.A*0.01D0) B=A*0.01D0
  36       IF (E.GE.3.D0) GO TO 10
  37       IF (B.LT.0.125D0) B=0.125D0
  38       EEXI=0.125D0+(E-0.125D0)*(B-0.125D0)/2.875D0
  39       RETURN
  40    10 CONTINUE
  41       EEXI=B
  42       RETURN
  43    20 CONTINUE
  44       EEXI=E
  45       RETURN
  46       END