GEANT321/fluka/xseneu.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 XSENEU.FOR
  13 *COPY XSENEU
  14 *
  15 *=== xseneu ===========================================================*
  16 *
  17       FUNCTION XSENEU ( 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 *
  35        IZ = NINT ( ZTAR )
  36        IF ( E .GT. 0.1D+00 ) THEN
  37           EMEV = E * 1.D+03
  38           XSENEU = SITSAO ( EMEV, IZ, ATAR )
  39        ELSE IF ( E .GT. 0.012D+00 ) THEN
  40           X  = 1.D+01 * E
  41           XSENEU = AANXSE (IZ) * ( X - ZZNXSE (IZ) ) / X * ( 1.D+00
  42      &           + BBNXSE (IZ) * X / ( 1.D+00 + (CCNXSE(IZ)*X)**2 )
  43      &           + DDNXSE (IZ) * X * X / ( 1.D+00 + (EENXSE(IZ)*X)**4 ))
  44        ELSE
  45           XSENEU = 0.D+00
  46        END IF
  47        XSENEU = MAX ( XSENEU, ZERZER )
  48 *=== End of function Xseneu ===========================================*
  49        RETURN
  50        END