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