This commit was generated by cvs2svn to compensate for changes in r2,
[u/mrichter/AliRoot.git] / 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