1 |
cafagna |
3.2 |
|
2 |
cafagna |
3.1 |
* |
3 |
pam-ba |
3.5 |
* $Id: gpspev.F,v 3.4 2005/12/05 12:15:21 pam-ba Exp $ |
4 |
cafagna |
3.2 |
* |
5 |
|
|
* $Log: gpspev.F,v $ |
6 |
pam-ba |
3.5 |
* Revision 3.4 2005/12/05 12:15:21 pam-ba |
7 |
|
|
* new spectrometer geometry and internal magnetic field |
8 |
|
|
* |
9 |
pam-ba |
3.4 |
* Revision 3.3 2005/07/25 11:53:21 cafagna |
10 |
|
|
* Several updates. See history for details |
11 |
|
|
* |
12 |
cafagna |
3.3 |
* Revision 3.2 2005/06/21 02:42:27 cafagna |
13 |
|
|
* Major modification to the geometry and to the random number chain |
14 |
|
|
* |
15 |
cafagna |
3.2 |
* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna |
16 |
|
|
* First GPAMELA release on CVS |
17 |
cafagna |
3.1 |
* |
18 |
|
|
* |
19 |
|
|
*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola |
20 |
|
|
*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna |
21 |
|
|
*CMZ : 1.02/00 07/03/97 14.48.41 by Francesco Cafagna |
22 |
|
|
*CMZ : 1.00/02 06/02/96 17.34.50 by Francesco Cafagna |
23 |
|
|
*-- Author : Francesco Cafagna 09/12/95 |
24 |
|
|
SUBROUTINE GPSPEV |
25 |
|
|
************************************************************************ |
26 |
|
|
* * |
27 |
|
|
* Volumes definition for the silicon tracking system and magnet * |
28 |
|
|
* Called by: GPGEM * |
29 |
|
|
* Author: Francesco Cafagna, 05/12/95 13.39.41 * |
30 |
|
|
* * |
31 |
|
|
************************************************************************ |
32 |
|
|
#include "gpgeo.inc" |
33 |
|
|
#include "gpmed.inc" |
34 |
|
|
* |
35 |
|
|
INTEGER IROT,IVOLU,N,NMED,NUM |
36 |
|
|
REAL X,Y,Z |
37 |
|
|
* |
38 |
|
|
* Define the SPEB volume |
39 |
|
|
* |
40 |
|
|
NMED= MN2 |
41 |
|
|
CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU) |
42 |
|
|
* |
43 |
cafagna |
3.3 |
* Define the MGFR volume |
44 |
|
|
* |
45 |
|
|
NMED= MAL |
46 |
|
|
CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU) |
47 |
|
|
* |
48 |
cafagna |
3.1 |
* Define the MGPL volume |
49 |
|
|
* |
50 |
|
|
NMED= MMAG |
51 |
|
|
CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU) |
52 |
|
|
* |
53 |
|
|
* Define the MGPI volume |
54 |
|
|
* |
55 |
|
|
NMED= MN2 |
56 |
|
|
CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU) |
57 |
|
|
* |
58 |
|
|
* Define the TRPB volume |
59 |
|
|
* |
60 |
|
|
NMED= MAL |
61 |
|
|
CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU) |
62 |
cafagna |
3.2 |
***** |
63 |
|
|
cC ml 20/04/05: |
64 |
cafagna |
3.1 |
* |
65 |
cafagna |
3.2 |
cc* Define the TRPL volume |
66 |
|
|
* Define the TPAS and TPAI volume |
67 |
cafagna |
3.1 |
* |
68 |
|
|
NMED= MN2 |
69 |
cafagna |
3.2 |
CC CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU) |
70 |
|
|
CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU) |
71 |
cafagna |
3.3 |
NMED= MAL |
72 |
cafagna |
3.2 |
CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU) |
73 |
|
|
CC end ml. |
74 |
|
|
**************** |
75 |
cafagna |
3.1 |
* |
76 |
|
|
* Define the TRSL volume |
77 |
|
|
* |
78 |
cafagna |
3.2 |
c ml: 26/11/04: |
79 |
|
|
c NMED= MN2 |
80 |
|
|
NMED=MSIT |
81 |
cafagna |
3.1 |
CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU) |
82 |
|
|
* |
83 |
cafagna |
3.2 |
* Define the THBP volume |
84 |
|
|
* |
85 |
cafagna |
3.3 |
NMED=MG10C |
86 |
cafagna |
3.2 |
CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU) |
87 |
|
|
C end ml. |
88 |
|
|
* |
89 |
cafagna |
3.1 |
* Define the TSPA volume |
90 |
|
|
* |
91 |
|
|
NMED= MSIT |
92 |
|
|
CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU) |
93 |
|
|
* |
94 |
|
|
* Define the TRCP volume |
95 |
|
|
* |
96 |
|
|
NMED= MCF |
97 |
|
|
CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU) |
98 |
|
|
* |
99 |
pam-ba |
3.4 |
* Define the MGPA volume |
100 |
|
|
* |
101 |
|
|
NMED= MAL |
102 |
|
|
CALL GSVOLU('MGPA','BOX ',NMED,MGPA, 3,IVOLU) |
103 |
|
|
* |
104 |
|
|
* Define the TPGA volume |
105 |
|
|
* |
106 |
|
|
NMED= MN2 |
107 |
|
|
CALL GSVOLU('TPGA','BOX ',NMED,TPGA, 3,IVOLU) |
108 |
|
|
* |
109 |
pam-ba |
3.5 |
* Define the TPGI volume |
110 |
|
|
* |
111 |
|
|
NMED= MN2 |
112 |
|
|
CALL GSVOLU('TPGI','BOX ',NMED,TPGI, 3,IVOLU) |
113 |
|
|
* |
114 |
pam-ba |
3.4 |
* Define the TPGU volume |
115 |
|
|
* |
116 |
|
|
NMED= MN2 |
117 |
|
|
CALL GSVOLU('TPGU','BOX ',NMED,TPGU, 3,IVOLU) |
118 |
|
|
* |
119 |
|
|
* Define the TPGD volume |
120 |
|
|
* |
121 |
|
|
NMED= MN2 |
122 |
|
|
CALL GSVOLU('TPGD','BOX ',NMED,TPGD, 3,IVOLU) |
123 |
|
|
|
124 |
|
|
* |
125 |
cafagna |
3.1 |
* Positioning volume TSPA into the mother TRSL |
126 |
|
|
* |
127 |
|
|
N= 1 |
128 |
|
|
X= 0. |
129 |
cafagna |
3.3 |
Y= 0. |
130 |
cafagna |
3.1 |
Z= 0. |
131 |
|
|
CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY') |
132 |
cafagna |
3.3 |
C # N= 2 |
133 |
|
|
C # C Y= Y - 2.*TSPA(2) |
134 |
|
|
C # Y=-TRSL(2)/2 |
135 |
|
|
C # CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY') |
136 |
cafagna |
3.2 |
c ml: 26/11/04: |
137 |
pam-ba |
3.4 |
* Positioning volume TRCP into mother TPAS |
138 |
cafagna |
3.3 |
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one |
139 |
|
|
C # Z=-TRPB(3)+0.11+TRCP(3) |
140 |
|
|
C # Y=-TRPB(2)+1.5-0.225+TRCP(2) |
141 |
cafagna |
3.2 |
CC end ml. |
142 |
pam-ba |
3.4 |
c Z=-TRPB(3)+2*TPAI(3)+TRCP(3) |
143 |
|
|
Z=TRPB(3)-2*TPAI(3)-TRCP(3) |
144 |
cafagna |
3.3 |
y=-0.225-TPAS(2)+TRCP(2) |
145 |
|
|
C # end CAF mod |
146 |
cafagna |
3.2 |
DO I=1,4 |
147 |
|
|
CC ml 20/04/05: |
148 |
|
|
CC X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1) |
149 |
cafagna |
3.3 |
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one |
150 |
|
|
C # X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1) |
151 |
pam-ba |
3.4 |
c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders. |
152 |
|
|
c for now the glue is a space of nitrogen of 0.015 mm |
153 |
|
|
c X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1) |
154 |
|
|
X=-TPAS(1)+2*(I-1)*0.0015+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1) |
155 |
cafagna |
3.2 |
CC end ml |
156 |
cafagna |
3.3 |
CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY') |
157 |
|
|
C # end CAF mod |
158 |
cafagna |
3.2 |
ENDDO |
159 |
|
|
cc ml 20/04/05: |
160 |
pam-ba |
3.4 |
* Positioning volume TRSL into the mother TPAS |
161 |
cafagna |
3.3 |
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one |
162 |
|
|
C # Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2) |
163 |
cafagna |
3.2 |
C Y=2.75 |
164 |
|
|
CC Z= 0. |
165 |
pam-ba |
3.4 |
C Z=-0.005 |
166 |
|
|
c Z=-TRPB(3)+2*TPAI(3)+TRCP(3) |
167 |
|
|
Z=TRPB(3)-2*TPAI(3)-TRCP(3) |
168 |
cafagna |
3.2 |
DO I=1, NTRSL |
169 |
pam-ba |
3.4 |
c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders. |
170 |
|
|
c X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+ |
171 |
|
|
c + (2*(MOD(I-1,3)+1)-1)*TRSL(1) |
172 |
cafagna |
3.3 |
X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+ |
173 |
pam-ba |
3.4 |
+ (2*(MOD(I-1,3)+1)-1)*0.0015+(2*(MOD(I-1,3)+1)-1)*TRSL(1) |
174 |
cafagna |
3.3 |
IF(I.LE.3) THEN |
175 |
pam-ba |
3.4 |
C ml: 22/9/05: |
176 |
|
|
c Y=TPAS(2)-TRSL(2) |
177 |
|
|
Y=TPAS(2)-0.1485-TRSL(2) |
178 |
cafagna |
3.3 |
ELSE |
179 |
pam-ba |
3.4 |
C Y=TPAS(2)-3*TRSL(2) |
180 |
|
|
Y=TPAS(2)-0.1485-0.0015-3*TRSL(2) |
181 |
cafagna |
3.3 |
ENDIF |
182 |
|
|
CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY') |
183 |
cafagna |
3.2 |
ENDDO |
184 |
cafagna |
3.3 |
C # end CAF mod |
185 |
cafagna |
3.2 |
********** |
186 |
|
|
cc ml: 20/04/05 |
187 |
pam-ba |
3.4 |
* Positioning volume THBP into the mother TPAS |
188 |
cafagna |
3.3 |
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one |
189 |
|
|
C # Y=-TRPB(2)+1.5+THBP(2) |
190 |
|
|
C # cc Z= 0. |
191 |
|
|
C # Z=-TRPB(3)+0.11+TRCP(3) |
192 |
pam-ba |
3.4 |
c ml: 22/9/05: |
193 |
|
|
c Y=-TPAS(2)+THBP(2) |
194 |
|
|
Y=-TPAS(2)+0.1485+THBP(2) |
195 |
|
|
C Z=-0.005 |
196 |
|
|
c Z=-TRPB(3)+2*TPAI(3)+TRCP(3) |
197 |
|
|
Z=TRPB(3)-2*TPAI(3)-TRCP(3) |
198 |
cafagna |
3.2 |
DO I=1, NTHBP |
199 |
|
|
CC X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1) |
200 |
cafagna |
3.3 |
C # X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1) |
201 |
pam-ba |
3.4 |
c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon |
202 |
|
|
c ladders |
203 |
|
|
c X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1) |
204 |
|
|
X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)+ |
205 |
|
|
+ (2*(MOD(I-1,3)+1)-1)*0.0015 |
206 |
cafagna |
3.3 |
CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY') |
207 |
cafagna |
3.2 |
ENDDO |
208 |
|
|
c end ml. |
209 |
|
|
cc ml: 20/04/05 |
210 |
cafagna |
3.1 |
* |
211 |
pam-ba |
3.4 |
* Positioning volume TPAI into the mother TPAS |
212 |
cafagna |
3.1 |
* |
213 |
cafagna |
3.3 |
C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS |
214 |
|
|
X=0 |
215 |
|
|
Y=-TPAS(2)+TPAI(2) |
216 |
pam-ba |
3.4 |
c Z=-TPAS(3)+TPAI(3) |
217 |
|
|
Z=TRPB(3)-TPAI(3) |
218 |
cafagna |
3.3 |
C # X=-TRPB(1)+2.29+TPAI(1) |
219 |
|
|
C # Y=-TRPB(2)+1.5+2.5+TPAI(2) |
220 |
|
|
C # Z=-TRPB(3)+TPAI(3) |
221 |
|
|
CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY') |
222 |
|
|
C # end CAF mod |
223 |
cafagna |
3.2 |
* |
224 |
|
|
* Positioning volume TPAS into the mother TRPB |
225 |
|
|
* |
226 |
cafagna |
3.3 |
C # Mod By Caf on 15th Jul 2005 |
227 |
|
|
X=0 |
228 |
cafagna |
3.2 |
Y=-TRPB(2)+1.5+TPAS(2) |
229 |
cafagna |
3.3 |
Z=0 |
230 |
|
|
C # X=-TRPB(1)+2.29+TPAS(1) |
231 |
|
|
C # Y=-TRPB(2)+1.5+TPAS(2) |
232 |
|
|
C # Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3) |
233 |
cafagna |
3.2 |
CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY') |
234 |
|
|
c end ml. |
235 |
pam-ba |
3.4 |
* |
236 |
|
|
* Positioning volume MGPA into the mother MGPL |
237 |
|
|
* |
238 |
|
|
N= 1 |
239 |
|
|
X= 0. |
240 |
|
|
Y= 0. |
241 |
|
|
Z= 0. |
242 |
|
|
CALL GSPOS('MGPA',N,'MGPL',X,Y,Z,0,'ONLY') |
243 |
|
|
* |
244 |
|
|
* Positioning volume TPGD into the mother MGFR |
245 |
|
|
* |
246 |
|
|
X=0. |
247 |
|
|
Z=-MGFR(3)+TPGD(3) |
248 |
|
|
Y=0. |
249 |
|
|
CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY') |
250 |
cafagna |
3.1 |
* |
251 |
pam-ba |
3.5 |
* Positioning volume TPGI into the mother MGFR under MGPL |
252 |
|
|
* |
253 |
|
|
N=1 |
254 |
|
|
X=0. |
255 |
|
|
Z=-MGFR(3)+2*TPGD(3)+TPGI(3) |
256 |
|
|
Y=0. |
257 |
|
|
CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY') |
258 |
|
|
* |
259 |
cafagna |
3.1 |
* Positioning volume MGPI into the mother MGPL |
260 |
|
|
* |
261 |
|
|
N= 1 |
262 |
|
|
X= 0. |
263 |
|
|
Y= 0. |
264 |
pam-ba |
3.5 |
CC ML 10/11/05: Positioning volume MGPI into the mother MGPA |
265 |
cafagna |
3.1 |
Z= 0. |
266 |
pam-ba |
3.4 |
CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY') |
267 |
|
|
C Z= -MGFR(3)+2*TPGD(3)+MGPI(3) |
268 |
|
|
C CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY') |
269 |
cafagna |
3.3 |
* |
270 |
|
|
* Positioning volume MGPL into the mother MGFR |
271 |
|
|
* |
272 |
pam-ba |
3.4 |
c ml: 10/11/05 |
273 |
pam-ba |
3.5 |
N=1 |
274 |
pam-ba |
3.4 |
c Z=-MGFR(3)+MGPL(3) |
275 |
pam-ba |
3.5 |
Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3) |
276 |
cafagna |
3.3 |
CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY') |
277 |
pam-ba |
3.4 |
* |
278 |
pam-ba |
3.5 |
* Positioning volume TPGI into the mother MGFR over MGPL |
279 |
|
|
* |
280 |
|
|
N=2 |
281 |
|
|
X=0. |
282 |
|
|
Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3))+TPGI(3) |
283 |
|
|
Y=0. |
284 |
|
|
CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY') |
285 |
|
|
* |
286 |
pam-ba |
3.4 |
* Positioning volume TPGA into the mother MGFR |
287 |
|
|
* |
288 |
|
|
C ML: 10/11/05 |
289 |
pam-ba |
3.5 |
N=1 |
290 |
pam-ba |
3.4 |
X=0. |
291 |
|
|
Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2) |
292 |
|
|
C Z=-MGFR(3)+2*MGPL(3)+TPGA(3) |
293 |
pam-ba |
3.5 |
C Z=MGFR(3)-2*TPGU(3)-2*TRPB(3)-TPGA(3) |
294 |
|
|
Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3))+TPGA(3) |
295 |
pam-ba |
3.4 |
C Y=-0.75 |
296 |
|
|
CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY') |
297 |
cafagna |
3.1 |
* |
298 |
|
|
* Positioning volumes TRPB & MGPL into the mother SPEB |
299 |
|
|
* |
300 |
|
|
N= 1 |
301 |
|
|
X= 0. |
302 |
pam-ba |
3.4 |
C Y= SPEB(2) - TRPB(2) |
303 |
|
|
Y= -MGFR(2) + TRPB(2) |
304 |
|
|
C Z= SPEB(3) - TRPB(3) |
305 |
|
|
c Z=-MGFR(3)+2*MGPL(3)+2*TPGA(3)+TRPB(3) |
306 |
pam-ba |
3.5 |
C Z=MGFR(3)-2*TPGU(3)-TRPB(3) |
307 |
|
|
Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3))+TRPB(3) |
308 |
pam-ba |
3.4 |
C CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY') |
309 |
|
|
CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY') |
310 |
|
|
c CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY') |
311 |
|
|
Z=SPEB(3) |
312 |
cafagna |
3.1 |
DO I=1, 5 |
313 |
cafagna |
3.3 |
X= 0. |
314 |
|
|
Y= 0. |
315 |
pam-ba |
3.4 |
C Z= Z - TRPB(3) - MGFR(3) |
316 |
|
|
Z=Z-MGFR(3) |
317 |
cafagna |
3.1 |
N= I |
318 |
cafagna |
3.3 |
CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY') |
319 |
pam-ba |
3.4 |
C Y= SPEB(2) - TRPB(2) |
320 |
|
|
C Z= Z - MGFR(3) - TRPB(3) |
321 |
|
|
C N= (I+1) |
322 |
|
|
C IF(I.EQ.5) THEN |
323 |
|
|
CC CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY') |
324 |
|
|
C CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,10,'ONLY') |
325 |
|
|
C ELSE |
326 |
|
|
CC CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY') |
327 |
|
|
C CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY') |
328 |
|
|
C ENDIF |
329 |
|
|
Z=Z-MGFR(3) |
330 |
cafagna |
3.1 |
ENDDO |
331 |
pam-ba |
3.4 |
N=2 |
332 |
|
|
X=0. |
333 |
|
|
Y=MGFR(2) - TRPB(2) |
334 |
|
|
Z=Z-TRPB(3) |
335 |
|
|
c CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY') |
336 |
|
|
CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY') |
337 |
|
|
* |
338 |
|
|
* Positioning volume TPGU into the mother MGFR |
339 |
|
|
* |
340 |
|
|
X=0. |
341 |
pam-ba |
3.5 |
C Z=MGFR(3)-TPGU(3) |
342 |
|
|
Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3)+TRPB(3))+ |
343 |
|
|
+ TPGU(3) |
344 |
pam-ba |
3.4 |
Y=-MGFR(2)+TRPB(2) |
345 |
|
|
CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY') |
346 |
cafagna |
3.1 |
* |
347 |
|
|
RETURN |
348 |
|
|
END |