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

Diff of /gpamela/gpspe/gpspev.F

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

revision 3.1.1.1 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.4 by pam-ba, Mon Dec 5 12:15:21 2005 UTC
# Line 1  Line 1 
1    
2  *  *
3  * $Id$  * $Id: gpspev.F,v 3.3 2005/07/25 11:53:21 cafagna Exp $
4    *
5    * $Log: gpspev.F,v $
6    * Revision 3.3  2005/07/25 11:53:21  cafagna
7    * Several updates. See history for details
8    *
9    * Revision 3.2  2005/06/21 02:42:27  cafagna
10    * Major modification to the geometry and to the random number chain
11    *
12    * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
13    * First GPAMELA release on CVS
14  *  *
 * $Log$  
15  *  *
16  *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
17  *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 27  Line 37 
37        NMED= MN2        NMED= MN2
38        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
39  *  *
40    *  Define the MGFR     volume
41    *
42          NMED= MAL
43          CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
44    *
45  *  Define the MGPL     volume  *  Define the MGPL     volume
46  *  *
47        NMED= MMAG        NMED= MMAG
# Line 41  Line 56 
56  *  *
57        NMED= MAL        NMED= MAL
58        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
59    *****
60    cC ml 20/04/05:
61  *  *
62  *  Define the TRPL     volume  cc*  Define the TRPL     volume
63    *  Define the TPAS and TPAI volume
64  *  *
65        NMED= MN2        NMED= MN2
66        CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)  CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
67          CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
68          NMED= MAL
69          CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
70    CC end ml.
71    ****************
72  *  *
73  *  Define the TRSL     volume  *  Define the TRSL     volume
74  *  *
75        NMED= MN2  c ml: 26/11/04:
76    c      NMED= MN2
77          NMED=MSIT
78        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
79  *  *
80    * Define the THBP volume
81    *
82          NMED=MG10C
83          CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
84    C end ml.
85    *
86  *  Define the TSPA     volume  *  Define the TSPA     volume
87  *  *
88        NMED= MSIT        NMED= MSIT
# Line 62  Line 93 
93        NMED= MCF        NMED= MCF
94        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)        CALL GSVOLU('TRCP','BOX ',NMED,TRCP, 3,IVOLU)
95  *  *
96  * Positioning volume TSPA     into the mother TRSL  *  Define the MGPA     volume
97  *  *
98        N= 1        NMED= MAL
99        X= 0.        CALL GSVOLU('MGPA','BOX ',NMED,MGPA, 3,IVOLU)
       Y= TRSL(2) - TSPA(2)  
       Z= 0.  
       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  
       N= 2  
       Y= Y - 2.*TSPA(2)  
       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  
100  *  *
101  * Positioning volume TRCP     into the mother TRSL  *  Define the TPGA     volume
102  *  *
103        N= 1        NMED= MN2
104        Z= 0.        CALL GSVOLU('TPGA','BOX ',NMED,TPGA, 3,IVOLU)
       Y= 0.  
       X= -TRSL(1) + TRCP(1)  
       CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')  
       N= 2  
       X= -X  
       CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')  
105  *  *
106  * Positioning volume TRSL     into the mother TRPL  *  Define the TPGU     volume
107  *  *
108          NMED= MN2
109          CALL GSVOLU('TPGU','BOX ',NMED,TPGU, 3,IVOLU)
110    *
111    *  Define the TPGD     volume
112    *
113          NMED= MN2
114          CALL GSVOLU('TPGD','BOX ',NMED,TPGD, 3,IVOLU)
115    
116    *
117    * Positioning volume TSPA     into the mother TRSL
118    *
119          N= 1
120          X= 0.
121        Y= 0.        Y= 0.
122        Z= 0.        Z= 0.
123          CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
124    C #       N= 2
125    C # C      Y= Y - 2.*TSPA(2)
126    C #       Y=-TRSL(2)/2
127    C #       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
128    c ml: 26/11/04:
129    * Positioning volume TRCP into mother TPAS
130    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
131    C #       Z=-TRPB(3)+0.11+TRCP(3)
132    C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
133    CC end ml.
134    c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
135          Z=TRPB(3)-2*TPAI(3)-TRCP(3)
136          y=-0.225-TPAS(2)+TRCP(2)
137    C # end CAF mod
138          DO I=1,4
139    CC ml 20/04/05:        
140    CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
141    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
142    C #          X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
143    c ml: 19/9/05: introducing the glue between TRCP and the silicon ladders.
144    c              for now the glue is a space of nitrogen of 0.015 mm
145    c         X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
146             X=-TPAS(1)+2*(I-1)*0.0015+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
147    CC end ml
148             CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
149    C # end CAF mod
150          ENDDO  
151    cc ml 20/04/05:
152    * Positioning volume TRSL     into the mother TPAS
153    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
154    C #       Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
155    C      Y=2.75
156    CC      Z= 0.
157    C      Z=-0.005
158    c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
159          Z=TRPB(3)-2*TPAI(3)-TRCP(3)
160        DO I=1, NTRSL        DO I=1, NTRSL
161           X= -TRPL(1) + TRSL(1) +2.*(I-1)*TRSL(1)  c ml: 20/9/05 introducing the glue between TRCP and the silicon ladders.
162           CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')  c         X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
163    c     +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
164             X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
165         +   (2*(MOD(I-1,3)+1)-1)*0.0015+(2*(MOD(I-1,3)+1)-1)*TRSL(1)        
166             IF(I.LE.3) THEN
167    C ml: 22/9/05:
168    c            Y=TPAS(2)-TRSL(2)
169                Y=TPAS(2)-0.1485-TRSL(2)
170             ELSE
171    C            Y=TPAS(2)-3*TRSL(2)
172                Y=TPAS(2)-0.1485-0.0015-3*TRSL(2)
173             ENDIF
174             CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
175        ENDDO        ENDDO
176    C # end CAF mod
177    **********
178    cc ml: 20/04/05
179    * Positioning volume THBP     into the mother TPAS
180    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
181    C #       Y=-TRPB(2)+1.5+THBP(2)
182    C # cc      Z= 0.
183    C #       Z=-TRPB(3)+0.11+TRCP(3)
184    c ml: 22/9/05:
185    c      Y=-TPAS(2)+THBP(2)
186          Y=-TPAS(2)+0.1485+THBP(2)
187    C      Z=-0.005
188    c      Z=-TRPB(3)+2*TPAI(3)+TRCP(3)
189          Z=TRPB(3)-2*TPAI(3)-TRCP(3)
190          DO I=1, NTHBP
191    CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)
192    C #          X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)
193    c ml: 21/9/05 introducing the glue between the carbon barrel and the silicon
194    c       ladders
195    c         X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)
196             X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)+
197         +      (2*(MOD(I-1,3)+1)-1)*0.0015
198             CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
199          ENDDO
200    c end ml.
201    cc ml: 20/04/05
202  *  *
203  * Positioning volume MGPI     into the mother MGPL  * Positioning volume TPAI into the mother TPAS
204    *
205    C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
206          X=0
207          Y=-TPAS(2)+TPAI(2)
208    c      Z=-TPAS(3)+TPAI(3)
209          Z=TRPB(3)-TPAI(3)
210    C #       X=-TRPB(1)+2.29+TPAI(1)
211    C #       Y=-TRPB(2)+1.5+2.5+TPAI(2)
212    C #       Z=-TRPB(3)+TPAI(3)
213          CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')
214    C # end CAF mod
215    *
216    * Positioning volume TPAS into the mother TRPB
217    *
218    C # Mod By Caf on 15th Jul 2005
219          X=0
220          Y=-TRPB(2)+1.5+TPAS(2)
221          Z=0
222    C #       X=-TRPB(1)+2.29+TPAS(1)
223    C #       Y=-TRPB(2)+1.5+TPAS(2)
224    C #       Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
225          CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
226    c end ml.
227    *
228    * Positioning volume MGPA     into the mother MGPL
229  *  *
230        N= 1        N= 1
231        X= 0.        X= 0.
232        Y= 0.        Y= 0.
233        Z= 0.        Z= 0.
234        CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')        CALL GSPOS('MGPA',N,'MGPL',X,Y,Z,0,'ONLY')
235    *
236    * Positioning volume TPGD     into the mother MGFR
237  *  *
238  * Positioning volume TRPL     into the mother TRPB        X=0.
239          Z=-MGFR(3)+TPGD(3)
240          Y=0.
241          CALL GSPOS('TPGD',N,'MGFR',X,Y,Z,0,'ONLY')
242    *
243    * Positioning volume MGPI     into the mother MGPL
244  *  *
245        N= 1        N= 1
246          X= 0.
247          Y= 0.
248    CC ML 10/11/05:  Positioning volume MGPI     into the mother MGFR
249        Z= 0.        Z= 0.
250        CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')        CALL GSPOS('MGPI',N,'MGPA',X,Y,Z,0,'ONLY')
251    C      Z= -MGFR(3)+2*TPGD(3)+MGPI(3)
252    C      CALL GSPOS('MGPI',N,'MGFR',X,Y,Z,0,'ONLY')
253    *
254    * Positioning volume MGPL     into the mother MGFR
255    *
256    c ml: 10/11/05
257    c      Z=-MGFR(3)+MGPL(3)
258          Z=-MGFR(3)+2*TPGD(3)+MGPI(3)
259          CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
260    *
261    * Positioning volume TPGA     into the mother MGFR
262    *
263    C ML: 10/11/05
264          X=0.
265          Y=MGFR(2)-2*(MGFR(2)-TRPB(2))-TPGA(2)
266    C      Z=-MGFR(3)+2*MGPL(3)+TPGA(3)
267          Z=MGFR(3)-2*TPGU(3)-2*TRPB(3)-TPGA(3)
268    C      Y=-0.75
269          CALL GSPOS('TPGA',N,'MGFR',X,Y,Z,0,'ONLY')
270  *  *
271  * Positioning volumes TRPB & MGPL into the mother SPEB  * Positioning volumes TRPB & MGPL into the mother SPEB
272  *  *
273        N= 1        N= 1
274        X= 0.        X= 0.
275        Y= 0.  C      Y= SPEB(2) - TRPB(2)
276        Z= SPEB(3) - TRPB(3)        Y= -MGFR(2) + TRPB(2)
277        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  C      Z= SPEB(3) - TRPB(3)
278    c      Z=-MGFR(3)+2*MGPL(3)+2*TPGA(3)+TRPB(3)
279          Z=MGFR(3)-2*TPGU(3)-TRPB(3)
280    C      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
281          CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
282    c      CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,9,'ONLY')
283          Z=SPEB(3)
284        DO I=1, 5        DO I=1, 5
285           Z= Z - TRPB(3) - MGPL(3)           X= 0.
286             Y= 0.
287    C         Z= Z - TRPB(3) - MGFR(3)
288             Z=Z-MGFR(3)
289           N= I           N= I
290           CALL GSPOS('MGPL',N,'SPEB',X,Y,Z,0,'ONLY')           CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
291           Z= Z - MGPL(3) - TRPB(3)  C         Y= SPEB(2) - TRPB(2)
292           N= (I+1)  C         Z= Z - MGFR(3) - TRPB(3)
293           CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')  C         N= (I+1)
294    C         IF(I.EQ.5) THEN
295    CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
296    C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,10,'ONLY')
297    C         ELSE
298    CC            CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
299    C            CALL GSPOS('TRPB',N,'MGFR',X,Y,Z,0,'ONLY')
300    C         ENDIF
301             Z=Z-MGFR(3)
302        ENDDO        ENDDO
303          N=2
304          X=0.
305          Y=MGFR(2) - TRPB(2)
306          Z=Z-TRPB(3)
307    c      CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,4,'ONLY')
308          CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
309    *
310    * Positioning volume TPGU     into the mother MGFR
311    *
312          X=0.
313          Z=MGFR(3)-TPGU(3)
314          Y=-MGFR(2)+TRPB(2)
315          CALL GSPOS('TPGU',N,'MGFR',X,Y,Z,0,'ONLY')
316  *  *
317        RETURN        RETURN
318        END        END
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.4
  ViewVC Help
Powered by ViewVC 1.1.23