GEANT321/gstrag/gsttab.F

   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
  26
  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