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

Annotation of /gpamela/gpspe/gpspev.F

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3.9 - (hide annotations) (download)
Fri Jun 30 15:38:16 2006 UTC (18 years, 5 months ago) by pam-ba
Branch: MAIN
CVS Tags: v4r7, v4r8
Changes since 3.8: +13 -12 lines
S22 and S12 heights positioned in GPAMELA at the nominal heights in PAMELA (see document: Main geometrical parameters of the PAMELA sub-detectors, 20 December 2005)

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

  ViewVC Help
Powered by ViewVC 1.1.23