1 *
2 * \$Id\$
3 *
4 * \$Log\$
5 * Revision 1.1.1.1  1995/10/24 10:21:39  cernlib
6 * Geant
7 *
8 *
9 #include "geant321/pilot.h"
10 *CMZ :  3.21/02 29/03/94  15.41.25  by  S.Giani
11 *-- Author :
12       SUBROUTINE GSTTAB(GAMMA,NT,EN,FN)
13 C   ********************************************************************
14 C   *                                                                  *
15 C   *  dN^2/dEdx table for a gamma factor                              *
16 C   *                                                                  *
17 C   *  ==> Called by : GSTINI                                          *
18 C   *      Authors   : P. Nevski                                       *
19 C   *                                                                  *
20 C   ********************************************************************
21 #include "geant321/gcstra.inc"
22       EXTERNAL GSTDN
23 #if !defined(CERNLIB_SINGLE)
24       DOUBLE PRECISION DGAMM, DGAM2, DBET2
25 #endif
27       DIMENSION EN(1),FN(1)
28       PARAMETER (RENER=0.05)
29 *
30 *     ------------------------------------------------------------------
31 *
32 C   ----  final tables for dN/dx integral vs Energy for a gamma   ------
33       DGAMM = GAMMA
34       DGAM2 = DGAMM**2
35       GAMMA2= DGAM2
36       DBET2 = 1-1/DGAM2
37       BETA2 = DBET2
38       S=0
39       NT=0
40       ES=0
41       EI=EMAX
42 C   Loop over I=NTAB,1,-1
43       DO 10 I=NTAB,1,-1
44          EO=EI
45          EI=EMM(I)
46          S=S+GXGINT(GSTDN,EI,EO,EPPS)
47 C      Skip Unless ABS(ES-EI)>RENER \ I==1
48          IF ((ABS(ES-EI).LE.RENER).AND.(I.NE.1))GO TO 10
49          ES=EI
50          NT=NT+1
51          EN(NT)=EI
52          FN(NT)=S
53    10 CONTINUE
54       DNDX=S
55       END