| 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 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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame