Some function moved to AliZDC
[u/mrichter/AliRoot.git] / GEANT321 / gstrag / gsttab.F
CommitLineData
fe4da5cc 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)
13C ********************************************************************
14C * *
15C * dN^2/dEdx table for a gamma factor *
16C * *
17C * ==> Called by : GSTINI *
18C * Authors : P. Nevski *
19C * *
20C ********************************************************************
21#include "geant321/gcstra.inc"
22 EXTERNAL GSTDN
23#if !defined(CERNLIB_SINGLE)
24 DOUBLE PRECISION DGAMM, DGAM2, DBET2
25#endif
26
27 DIMENSION EN(1),FN(1)
28 PARAMETER (RENER=0.05)
29*
30* ------------------------------------------------------------------
31*
32C ---- 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
42C 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)
47C 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