5 * Revision 1.1.1.1 1995/10/24 10:21:38 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.25 by S.Giani
12 FUNCTION GOSCIN(EIN1EV,EIN2EV)
14 C. ******************************************************************
16 C. * Input energy in eV, Sandia tables in keV *
19 C. ******************************************************************
21 #include "geant321/gcbank.inc"
22 #include "geant321/gcjloc.inc"
23 #include "geant321/gconsp.inc"
24 #include "geant321/gcmate.inc"
25 #include "geant321/gc10ev.inc"
26 #if !defined(CERNLIB_SINGLE)
27 DOUBLE PRECISION ONE,RES,ZERO,EIN1,EIN2,REST,EBEG,EEND
28 DOUBLE PRECISION E1IN,E2IN,TWO,THREE,C1,C2,C3,C4
30 PARAMETER (ZERO=0,ONE=1,TWO=2,THREE=3)
32 C. ------------------------------------------------------------------
52 IF(EIN2.GT.EBEG.AND.EIN1.LT.EEND) THEN
53 E1IN = ONE/MAX(EBEG,EIN1)
54 E2IN = ONE/MIN(EEND,EIN2)
61 + +C1*LOG(E2IN)+E2IN*(C2+E2IN*(C3+E2IN*C4))
62 + -C1*LOG(E1IN)-E1IN*(C2+E1IN*(C3+E1IN*C4))
66 20 REST = REST+WEIGHT*RES
68 C RES value is in cm**2/(g keV)
69 GOSCIN = REST*1E-3*A/AVO