/[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.10 by cafagna, Fri Dec 1 12:25:51 2006 UTC revision 3.11 by pam-fi, Fri Apr 23 10:04:13 2010 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpspev.F,v 3.9 2006/06/30 15:38:16 pam-ba Exp $  * $Id: gpspev.F,v 3.10 2006/12/01 12:25:51 cafagna Exp $
3  *  *
4  * $Log: gpspev.F,v $  * $Log: gpspev.F,v $
5    * Revision 3.10  2006/12/01 12:25:51  cafagna
6    *  New Antiproton annihilation cross section added
7    *
8  * Revision 3.9  2006/06/30 15:38:16  pam-ba  * Revision 3.9  2006/06/30 15:38:16  pam-ba
9  * 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)  * 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)
10  *  *
# Line 48  Line 51 
51  *  *
52        INTEGER IROT,IVOLU,N,NMED,NUM        INTEGER IROT,IVOLU,N,NMED,NUM
53        REAL X,Y,Z        REAL X,Y,Z
54    C      
55          REAL ROWMAT(3)
56          REAL OMEGA(36),BETA(36),GAMMA(36)
57          REAL DX(36),DY(36),DZ(36)
58          REAL DXFLO(6),DYFLO(6),DZFLO(6)
59          CHARACTER*19 FILEALLI
60          LOGICAL ROTATE
61    c
62    c
63    
64    cv
65    cv
66          DO IK=1,6
67    
68          IF(IK.EQ.1)  FILEALLI='parameters_l1s2.dat'
69          IF(IK.EQ.2)  FILEALLI='parameters_l2s2.dat'
70          IF(IK.EQ.3)  FILEALLI='parameters_l3s2.dat'
71          IF(IK.EQ.4)  FILEALLI='parameters_l1s1.dat'
72          IF(IK.EQ.5)  FILEALLI='parameters_l2s1.dat'
73          IF(IK.EQ.6)  FILEALLI='parameters_l3s1.dat'
74          
75          write(6,*) 'OPEN ',FILEALLI
76          OPEN(91,file=FILEALLI,status='old')
77          DO IH=1,6
78            II=(6-IH)*6+IK
79            READ(91,*) OMEGA(II)
80            READ(91,*) BETA(II)
81            READ(91,*) GAMMA(II)
82            READ(91,*) DX(II)
83            READ(91,*) DY(II)
84            READ(91,*) DZ(II)
85          ENDDO
86          
87          CLOSE(91)
88          
89          ENDDO
90        
91          DO IJ=1,6
92           DXFLO(IJ)=0.
93           DYFLO(IJ)=0.
94           DO IL=1,6
95            II=(IJ-1)*6+IL
96            DXFLO(IJ)=DXFLO(IJ)+DX(II)
97            DYFLO(IJ)=DYFLO(IJ)+DY(II)
98           ENDDO
99           DXFLO(IJ)=DXFLO(IJ)/6.
100           DYFLO(IJ)=DYFLO(IJ)/6.
101          ENDDO
102          
103          DO IJ=1,6
104           DO IL=1,6
105            II=(IJ-1)*6+IL
106            DX(II)=( DX(II)-DXFLO(IJ) )*1.E-4
107            DY(II)=( DY(II)-DYFLO(IJ) )*1.E-4
108            DZ(II)=(DZ(II)-100.)*1.E-4
109           ENDDO
110           DXFLO(IJ)=DXFLO(IJ)*1.E-4
111           DYFLO(IJ)=DYFLO(IJ)*1.E-4
112          ENDDO
113    
114          
115          DO IJ=1,6
116           DO IL=1,6
117           II=(IJ-1)*6+IL      
118           OMEGA_1=OMEGA(II)*1.E-6
119           BETA_1=BETA(II)*1.E-6
120           GAMMA_1=GAMMA(II)*1.E-6
121    c
122    c
123           IF(IJ.EQ.6) THEN
124            OMEGA_1=-OMEGA_1
125            BETA_1=BETA_1
126            GAMMA_1=-GAMMA_1
127           ENDIF
128    
129    
130    
131           ROWMAT(1)=1.
132           ROWMAT(2)=OMEGA_1
133           ROWMAT(3)=-GAMMA_1
134           CALL GFANG(ROWMAT,COSTH,SINTH,COSPH,SINPH,ROTATE)
135           THETA1=ATAN2(SINTH,COSTH)/acos(-1.)*180.
136           PHI1=ATAN2(SINPH,COSPH)/acos(-1.)*180.
137        
138           ROWMAT(1)=-OMEGA_1
139           ROWMAT(2)=1.
140           ROWMAT(3)=BETA_1
141    
142           CALL GFANG(ROWMAT,COSTH,SINTH,COSPH,SINPH,ROTATE)
143           THETA2=ATAN2(SINTH,COSTH)/acos(-1.)*180.
144           PHI2=ATAN2(SINPH,COSPH)/acos(-1.)*180.
145    
146           ROWMAT(1)=GAMMA_1
147           ROWMAT(2)=-BETA_1
148           ROWMAT(3)=1.
149    
150           CALL GFANG(ROWMAT,COSTH,SINTH,COSPH,SINPH,ROTATE)
151           THETA3=ATAN2(SINTH,COSTH)/acos(-1.)*180.
152           PHI3=ATAN2(SINPH,COSPH)/acos(-1.)*180.
153    
154           write(6,*) 'GSROTM(50+II)',50+II,THETA1,PHI1,THETA2,PHI2,THETA3,PHI3
155           IIKROT=50+II
156           CALL GSROTM(IIKROT,THETA1,PHI1,THETA2,PHI2,THETA3,PHI3)
157          
158           ENDDO
159          ENDDO
160    
161  *  *
162  *  Define the SPEB     volume  *  Define the SPEB     volume
163  *  *
# Line 57  Line 167 
167  *  Define the MGFR     volume  *  Define the MGFR     volume
168  *  *
169        NMED= MAL        NMED= MAL
170        CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)  cv      CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
171          CALL GSVOLU('MGF1','BOX ',NMED,MGFR, 3,IVOLU)
172          CALL GSVOLU('MGF2','BOX ',NMED,MGFR, 3,IVOLU)
173          CALL GSVOLU('MGF3','BOX ',NMED,MGFR, 3,IVOLU)
174          CALL GSVOLU('MGF4','BOX ',NMED,MGFR, 3,IVOLU)
175          CALL GSVOLU('MGF5','BOX ',NMED,MGFR, 3,IVOLU)
176        
177  *  *
178  *  Define the MGPL     volume  *  Define the MGPL     volume
179  *  *
# Line 72  Line 188 
188  *  Define the TRPB     volume  *  Define the TRPB     volume
189  *  *
190        NMED= MAL        NMED= MAL
191        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)  cv      CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
192          CALL GSVOLU('TRP1','BOX ',NMED,TRPB, 3,IVOLU)
193          CALL GSVOLU('TRP2','BOX ',NMED,TRPB, 3,IVOLU)
194          CALL GSVOLU('TRP3','BOX ',NMED,TRPB, 3,IVOLU)
195          CALL GSVOLU('TRP4','BOX ',NMED,TRPB, 3,IVOLU)
196          CALL GSVOLU('TRP5','BOX ',NMED,TRPB, 3,IVOLU)
197          CALL GSVOLU('TRP6','BOX ',NMED,TRPB, 3,IVOLU)            
198  *****  *****
199  cC ml 20/04/05:  cC ml 20/04/05:
200  *  *
# Line 80  cc*  Define the TRPL     volume Line 202  cc*  Define the TRPL     volume
202  *  Define the TPAS and TPAI volume  *  Define the TPAS and TPAI volume
203  *  *
204        NMED= MN2        NMED= MN2
205  CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)  cv      CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
206        CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)        CALL GSVOLU('TPA1','BOX ',NMED,TPAS, 3,IVOLU)
207          CALL GSVOLU('TPA2','BOX ',NMED,TPAS, 3,IVOLU)
208          CALL GSVOLU('TPA3','BOX ',NMED,TPAS, 3,IVOLU)
209          CALL GSVOLU('TPA4','BOX ',NMED,TPAS, 3,IVOLU)
210          CALL GSVOLU('TPA5','BOX ',NMED,TPAS, 3,IVOLU)
211          CALL GSVOLU('TPA6','BOX ',NMED,TPAS, 3,IVOLU)
212        NMED= MAL        NMED= MAL
213        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
214  CC end ml.  CC end ml.
# Line 150  C end ml. Line 277  C end ml.
277        Y= 0.        Y= 0.
278        Z= 0.        Z= 0.
279        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
280  C #       N= 2  
 C # C      Y= Y - 2.*TSPA(2)  
 C #       Y=-TRSL(2)/2  
 C #       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  
281  c ml: 26/11/04:  c ml: 26/11/04:
282  * Positioning volume TRCP into mother TPAS  * Positioning volume TRCP into mother TPAS
283  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
284  C #       Z=-TRPB(3)+0.11+TRCP(3)  
 C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)  
285  CC end ml.  CC end ml.
286  c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)  
287        Z=TRPB(3)-2*TPAI(3)-TRCP(3)        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
288        y=-TPAS(2)+TRCP(2)        y=-TPAS(2)+TRCP(2)
289  C # end CAF mod  C # end CAF mod
290        TPASSPACE=TPAS(1)-(3*TRSL(1)+4*TRCP(1)+3*GLUEX)        TPASSPACE=TPAS(1)-(3*TRSL(1)+4*TRCP(1)+3*GLUEX)
291        DO I=1,4        DO I=1,4
292  CC ml 20/04/05:          CC ml 20/04/05:        
 CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  
 C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one  
 C #          X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  
 c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.  
 c              for now the glue is a space of nitrogen of 0.015 mm  
 c         X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  
293  C Ml. 15/12/05: introducing a space (along x) between the inner part of the  C Ml. 15/12/05: introducing a space (along x) between the inner part of the
294  c  aluminum frame and respectively the first and the last carbon fiber barrel.  c  aluminum frame and respectively the first and the last carbon fiber barrel.
295  c It is has been done increasing the dimension of TPAS along x:  c It is has been done increasing the dimension of TPAS along x:
296  C         X=-TPAS(1)+2*(I-1)*XGLUE+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  
297           X=-TPAS(1)+TPASSPACE+2*(I-1)*GLUEX+(2*I-1)*TRCP(1)+           X=-TPAS(1)+TPASSPACE+2*(I-1)*GLUEX+(2*I-1)*TRCP(1)+
298       +      2*(I-1)*TRSL(1)       +      2*(I-1)*TRSL(1)
299  CC end ml  CC end ml
300           CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')  cv         CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'MANY')
301               CALL GSPOS('TRCP',I,'TPA1',X,Y,Z,0,'MANY')
302               CALL GSPOS('TRCP',I,'TPA2',X,Y,Z,0,'MANY')
303               CALL GSPOS('TRCP',I,'TPA3',X,Y,Z,0,'MANY')
304               CALL GSPOS('TRCP',I,'TPA4',X,Y,Z,0,'MANY')
305               CALL GSPOS('TRCP',I,'TPA5',X,Y,Z,0,'MANY')
306               CALL GSPOS('TRCP',I,'TPA6',X,Y,Z,0,'MANY')
307  C # end CAF mod  C # end CAF mod
308        ENDDO          ENDDO  
309  cc ml 20/04/05:  cc ml 20/04/05:
310  * Positioning volume TRSL     into the mother TPAS  * Positioning volume TRSL     into the mother TPAS
311  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
312  C #       Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)  
 C      Y=2.75  
 CC      Z= 0.  
 C      Z=-0.005  
 c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)  
313        Z=TRPB(3)-2*TPAI(3)-TRCP(3)        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
314        DO I=1, NTRSL        DO I=1, NTRSL
315  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.
316  c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+  
 c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)  
317           X=-TPAS(1)+TPASSPACE+2*(MOD(I-1,3)+1)*TRCP(1)+           X=-TPAS(1)+TPASSPACE+2*(MOD(I-1,3)+1)*TRCP(1)+
318       +   (2*(MOD(I-1,3)+1)-1)*GLUEX+(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)
319           IF(I.LE.3) THEN           IF(I.LE.3) THEN
320  C ml: 22/9/05:  C ml: 22/9/05:
321  c            Y=TPAS(2)-TRSL(2)  
322              Y=TPAS(2)-0.1485+0.0005-TRSL(2)              Y=TPAS(2)-0.1485+0.0005-TRSL(2)
323           ELSE           ELSE
324  C            Y=TPAS(2)-3*TRSL(2)  
325              Y=TPAS(2)-0.1485-GLUEY-3*TRSL(2)              Y=TPAS(2)-0.1485-GLUEY-3*TRSL(2)
326           ENDIF           ENDIF
327           CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')  
328    cv        X=X+DX_1
329    cv        Y=Y+DY_1
330    cv        Z=Z+DZ_1
331    cv          CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,51,'ONLY')
332    cv          CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
333               X1=X+DX(I)
334               Y1=Y+DY(I)
335               Z1=Z+DZ(I)
336               IIRR=50+I
337               CALL GSPOS('TRSL',I,'TPA1',X1,Y1,Z1,IIRR,'ONLY')
338    
339               X2=X+DX(NTRSL+I)
340               Y2=Y+DY(NTRSL+I)
341               Z2=Z+DZ(NTRSL+I)
342               IIRR=50+NTRSL+I
343               CALL GSPOS('TRSL',NTRSL+I,'TPA2',X2,Y2,Z2,IIRR,'ONLY')
344    
345               X3=X+DX(2*NTRSL+I)
346               Y3=Y+DY(2*NTRSL+I)
347               Z3=Z+DZ(2*NTRSL+I)
348               IIRR=50+2*NTRSL+I
349               CALL GSPOS('TRSL',2*NTRSL+I,'TPA3',X3,Y3,Z3,IIRR,'ONLY')
350    
351    
352               X4=X+DX(3*NTRSL+I)
353               Y4=Y+DY(3*NTRSL+I)
354               Z4=Z+DZ(3*NTRSL+I)
355               IIRR=50+3*NTRSL+I
356               CALL GSPOS('TRSL',3*NTRSL+I,'TPA4',X4,Y4,Z4,IIRR,'ONLY')
357    
358               X5=X+DX(4*NTRSL+I)
359               Y5=Y+DY(4*NTRSL+I)
360               Z5=Z+DZ(4*NTRSL+I)
361               IIRR=50+4*NTRSL+I
362               CALL GSPOS('TRSL',4*NTRSL+I,'TPA5',X5,Y5,Z5,IIRR,'ONLY')
363    
364               IF(I.LE.3) KK=3+I
365               IF(I.GE.4) KK=I-3
366               X6=X+DX(5*NTRSL+KK)
367               Y6=Y-DY(5*NTRSL+KK)
368               Z6=Z-DZ(5*NTRSL+KK)    
369               IIRR=50+5*NTRSL+KK
370               CALL GSPOS('TRSL',5*NTRSL+I,'TPA6',X6,Y6,Z6,IIRR,'ONLY')
371        ENDDO        ENDDO
372  C # end CAF mod  C # end CAF mod
373  **********  **********
374  cc ml: 20/04/05  cc ml: 20/04/05
375  * Positioning volume THBP     into the mother TPAS  * Positioning volume THBP     into the mother TPAS
376  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
377  C #       Y=-TRPB(2)+1.5+THBP(2)  
 C # cc      Z= 0.  
 C #       Z=-TRPB(3)+0.11+TRCP(3)  
378  c ml: 22/9/05:  c ml: 22/9/05:
379  c      Y=-TPAS(2)+THBP(2)  
380        Y=-TPAS(2)+0.1485+THBP(2)        Y=-TPAS(2)+0.1485+THBP(2)
381  C      Z=-0.005  
 c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)  
382        Z=TRPB(3)-2*TPAI(3)-TRCP(3)        Z=TRPB(3)-2*TPAI(3)-TRCP(3)
383        DO I=1, NTHBP        DO I=1, NTHBP
384  CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)  
 C #          X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)  
385  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
386  c       ladders  c       ladders
387  c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)  
388           X=-TPAS(1)+TPASSPACE+2*I*TRCP(1)+(2*I-1)*THBP(1)+           X=-TPAS(1)+TPASSPACE+2*I*TRCP(1)+(2*I-1)*THBP(1)+
389       +      (2*(MOD(I-1,3)+1)-1)*GLUEX       +      (2*(MOD(I-1,3)+1)-1)*GLUEX
390           CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')  cv         CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'MANY')
391            CALL GSPOS('THBP',I,'TPA1',X,Y,Z,0,'MANY')
392            CALL GSPOS('THBP',I,'TPA2',X,Y,Z,0,'MANY')
393            CALL GSPOS('THBP',I,'TPA3',X,Y,Z,0,'MANY')
394            CALL GSPOS('THBP',I,'TPA4',X,Y,Z,0,'MANY')
395            CALL GSPOS('THBP',I,'TPA5',X,Y,Z,0,'MANY')
396            CALL GSPOS('THBP',I,'TPA6',X,Y,Z,0,'MANY')
397        ENDDO        ENDDO
398  c end ml.  c end ml.
399  cc ml: 20/04/05  cc ml: 20/04/05
# Line 240  cc ml: 20/04/05 Line 403  cc ml: 20/04/05
403  C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS  C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
404        X=0        X=0
405        Y=-TPAS(2)+TPAI(2)        Y=-TPAS(2)+TPAI(2)
 c      Z=-TPAS(3)+TPAI(3)  
406        Z=TRPB(3)-TPAI(3)        Z=TRPB(3)-TPAI(3)
407  C #       X=-TRPB(1)+2.29+TPAI(1)  
408  C #       Y=-TRPB(2)+1.5+2.5+TPAI(2)  cv      CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'MANY')
409  C #       Z=-TRPB(3)+TPAI(3)        CALL GSPOS('TPAI',1,'TPA1',X,Y,Z,0,'MANY')
410        CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')        CALL GSPOS('TPAI',1,'TPA2',X,Y,Z,0,'MANY')
411          CALL GSPOS('TPAI',1,'TPA3',X,Y,Z,0,'MANY')
412          CALL GSPOS('TPAI',1,'TPA4',X,Y,Z,0,'MANY')
413          CALL GSPOS('TPAI',1,'TPA5',X,Y,Z,0,'MANY')
414          CALL GSPOS('TPAI',1,'TPA6',X,Y,Z,0,'MANY')
415  C # end CAF mod  C # end CAF mod
416  *  *
417  * Positioning volume TPAS into the mother TRPB  * Positioning volume TPAS into the mother TRPB
# Line 254  C # Mod By Caf on 15th Jul 2005 Line 420  C # Mod By Caf on 15th Jul 2005
420        X=0        X=0
421        Y=-TRPB(2)+1.5+TPAS(2)        Y=-TRPB(2)+1.5+TPAS(2)
422        Z=0        Z=0
423  C #       X=-TRPB(1)+2.29+TPAS(1)  
424  C #       Y=-TRPB(2)+1.5+TPAS(2)  cv      CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
425  C #       Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)        CALL GSPOS('TPA1',1,'TRP1',X,Y,Z,0,'ONLY')
426        CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')        CALL GSPOS('TPA2',1,'TRP2',X,Y,Z,0,'ONLY')
427          CALL GSPOS('TPA3',1,'TRP3',X,Y,Z,0,'ONLY')
428          CALL GSPOS('TPA4',1,'TRP4',X,Y,Z,0,'ONLY')
429          CALL GSPOS('TPA5',1,'TRP5',X,Y,Z,0,'ONLY')
430          CALL GSPOS('TPA6',1,'TRP6',X,Y,Z,0,'ONLY')
431  c end ml.  c end ml.
432  *  *
433  * Positioning volume MGPA     into the mother MGPL  * Positioning volume MGPA     into the mother MGPL
# Line 273  c end ml. Line 443  c end ml.
443        X=0.        X=0.
444        Z=-MGFR(3)+TPGD(3)        Z=-MGFR(3)+TPGD(3)
445        Y=0.        Y=0.
446        CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
447          CALL GSPOS('TPGD',1,'MGF1',X,Y,Z,0,'ONLY')
448          CALL GSPOS('TPGD',1,'MGF2',X,Y,Z,0,'ONLY')
449          CALL GSPOS('TPGD',1,'MGF3',X,Y,Z,0,'ONLY')
450          CALL GSPOS('TPGD',1,'MGF4',X,Y,Z,0,'ONLY')
451          CALL GSPOS('TPGD',1,'MGF5',X,Y,Z,0,'ONLY')
452  *  *
453  * Positioning volume TPGI     into the mother MGFR under MGPL  * Positioning volume TPGI     into the mother MGFR under MGPL
454  *  *
# Line 281  c end ml. Line 456  c end ml.
456        X=0.        X=0.
457        Z=-MGFR(3)+2*TPGD(3)+TPGI(3)        Z=-MGFR(3)+2*TPGD(3)+TPGI(3)
458        Y=0.        Y=0.
459        CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')
460          CALL GSPOS('TPGI',1,'MGF1',X,Y,Z,0,'ONLY')
461          CALL GSPOS('TPGI',1,'MGF2',X,Y,Z,0,'ONLY')
462          CALL GSPOS('TPGI',1,'MGF3',X,Y,Z,0,'ONLY')
463          CALL GSPOS('TPGI',1,'MGF4',X,Y,Z,0,'ONLY')
464          CALL GSPOS('TPGI',1,'MGF5',X,Y,Z,0,'ONLY')
465  *  *
466  * Positioning volume MGPI     into the mother MGPL  * Positioning volume MGPI     into the mother MGPL
467  *  *
# Line 291  c end ml. Line 471  c end ml.
471  CC ML 10/11/05:  Positioning volume MGPI     into the mother MGPA  CC ML 10/11/05:  Positioning volume MGPI     into the mother MGPA
472        Z= 0.        Z= 0.
473        CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')        CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')
474  C      Z= -MGFR(3)+2*TPGD(3)+MGPI(3)  
 C      CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY')  
475  *  *
476  * Positioning volume MGPL     into the mother MGFR  * Positioning volume MGPL     into the mother MGFR
477  *  *
478  c ml: 10/11/05  c ml: 10/11/05
479        N=1        N=1
 c      Z=-MGFR(3)+MGPL(3)  
480        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3)        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3))+MGPL(3)
481        CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
482          CALL GSPOS('MGPL',1,'MGF1',X,Y,Z,0,'ONLY')    
483          CALL GSPOS('MGPL',1,'MGF2',X,Y,Z,0,'ONLY')
484          CALL GSPOS('MGPL',1,'MGF3',X,Y,Z,0,'ONLY')
485          CALL GSPOS('MGPL',1,'MGF4',X,Y,Z,0,'ONLY')
486          CALL GSPOS('MGPL',1,'MGF5',X,Y,Z,0,'ONLY')
487  *  *
488  * Positioning volume TPGI     into the mother MGFR above MGPL  * Positioning volume TPGI     into the mother MGFR above MGPL
489  *  *
# Line 308  c      Z=-MGFR(3)+MGPL(3) Line 491  c      Z=-MGFR(3)+MGPL(3)
491        X=0.        X=0.
492        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3))+TPGI(3)        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3))+TPGI(3)
493        Y=0.        Y=0.
494        CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('TPGI',N,'MGFR',X,Y,Z,0,'ONLY')
495          CALL GSPOS('TPGI',2,'MGF1',X,Y,Z,0,'ONLY')    
496          CALL GSPOS('TPGI',2,'MGF2',X,Y,Z,0,'ONLY')
497          CALL GSPOS('TPGI',2,'MGF3',X,Y,Z,0,'ONLY')
498          CALL GSPOS('TPGI',2,'MGF4',X,Y,Z,0,'ONLY')
499          CALL GSPOS('TPGI',2,'MGF5',X,Y,Z,0,'ONLY')
500  *  *
501  * Positioning volume TPGA     into the mother MGFR  * Positioning volume TPGA     into the mother MGFR
502  *  *
# Line 316  C ML: 10/11/05 Line 504  C ML: 10/11/05
504        N=1        N=1
505        X=0.        X=0.
506        Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)        Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)
507  C      Z=-MGFR(3)+2*MGPL(3)+TPGA(3)  
 C      Z=MGFR(3)-2*TPGU(3)-2*TRPB(3)-TPGA(3)  
508        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3))+TPGA(3)        Z=-MGFR(3)+2*(TPGD(3)+TPGI(3)+MGPL(3)+TPGI(3))+TPGA(3)
509  C      Y=-0.75  
510        CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')
511          CALL GSPOS('TPGA',1,'MGF1',X,Y,Z,0,'ONLY')
512          CALL GSPOS('TPGA',1,'MGF2',X,Y,Z,0,'ONLY')
513          CALL GSPOS('TPGA',1,'MGF3',X,Y,Z,0,'ONLY')
514          CALL GSPOS('TPGA',1,'MGF4',X,Y,Z,0,'ONLY')
515          CALL GSPOS('TPGA',1,'MGF5',X,Y,Z,0,'ONLY')
516  *  *
517  * Positioning volumes TRPB & MGPL into the mother SPEB  * Positioning volumes TRPB & MGPL into the mother SPEB
518  *  *
519        N= 1        N= 1
520        X= 0.        X= 0.
521  C      Y= SPEB(2) - TRPB(2)  
522        Y= -MGFR(2) + TRPB(2)        Y= -MGFR(2) + TRPB(2)
523  C      Z= SPEB(3) - TRPB(3)  
 c      Z=-MGFR(3)+2*MGPL(3)+2*TPGA(3)+TRPB(3)  
 C      Z=MGFR(3)-2*TPGU(3)-TRPB(3)  
524        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)
525  C      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  
526        CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')  cv      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
527  c      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY')        XX1=X+DXFLO(1)
528          YY1=Y+DYFLO(1)
529          CALL GSPOS('TRP1',1,'MGF1',XX1,YY1,Z,0,'ONLY')
530          XX2=X+DXFLO(2)
531          YY2=Y+DYFLO(2)
532          CALL GSPOS('TRP2',1,'MGF2',XX2,YY2,Z,0,'ONLY')
533          XX3=X+DXFLO(3)
534          YY3=Y+DYFLO(3)
535          CALL GSPOS('TRP3',1,'MGF3',XX3,YY3,Z,0,'ONLY')
536          XX4=X+DXFLO(4)
537          YY4=Y+DYFLO(4)
538          CALL GSPOS('TRP4',1,'MGF4',XX4,YY4,Z,0,'ONLY')
539          XX5=X+DXFLO(5)
540          YY5=Y+DYFLO(5)
541          CALL GSPOS('TRP5',1,'MGF5',XX5,YY5,Z,0,'ONLY')
542        Z=SPEB(3)        Z=SPEB(3)
543        DO I=1, 5        DO I=1, 5
544           X= 0.           X= 0.
545           Y= 0.           Y= 0.
546  C         Z= Z - TRPB(3) - MGFR(3)  
547           Z=Z-MGFR(3)           Z=Z-MGFR(3)
548           N= I           N= I
549           CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')  CV         CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
550  C         Y= SPEB(2) - TRPB(2)           IF(N.EQ.1) CALL GSPOS('MGF1',1,'SPEB',X,Y,Z,0,'ONLY')
551  C         Z= Z - MGFR(3) - TRPB(3)           IF(N.EQ.2) CALL GSPOS('MGF2',1,'SPEB',X,Y,Z,0,'ONLY')
552  C         N= (I+1)           IF(N.EQ.3) CALL GSPOS('MGF3',1,'SPEB',X,Y,Z,0,'ONLY')
553  C         IF(I.EQ.5) THEN           IF(N.EQ.4) CALL GSPOS('MGF4',1,'SPEB',X,Y,Z,0,'ONLY')
554  CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')           IF(N.EQ.5) CALL GSPOS('MGF5',1,'SPEB',X,Y,Z,0,'ONLY')
555  C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,10,'ONLY')  
 C         ELSE  
 CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  
 C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')  
 C         ENDIF  
556           Z=Z-MGFR(3)           Z=Z-MGFR(3)
557        ENDDO        ENDDO
558        N=2        N=2
# Line 360  C         ENDIF Line 560  C         ENDIF
560        Y=MGFR(2) - TRPB(2)        Y=MGFR(2) - TRPB(2)
561        Z=Z-TRPB(3)        Z=Z-TRPB(3)
562        ZTRPB=Z        ZTRPB=Z
563  c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')  
564        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')  cv      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
565          XX6=X+DXFLO(6)
566          YY6=Y+DYFLO(6)
567          CALL GSPOS('TRP6',1,'SPEB',XX6,YY6,Z,10,'ONLY')
568  *  *
569  * Positioning volume TBAL into the mother SPEB  * Positioning volume TBAL into the mother SPEB
570  *  *
# Line 382  C DTBAL is the distance along x of the v Line 585  C DTBAL is the distance along x of the v
585  * Positioning volume TPGU     into the mother MGFR  * Positioning volume TPGU     into the mother MGFR
586  *  *
587        X=0.        X=0.
588  C      Z=MGFR(3)-TPGU(3)  
589        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))+
590       +   TPGU(3)             +   TPGU(3)      
591        Y=-MGFR(2)+TRPB(2)        Y=-MGFR(2)+TRPB(2)
592        CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')  CV      CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
593          CALL GSPOS('TPGU',1,'MGF1',X,Y,Z,0,'ONLY')
594          CALL GSPOS('TPGU',1,'MGF2',X,Y,Z,0,'ONLY')
595          CALL GSPOS('TPGU',1,'MGF3',X,Y,Z,0,'ONLY')
596          CALL GSPOS('TPGU',1,'MGF4',X,Y,Z,0,'ONLY')
597          CALL GSPOS('TPGU',1,'MGF5',X,Y,Z,0,'ONLY')
598  *  *
599        RETURN        RETURN
600        END        END
Legend:
Removed from v.3.10  
changed lines
  Added in v.3.11
  ViewVC Help
Powered by ViewVC 1.1.23