]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:21:36 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | *CMZ : 3.21/02 29/03/94 15.41.23 by S.Giani | |
11 | *-- Author : | |
12 | FUNCTION GHSIG(P,ITYPE,A) | |
13 | C. | |
14 | C. ****************************************************************** | |
15 | C. * * | |
16 | C. * Returns absorption cross-section in millibarns * | |
17 | C. * for a particle with * | |
18 | C. * momentum P (GEV/C) * | |
19 | C. * type ITYPE * | |
20 | C. * on a nucleus of atomic number A * | |
21 | C. * * | |
22 | C. * ==>Called by : GHSIGM * | |
23 | C. * Author R.Barlow ********* * | |
24 | C. * * | |
25 | C. ****************************************************************** | |
26 | C | |
27 | C The internal particle types are as follows... | |
28 | C 1 proton | |
29 | C 2 antiproton | |
30 | C 3 K+ | |
31 | C 4 K- | |
32 | C 5 pion | |
33 | C | |
34 | DIMENSION SIGB(100),ALPHA(100),SIG(162,5) | |
35 | DIMENSION SIGPR(162,1),SIGPB(162,1),SIGKP(162,1),SIGKM(162,1) | |
36 | DIMENSION SIGPI(162,1) | |
37 | EQUIVALENCE (SIGPR(1,1),SIG(1,1)),(SIGPB(1,1),SIG(1,2)) | |
38 | EQUIVALENCE (SIGKP(1,1),SIG(1,3)),(SIGKM(1,1),SIG(1,4)) | |
39 | EQUIVALENCE (SIGPI(1,1),SIG(1,5)) | |
40 | C | |
41 | DATA ALPHA/ | |
42 | + 0.9826,0.9659,0.9500,0.9348,0.9203,0.9064,0.8932,0.8807,0.8687, | |
43 | + 0.8574,0.8466,0.8363,0.8265,0.8172,0.8084,0.8000,0.7921,0.7845, | |
44 | + 0.7773,0.7705,0.7640,0.7579,0.7520,0.7464,0.7411,0.7361,0.7313, | |
45 | + 0.7267,0.7224,0.7182,0.7142,0.7105,0.7069,0.7034,0.7001,0.6970, | |
46 | + 0.6939,0.6911,0.6883,0.6856,0.6831,0.6807,0.6783,0.6761,0.6739, | |
47 | + 0.6718,0.6698,0.6679,0.6660,0.6642,0.6625,0.6609,0.6592,0.6577, | |
48 | + 0.6562,0.6547,0.6533,0.6519,0.6506,0.6493,0.6481,0.6469,0.6457, | |
49 | + 0.6446,0.6435,0.6424,0.6413,0.6403,0.6393,0.6384,0.6374,0.6365, | |
50 | + 0.6356,0.6347,0.6339,0.6330,0.6322,0.6314,0.6307,0.6299,0.6292, | |
51 | + 0.6284,0.6277,0.6270,0.6263,0.6257,0.6250,0.6244,0.6237,0.6231, | |
52 | + 0.6225,0.6219,0.6213,0.6208,0.6202,0.6197,0.6191,0.6186,0.6181, | |
53 | + 0.6175/ | |
54 | DATA SIGB/ | |
55 | + 1.02, 2.10, 3.21, 4.36, 5.54, 6.75, 7.99, 9.24, 10.51, | |
56 | + 11.79, 13.08, 14.37, 15.66, 16.95, 18.22, 19.49, 20.75, 21.99, | |
57 | + 23.22, 24.43, 25.62, 26.79, 27.93, 29.06, 30.16, 31.24, 32.30, | |
58 | + 33.33, 34.34, 35.32, 36.29, 37.22, 38.14, 39.03, 39.91, 40.75, | |
59 | + 41.58, 42.39, 43.18, 43.95, 44.70, 45.43, 46.14, 46.83, 47.51, | |
60 | + 48.17, 48.82, 49.45, 50.06, 50.66, 51.25, 51.83, 52.39, 52.93, | |
61 | + 53.47, 53.99, 54.51, 55.01, 55.50, 55.98, 56.45, 56.92, 57.37, | |
62 | + 57.81, 58.25, 58.68, 59.10, 59.51, 59.91, 60.31, 60.70, 61.08, | |
63 | + 61.46, 61.83, 62.19, 62.55, 62.90, 63.25, 63.59, 63.93, 64.26, | |
64 | + 64.58, 64.90, 65.22, 65.53, 65.84, 66.14, 66.44, 66.74, 67.03, | |
65 | + 67.31, 67.60, 67.88, 68.15, 68.42, 68.69, 68.96, 69.22, 69.48, | |
66 | + 69.74/ | |
67 | DATA SIGPR / | |
68 | + 79.01, 76.05, 73.81, 71.75, 69.67, 66.70, 64.89, 63.09, 60.58, | |
69 | + 58.99, 57.49, 56.09, 54.02, 52.75, 51.54, 50.39, 48.78, 47.79, | |
70 | + 46.86, 45.99, 44.73, 43.97, 43.25, 42.59, 41.59, 40.99, 40.41, | |
71 | + 39.85, 38.88, 38.26, 37.67, 37.06, 36.31, 35.82, 35.37, 35.00, | |
72 | + 34.40, 34.03, 33.71, 33.45, 33.00, 32.72, 32.50, 32.31, 32.01, | |
73 | + 31.84, 31.69, 31.56, 31.38, 31.29, 31.17, 31.09, 31.17, 31.34, | |
74 | + 31.61, 31.87, 32.62, 33.16, 33.84, 34.62, 35.85, 36.67, 37.51, | |
75 | + 38.72, 39.52, 40.29, 41.00, 41.97, 42.54, 43.01, 43.38, 43.82, | |
76 | + 44.01, 44.14, 44.24, 44.34, 44.38, 44.40, 44.40, 44.36, 44.33, | |
77 | + 44.28, 44.21, 44.07, 43.96, 43.84, 43.71, 43.48, 43.33, 43.17, | |
78 | + 43.02, 42.79, 42.64, 42.49, 42.36, 42.16, 42.05, 41.94, 41.85, | |
79 | + 41.71, 41.63, 41.54, 41.46, 41.35, 41.27, 41.20, 41.13, 41.02, | |
80 | + 40.95, 40.88, 40.81, 40.72, 40.66, 40.61, 40.57, 40.49, 40.44, | |
81 | + 40.37, 40.30, 40.19, 40.11, 40.04, 39.92, 39.84, 39.76, 39.69, | |
82 | + 39.57, 39.51, 39.45, 39.40, 39.32, 39.26, 39.21, 39.16, 39.09, | |
83 | + 39.05, 39.01, 38.98, 38.93, 38.90, 38.88, 38.85, 38.82, 38.81, | |
84 | + 38.79, 38.77, 38.75, 38.73, 38.72, 38.70, 38.68, 38.66, 38.65, | |
85 | + 38.64, 38.62, 38.60, 38.59, 38.57, 38.55, 38.53, 38.51, 38.49/ | |
86 | DATA SIGPB / | |
87 | + 505.81,485.41,470.03,455.70,441.22,420.60,407.94,395.19,377.36, | |
88 | + 366.03,355.16,345.07,329.94,320.60,311.79,303.19,291.17,283.77, | |
89 | + 276.68,270.09,260.15,254.24,248.43,243.19,234.87,229.91,225.10, | |
90 | + 220.26,213.46,209.12,204.98,200.98,195.36,191.56,187.98,184.49, | |
91 | + 179.41,176.04,172.86,169.81,165.19,162.28,159.35,156.55,152.41, | |
92 | + 149.95,147.59,145.11,141.66,139.59,137.45,135.41,132.43,130.53, | |
93 | + 128.59,126.81,124.18,122.54,120.90,119.34,117.13,115.67,114.22, | |
94 | + 112.04,110.54,109.13,107.70,105.63,104.23,102.97,101.64, 99.62, | |
95 | + 98.30, 96.99, 95.81, 93.97, 92.88, 91.80, 90.69, 89.03, 87.97, | |
96 | + 86.92, 85.89, 84.36, 83.38, 82.39, 81.40, 79.89, 78.92, 77.95, | |
97 | + 76.96, 75.50, 74.59, 73.68, 72.81, 71.51, 70.73, 69.95, 69.19, | |
98 | + 68.07, 67.39, 66.71, 66.04, 65.13, 64.54, 63.99, 63.42, 62.55, | |
99 | + 62.05, 61.52, 61.02, 60.23, 59.74, 59.27, 58.79, 58.07, 57.64, | |
100 | + 57.20, 56.77, 56.10, 55.66, 55.23, 54.41, 53.90, 53.42, 52.94, | |
101 | + 52.20, 51.73, 51.30, 50.85, 50.20, 49.79, 49.47, 49.16, 48.72, | |
102 | + 48.46, 48.15, 47.84, 47.42, 47.17, 46.92, 46.68, 46.31, 46.12, | |
103 | + 45.90, 45.70, 45.42, 45.28, 45.13, 44.99, 44.80, 44.68, 44.58, | |
104 | + 44.48, 44.35, 44.28, 44.21, 44.15, 44.06, 44.00, 43.95, 43.90/ | |
105 | DATA SIGKP / | |
106 | + 9.05, 8.86, 8.71, 8.58, 8.45, 8.28, 8.19, 8.11, 8.02, | |
107 | + 7.97, 7.94, 7.92, 7.90, 7.91, 7.92, 7.95, 8.01, 8.05, | |
108 | + 8.10, 8.17, 8.28, 8.37, 8.47, 8.56, 8.74, 8.86, 8.98, | |
109 | + 9.13, 9.35, 9.51, 9.66, 9.83, 10.10, 10.27, 10.45, 10.63, | |
110 | + 10.89, 11.07, 11.24, 11.40, 11.65, 11.81, 11.96, 12.12, 12.37, | |
111 | + 12.54, 12.71, 12.87, 13.19, 13.43, 13.70, 13.99, 14.49, 14.87, | |
112 | + 15.28, 15.71, 16.37, 16.79, 17.19, 17.57, 18.05, 18.30, 18.49, | |
113 | + 18.68, 18.74, 18.76, 18.75, 18.69, 18.63, 18.57, 18.50, 18.39, | |
114 | + 18.34, 18.29, 18.25, 18.18, 18.14, 18.09, 18.05, 17.98, 17.92, | |
115 | + 17.87, 17.81, 17.72, 17.66, 17.62, 17.59, 17.55, 17.52, 17.49, | |
116 | + 17.46, 17.44, 17.43, 17.42, 17.41, 17.41, 17.40, 17.40, 17.40, | |
117 | + 17.40, 17.40, 17.40, 17.40, 17.41, 17.41, 17.42, 17.42, 17.43, | |
118 | + 17.43, 17.44, 17.44, 17.45, 17.45, 17.46, 17.47, 17.48, 17.48, | |
119 | + 17.48, 17.49, 17.49, 17.50, 17.50, 17.50, 17.49, 17.49, 17.49, | |
120 | + 17.49, 17.50, 17.50, 17.50, 17.51, 17.52, 17.54, 17.55, 17.58, | |
121 | + 17.60, 17.62, 17.63, 17.66, 17.68, 17.70, 17.73, 17.76, 17.79, | |
122 | + 17.81, 17.83, 17.87, 17.90, 17.92, 17.95, 17.99, 18.01, 18.03, | |
123 | + 18.06, 18.10, 18.13, 18.16, 18.19, 18.23, 18.26, 18.29, 18.33/ | |
124 | DATA SIGKM / | |
125 | + 76.56, 74.23, 72.44, 70.86, 69.28, 66.99, 65.66, 64.35, 62.47, | |
126 | + 61.33, 60.23, 59.18, 57.65, 56.68, 55.73, 54.78, 53.21, 52.17, | |
127 | + 51.31, 50.48, 49.27, 48.55, 47.89, 47.30, 46.48, 45.91, 45.31, | |
128 | + 44.74, 43.78, 43.13, 42.48, 41.77, 40.57, 39.85, 39.01, 38.22, | |
129 | + 37.05, 36.34, 35.66, 35.01, 34.23, 33.98, 33.73, 33.53, 33.26, | |
130 | + 33.24, 33.41, 33.68, 34.13, 34.53, 34.94, 35.52, 36.37, 37.02, | |
131 | + 37.58, 38.00, 38.11, 37.88, 37.44, 36.88, 35.77, 35.12, 34.34, | |
132 | + 33.26, 32.60, 31.91, 31.29, 30.54, 30.13, 29.77, 29.45, 29.03, | |
133 | + 28.80, 28.58, 28.35, 27.97, 27.70, 27.40, 27.10, 26.70, 26.48, | |
134 | + 26.29, 26.07, 25.73, 25.52, 25.31, 25.13, 24.86, 24.72, 24.57, | |
135 | + 24.42, 24.17, 24.01, 23.89, 23.77, 23.60, 23.49, 23.38, 23.29, | |
136 | + 23.14, 23.05, 22.97, 22.90, 22.79, 22.71, 22.64, 22.56, 22.45, | |
137 | + 22.37, 22.28, 22.20, 22.07, 21.99, 21.91, 21.82, 21.69, 21.61, | |
138 | + 21.53, 21.45, 21.34, 21.27, 21.20, 21.11, 21.06, 21.00, 20.95, | |
139 | + 20.87, 20.83, 20.80, 20.77, 20.71, 20.68, 20.65, 20.62, 20.57, | |
140 | + 20.54, 20.52, 20.50, 20.46, 20.43, 20.40, 20.38, 20.35, 20.32, | |
141 | + 20.30, 20.28, 20.25, 20.24, 20.23, 20.21, 20.20, 20.20, 20.19, | |
142 | + 20.19, 20.18, 20.18, 20.18, 20.17, 20.16, 20.15, 20.14, 20.12/ | |
143 | DATA SIGPI / | |
144 | + 5.02, 5.69, 6.20, 6.81, 7.56, 9.29, 10.54, 12.24, 15.11, | |
145 | + 17.54, 20.45, 23.24, 28.67, 32.75, 37.60, 42.78, 52.32, 59.55, | |
146 | + 67.49, 75.47, 87.11, 94.75,101.54,106.41,111.56,113.00,112.73, | |
147 | + 110.52,105.96,101.22, 95.06, 88.38, 77.77, 70.51, 63.86, 57.83, | |
148 | + 49.65, 45.30, 41.45, 38.23, 34.33, 32.40, 30.86, 29.68, 28.76, | |
149 | + 28.61, 28.53, 28.41, 28.04, 27.81, 27.89, 28.24, 29.65, 31.04, | |
150 | + 32.53, 33.91, 34.81, 34.92, 35.12, 35.32, 35.83, 36.30, 36.71, | |
151 | + 36.93, 36.79, 36.49, 36.18, 35.80, 35.59, 35.35, 35.12, 34.65, | |
152 | + 34.26, 33.85, 33.46, 32.97, 32.75, 32.61, 32.51, 32.37, 32.25, | |
153 | + 32.12, 31.96, 31.67, 31.47, 31.25, 31.02, 30.67, 30.44, 30.21, | |
154 | + 29.98, 29.65, 29.44, 29.23, 29.03, 28.72, 28.52, 28.33, 28.14, | |
155 | + 27.85, 27.66, 27.47, 27.30, 27.04, 26.90, 26.76, 26.65, 26.51, | |
156 | + 26.43, 26.35, 26.28, 26.19, 26.13, 26.07, 26.02, 25.92, 25.85, | |
157 | + 25.78, 25.70, 25.58, 25.51, 25.43, 25.32, 25.25, 25.18, 25.11, | |
158 | + 25.01, 24.94, 24.87, 24.81, 24.72, 24.65, 24.59, 24.54, 24.46, | |
159 | + 24.41, 24.36, 24.32, 24.25, 24.21, 24.18, 24.14, 24.09, 24.06, | |
160 | + 24.03, 24.00, 23.96, 23.93, 23.91, 23.88, 23.84, 23.82, 23.80, | |
161 | + 23.78, 23.76, 23.74, 23.73, 23.73, 23.72, 23.72, 23.72, 23.73/ | |
162 | C | |
163 | C 25.49675=1./LOG(1.04) | |
164 | C | |
165 | X = 25.49675 * LOG(P * 10.) | |
166 | N = X | |
167 | DIFF = X - N | |
168 | N = N + 1 | |
169 | IF (N .LE. 0) N = 1 | |
170 | IF (N .GT. 162) N = 162 | |
171 | C | |
172 | C Interpolate in SIG tables. | |
173 | C | |
174 | S = SIG(N,ITYPE) | |
175 | IF (N.EQ.1 .OR. N.EQ.162) GO TO 50 | |
176 | S = S + DIFF * (SIG(N + 1,ITYPE) - S) | |
177 | 50 CONTINUE | |
178 | C | |
179 | C NOW FIND ABUL-MAGD PARAMETERS | |
180 | C | |
181 | IS = S | |
182 | IF (IS.LE.0) IS = 1 | |
183 | IF (IS.GE.100) IS = 99 | |
184 | SB = SIGB(IS) + (S - IS) * (SIGB(IS + 1) - SIGB(IS)) | |
185 | AL = ALPHA(IS) + (S - IS) * (ALPHA(IS + 1) - ALPHA(IS)) | |
186 | GHSIG = SB*A**AL | |
187 | END |