]> git.uio.no Git - u/mrichter/AliRoot.git/blame - GEANT321/fluka/xsepro.F
Default compile option changed to -g (Alpha)
[u/mrichter/AliRoot.git] / GEANT321 / fluka / xsepro.F
CommitLineData
fe4da5cc 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