gpamela/gpspe/gpspev.F


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