5 * Revision 1.1.1.1 1995/10/24 10:21:37 cernlib
9 #include "geant321/pilot.h"
10 *CMZ : 3.21/02 29/03/94 15.41.25 by S.Giani
12 FUNCTION GKOKRI(E,EMINEV,EMAXEV)
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 EINV,ECUR,ONE,RES,ZERO,EBEG,EEND
28 DOUBLE PRECISION ALPH,BETA,GAMM,E1,E2,GPSCIN,POLE
29 DOUBLE PRECISION EMAX,EMIN,REST
31 PARAMETER (ONE=1,ZERO=0)
33 C. ------------------------------------------------------------------
59 IF(ECUR.LE.E2.AND.ECUR.GE.E1) THEN
61 * *** The pole of the integration is in this interval
62 EPS1 = (ECUR-EBEG)/ECUR
63 EPS2 = (EEND-ECUR)/EEND
66 * *** First the pole and then a simple integration
72 * *** First a simple integration and then the pole
79 POLE = EINV*(Q(J)*BETA+EINV*( Q(J+1)*(GAMM-ALPH)+
80 + EINV*( Q(J+2)*(0.5*GAMM**2+BETA-0.5*ALPH**2)+EINV*
81 + Q(J+3)*(GAMM**3/3+GAMM-ALPH-ALPH**3/3))))
85 * *** This is a normal integration
86 RES = RES + GPSCIN(E1,E2,ECUR,Q(J))
89 IF(EBEG.GE.EMAX) GO TO 20
91 20 REST = REST+WEIGHT*RES
93 C RES value is in cm**2/(g keV)
94 GKOKRI = REST*1E-6*A/AVO