gpamela/gpspe/gpspev.F


*
* $Id: gpspev.F,v 3.3 2005/07/25 11:53:21 cafagna Exp $
*
* $Log: gpspev.F,v $
* Revision 3.3  2005/07/25 11:53:21  cafagna
* Several updates. See history for details
*
* Revision 3.2  2005/06/21 02:42:27  cafagna
* Major modification to the geometry and to the random number chain
*
* Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
* First GPAMELA release on CVS
*
*
*CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
*CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna
*CMZ :  1.02/00 07/03/97  14.48.41  by  Francesco Cafagna
*CMZ :  1.00/02 06/02/96  17.34.50  by  Francesco Cafagna
*-- Author :    Francesco Cafagna   09/12/95
      SUBROUTINE GPSPEV
************************************************************************
*                                                                      *
* Volumes definition for the silicon tracking system and magnet        *
* Called by: GPGEM                                                     *
* Author: Francesco Cafagna, 05/12/95 13.39.41                         *
*                                                                      *
************************************************************************
#include "gpgeo.inc"
#include "gpmed.inc"
*
      INTEGER IROT,IVOLU,N,NMED,NUM
      REAL X,Y,Z
*
*  Define the SPEB     volume
*
      NMED= MN2
      CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
*
*  Define the MGFR     volume
*
      NMED= MAL
      CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
*
*  Define the MGPL     volume
*
      NMED= MMAG
      CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU)
*
*  Define the MGPI     volume
*
      NMED= MN2
      CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU)
*
*  Define the TRPB     volume
*
      NMED= MAL
      CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
*****
cC ml 20/04/05:
*
cc*  Define the TRPL     volume
*  Define the TPAS and TPAI volume
*
      NMED= MN2
CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
      CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
      NMED= MAL
      CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
CC end ml.
****************
*
*  Define the TRSL     volume
*
c ml: 26/11/04:
c      NMED= MN2
      NMED=MSIT
      CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
*
* Define the THBP volume
*
      NMED=MG10C
      CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
C end ml.
*
*  Define the TSPA     volume
*
      NMED= MSIT
      CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU)
*
*  Define the TRCP     volume
*
      NMED= MCF
      CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
*
*  Define the MGPA     volume
*
      NMED= MAL
      CALL GSVOLU('MGPA','BOX ',NMED,MGPA, 3,IVOLU)
*
*  Define the TPGA     volume
*
      NMED= MN2
      CALL GSVOLU('TPGA','BOX ',NMED,TPGA, 3,IVOLU)
*
*  Define the TPGU     volume
*
      NMED= MN2
      CALL GSVOLU('TPGU','BOX ',NMED,TPGU, 3,IVOLU)
*
*  Define the TPGD     volume
*
      NMED= MN2
      CALL GSVOLU('TPGD','BOX ',NMED,TPGD, 3,IVOLU)

*
* Positioning volume TSPA     into the mother TRSL
*
      N= 1
      X= 0.
      Y= 0.
      Z= 0.
      CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
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')
c ml: 26/11/04:
* Positioning volume TRCP into mother TPAS
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
C #       Z=-TRPB(3)+0.11+TRCP(3)
C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
CC end ml.
c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
      Z=TRPB(3)-2*TPAI(3)-TRCP(3)
      y=-0.225-TPAS(2)+TRCP(2)
C # end CAF mod
      DO I=1,4
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)
         X=-TPAS(1)+2*(I-1)*0.0015+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
CC end ml
         CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
C # end CAF mod
      ENDDO  
cc ml 20/04/05:
* Positioning volume TRSL     into the mother TPAS
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
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)
      Z=TRPB(3)-2*TPAI(3)-TRCP(3)
      DO I=1, NTRSL
c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.
c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
     +   (2*(MOD(I-1,3)+1)-1)*0.0015+(2*(MOD(I-1,3)+1)-1)*TRSL(1)         
         IF(I.LE.3) THEN
C ml: 22/9/05:
c            Y=TPAS(2)-TRSL(2)
            Y=TPAS(2)-0.1485-TRSL(2)
         ELSE
C            Y=TPAS(2)-3*TRSL(2)
            Y=TPAS(2)-0.1485-0.0015-3*TRSL(2)
         ENDIF
         CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
      ENDDO
C # end CAF mod
**********
cc ml: 20/04/05
* Positioning volume THBP     into the mother TPAS
C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
C #       Y=-TRPB(2)+1.5+THBP(2)
C # cc      Z= 0.
C #       Z=-TRPB(3)+0.11+TRCP(3)
c ml: 22/9/05:
c      Y=-TPAS(2)+THBP(2)
      Y=-TPAS(2)+0.1485+THBP(2)
C      Z=-0.005
c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
      Z=TRPB(3)-2*TPAI(3)-TRCP(3)
      DO I=1, NTHBP
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)
c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon
c       ladders 
c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)
         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)+
     +      (2*(MOD(I-1,3)+1)-1)*0.0015
         CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
      ENDDO
c end ml.
cc ml: 20/04/05
*
* Positioning volume TPAI into the mother TPAS
*
C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
      X=0
      Y=-TPAS(2)+TPAI(2)
c      Z=-TPAS(3)+TPAI(3)
      Z=TRPB(3)-TPAI(3)
C #       X=-TRPB(1)+2.29+TPAI(1)
C #       Y=-TRPB(2)+1.5+2.5+TPAI(2)
C #       Z=-TRPB(3)+TPAI(3)
      CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')
C # end CAF mod
*
* Positioning volume TPAS into the mother TRPB
*
C # Mod By Caf on 15th Jul 2005
      X=0
      Y=-TRPB(2)+1.5+TPAS(2)
      Z=0
C #       X=-TRPB(1)+2.29+TPAS(1)
C #       Y=-TRPB(2)+1.5+TPAS(2)
C #       Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
      CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
c end ml.
*
* Positioning volume MGPA     into the mother MGPL
*
      N= 1
      X= 0.
      Y= 0.
      Z= 0.
      CALL GSPOS('MGPA',N,'MGPL',X,Y,Z,0,'ONLY')
*
* Positioning volume TPGD     into the mother MGFR
*
      X=0.
      Z=-MGFR(3)+TPGD(3)
      Y=0.
      CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
*
* Positioning volume MGPI     into the mother MGPL
*
      N= 1
      X= 0.
      Y= 0.
CC ML 10/11/05:  Positioning volume MGPI     into the mother MGFR
      Z= 0.
      CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')
C      Z= -MGFR(3)+2*TPGD(3)+MGPI(3)
C      CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY')
*
* Positioning volume MGPL     into the mother MGFR
*
c ml: 10/11/05
c      Z=-MGFR(3)+MGPL(3)
      Z=-MGFR(3)+2*TPGD(3)+MGPI(3)
      CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
*
* Positioning volume TPGA     into the mother MGFR
*
C ML: 10/11/05
      X=0.
      Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)
C      Z=-MGFR(3)+2*MGPL(3)+TPGA(3)
      Z=MGFR(3)-2*TPGU(3)-2*TRPB(3)-TPGA(3)
C      Y=-0.75
      CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')
*
* Positioning volumes TRPB & MGPL into the mother SPEB
*
      N= 1
      X= 0.
C      Y= SPEB(2) - TRPB(2)
      Y= -MGFR(2) + TRPB(2)
C      Z= SPEB(3) - TRPB(3)
c      Z=-MGFR(3)+2*MGPL(3)+2*TPGA(3)+TRPB(3)
      Z=MGFR(3)-2*TPGU(3)-TRPB(3)
C      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
c      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY')
      Z=SPEB(3)
      DO I=1, 5
         X= 0.
         Y= 0.
C         Z= Z - TRPB(3) - MGFR(3)
         Z=Z-MGFR(3)
         N= I
         CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
C         Y= SPEB(2) - TRPB(2)
C         Z= Z - MGFR(3) - TRPB(3)
C         N= (I+1)
C         IF(I.EQ.5) THEN
CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
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
         Z=Z-MGFR(3)
      ENDDO
      N=2
      X=0.
      Y=MGFR(2) - TRPB(2)
      Z=Z-TRPB(3)
c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')
      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
*
* Positioning volume TPGU     into the mother MGFR
*
      X=0.
      Z=MGFR(3)-TPGU(3)
      Y=-MGFR(2)+TRPB(2)
      CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
*
      RETURN
      END
1	cafagna	3.2
2	cafagna	3.1	*
3	pam-ba	3.4	* $Id: gpspev.F,v 3.3 2005/07/25 11:53:21 cafagna Exp $
4	cafagna	3.2	*
5			* $Log: gpspev.F,v $
6	pam-ba	3.4	* Revision 3.3 2005/07/25 11:53:21 cafagna
7			* Several updates. See history for details
8			*
9	cafagna	3.3	* Revision 3.2 2005/06/21 02:42:27 cafagna
10			* Major modification to the geometry and to the random number chain
11			*
12	cafagna	3.2	* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna
13			* First GPAMELA release on CVS
14	cafagna	3.1	*
15			*
16			*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
17			*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna
18			*CMZ : 1.02/00 07/03/97 14.48.41 by Francesco Cafagna
19			*CMZ : 1.00/02 06/02/96 17.34.50 by Francesco Cafagna
20			*-- Author : Francesco Cafagna 09/12/95
21			SUBROUTINE GPSPEV
22			************************************************************************
23			* *
24			* Volumes definition for the silicon tracking system and magnet *
25			* Called by: GPGEM *
26			* Author: Francesco Cafagna, 05/12/95 13.39.41 *
27			* *
28			************************************************************************
29			#include "gpgeo.inc"
30			#include "gpmed.inc"
31			*
32			INTEGER IROT,IVOLU,N,NMED,NUM
33			REAL X,Y,Z
34			*
35			* Define the SPEB volume
36			*
37			NMED= MN2
38			CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
39			*
40	cafagna	3.3	* Define the MGFR volume
41			*
42			NMED= MAL
43			CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
44			*
45	cafagna	3.1	* Define the MGPL volume
46			*
47			NMED= MMAG
48			CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU)
49			*
50			* Define the MGPI volume
51			*
52			NMED= MN2
53			CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU)
54			*
55			* Define the TRPB volume
56			*
57			NMED= MAL
58			CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
59	cafagna	3.2	*****
60			cC ml 20/04/05:
61	cafagna	3.1	*
62	cafagna	3.2	cc* Define the TRPL volume
63			* Define the TPAS and TPAI volume
64	cafagna	3.1	*
65			NMED= MN2
66	cafagna	3.2	CC CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
67			CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
68	cafagna	3.3	NMED= MAL
69	cafagna	3.2	CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
70			CC end ml.
71			****************
72	cafagna	3.1	*
73			* Define the TRSL volume
74			*
75	cafagna	3.2	c ml: 26/11/04:
76			c NMED= MN2
77			NMED=MSIT
78	cafagna	3.1	CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
79			*
80	cafagna	3.2	* Define the THBP volume
81			*
82	cafagna	3.3	NMED=MG10C
83	cafagna	3.2	CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
84			C end ml.
85			*
86	cafagna	3.1	* Define the TSPA volume
87			*
88			NMED= MSIT
89			CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU)
90			*
91			* Define the TRCP volume
92			*
93			NMED= MCF
94			CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
95			*
96	pam-ba	3.4	* Define the MGPA volume
97			*
98			NMED= MAL
99			CALL GSVOLU('MGPA','BOX ',NMED,MGPA, 3,IVOLU)
100			*
101			* Define the TPGA volume
102			*
103			NMED= MN2
104			CALL GSVOLU('TPGA','BOX ',NMED,TPGA, 3,IVOLU)
105			*
106			* Define the TPGU volume
107			*
108			NMED= MN2
109			CALL GSVOLU('TPGU','BOX ',NMED,TPGU, 3,IVOLU)
110			*
111			* Define the TPGD volume
112			*
113			NMED= MN2
114			CALL GSVOLU('TPGD','BOX ',NMED,TPGD, 3,IVOLU)
115
116			*
117	cafagna	3.1	* Positioning volume TSPA into the mother TRSL
118			*
119			N= 1
120			X= 0.
121	cafagna	3.3	Y= 0.
122	cafagna	3.1	Z= 0.
123			CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
124	cafagna	3.3	C # N= 2
125			C # C Y= Y - 2.*TSPA(2)
126			C # Y=-TRSL(2)/2
127			C # CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
128	cafagna	3.2	c ml: 26/11/04:
129	pam-ba	3.4	* Positioning volume TRCP into mother TPAS
130	cafagna	3.3	C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
131			C # Z=-TRPB(3)+0.11+TRCP(3)
132			C # Y=-TRPB(2)+1.5-0.225+TRCP(2)
133	cafagna	3.2	CC end ml.
134	pam-ba	3.4	c Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
135			Z=TRPB(3)-2*TPAI(3)-TRCP(3)
136	cafagna	3.3	y=-0.225-TPAS(2)+TRCP(2)
137			C # end CAF mod
138	cafagna	3.2	DO I=1,4
139			CC ml 20/04/05:
140			CC X=-TRPL(1)+(2I-1)TRCP(1)+2(I-1)TRSL(1)
141	cafagna	3.3	C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
142			C # X=-TRPB(1)+2.29+(2I-1)TRCP(1)+2(I-1)TRSL(1)
143	pam-ba	3.4	c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.
144			c for now the glue is a space of nitrogen of 0.015 mm
145			c X=-TPAS(1)+(2I-1)TRCP(1)+2(I-1)TRSL(1)
146			X=-TPAS(1)+2(I-1)0.0015+(2I-1)TRCP(1)+2(I-1)TRSL(1)
147	cafagna	3.2	CC end ml
148	cafagna	3.3	CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
149			C # end CAF mod
150	cafagna	3.2	ENDDO
151			cc ml 20/04/05:
152	pam-ba	3.4	* Positioning volume TRSL into the mother TPAS
153	cafagna	3.3	C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
154			C # Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
155	cafagna	3.2	C Y=2.75
156			CC Z= 0.
157	pam-ba	3.4	C Z=-0.005
158			c Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
159			Z=TRPB(3)-2*TPAI(3)-TRCP(3)
160	cafagna	3.2	DO I=1, NTRSL
161	pam-ba	3.4	c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.
162			c X=-TPAS(1)+2(MOD(I-1,3)+1)TRCP(1)+
163			c + (2(MOD(I-1,3)+1)-1)TRSL(1)
164	cafagna	3.3	X=-TPAS(1)+2(MOD(I-1,3)+1)TRCP(1)+
165	pam-ba	3.4	+ (2(MOD(I-1,3)+1)-1)0.0015+(2(MOD(I-1,3)+1)-1)TRSL(1)
166	cafagna	3.3	IF(I.LE.3) THEN
167	pam-ba	3.4	C ml: 22/9/05:
168			c Y=TPAS(2)-TRSL(2)
169			Y=TPAS(2)-0.1485-TRSL(2)
170	cafagna	3.3	ELSE
171	pam-ba	3.4	C Y=TPAS(2)-3*TRSL(2)
172			Y=TPAS(2)-0.1485-0.0015-3*TRSL(2)
173	cafagna	3.3	ENDIF
174			CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
175	cafagna	3.2	ENDDO
176	cafagna	3.3	C # end CAF mod
177	cafagna	3.2	**********
178			cc ml: 20/04/05
179	pam-ba	3.4	* Positioning volume THBP into the mother TPAS
180	cafagna	3.3	C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
181			C # Y=-TRPB(2)+1.5+THBP(2)
182			C # cc Z= 0.
183			C # Z=-TRPB(3)+0.11+TRCP(3)
184	pam-ba	3.4	c ml: 22/9/05:
185			c Y=-TPAS(2)+THBP(2)
186			Y=-TPAS(2)+0.1485+THBP(2)
187			C Z=-0.005
188			c Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
189			Z=TRPB(3)-2*TPAI(3)-TRCP(3)
190	cafagna	3.2	DO I=1, NTHBP
191			CC X= -TRPL(1) + 2ITRCP(1)+(2I-1)THBP(1)
192	cafagna	3.3	C # X=-TRPB(1)+2.29+2ITRCP(1)+(2I-1)THBP(1)
193	pam-ba	3.4	c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon
194			c ladders
195			c X=-TPAS(1)+2ITRCP(1)+(2I-1)THBP(1)
196			X=-TPAS(1)+2ITRCP(1)+(2I-1)THBP(1)+
197			+ (2(MOD(I-1,3)+1)-1)0.0015
198	cafagna	3.3	CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
199	cafagna	3.2	ENDDO
200			c end ml.
201			cc ml: 20/04/05
202	cafagna	3.1	*
203	pam-ba	3.4	* Positioning volume TPAI into the mother TPAS
204	cafagna	3.1	*
205	cafagna	3.3	C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
206			X=0
207			Y=-TPAS(2)+TPAI(2)
208	pam-ba	3.4	c Z=-TPAS(3)+TPAI(3)
209			Z=TRPB(3)-TPAI(3)
210	cafagna	3.3	C # X=-TRPB(1)+2.29+TPAI(1)
211			C # Y=-TRPB(2)+1.5+2.5+TPAI(2)
212			C # Z=-TRPB(3)+TPAI(3)
213			CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')
214			C # end CAF mod
215	cafagna	3.2	*
216			* Positioning volume TPAS into the mother TRPB
217			*
218	cafagna	3.3	C # Mod By Caf on 15th Jul 2005
219			X=0
220	cafagna	3.2	Y=-TRPB(2)+1.5+TPAS(2)
221	cafagna	3.3	Z=0
222			C # X=-TRPB(1)+2.29+TPAS(1)
223			C # Y=-TRPB(2)+1.5+TPAS(2)
224			C # Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
225	cafagna	3.2	CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
226			c end ml.
227	pam-ba	3.4	*
228			* Positioning volume MGPA into the mother MGPL
229			*
230			N= 1
231			X= 0.
232			Y= 0.
233			Z= 0.
234			CALL GSPOS('MGPA',N,'MGPL',X,Y,Z,0,'ONLY')
235			*
236			* Positioning volume TPGD into the mother MGFR
237			*
238			X=0.
239			Z=-MGFR(3)+TPGD(3)
240			Y=0.
241			CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
242	cafagna	3.1	*
243			* Positioning volume MGPI into the mother MGPL
244			*
245			N= 1
246			X= 0.
247			Y= 0.
248	pam-ba	3.4	CC ML 10/11/05: Positioning volume MGPI into the mother MGFR
249	cafagna	3.1	Z= 0.
250	pam-ba	3.4	CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')
251			C Z= -MGFR(3)+2*TPGD(3)+MGPI(3)
252			C CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY')
253	cafagna	3.3	*
254			* Positioning volume MGPL into the mother MGFR
255			*
256	pam-ba	3.4	c ml: 10/11/05
257			c Z=-MGFR(3)+MGPL(3)
258			Z=-MGFR(3)+2*TPGD(3)+MGPI(3)
259	cafagna	3.3	CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
260	pam-ba	3.4	*
261			* Positioning volume TPGA into the mother MGFR
262			*
263			C ML: 10/11/05
264			X=0.
265			Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)
266			C Z=-MGFR(3)+2*MGPL(3)+TPGA(3)
267			Z=MGFR(3)-2TPGU(3)-2TRPB(3)-TPGA(3)
268			C Y=-0.75
269			CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')
270	cafagna	3.1	*
271			* Positioning volumes TRPB & MGPL into the mother SPEB
272			*
273			N= 1
274			X= 0.
275	pam-ba	3.4	C Y= SPEB(2) - TRPB(2)
276			Y= -MGFR(2) + TRPB(2)
277			C Z= SPEB(3) - TRPB(3)
278			c Z=-MGFR(3)+2MGPL(3)+2TPGA(3)+TRPB(3)
279			Z=MGFR(3)-2*TPGU(3)-TRPB(3)
280			C CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
281			CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
282			c CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY')
283			Z=SPEB(3)
284	cafagna	3.1	DO I=1, 5
285	cafagna	3.3	X= 0.
286			Y= 0.
287	pam-ba	3.4	C Z= Z - TRPB(3) - MGFR(3)
288			Z=Z-MGFR(3)
289	cafagna	3.1	N= I
290	cafagna	3.3	CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
291	pam-ba	3.4	C Y= SPEB(2) - TRPB(2)
292			C Z= Z - MGFR(3) - TRPB(3)
293			C N= (I+1)
294			C IF(I.EQ.5) THEN
295			CC CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
296			C CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,10,'ONLY')
297			C ELSE
298			CC CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
299			C CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
300			C ENDIF
301			Z=Z-MGFR(3)
302	cafagna	3.1	ENDDO
303	pam-ba	3.4	N=2
304			X=0.
305			Y=MGFR(2) - TRPB(2)
306			Z=Z-TRPB(3)
307			c CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')
308			CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
309			*
310			* Positioning volume TPGU into the mother MGFR
311			*
312			X=0.
313			Z=MGFR(3)-TPGU(3)
314			Y=-MGFR(2)+TRPB(2)
315			CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
316	cafagna	3.1	*
317			RETURN
318			END