GEANT321/peanut/xinpro.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:22:04  cernlib
   6 * Geant
   7 *
   8 *
   9 #include "geant321/pilot.h"
  10 *CMZ :  3.21/02 29/03/94  15.41.46  by  S.Giani
  11 *-- Author :
  12 *$ CREATE XINPRO.FOR
  13 *COPY XINPRO
  14 *
  15 *=== xinpro ===========================================================*
  16 *
  17       FUNCTION XINPRO ( E, ZTAR, VCOUL )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *----------------------------------------------------------------------*
  25 *
  26 #include "geant321/isotop.inc"
  27 #include "geant321/xsepar.inc"
  28       COMMON / FKNEGX / XSENEG
  29 *
  30       IZ = NINT ( ZTAR )
  31       IF ( E .GT. 0.1D+00 ) THEN
  32          ATAR = 0.D+00
  33          DO 25 IS = ISONDX (1,IZ), ISONDX (2,IZ)
  34              ATAR = ATAR + ISOMNM (IS) * ABUISO (IS)
  35    25    CONTINUE
  36          EMEV = E * 1.D+03
  37          XINPRO = SITSAO ( EMEV, IZ, ATAR )
  38       ELSE
  39          ZZCOUL = 1.D+01 * VCOUL
  40          X  = 1.D+01 * E
  41          X2 = X * X
  42          XINPRO = AAPXSE (IZ) * ( X - ZZCOUL ) / X * ( 1.D+00
  43      &          + BBPXSE (IZ) * X + CCPXSE (IZ) * X2 + DDPXSE (IZ)
  44      &          * X * X2 + EEPXSE (IZ) * X2 * X2 + FFPXSE (IZ) * X2
  45      &          * X * X2 )
  46          IF ( XINPRO .LT. 0.D+00 ) THEN
  47             XSENEG = XINPRO
  48             XINPRO = 0.D+00
  49          END IF
  50       END IF
  51 *=== End of function XINPRO ===========================================*
  52        RETURN
  53        END