* * $Id: xseneu.F,v 1.1.1.1 1995/10/24 10:20:05 cernlib Exp $ * * $Log: xseneu.F,v $ * Revision 1.1.1.1 1995/10/24 10:20:05 cernlib * Geant * * #include "geant321/pilot.h" *CMZ : 3.21/02 29/03/94 15.41.45 by S.Giani *-- Author : *$ CREATE XSENEU.FOR *COPY XSENEU * *=== xseneu ===========================================================* * FUNCTION XSENEU ( E, ZTAR, ATAR ) #include "geant321/dblprc.inc" #include "geant321/dimpar.inc" #include "geant321/iounit.inc" * *----------------------------------------------------------------------* * * * Created on 20 september 1991 by Alfredo Ferrari & Paola Sala * * Infn - Milan * * * * Last change on 20-sep-91 by Alfredo Ferrari * * * * * *----------------------------------------------------------------------* * #include "geant321/xsepar.inc" * IZ = NINT ( ZTAR ) IF ( E .GT. 0.1D+00 ) THEN EMEV = E * 1.D+03 XSENEU = SITSAO ( EMEV, IZ, ATAR ) ELSE IF ( E .GT. 0.012D+00 ) THEN X = 1.D+01 * E XSENEU = AANXSE (IZ) * ( X - ZZNXSE (IZ) ) / X * ( 1.D+00 & + BBNXSE (IZ) * X / ( 1.D+00 + (CCNXSE(IZ)*X)**2 ) & + DDNXSE (IZ) * X * X / ( 1.D+00 + (EENXSE(IZ)*X)**4 )) ELSE XSENEU = 0.D+00 END IF XSENEU = MAX ( XSENEU, ZERZER ) *=== End of function Xseneu ===========================================* RETURN END