gpamela/gpspe/gpspev.F


*
* $Id: gpspev.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
*
* $Log: gpspev.F,v $
* 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 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)
      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=MAL
      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)
*
* Positioning volume TSPA     into the mother TRSL
*
      N= 1
      X= 0.
c ml: 26/11/66:
c      Y= TRSL(2) - TSPA(2)
      Y=+TRSL(2)/2
      Z= 0.
      CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
      N= 2
C      Y= Y - 2.*TSPA(2)
      Y=-TRSL(2)/2
      CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
c ml: 26/11/04:
*
c* Positioning volume TRCP     into the mother TRSL
CC ml 20/04/05: 
CC* Positioning volume TRCP into mother TRPL
* Positioning volume TRCP into mother TRPB
*
c      N= 1
CC ml 20/04/05: 
cc      Z= 0.
cc      Y= 0.
      Z=-TRPB(3)+0.11+TRCP(3)
      Y=-TRPB(2)+1.5-0.225+TRCP(2)
CC end ml.
      DO I=1,4
CC ml 20/04/05:         
CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
         X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
CC end ml
c         CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')
         CALL GSPOS('TRCP',I,'TRPB',X,Y,Z,0,'ONLY')
      ENDDO  
c      X= -TRSL(1) + TRCP(1)
c      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
C      N= 2
C      X= -X
C      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
*****
cc ml 20/04/05:
*
*** Positioning volume TRSL     into the mother TRPL
* Positioning volume TRSL     into the mother TRPB
*
c ml: 26/11/04:
C      Y= 0.
cc      Y=THBP(2)
      Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
C      Y=2.75
CC      Z= 0.
      Z=-TRPB(3)+0.11+TRCP(3)
      DO I=1, NTRSL
c         X= -TRPL(1) + TRSL(1) +2.*(I-1)*TRSL(1)
cc         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*TRSL(1)
         X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*TRSL(1)
cc         CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')
         CALL GSPOS('TRSL',I,'TRPB',X,Y,Z,0,'ONLY')
      ENDDO
**********
cc ml: 20/04/05
*
********* Positioning volume THBP     into the mother TRPL
* Positioning volume THBP     into the mother TRPB
*
cc      Y=-TRSL(2)
      Y=-TRPB(2)+1.5+THBP(2)
cc      Z= 0.
      Z=-TRPB(3)+0.11+TRCP(3)
      DO I=1, NTHBP
CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)
         X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)
cc         CALL GSPOS('THBP',I,'TRPL',X,Y,Z,0,'ONLY')
         CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')
      ENDDO
c end ml.
cc ml: 20/04/05
*
* Positioning volume TPAI into the mother TRPB
*
      X=-TRPB(1)+2.29+TPAI(1)
      Y=-TRPB(2)+1.5+2.5+TPAI(2)
      Z=-TRPB(3)+TPAI(3)
      CALL GSPOS('TPAI',I,'TRPB',X,Y,Z,0,'ONLY')
*
* Positioning volume TPAS into the mother TRPB
*
      X=-TRPB(1)+2.29+TPAS(1)
      Y=-TRPB(2)+1.5+TPAS(2)
      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 MGPI     into the mother MGPL
*
      N= 1
      X= 0.
      Y= 0.
      Z= 0.
      CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')
C Ml: 20/04/05:
c*
c* Positioning volume TRPL     into the mother TRPB
c*
c      N= 1
c      Z= 0.
c      CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')
c end ml.
*
* Positioning volumes TRPB & MGPL into the mother SPEB
*
      N= 1
      X= 0.
      Y= 0.
      Z= SPEB(3) - TRPB(3)
      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
      DO I=1, 5
         Z= Z - TRPB(3) - MGPL(3)
         N= I
         CALL GSPOS('MGPL',N,'SPEB',X,Y,Z,0,'ONLY')
         Z= Z - MGPL(3) - TRPB(3)
         N= (I+1)
         CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
      ENDDO
*
      RETURN
      END
1	cafagna	3.2
2	cafagna	3.1	*
3	cafagna	3.2	* $Id: gpspev.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
4			*
5			* $Log: gpspev.F,v $
6			* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna
7			* First GPAMELA release on CVS
8	cafagna	3.1	*
9			*
10			*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
11			*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna
12			*CMZ : 1.02/00 07/03/97 14.48.41 by Francesco Cafagna
13			*CMZ : 1.00/02 06/02/96 17.34.50 by Francesco Cafagna
14			*-- Author : Francesco Cafagna 09/12/95
15			SUBROUTINE GPSPEV
16			************************************************************************
17			* *
18			* Volumes definition for the silicon tracking system and magnet *
19			* Called by: GPGEM *
20			* Author: Francesco Cafagna, 05/12/95 13.39.41 *
21			* *
22			************************************************************************
23			#include "gpgeo.inc"
24			#include "gpmed.inc"
25			*
26			INTEGER IROT,IVOLU,N,NMED,NUM
27			REAL X,Y,Z
28			*
29			* Define the SPEB volume
30			*
31			NMED= MN2
32			CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
33			*
34			* Define the MGPL volume
35			*
36			NMED= MMAG
37			CALL GSVOLU('MGPL','BOX ',NMED,MGPL, 3,IVOLU)
38			*
39			* Define the MGPI volume
40			*
41			NMED= MN2
42			CALL GSVOLU('MGPI','BOX ',NMED,MGPI, 3,IVOLU)
43			*
44			* Define the TRPB volume
45			*
46			NMED= MAL
47			CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
48	cafagna	3.2	*****
49			cC ml 20/04/05:
50	cafagna	3.1	*
51	cafagna	3.2	cc* Define the TRPL volume
52			* Define the TPAS and TPAI volume
53	cafagna	3.1	*
54			NMED= MN2
55	cafagna	3.2	CC CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
56			CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
57			CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
58			CC end ml.
59			****************
60	cafagna	3.1	*
61			* Define the TRSL volume
62			*
63	cafagna	3.2	c ml: 26/11/04:
64			c NMED= MN2
65			NMED=MSIT
66	cafagna	3.1	CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
67			*
68	cafagna	3.2	* Define the THBP volume
69			*
70			NMED=MAL
71			CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
72			C end ml.
73			*
74	cafagna	3.1	* Define the TSPA volume
75			*
76			NMED= MSIT
77			CALL GSVOLU('TSPA','BOX ',NMED,TSPA, 3,IVOLU)
78			*
79			* Define the TRCP volume
80			*
81			NMED= MCF
82			CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
83			*
84			* Positioning volume TSPA into the mother TRSL
85			*
86			N= 1
87			X= 0.
88	cafagna	3.2	c ml: 26/11/66:
89			c Y= TRSL(2) - TSPA(2)
90			Y=+TRSL(2)/2
91	cafagna	3.1	Z= 0.
92			CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
93			N= 2
94	cafagna	3.2	C Y= Y - 2.*TSPA(2)
95			Y=-TRSL(2)/2
96	cafagna	3.1	CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
97	cafagna	3.2	c ml: 26/11/04:
98			*
99			c* Positioning volume TRCP into the mother TRSL
100			CC ml 20/04/05:
101			CC* Positioning volume TRCP into mother TRPL
102			* Positioning volume TRCP into mother TRPB
103			*
104			c N= 1
105			CC ml 20/04/05:
106			cc Z= 0.
107			cc Y= 0.
108			Z=-TRPB(3)+0.11+TRCP(3)
109			Y=-TRPB(2)+1.5-0.225+TRCP(2)
110			CC end ml.
111			DO I=1,4
112			CC ml 20/04/05:
113			CC X=-TRPL(1)+(2I-1)TRCP(1)+2(I-1)TRSL(1)
114			X=-TRPB(1)+2.29+(2I-1)TRCP(1)+2(I-1)TRSL(1)
115			CC end ml
116			c CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')
117			CALL GSPOS('TRCP',I,'TRPB',X,Y,Z,0,'ONLY')
118			ENDDO
119			c X= -TRSL(1) + TRCP(1)
120			c CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
121			C N= 2
122			C X= -X
123			C CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
124			*****
125			cc ml 20/04/05:
126			*
127			*** Positioning volume TRSL into the mother TRPL
128			* Positioning volume TRSL into the mother TRPB
129			*
130			c ml: 26/11/04:
131			C Y= 0.
132			cc Y=THBP(2)
133			Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
134			C Y=2.75
135			CC Z= 0.
136			Z=-TRPB(3)+0.11+TRCP(3)
137			DO I=1, NTRSL
138			c X= -TRPL(1) + TRSL(1) +2.(I-1)TRSL(1)
139			cc X= -TRPL(1) + 2ITRCP(1)+(2I-1)TRSL(1)
140			X=-TRPB(1)+2.29+2ITRCP(1)+(2I-1)TRSL(1)
141			cc CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')
142			CALL GSPOS('TRSL',I,'TRPB',X,Y,Z,0,'ONLY')
143			ENDDO
144			**********
145			cc ml: 20/04/05
146	cafagna	3.1	*
147	cafagna	3.2	********* Positioning volume THBP into the mother TRPL
148			* Positioning volume THBP into the mother TRPB
149	cafagna	3.1	*
150	cafagna	3.2	cc Y=-TRSL(2)
151			Y=-TRPB(2)+1.5+THBP(2)
152			cc Z= 0.
153			Z=-TRPB(3)+0.11+TRCP(3)
154			DO I=1, NTHBP
155			CC X= -TRPL(1) + 2ITRCP(1)+(2I-1)THBP(1)
156			X=-TRPB(1)+2.29+2ITRCP(1)+(2I-1)THBP(1)
157			cc CALL GSPOS('THBP',I,'TRPL',X,Y,Z,0,'ONLY')
158			CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')
159			ENDDO
160			c end ml.
161			cc ml: 20/04/05
162	cafagna	3.1	*
163	cafagna	3.2	* Positioning volume TPAI into the mother TRPB
164	cafagna	3.1	*
165	cafagna	3.2	X=-TRPB(1)+2.29+TPAI(1)
166			Y=-TRPB(2)+1.5+2.5+TPAI(2)
167			Z=-TRPB(3)+TPAI(3)
168			CALL GSPOS('TPAI',I,'TRPB',X,Y,Z,0,'ONLY')
169			*
170			* Positioning volume TPAS into the mother TRPB
171			*
172			X=-TRPB(1)+2.29+TPAS(1)
173			Y=-TRPB(2)+1.5+TPAS(2)
174			Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
175			CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
176			c end ml.
177
178	cafagna	3.1	*
179			* Positioning volume MGPI into the mother MGPL
180			*
181			N= 1
182			X= 0.
183			Y= 0.
184			Z= 0.
185			CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')
186	cafagna	3.2	C Ml: 20/04/05:
187			c*
188			c* Positioning volume TRPL into the mother TRPB
189			c*
190			c N= 1
191			c Z= 0.
192			c CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')
193			c end ml.
194	cafagna	3.1	*
195			* Positioning volumes TRPB & MGPL into the mother SPEB
196			*
197			N= 1
198			X= 0.
199			Y= 0.
200			Z= SPEB(3) - TRPB(3)
201			CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
202			DO I=1, 5
203			Z= Z - TRPB(3) - MGPL(3)
204			N= I
205			CALL GSPOS('MGPL',N,'SPEB',X,Y,Z,0,'ONLY')
206			Z= Z - MGPL(3) - TRPB(3)
207			N= (I+1)
208			CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
209			ENDDO
210			*
211			RETURN
212			END