GEANT321/gstrag/gstcor.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 GSTCOR(CC,POT,X,DENS,DELTA)
  13
  14 C.
  15 C.    *****************************************************************
  16 C.    *                                                               *
  17 C.    *      Sternheimer correction as computed in GDRELX             *
  18 C.    *                                                               *
  19 C.    *    ==>Called by :                                             *
  20 C.    *       Author    L.Urban    *********                          *
  21 C.    *                                                               *
  22 C.    *****************************************************************
  23 *
  24 *         condensed material ? ( dens .gt. 0.05 ? )
  25       IF(DENS.GT.0.05) THEN
  26 *        IF(POT.LT.1.E-7) THEN   (POT was in GeV, now in keV)
  27          IF(POT.LT.0.1) THEN
  28             IF(CC.LT.3.681) THEN
  29                X0=0.2
  30             ELSE
  31                X0=0.326*CC-1.
  32             ENDIF
  33             X1=2.
  34          ELSE
  35             IF(CC.LT.5.215) THEN
  36                X0=0.2
  37             ELSE
  38                X0=0.326*CC-1.5
  39             ENDIF
  40             X1=3.
  41          ENDIF
  42 *         gas ?   ( dens . le . 0.05 ? )
  43       ELSE
  44          IF(CC.LE.12.25) THEN
  45             IP=INT((CC-10.)/0.5)+1
  46             IF(IP.LT.0) IP=0
  47             IF(IP.GT.4) IP=4
  48             X0=1.6+0.1*FLOAT(IP)
  49             X1=4.
  50          ELSE
  51             IF(CC.LE.13.804) THEN
  52                X0=2.
  53                X1=5.
  54             ELSE
  55                X0=0.326*CC-2.5
  56                X1=5.
  57             ENDIF
  58          ENDIF
  59       ENDIF
  60 *
  61       XA=CC/4.606
  62       XM=3.
  63       AA=4.606*(XA-X0)/(X1-X0)**XM
  64 *
  65       DELTA=0.
  66       IF(X.GT.X0) THEN
  67          DELTA=4.606*X-CC
  68          IF(X.LT.X1) DELTA=DELTA+AA*(X1-X)**XM
  69       ENDIF
  70
  71       END