GEANT321/fluka/xsepro.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 XSEPRO.FOR
  13 *COPY XSEPRO
  14 *
  15 *=== xsepro ===========================================================*
  16 *
  17       FUNCTION XSEPRO ( E, ZTAR, ATAR )
  18
  19 #include "geant321/dblprc.inc"
  20 #include "geant321/dimpar.inc"
  21 #include "geant321/iounit.inc"
  22 *
  23 *----------------------------------------------------------------------*
  24 *                                                                      *
  25 *     Created on 20 september 1991 by    Alfredo Ferrari & Paola Sala  *
  26 *                                                   Infn - Milan       *
  27 *                                                                      *
  28 *     Last change on 20-sep-91     by    Alfredo Ferrari               *
  29 *                                                                      *
  30 *                                                                      *
  31 *----------------------------------------------------------------------*
  32 *
  33 #include "geant321/xsepar.inc"
  34        COMMON / FKNEGX / XSENEG
  35 *
  36        IZ = NINT ( ZTAR )
  37        IF ( E .GT. 0.1D+00 ) THEN
  38           EMEV = E * 1.D+03
  39           XSEPRO = SITSAO ( EMEV, IZ, ATAR )
  40        ELSE
  41           X  = 1.D+01 * E
  42           X2 = X * X
  43           XSEPRO = AAPXSE (IZ) * ( X - ZZPXSE (IZ) ) / X * ( 1.D+00
  44      &           + BBPXSE (IZ) * X + CCPXSE (IZ) * X2 + DDPXSE (IZ)
  45      &           * X * X2 + EEPXSE (IZ) * X2 * X2 + FFPXSE (IZ) * X2
  46      &           * X * X2 )
  47           IF ( XSEPRO .LT. 0.D+00 ) THEN
  48              XSENEG = XSEPRO
  49              XSEPRO = 0.D+00
  50           END IF
  51        END IF
  52 *=== End of function XSEPRO ===========================================*
  53        RETURN
  54        END