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

Diff of /gpamela/gpspe/gpspev.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.5 by pam-ba, Tue Dec 13 10:31:29 2005 UTC revision 3.10 by cafagna, Fri Dec 1 12:25:51 2006 UTC
# Line 1  Line 1 
   
1  *  *
2  * $Id: gpspev.F,v 3.4 2005/12/05 12:15:21 pam-ba Exp $  * $Id: gpspev.F,v 3.9 2006/06/30 15:38:16 pam-ba Exp $
3  *  *
4  * $Log: gpspev.F,v $  * $Log: gpspev.F,v $
5    * Revision 3.9  2006/06/30 15:38:16  pam-ba
6    * 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)
7    *
8    * Revision 3.8  2006/05/02 12:02:43  bottai
9    * micrometric adjustment of spe geometry
10    *
11    * Revision 3.7  2005/12/16 10:20:23  cafagna
12    * Bug fixed in the new update
13    *
14    * Revision 3.6  2005/12/16 09:13:59  cafagna
15    * New small valume added to the tracker frame
16    *
17    * Revision 3.5  2005/12/13 10:31:29  pam-ba
18    * added a new volume, TPGI, in the spectrometer
19    *
20  * Revision 3.4  2005/12/05 12:15:21  pam-ba  * Revision 3.4  2005/12/05 12:15:21  pam-ba
21  * new spectrometer geometry and internal magnetic field  * new spectrometer geometry and internal magnetic field
22  *  *
# Line 96  C end ml. Line 110  C end ml.
110        NMED= MCF        NMED= MCF
111        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
112  *  *
113    *  Define the TBAL     volume
114    *
115          NMED= MAL
116          CALL GSVOLU('TBAL','BOX ',NMED,TBAL, 3,IVOLU)
117    *
118  *  Define the MGPA     volume  *  Define the MGPA     volume
119  *  *
120        NMED= MAL        NMED= MAL
# Line 124  C end ml. Line 143  C end ml.
143  *  *
144  * Positioning volume TSPA     into the mother TRSL  * Positioning volume TSPA     into the mother TRSL
145  *  *
146          GLUEX=XGLUE-0.0005
147          GLUEY=XGLUE
148        N= 1        N= 1
149        X= 0.        X= 0.
150        Y= 0.        Y= 0.
# Line 141  C #       Y=-TRPB(2)+1.5-0.225+TRCP(2) Line 162  C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
162  CC end ml.  CC end ml.
163  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
164        Z=TRPB(3)-2*TPAI(3)-TRCP(3)        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
165        y=-0.225-TPAS(2)+TRCP(2)        y=-TPAS(2)+TRCP(2)
166  C # end CAF mod  C # end CAF mod
167          TPASSPACE=TPAS(1)-(3*TRSL(1)+4*TRCP(1)+3*GLUEX)
168        DO I=1,4        DO I=1,4
169  CC ml 20/04/05:          CC ml 20/04/05:        
170  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
# Line 151  C #          X=-TRPB(1)+2.29+(2*I-1)*TRC Line 173  C #          X=-TRPB(1)+2.29+(2*I-1)*TRC
173  c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.  c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.
174  c              for now the glue is a space of nitrogen of 0.015 mm  c              for now the glue is a space of nitrogen of 0.015 mm
175  c         X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  c         X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
176           X=-TPAS(1)+2*(I-1)*0.0015+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  C Ml. 15/12/05: introducing a space (along x) between the inner part of the
177    c  aluminum frame and respectively the first and the last carbon fiber barrel.
178    c It is has been done increasing the dimension of TPAS along x:
179    C         X=-TPAS(1)+2*(I-1)*XGLUE+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
180             X=-TPAS(1)+TPASSPACE+2*(I-1)*GLUEX+(2*I-1)*TRCP(1)+
181         +      2*(I-1)*TRSL(1)
182  CC end ml  CC end ml
183           CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')           CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
184  C # end CAF mod  C # end CAF mod
# Line 169  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3) Line 196  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
196  c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.  c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.
197  c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+  c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
198  c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)  c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
199           X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+           X=-TPAS(1)+TPASSPACE+2*(MOD(I-1,3)+1)*TRCP(1)+
200       +   (2*(MOD(I-1,3)+1)-1)*0.0015+(2*(MOD(I-1,3)+1)-1)*TRSL(1)               +   (2*(MOD(I-1,3)+1)-1)*GLUEX+(2*(MOD(I-1,3)+1)-1)*TRSL(1)
201           IF(I.LE.3) THEN           IF(I.LE.3) THEN
202  C ml: 22/9/05:  C ml: 22/9/05:
203  c            Y=TPAS(2)-TRSL(2)  c            Y=TPAS(2)-TRSL(2)
204              Y=TPAS(2)-0.1485-TRSL(2)              Y=TPAS(2)-0.1485+0.0005-TRSL(2)
205           ELSE           ELSE
206  C            Y=TPAS(2)-3*TRSL(2)  C            Y=TPAS(2)-3*TRSL(2)
207              Y=TPAS(2)-0.1485-0.0015-3*TRSL(2)              Y=TPAS(2)-0.1485-GLUEY-3*TRSL(2)
208           ENDIF           ENDIF
209           CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')           CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
210        ENDDO        ENDDO
# Line 201  C #          X=-TRPB(1)+2.29+2*I*TRCP(1) Line 228  C #          X=-TRPB(1)+2.29+2*I*TRCP(1)
228  c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon  c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon
229  c       ladders  c       ladders
230  c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)  c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)
231           X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)+           X=-TPAS(1)+TPASSPACE+2*I*TRCP(1)+(2*I-1)*THBP(1)+
232       +      (2*(MOD(I-1,3)+1)-1)*0.0015       +      (2*(MOD(I-1,3)+1)-1)*GLUEX
233           CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')           CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
234        ENDDO        ENDDO
235  c end ml.  c end ml.
# Line 248  c end ml. Line 275  c end ml.
275        Y=0.        Y=0.
276        CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')        CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
277  *  *
278  * Positioning volume TPGI     into the mother MGFR under MGPL  * Positioning volume TPGI     into the mother MGFR under MGPL
279  *  *
280        N=1        N=1
281        X=0.        X=0.
# Line 275  c      Z=-MGFR(3)+MGPL(3) Line 302  c      Z=-MGFR(3)+MGPL(3)
302        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3)        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3)
303        CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')        CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
304  *  *
305  * Positioning volume TPGI     into the mother MGFR over MGPL  * Positioning volume TPGI     into the mother MGFR above MGPL
306  *  *
307        N=2        N=2
308        X=0.        X=0.
# Line 332  C         ENDIF Line 359  C         ENDIF
359        X=0.        X=0.
360        Y=MGFR(2) - TRPB(2)        Y=MGFR(2) - TRPB(2)
361        Z=Z-TRPB(3)        Z=Z-TRPB(3)
362          ZTRPB=Z
363  c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')  c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')
364        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
365  *  *
366    * Positioning volume TBAL into the mother SPEB
367    *
368          DO I=1,2
369    C DTBAL is the distance along x of the volumes TBAL from the end of TRPB
370             DTBAL=3.6
371             N=I
372             IF(I.EQ.1) THEN
373                X=TRPB(1)-3.6-TBAL(1)
374             ELSE
375                X=-(TRPB(1)-3.6-TBAL(1))
376             ENDIF
377             Y=MGFR(2) - 2*TRPB(2)-TBAL(2)
378             Z=ZTRPB
379             CALL GSPOS('TBAL',N,'SPEB',X,Y,Z,0,'ONLY')
380          ENDDO
381    *
382  * Positioning volume TPGU     into the mother MGFR  * Positioning volume TPGU     into the mother MGFR
383  *  *
384        X=0.        X=0.
385  C      Z=MGFR(3)-TPGU(3)  C      Z=MGFR(3)-TPGU(3)
386        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3)+TRPB(3))+        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3)+TPGA(3)+TRPB(3))+
387       +   TPGU(3)       +   TPGU(3)      
388        Y=-MGFR(2)+TRPB(2)        Y=-MGFR(2)+TRPB(2)
389        CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')        CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
390  *  *
Legend:
Removed from v.3.5  
changed lines
  Added in v.3.10
  ViewVC Help
Powered by ViewVC 1.1.23