/[PAMELA software]/gpamela/gpcas/gpcasv.F
ViewVC logotype

Diff of /gpamela/gpcas/gpcasv.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.1 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.2 by pamela, Thu Dec 5 10:17:41 2002 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpcasv.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $
3  *  *
4  * $Log$  * $Log: gpcasv.F,v $
5    * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
6    * First GPAMELA release on CVS
7  *  *
8    *
9    *               17/10/2002  11.58.25  by  Jens Lund
10  *CMZ :  3.00/00 13/11/2000  09.08.20  by  Emiliano Mocchiutti  *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  *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  *CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna
# Line 23  Line 27 
27  *  *
28        INTEGER IROT,IVOLU,N,NMED,NUM        INTEGER IROT,IVOLU,N,NMED,NUM
29        REAL X,Y,Z        REAL X,Y,Z
30    *JeL:
31          REAL SBOX(3),STRAP(11)
32  *  *
33        NMED=MAL  C      NMED=MAL
34          NMED=MN2
35    *
36    *  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  *  *
41  *  Define the CAS cradle volume  * Side AC scintillator (CAS), 12 pieces of scint., N2 and Aluminum
42    *                            to define the right shape
43  *  *
44        CALL GSVOLU('CSSX','BOX ',NMED,CSSX, 3,IVOLU)  * Aluminum (cover) box
45        CALL GSVOLU('CSSY','BOX ',NMED,CSSY, 3,IVOLU)        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  *EM:  *EM:
177  * Define the CAS scintillator volume  * Define the CAS scintillator volume
178  *  *
179        NMED=MSCIN  C      NMED=MSCIN
180        CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)  C      CALL GSVOLU('CASX','BOX ',NMED,CASX,3,IVOLU)
181        CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)  C      CALL GSVOLU('CASY','BOX ',NMED,CASY,3,IVOLU)
182  *  *
183  * Positioning Scintillators in cradles  * Positioning Scintillators in cradles
184  *  *
# Line 43  Line 186 
186        X = 0.        X = 0.
187        Y = 0.        Y = 0.
188        Z = 0.        Z = 0.
189        CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')  C      CALL GSPOS('CASX',N,'CSSX',X,Y,Z,0,'ONLY')
190        CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')  C      CALL GSPOS('CASY',N,'CSSY',X,Y,Z,0,'ONLY')
191  *END: EM.  *END: EM.
192  *  *
193    *
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        RETURN        RETURN
284        END        END
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.2
  ViewVC Help
Powered by ViewVC 1.1.23