Some function moved to AliZDC
[u/mrichter/AliRoot.git] / GEANT321 / gstrag / gstcor.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 GSTCOR(CC,POT,X,DENS,DELTA)
13
14C.
15C. *****************************************************************
16C. * *
17C. * Sternheimer correction as computed in GDRELX *
18C. * *
19C. * ==>Called by : *
20C. * Author L.Urban ********* *
21C. * *
22C. *****************************************************************
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