/[PAMELA software]/gpamela/gpcat/gpcatv.F
ViewVC logotype

Annotation of /gpamela/gpcat/gpcatv.F

Parent Directory Parent Directory | Revision Log Revision Log

Revision 3.2 - (hide annotations) (download)
Thu Dec 5 10:17:42 2002 UTC (23 years, 3 months ago) by pamela
Branch: MAIN
CVS Tags: v4r4, v4r5, v4r6, v4r7, v4r0, v4r1, v4r2, v4r3, v4r8, v4r9, v4r14, v4r12, v4r13, v4r10, v4r11, v3r3, v3r1, v3r2, HEAD
Changes since 3.1: +420 -75 lines 
Update CAS and CALO geometries and positions. Makefile updated as well
1 cafagna 3.1 *
2 pamela 3.2 * $Id: gpcatv.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
3 cafagna 3.1 *
4 pamela 3.2 * $Log: gpcatv.F,v $
5     * Revision 3.1.1.1 2002/07/11 16:02:14 cafagna
6     * First GPAMELA release on CVS
7 cafagna 3.1 *
8 pamela 3.2 *
9     * 17/10/2002 22.28.20 by Jens Lund
10 cafagna 3.1 *CMZ : 3.00/00 14/11/2000 09.13.05 by Emiliano Mocchiutti
11     *CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
12     *CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna
13     *CMZ : 1.00/02 08/02/96 17.12.36 by Francesco Cafagna
14     *-- Author : Francesco Cafagna 08/02/96
15     SUBROUTINE GPCATV
16     ************************************************************************
17     * *
18     * Volumes definition for the TOP anticoincidence counters *
19     * Called by: GPGEO *
20     * *
21     * Author: Francesco Cafagna, 08/02/96 16.58.49 *
22     * *
23     ************************************************************************
24     #include "gpgeo.inc"
25     #include "gpmed.inc"
26     *
27     INTEGER IROT,IVOLU,N,NMED,NUM
28     REAL X,Y,Z
29 pamela 3.2 *JeL:
30     REAL TBOX(3),TTRD1(4)
31     NMED=MN2
32 cafagna 3.1 *
33 pamela 3.2 * Define the CAT 'virtual' box
34     CALL GSVOLU('CATA','BOX ',NMED,CATA, 3,IVOLU)
35 cafagna 3.1 *
36 pamela 3.2 * Define the CAT (cover) Aluminum box
37     NMED=MAL
38     TBOX(1) = 0.5*56.0
39     TBOX(2) = 0.5*1.2
40     TBOX(3) = 0.5*54.33
41     CALL GSVOLU('ATBB','BOX ',NMED,TBOX, 3,IVOLU)
42     *
43     * Define the CAT rims
44     NMED=MPLAS
45     TBOX(1) = 0.5*55.3
46     TBOX(2) = 0.5*0.8
47     TBOX(3) = 0.5*53.63
48     CALL GSVOLU('PLA2','BOX ',NMED,TBOX, 3,IVOLU)
49 cafagna 3.1 *
50 pamela 3.2 * Define the CAT scintillator
51     NMED=MSCIN
52     TBOX(1) = 0.5*53.7
53     TBOX(2) = 0.5*0.8
54     TBOX(3) = 0.5*52.03
55     CALL GSVOLU('TOP1','BOX ',NMED,TBOX, 3,IVOLU)
56     *
57     * Define middle Pl. box
58     NMED=MPLAS
59     TBOX(1) = 0.5*20.22
60     TBOX(2) = 0.5*0.8
61     TBOX(3) = 0.5*16.88
62     CALL GSVOLU('PTBM','BOX ',NMED,TBOX, 3,IVOLU)
63 cafagna 3.1 *
64 pamela 3.2 * Define middle Al. box
65     NMED=MAL
66     TBOX(1) = 0.5*18.82
67     TBOX(2) = 0.5*0.8
68     TBOX(3) = 0.5*15.48
69     CALL GSVOLU('ATBM','BOX ',NMED,TBOX, 3,IVOLU)
70 cafagna 3.1 *
71 pamela 3.2 * Define middle 'air' box
72 cafagna 3.1 NMED=MN2
73 pamela 3.2 TBOX(1) = 0.5*18.22
74     TBOX(2) = 0.5*1.2
75     TBOX(3) = 0.5*14.88
76     CALL GSVOLU('LTBM','BOX ',NMED,TBOX, 3,IVOLU)
77     *
78     * Define corner Pl. box
79     NMED=MPLAS
80     TBOX(1) = 0.5*9.95
81     TBOX(2) = 0.5*0.8
82     TBOX(3) = 0.5*11.115
83     CALL GSVOLU('PTBC','BOX ',NMED,TBOX, 3,IVOLU)
84 cafagna 3.1 *
85 pamela 3.2 * Define corner Al. box
86     NMED=MAL
87     TBOX(1) = 0.5*9.95
88     TBOX(2) = 0.5*0.8
89     TBOX(3) = 0.5*11.115
90     CALL GSVOLU('ATBC','BOX ',NMED,TBOX, 3,IVOLU)
91 cafagna 3.1 *
92 pamela 3.2 * Define corner 'air' box
93     NMED=MN2
94     TBOX(1) = 0.5*9.95
95     TBOX(2) = 0.5*1.2
96     TBOX(3) = 0.5*11.115
97     CALL GSVOLU('LTBC','BOX ',NMED,TBOX, 3,IVOLU)
98     *
99     * Define Pl. trapezoid (TRD1) 'top'
100     NMED=MPLAS
101     TTRD1(1) = 0.
102     TTRD1(2) = 16.9
103     TTRD1(3) = 0.5*0.8
104     TTRD1(4) = 0.5*4.81
105     CALL GSVOLU('PTRT','TRD1',NMED,TTRD1, 4,IVOLU)
106 cafagna 3.1 *
107 pamela 3.2 * Define Al. trapezoid (TRD1) 'top'
108     NMED=MAL
109     TTRD1(1) = 0.
110     TTRD1(2) = 17.7
111     TTRD1(3) = 0.5*0.8
112     TTRD1(4) = 0.5*4.81
113     CALL GSVOLU('ATRT','TRD1',NMED,TTRD1, 4,IVOLU)
114 cafagna 3.1 *
115 pamela 3.2 * Define 'air' trapezoid (TRD1) 'top'
116     NMED=MN2
117     TTRD1(1) = 0.
118     TTRD1(2) = 18.05
119     TTRD1(3) = 0.5*1.2
120     TTRD1(4) = 0.5*4.81
121     CALL GSVOLU('LTRT','TRD1',NMED,TTRD1, 4,IVOLU)
122     *
123     * Define Pl. trapezoid (TRD1) 'bottom'
124     NMED=MPLAS
125     TTRD1(1) = 16.9
126     TTRD1(2) = 0.
127     TTRD1(3) = 0.5*0.8
128     TTRD1(4) = 0.5*4.81
129     CALL GSVOLU('PTRB','TRD1',NMED,TTRD1, 4,IVOLU)
130 cafagna 3.1 *
131 pamela 3.2 * Define Al. trapezoid (TRD1) 'bottom'
132 cafagna 3.1 NMED=MAL
133 pamela 3.2 TTRD1(1) = 17.7
134     TTRD1(2) = 0.
135     TTRD1(3) = 0.5*0.8
136     TTRD1(4) = 0.5*4.81
137     CALL GSVOLU('ATRB','TRD1',NMED,TTRD1, 4,IVOLU)
138 cafagna 3.1 *
139 pamela 3.2 * Define 'air' trapezoid (TRD1) 'bottom'
140     NMED=MN2
141     TTRD1(1) = 18.05
142     TTRD1(2) = 0.
143     TTRD1(3) = 0.5*1.2
144     TTRD1(4) = 0.5*4.81
145     CALL GSVOLU('LTRB','TRD1',NMED,TTRD1, 4,IVOLU)
146     *
147     * Define Pl. trapezoid (TRD1) 'left&right'
148     NMED=MPLAS
149     TTRD1(1) = 0.
150     TTRD1(2) = 14.9
151     TTRD1(3) = 0.5*0.8
152     TTRD1(4) = 0.5*4.75
153     CALL GSVOLU('PTRL','TRD1',NMED,TTRD1, 4,IVOLU)
154 cafagna 3.1 *
155 pamela 3.2 * Define Al. trapezoid (TRD1) 'left&right'
156     NMED=MAL
157     TTRD1(1) = 0.
158     TTRD1(2) = 15.7
159     TTRD1(3) = 0.5*0.8
160     TTRD1(4) = 0.5*4.75
161     CALL GSVOLU('ATRL','TRD1',NMED,TTRD1, 4,IVOLU)
162 cafagna 3.1 *
163 pamela 3.2 * Define 'air' trapezoid (TRD1) 'left&right'
164     NMED=MN2
165     TTRD1(1) = 0.
166     TTRD1(2) = 16.05
167     TTRD1(3) = 0.5*1.2
168     TTRD1(4) = 0.5*4.75
169     CALL GSVOLU('LTRL','TRD1',NMED,TTRD1, 4,IVOLU)
170 cafagna 3.1 *
171     *
172 pamela 3.2 * POSITIONING OF THE CAT VOLUMES
173 cafagna 3.1 *
174     *
175 pamela 3.2 * Position the aluminum box (ATBB) in CATA
176     N = 1
177     X = 0.
178     Y = 0.
179     Z = 0.
180     CALL GSPOS('ATBB',N,'CATA',X,Y,Z,0,'ONLY')
181 cafagna 3.1 *
182 pamela 3.2 * Position the plastic box (PLA2) in ATBB
183 cafagna 3.1 N = 1
184 pamela 3.2 X = 0.
185     Y = 0.
186     Z = 0.
187     CALL GSPOS('PLA2',N,'ATBB',X,Y,Z,0,'ONLY')
188     *
189     * Position the CAT scintillator (TOP1) in PLA2
190     X = 0.
191 cafagna 3.1 Y = 0.
192     Z = 0.
193 pamela 3.2 CALL GSPOS('TOP1',N,'PLA2',X,Y,Z,0,'ONLY')
194     *
195     * Position the middle plastic box (PTBM) in TOP1
196     X = 0.
197     Y = 0.
198 cafagna 3.1 Z = 0.
199 pamela 3.2 CALL GSPOS('PTBM',N,'TOP1',X,Y,Z,0,'ONLY')
200 cafagna 3.1 *
201 pamela 3.2 * Position the middle aluminum box (ATBM) in PLA2
202     X = 0.
203     Y = 0.
204     Z = 0.
205     CALL GSPOS('ATBM',N,'PLA2',X,Y,Z,0,'ONLY')
206 cafagna 3.1 *
207 pamela 3.2 * Position the middle 'air' box (LTBM) in ATBB
208     X = 0.
209 cafagna 3.1 Y = 0.
210     Z = 0.
211 pamela 3.2 CALL GSPOS('LTBM',N,'ATBB',X,Y,Z,0,'ONLY')
212     *
213     * Position the plastic U.L corner box (PTBC) in TOP1
214     X = -21.875
215     Y = 0.
216     Z = 20.4575
217     CALL GSPOS('PTBC',N,'TOP1',X,Y,Z,0,'ONLY')
218     *
219     * Position the aluminum U.L corner box (ATBC) in PLA2
220     X = -22.675
221     Y = 0.
222     Z = 21.2575
223     CALL GSPOS('ATBC',N,'PLA2',X,Y,Z,0,'ONLY')
224     *
225     * Position the 'air' U.L corner box (LTBC) in ATBB
226     X = -23.025
227     Y = 0.
228     Z = 21.6075
229     CALL GSPOS('LTBC',N,'ATBB',X,Y,Z,0,'ONLY')
230     *
231     * Position the plastic U.R corner box (PTBC) in TOP1
232     N = 2
233     X = 21.875
234     Y = 0.
235     Z = 20.4575
236     CALL GSPOS('PTBC',N,'TOP1',X,Y,Z,0,'ONLY')
237     *
238     * Position the aluminum U.R corner box (ATBC) in PLA2
239 cafagna 3.1 N = 2
240 pamela 3.2 X = 22.675
241     Y = 0.
242     Z = 21.2575
243     CALL GSPOS('ATBC',N,'PLA2',X,Y,Z,0,'ONLY')
244     *
245     * Position the 'air' U.R corner box (LTBC) in ATBB
246     X = 23.025
247     Y = 0.
248     Z = 21.6075
249     CALL GSPOS('LTBC',N,'ATBB',X,Y,Z,0,'ONLY')
250     *
251     * Position the plastic L.L corner box (PTBC) in TOP1
252     N = 3
253     X = -21.875
254     Y = 0.
255     Z = -20.4575
256     CALL GSPOS('PTBC',N,'TOP1',X,Y,Z,0,'ONLY')
257     *
258     * Position the aluminum L.L corner box (ATBC) in PLA2
259     N = 3
260     X = -22.675
261     Y = 0.
262     Z = -21.2575
263     CALL GSPOS('ATBC',N,'PLA2',X,Y,Z,0,'ONLY')
264     *
265     * Position the 'air' L.L corner box (LTBC) in ATBB
266     X = -23.025
267     Y = 0.
268     Z = -21.6075
269     CALL GSPOS('LTBC',N,'ATBB',X,Y,Z,0,'ONLY')
270     *
271     * Position the plastic L.R corner box (PTBC) in TOP1
272     N = 4
273     X = 21.875
274     Y = 0.
275     Z = -20.4575
276     CALL GSPOS('PTBC',N,'TOP1',X,Y,Z,0,'ONLY')
277     *
278     * Position the aluminum L.R corner box (ATBC) in PLA2
279     N = 4
280     X = 22.675
281     Y = 0.
282     Z = -21.2575
283     CALL GSPOS('ATBC',N,'PLA2',X,Y,Z,0,'ONLY')
284     *
285     * Position the 'air' L.R corner box (LTBC) in ATBB
286     X = 23.025
287     Y = 0.
288     Z = -21.6075
289     CALL GSPOS('LTBC',N,'ATBB',X,Y,Z,0,'ONLY')
290 cafagna 3.1 *
291 pamela 3.2 * Position the plastic 'top' TRD1 (PTRT) in TOP1
292     N = 1
293     X = 0.
294     Y = 0.
295     Z = 23.61
296     CALL GSPOS('PTRT',N,'TOP1',X,Y,Z,0,'ONLY')
297 cafagna 3.1 *
298 pamela 3.2 * Position the aluminum 'top' TRD1 (ATRT) in PLA2
299 cafagna 3.1 N = 1
300     X = 0.
301     Y = 0.
302 pamela 3.2 Z = 24.41
303     CALL GSPOS('ATRT',N,'PLA2',X,Y,Z,0,'ONLY')
304     *
305     * Position the 'air' 'top' TRD1 (LTRT) in ATBB
306     X = 0.
307     Y = 0.
308     Z = 24.76
309     CALL GSPOS('LTRT',N,'ATBB',X,Y,Z,0,'ONLY')
310     *
311     * Position the plastic 'bottom' TRD1 (PTRB) in TOP1
312     X = 0.
313     Y = 0.
314     Z = -23.61
315     CALL GSPOS('PTRB',N,'TOP1',X,Y,Z,0,'ONLY')
316     *
317     * Position the aluminum 'bottom' TRD1 (ATRB) in PLA2
318     X = 0.
319     Y = 0.
320     Z = -24.41
321     CALL GSPOS('ATRB',N,'PLA2',X,Y,Z,0,'ONLY')
322     *
323     * Position the 'air' 'bottom' TRD1 (LTRB) in ATBB
324     X = 0.
325     Y = 0.
326     Z = -24.76
327     CALL GSPOS('LTRB',N,'ATBB',X,Y,Z,0,'ONLY')
328     *
329     * Position the plastic 'left' TRD1 (PTRL) in TOP1 and rotate
330     N = 1
331     X = -24.475
332     Y = 0.
333 cafagna 3.1 Z = 0.
334 pamela 3.2 CALL GSPOS('PTRL',N,'TOP1',X,Y,Z,7,'ONLY')
335 cafagna 3.1 *
336 pamela 3.2 * Position the aluminum 'left' TRD1 (ATRL) in PLA2 and rotate
337     N = 1
338     X = -25.275
339     Y = 0.
340     Z = 0.
341     CALL GSPOS('ATRL',N,'PLA2',X,Y,Z,7,'ONLY')
342 cafagna 3.1 *
343 pamela 3.2 * Position the 'air' 'left' TRD1 (LTRL) in ATBB and rotate
344     X = -25.625
345     Y = 0.
346 cafagna 3.1 Z = 0.
347 pamela 3.2 CALL GSPOS('LTRL',N,'ATBB',X,Y,Z,7,'ONLY')
348     *
349     * Position the plastic 'right' TRD1 (PTRL) in TOP1 and rotate
350 cafagna 3.1 N = 2
351 pamela 3.2 X = 24.475
352     Y = 0.
353 cafagna 3.1 Z = 0.
354 pamela 3.2 CALL GSPOS('PTRL',N,'TOP1',X,Y,Z,8,'ONLY')
355     *
356     * Position the aluminum 'right' TRD1 (ATRL) in PLA2 and rotate
357     N = 2
358     X = 25.275
359     Y = 0.
360 cafagna 3.1 Z = 0.
361 pamela 3.2 CALL GSPOS('ATRL',N,'PLA2',X,Y,Z,8,'ONLY')
362     *
363     * Position the 'air' 'right' TRD1 (LTRL) in ATBB and rotate
364     X = 25.625
365     Y = 0.
366 cafagna 3.1 Z = 0.
367 pamela 3.2 CALL GSPOS('LTRL',N,'ATBB',X,Y,Z,8,'ONLY')
368     *END: JeL.
369    
370     *EM:
371     C NMED=MSCIN
372     *
373     * Define the CATL volume
374     *
375     C CALL GSVOLU('CATL','BOX ',NMED,CATL, 3,IVOLU)
376     *
377     * Define the CATT volume
378     *
379     C CALL GSVOLU('CATT','BOX ',NMED,CATT, 3,IVOLU)
380     *
381     C NMED=MN2
382     *
383     * Define the CATH volume
384     *
385     C CALL GSVOLU('CATH','BOX ',NMED,CATH, 3,IVOLU)
386     *
387     * Define the CATF volume
388     *
389     C CALL GSVOLU('CATF','BOX ',NMED,CATF, 3,IVOLU)
390     *
391     C NMED=MAL
392     *
393     * Define the CATB volume
394     *
395     C CALL GSVOLU('CATB','BOX ',NMED,CATB, 3,IVOLU)
396     *
397     * Define the CATO volume
398     *
399     C CALL GSVOLU('CATO','BOX ',NMED,CATO, 3,IVOLU)
400     *
401     * Define the CATP volume
402     *
403     C CALL GSVOLU('CATP','BOX ',NMED,CATP, 3,IVOLU)
404     *
405     * Positioning volume CATT into the mother CATB
406     *
407     C N = 1
408     C X= CATL(1)-CATB(1)
409     C Y= CATL(2)+(2.*CATT(2)-CATB(2))
410     C X = -CATT(1)
411     C Y = 0.
412     C Z = 0.
413     C CALL GSPOS('CATT',N,'CATB',X,Y,Z,0,'ONLY')
414     C N = 2
415     C X = -X
416     C Y = -Y
417     C Z = 0.
418     C CALL GSPOS('CATT',N,'CATB',X,Y,Z,0,'ONLY')
419     *
420     * Positioning volume CATL into the mother CATB
421     *
422     C N = 1
423     C X = -CATT(1)-(2.*CATL(1)-CATB(1))
424     C Y = CATT(2)-CATB(2)
425     C X = -(2*CATT(1)+CATL(1))
426     C Y = 0.
427     C Z = 0.
428     C CALL GSPOS('CATL',N,'CATB',X,Y,Z,0,'ONLY')
429     C N = 2
430     C X = -X
431     C Y = -Y
432     C Z = 0.
433     C CALL GSPOS('CATL',N,'CATB',X,Y,Z,0,'ONLY')
434     *
435     * Positioning volume CATH into the mother CATB
436     *
437     C N = 1
438     C X = 0.
439     C Y = 0.
440     C Z = 0.
441     C CALL GSPOS('CATH',N,'CATB',X,Y,Z,0,'ONLY')
442     *
443     * Positioning volume CATF into the mother CATB
444     *
445     C N = 1
446     C X = -(2*CATT(1)+CATL(1))
447     C Y = CATL(2)+CATF(2)
448     C Z = 0.
449     C CALL GSPOS('CATF',N,'CATB',X,Y,Z,0,'ONLY')
450     C N = 2
451     C X = X
452     C Y = -Y
453     C Z = 0.
454     C CALL GSPOS('CATF',N,'CATB',X,Y,Z,0,'ONLY')
455     C N = 3
456     C X = -X
457     C Y = Y
458     C Z = 0.
459     C CALL GSPOS('CATF',N,'CATB',X,Y,Z,0,'ONLY')
460     C N = 4
461     C X = X
462     C Y = -Y
463     C Z = 0.
464     C CALL GSPOS('CATF',N,'CATB',X,Y,Z,0,'ONLY')
465 cafagna 3.1
466     *
467     * Positioning volume CATP into the mother CATB
468     *
469 pamela 3.2 C N = 1
470     C X = -CATH(1)+CATP(1)
471     C Y = 0.
472     C Z = 0.
473     C CALL GSPOS('CATP',N,'CATB',X,Y,Z,0,'ONLY')
474     C N = 2
475     C X = -X
476     C Y = Y
477     C Z = 0.
478     C CALL GSPOS('CATP',N,'CATB',X,Y,Z,0,'ONLY')
479 cafagna 3.1 *
480     * Positioning volume CATO into the mother CATB
481     *
482 pamela 3.2 C N = 1
483     C X = 0.
484     C Y = -CATH(2)+CATO(2)
485     C Z = 0.
486     C CALL GSPOS('CATO',N,'CATB',X,Y,Z,0,'ONLY')
487     C N = 2
488     C X = X
489     C Y = -Y
490     C Z = 0.
491     C CALL GSPOS('CATO',N,'CATB',X,Y,Z,0,'ONLY')
492 cafagna 3.1 *
493     RETURN
494     END
  ViewVC Help
Powered by ViewVC 1.1.23