/[PAMELA software]/gpamela/gpspe/gpspev.F
ViewVC logotype

Contents of /gpamela/gpspe/gpspev.F

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3.2 - (show annotations) (download)
Tue Jun 21 02:42:27 2005 UTC (19 years, 5 months ago) by cafagna
Branch: MAIN
CVS Tags: v4r2
Changes since 3.1: +114 -29 lines
Major modification to the geometry and to the random number chain

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

  ViewVC Help
Powered by ViewVC 1.1.23