/[PAMELA software]/gpamela/gpspe/gpspev.F
ViewVC logotype

Annotation of /gpamela/gpspe/gpspev.F

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3.5 - (hide annotations) (download)
Tue Dec 13 10:31:29 2005 UTC (18 years, 11 months ago) by pam-ba
Branch: MAIN
Changes since 3.4: +36 -6 lines
added a new volume, TPGI, in the spectrometer

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

  ViewVC Help
Powered by ViewVC 1.1.23