1 |
|
2 |
* |
3 |
* $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 |
* |
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 |
***** |
49 |
cC ml 20/04/05: |
50 |
* |
51 |
cc* Define the TRPL volume |
52 |
* Define the TPAS and TPAI volume |
53 |
* |
54 |
NMED= MN2 |
55 |
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 |
* |
61 |
* Define the TRSL volume |
62 |
* |
63 |
c ml: 26/11/04: |
64 |
c NMED= MN2 |
65 |
NMED=MSIT |
66 |
CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU) |
67 |
* |
68 |
* Define the THBP volume |
69 |
* |
70 |
NMED=MAL |
71 |
CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU) |
72 |
C end ml. |
73 |
* |
74 |
* 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 |
c ml: 26/11/66: |
89 |
c Y= TRSL(2) - TSPA(2) |
90 |
Y=+TRSL(2)/2 |
91 |
Z= 0. |
92 |
CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY') |
93 |
N= 2 |
94 |
C Y= Y - 2.*TSPA(2) |
95 |
Y=-TRSL(2)/2 |
96 |
CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY') |
97 |
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)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1) |
114 |
X=-TRPB(1)+2.29+(2*I-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) + 2*I*TRCP(1)+(2*I-1)*TRSL(1) |
140 |
X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-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 |
* |
147 |
********* Positioning volume THBP into the mother TRPL |
148 |
* Positioning volume THBP into the mother TRPB |
149 |
* |
150 |
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) + 2*I*TRCP(1)+(2*I-1)*THBP(1) |
156 |
X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-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 |
* |
163 |
* Positioning volume TPAI into the mother TRPB |
164 |
* |
165 |
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 |
* |
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 |
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 |
* |
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 |