gpamela/gpcas/gpcasv.F

*
* $Id: gpcasv.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
*
* $Log: gpcasv.F,v $
* Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
* First GPAMELA release on CVS
*
*
*               17/10/2002  11.58.25  by  Jens Lund
*CMZ :  3.00/00 13/11/2000  09.08.20  by  Emiliano Mocchiutti
*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  21.47.47  by  Francesco Cafagna
*CMZ :  1.00/02 08/02/96  17.18.05  by  Francesco Cafagna
*-- Author :    Francesco Cafagna   08/02/96
      SUBROUTINE GPCASV
************************************************************************
*                                                                      *
* Volumes definition for the side anticoincidence counters             *
* Called by: GPGEO                                                     *
*                                                                      *
* Author: Francesco Cafagna, 08/02/96 16.58.49                         *
*                                                                      *
************************************************************************
#include "gpgeo.inc"
#include "gpmed.inc"
*
      INTEGER IROT,IVOLU,N,NMED,NUM
      REAL X,Y,Z
*JeL:
      REAL SBOX(3),STRAP(11)
*
C      NMED=MAL
      NMED=MN2
*
*  Define the CAS 'virtual' volume
*
      CALL GSVOLU('CASA','BOX ',NMED,CASA, 3,IVOLU)
C      CALL GSVOLU('CSSY','BOX ',NMED,CSSY, 3,IVOLU)
*
* Side AC scintillator (CAS), 12 pieces of scint., N2 and Aluminum
*                            to define the right shape
*
* Aluminum (cover) box
      NMED=MAL
      SBOX(1) = 0.5*42.9
      SBOX(2) = 0.5*1.2
      SBOX(3) = 0.5*41.1
      CALL GSVOLU('STRP','BOX ',NMED,SBOX,3,IVOLU)
* Scintillator sheet
      NMED=MSCIN
      SBOX(1) = 0.5*40.6
      SBOX(2) = 0.5*0.8
      SBOX(3) = 0.5*38.8
      CALL GSVOLU('SID1','BOX ',NMED,SBOX,3,IVOLU)
* Plastic sheet to create rims
      NMED=MPLAS
      SBOX(1) = 0.5*42.2
      SBOX(2) = 0.5*0.8
      SBOX(3) = 0.5*40.4
      CALL GSVOLU('PLA1','BOX ',NMED,SBOX,3,IVOLU)
* Plastic trapezoid, TRAP 1,2 & 3
      NMED=MPLAS
      STRAP(1) = 0.5*32.2
      STRAP(2) = 7.361412634
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*8.32
      STRAP(10) = 0.5*8.32
      STRAP(11) = 0.
      CALL GSVOLU('STP1','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*4.8
      STRAP(2) = 20.55604522
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*3.6
      STRAP(10) = 0.5*3.6
      STRAP(11) = 0.
      CALL GSVOLU('STP2','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*1.8
      STRAP(2) = 58.17255342
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.5*3.6
      STRAP(6) = 0.5*3.6
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*9.4
      STRAP(10) = 0.5*9.4
      STRAP(11) = 0.
      CALL GSVOLU('STP3','TRAP',NMED,STRAP,11,IVOLU)
* Aluminum trapezoid, TRAP 1,2 & 3
      NMED=MAL
      STRAP(1) = 0.5*33.0
      STRAP(2) = 7.355684758
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*8.52
      STRAP(10) = 0.5*8.52
      STRAP(11) = 0.
      CALL GSVOLU('STA1','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*5.6
      STRAP(2) = 21.44773633
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*4.4
      STRAP(10) = 0.5*4.4
      STRAP(11) = 0.
      CALL GSVOLU('STA2','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*1.8
      STRAP(2) = 58.17255342
      STRAP(3) = 0.
      STRAP(4) = 0.5*0.8
      STRAP(5) = 0.5*4.4
      STRAP(6) = 0.5*4.4
      STRAP(7) = 0.
      STRAP(8) = 0.5*0.8
      STRAP(9) = 0.5*10.2
      STRAP(10) = 0.5*10.2
      STRAP(11) = 0.
      CALL GSVOLU('STA3','TRAP',NMED,STRAP,11,IVOLU)
* 'air' trapezoid, TRAP 1,2 & 3
      NMED=MN2
      STRAP(1) = 0.5*33.35
      STRAP(2) = 7.355377512
      STRAP(3) = 0.
      STRAP(4) = 0.5*1.2
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*1.2
      STRAP(9) = 0.5*8.61
      STRAP(10) = 0.5*8.61
      STRAP(11) = 0.
      CALL GSVOLU('STL1','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*5.6
      STRAP(2) = 21.44773633
      STRAP(3) = 0.
      STRAP(4) = 0.5*1.2
      STRAP(5) = 0.
      STRAP(6) = 0.
      STRAP(7) = 0.
      STRAP(8) = 0.5*1.2
      STRAP(9) = 0.5*4.4
      STRAP(10) = 0.5*4.4
      STRAP(11) = 0.
      CALL GSVOLU('STL2','TRAP',NMED,STRAP,11,IVOLU)
      STRAP(1) = 0.5*2.15
      STRAP(2) = 58.18078646
      STRAP(3) = 0.
      STRAP(4) = 0.5*1.2
      STRAP(5) = 0.5*4.3
      STRAP(6) = 0.5*4.3
      STRAP(7) = 0.
      STRAP(8) = 0.5*1.2
      STRAP(9) = 0.5*11.23
      STRAP(10) = 0.5*11.23
      STRAP(11) = 0.
      CALL GSVOLU('STL3','TRAP',NMED,STRAP,11,IVOLU)
*EM:
* Define the CAS scintillator volume
*
C      NMED=MSCIN
C      CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)
C      CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)
*
* Positioning Scintillators in cradles
*
      N = 1
      X = 0.
      Y = 0.
      Z = 0.
C      CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')
C      CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')
*END: EM.
*
*
* Position Al. box (STRP) in the CAS virt. box (CASA)
      CALL GSPOS('STRP',N,'CASA',X,Y,Z,0,'ONLY')
*
* Position the plastic sheet (PLA1) in STRP
      CALL GSPOS('PLA1',N,'STRP',X,Y,Z,0,'ONLY')
*
* Position the scintillator (SID1) in PLA1
      CALL GSPOS('SID1',N,'PLA1',X,Y,Z,0,'ONLY')

*
* Position Pl. TRAP (STP1) in SID1 and rotate
      X = 18.22
      Y = 0.
      Z = -3.3
      CALL GSPOS('STP1',N,'SID1',X,Y,Z,5,'ONLY')
*
* Position Pl. TRAP (STP2) in SID1 and rotate
      X = 19.4
      Y = 0.
      Z = 15.2
      CALL GSPOS('STP2',N,'SID1',X,Y,Z,4,'ONLY')
*
* Position Pl. TRAP (STP3) in SID1 and rotate
      X = 17.05
      Y = 0.
      Z = 18.5
      CALL GSPOS('STP3',N,'SID1',X,Y,Z,4,'ONLY')

*
* Position Al. TRAP (STA1) in PLA1 and rotate
      X = 18.97
      Y = 0.
      Z = -3.7
      CALL GSPOS('STA1',N,'PLA1',X,Y,Z,5,'ONLY')
*
* Position Al. TRAP (STA2) in PLA1 and rotate
      X = 20.
      Y = 0.
      Z = 15.6
      CALL GSPOS('STA2',N,'PLA1',X,Y,Z,4,'ONLY')
*
* Position Al. TRAP (STA3) in PLA1 and rotate
      X = 17.45
      Y = 0.
      Z = 19.3
      CALL GSPOS('STA3',N,'PLA1',X,Y,Z,4,'ONLY')
*
* Position 'air' TRAP (STL1) in STRP and rotate
      X = 19.2975
      Y = 0.
      Z = -3.875
      CALL GSPOS('STL1',N,'STRP',X,Y,Z,5,'ONLY')
*
* Position 'air' TRAP (STL2) in STRP and rotate
      X = 20.35
      Y = 0.
      Z = 15.6
      CALL GSPOS('STL2',N,'STRP',X,Y,Z,4,'ONLY')
*
* Position 'air' TRAP (STL3) in STRP and rotate
      X = 17.5675
      Y = 0.
      Z = 19.475
      CALL GSPOS('STL3',N,'STRP',X,Y,Z,4,'ONLY')
*END: JeL.
C*
C      NMED=MAL
C*
C*  Define the CAS cradle volume
C*
C      CALL GSVOLU('CSSX','BOX ',NMED,CSSX, 3,IVOLU)
C      CALL GSVOLU('CSSY','BOX ',NMED,CSSY, 3,IVOLU)
C*EM:
C* Define the CAS scintillator volume
C*
C      NMED=MSCIN
C      CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)
C      CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)
C*
C* Positioning Scintillators in cradles
C*
C      N = 1
C      X = 0.
C      Y = 0.
C      Z = 0.
C      CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')
C      CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')
C*END: EM.
C*
      RETURN
      END
1	cafagna	3.1	*
2	pamela	3.2	* $Id: gpcasv.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
3	cafagna	3.1	*
4	pamela	3.2	* $Log: gpcasv.F,v $
5			* Revision 3.1.1.1 2002/07/11 16:02:14 cafagna
6			* First GPAMELA release on CVS
7	cafagna	3.1	*
8	pamela	3.2	*
9			* 17/10/2002 11.58.25 by Jens Lund
10	cafagna	3.1	*CMZ : 3.00/00 13/11/2000 09.08.20 by Emiliano Mocchiutti
11			*CMZ : 2.01/00 05/04/2000 14.35.18 by Marialuigia Ambriola
12			*CMZ : 2.00/00 03/03/2000 15.39.06 by Francesco Cafagna
13			*CMZ : 1.02/00 07/03/97 21.47.47 by Francesco Cafagna
14			*CMZ : 1.00/02 08/02/96 17.18.05 by Francesco Cafagna
15			*-- Author : Francesco Cafagna 08/02/96
16			SUBROUTINE GPCASV
17			************************************************************************
18			* *
19			* Volumes definition for the side anticoincidence counters *
20			* Called by: GPGEO *
21			* *
22			* Author: Francesco Cafagna, 08/02/96 16.58.49 *
23			* *
24			************************************************************************
25			#include "gpgeo.inc"
26			#include "gpmed.inc"
27			*
28			INTEGER IROT,IVOLU,N,NMED,NUM
29			REAL X,Y,Z
30	pamela	3.2	*JeL:
31			REAL SBOX(3),STRAP(11)
32			*
33			C NMED=MAL
34			NMED=MN2
35	cafagna	3.1	*
36	pamela	3.2	* Define the CAS 'virtual' volume
37			*
38			CALL GSVOLU('CASA','BOX ',NMED,CASA, 3,IVOLU)
39			C CALL GSVOLU('CSSY','BOX ',NMED,CSSY, 3,IVOLU)
40	cafagna	3.1	*
41	pamela	3.2	* Side AC scintillator (CAS), 12 pieces of scint., N2 and Aluminum
42			* to define the right shape
43	cafagna	3.1	*
44	pamela	3.2	* Aluminum (cover) box
45			NMED=MAL
46			SBOX(1) = 0.5*42.9
47			SBOX(2) = 0.5*1.2
48			SBOX(3) = 0.5*41.1
49			CALL GSVOLU('STRP','BOX ',NMED,SBOX,3,IVOLU)
50			* Scintillator sheet
51			NMED=MSCIN
52			SBOX(1) = 0.5*40.6
53			SBOX(2) = 0.5*0.8
54			SBOX(3) = 0.5*38.8
55			CALL GSVOLU('SID1','BOX ',NMED,SBOX,3,IVOLU)
56			* Plastic sheet to create rims
57			NMED=MPLAS
58			SBOX(1) = 0.5*42.2
59			SBOX(2) = 0.5*0.8
60			SBOX(3) = 0.5*40.4
61			CALL GSVOLU('PLA1','BOX ',NMED,SBOX,3,IVOLU)
62			* Plastic trapezoid, TRAP 1,2 & 3
63			NMED=MPLAS
64			STRAP(1) = 0.5*32.2
65			STRAP(2) = 7.361412634
66			STRAP(3) = 0.
67			STRAP(4) = 0.5*0.8
68			STRAP(5) = 0.
69			STRAP(6) = 0.
70			STRAP(7) = 0.
71			STRAP(8) = 0.5*0.8
72			STRAP(9) = 0.5*8.32
73			STRAP(10) = 0.5*8.32
74			STRAP(11) = 0.
75			CALL GSVOLU('STP1','TRAP',NMED,STRAP,11,IVOLU)
76			STRAP(1) = 0.5*4.8
77			STRAP(2) = 20.55604522
78			STRAP(3) = 0.
79			STRAP(4) = 0.5*0.8
80			STRAP(5) = 0.
81			STRAP(6) = 0.
82			STRAP(7) = 0.
83			STRAP(8) = 0.5*0.8
84			STRAP(9) = 0.5*3.6
85			STRAP(10) = 0.5*3.6
86			STRAP(11) = 0.
87			CALL GSVOLU('STP2','TRAP',NMED,STRAP,11,IVOLU)
88			STRAP(1) = 0.5*1.8
89			STRAP(2) = 58.17255342
90			STRAP(3) = 0.
91			STRAP(4) = 0.5*0.8
92			STRAP(5) = 0.5*3.6
93			STRAP(6) = 0.5*3.6
94			STRAP(7) = 0.
95			STRAP(8) = 0.5*0.8
96			STRAP(9) = 0.5*9.4
97			STRAP(10) = 0.5*9.4
98			STRAP(11) = 0.
99			CALL GSVOLU('STP3','TRAP',NMED,STRAP,11,IVOLU)
100			* Aluminum trapezoid, TRAP 1,2 & 3
101			NMED=MAL
102			STRAP(1) = 0.5*33.0
103			STRAP(2) = 7.355684758
104			STRAP(3) = 0.
105			STRAP(4) = 0.5*0.8
106			STRAP(5) = 0.
107			STRAP(6) = 0.
108			STRAP(7) = 0.
109			STRAP(8) = 0.5*0.8
110			STRAP(9) = 0.5*8.52
111			STRAP(10) = 0.5*8.52
112			STRAP(11) = 0.
113			CALL GSVOLU('STA1','TRAP',NMED,STRAP,11,IVOLU)
114			STRAP(1) = 0.5*5.6
115			STRAP(2) = 21.44773633
116			STRAP(3) = 0.
117			STRAP(4) = 0.5*0.8
118			STRAP(5) = 0.
119			STRAP(6) = 0.
120			STRAP(7) = 0.
121			STRAP(8) = 0.5*0.8
122			STRAP(9) = 0.5*4.4
123			STRAP(10) = 0.5*4.4
124			STRAP(11) = 0.
125			CALL GSVOLU('STA2','TRAP',NMED,STRAP,11,IVOLU)
126			STRAP(1) = 0.5*1.8
127			STRAP(2) = 58.17255342
128			STRAP(3) = 0.
129			STRAP(4) = 0.5*0.8
130			STRAP(5) = 0.5*4.4
131			STRAP(6) = 0.5*4.4
132			STRAP(7) = 0.
133			STRAP(8) = 0.5*0.8
134			STRAP(9) = 0.5*10.2
135			STRAP(10) = 0.5*10.2
136			STRAP(11) = 0.
137			CALL GSVOLU('STA3','TRAP',NMED,STRAP,11,IVOLU)
138			* 'air' trapezoid, TRAP 1,2 & 3
139			NMED=MN2
140			STRAP(1) = 0.5*33.35
141			STRAP(2) = 7.355377512
142			STRAP(3) = 0.
143			STRAP(4) = 0.5*1.2
144			STRAP(5) = 0.
145			STRAP(6) = 0.
146			STRAP(7) = 0.
147			STRAP(8) = 0.5*1.2
148			STRAP(9) = 0.5*8.61
149			STRAP(10) = 0.5*8.61
150			STRAP(11) = 0.
151			CALL GSVOLU('STL1','TRAP',NMED,STRAP,11,IVOLU)
152			STRAP(1) = 0.5*5.6
153			STRAP(2) = 21.44773633
154			STRAP(3) = 0.
155			STRAP(4) = 0.5*1.2
156			STRAP(5) = 0.
157			STRAP(6) = 0.
158			STRAP(7) = 0.
159			STRAP(8) = 0.5*1.2
160			STRAP(9) = 0.5*4.4
161			STRAP(10) = 0.5*4.4
162			STRAP(11) = 0.
163			CALL GSVOLU('STL2','TRAP',NMED,STRAP,11,IVOLU)
164			STRAP(1) = 0.5*2.15
165			STRAP(2) = 58.18078646
166			STRAP(3) = 0.
167			STRAP(4) = 0.5*1.2
168			STRAP(5) = 0.5*4.3
169			STRAP(6) = 0.5*4.3
170			STRAP(7) = 0.
171			STRAP(8) = 0.5*1.2
172			STRAP(9) = 0.5*11.23
173			STRAP(10) = 0.5*11.23
174			STRAP(11) = 0.
175			CALL GSVOLU('STL3','TRAP',NMED,STRAP,11,IVOLU)
176	cafagna	3.1	*EM:
177			* Define the CAS scintillator volume
178			*
179	pamela	3.2	C NMED=MSCIN
180			C CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)
181			C CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)
182	cafagna	3.1	*
183			* Positioning Scintillators in cradles
184			*
185			N = 1
186			X = 0.
187			Y = 0.
188			Z = 0.
189	pamela	3.2	C CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')
190			C CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')
191	cafagna	3.1	*END: EM.
192			*
193	pamela	3.2	*
194			* Position Al. box (STRP) in the CAS virt. box (CASA)
195			CALL GSPOS('STRP',N,'CASA',X,Y,Z,0,'ONLY')
196			*
197			* Position the plastic sheet (PLA1) in STRP
198			CALL GSPOS('PLA1',N,'STRP',X,Y,Z,0,'ONLY')
199			*
200			* Position the scintillator (SID1) in PLA1
201			CALL GSPOS('SID1',N,'PLA1',X,Y,Z,0,'ONLY')
202
203			*
204			* Position Pl. TRAP (STP1) in SID1 and rotate
205			X = 18.22
206			Y = 0.
207			Z = -3.3
208			CALL GSPOS('STP1',N,'SID1',X,Y,Z,5,'ONLY')
209			*
210			* Position Pl. TRAP (STP2) in SID1 and rotate
211			X = 19.4
212			Y = 0.
213			Z = 15.2
214			CALL GSPOS('STP2',N,'SID1',X,Y,Z,4,'ONLY')
215			*
216			* Position Pl. TRAP (STP3) in SID1 and rotate
217			X = 17.05
218			Y = 0.
219			Z = 18.5
220			CALL GSPOS('STP3',N,'SID1',X,Y,Z,4,'ONLY')
221
222			*
223			* Position Al. TRAP (STA1) in PLA1 and rotate
224			X = 18.97
225			Y = 0.
226			Z = -3.7
227			CALL GSPOS('STA1',N,'PLA1',X,Y,Z,5,'ONLY')
228			*
229			* Position Al. TRAP (STA2) in PLA1 and rotate
230			X = 20.
231			Y = 0.
232			Z = 15.6
233			CALL GSPOS('STA2',N,'PLA1',X,Y,Z,4,'ONLY')
234			*
235			* Position Al. TRAP (STA3) in PLA1 and rotate
236			X = 17.45
237			Y = 0.
238			Z = 19.3
239			CALL GSPOS('STA3',N,'PLA1',X,Y,Z,4,'ONLY')
240			*
241			* Position 'air' TRAP (STL1) in STRP and rotate
242			X = 19.2975
243			Y = 0.
244			Z = -3.875
245			CALL GSPOS('STL1',N,'STRP',X,Y,Z,5,'ONLY')
246			*
247			* Position 'air' TRAP (STL2) in STRP and rotate
248			X = 20.35
249			Y = 0.
250			Z = 15.6
251			CALL GSPOS('STL2',N,'STRP',X,Y,Z,4,'ONLY')
252			*
253			* Position 'air' TRAP (STL3) in STRP and rotate
254			X = 17.5675
255			Y = 0.
256			Z = 19.475
257			CALL GSPOS('STL3',N,'STRP',X,Y,Z,4,'ONLY')
258			*END: JeL.
259			C*
260			C NMED=MAL
261			C*
262			C* Define the CAS cradle volume
263			C*
264			C CALL GSVOLU('CSSX','BOX ',NMED,CSSX, 3,IVOLU)
265			C CALL GSVOLU('CSSY','BOX ',NMED,CSSY, 3,IVOLU)
266			C*EM:
267			C* Define the CAS scintillator volume
268			C*
269			C NMED=MSCIN
270			C CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)
271			C CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)
272			C*
273			C* Positioning Scintillators in cradles
274			C*
275			C N = 1
276			C X = 0.
277			C Y = 0.
278			C Z = 0.
279			C CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')
280			C CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')
281			C*END: EM.
282			C*
283	cafagna	3.1	RETURN
284			END