/[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 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.3 by cafagna, Mon Jul 25 11:53:21 2005 UTC
# Line 1  Line 1 
1    
2  *  *
3  * $Id$  * $Id: gpspev.F,v 3.2 2005/06/21 02:42:27 cafagna Exp $
4    *
5    * $Log: gpspev.F,v $
6    * Revision 3.2  2005/06/21 02:42:27  cafagna
7    * Major modification to the geometry and to the random number chain
8    *
9    * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
10    * First GPAMELA release on CVS
11  *  *
 * $Log$  
12  *  *
13  *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
14  *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 34 
34        NMED= MN2        NMED= MN2
35        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)        CALL GSVOLU('SPEB','BOX ',NMED,SPEB, 3,IVOLU)
36  *  *
37    *  Define the MGFR     volume
38    *
39          NMED= MAL
40          CALL GSVOLU('MGFR','BOX ',NMED,MGFR, 3,IVOLU)
41    *
42  *  Define the MGPL     volume  *  Define the MGPL     volume
43  *  *
44        NMED= MMAG        NMED= MMAG
# Line 41  Line 53 
53  *  *
54        NMED= MAL        NMED= MAL
55        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)        CALL GSVOLU('TRPB','BOX ',NMED,TRPB, 3,IVOLU)
56    *****
57    cC ml 20/04/05:
58  *  *
59  *  Define the TRPL     volume  cc*  Define the TRPL     volume
60    *  Define the TPAS and TPAI volume
61  *  *
62        NMED= MN2        NMED= MN2
63        CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)  CC      CALL GSVOLU('TRPL','BOX ',NMED,TRPL, 3,IVOLU)
64          CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
65          NMED= MAL
66          CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
67    CC end ml.
68    ****************
69  *  *
70  *  Define the TRSL     volume  *  Define the TRSL     volume
71  *  *
72        NMED= MN2  c ml: 26/11/04:
73    c      NMED= MN2
74          NMED=MSIT
75        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)        CALL GSVOLU('TRSL','BOX ',NMED,TRSL, 3,IVOLU)
76  *  *
77    * Define the THBP volume
78    *
79          NMED=MG10C
80          CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
81    C end ml.
82    *
83  *  Define the TSPA     volume  *  Define the TSPA     volume
84  *  *
85        NMED= MSIT        NMED= MSIT
# Line 66  Line 94 
94  *  *
95        N= 1        N= 1
96        X= 0.        X= 0.
97        Y= TRSL(2) - TSPA(2)  c ml: 26/11/66:
98        Z= 0.  c      Y= TRSL(2) - TSPA(2)
99        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  C #       Y=+TRSL(2)/2
       N= 2  
       Y= Y - 2.*TSPA(2)  
       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  
 *  
 * Positioning volume TRCP     into the mother TRSL  
 *  
       N= 1  
       Z= 0.  
       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')  
 *  
 * Positioning volume TRSL     into the mother TRPL  
 *  
100        Y= 0.        Y= 0.
101        Z= 0.        Z= 0.
102          CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
103    C #       N= 2
104    C # C      Y= Y - 2.*TSPA(2)
105    C #       Y=-TRSL(2)/2
106    C #       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
107    c ml: 26/11/04:
108    *
109    c* Positioning volume TRCP     into the mother TRSL
110    CC ml 20/04/05:
111    CC* Positioning volume TRCP into mother TRPL
112    * Positioning volume TRCP into mother TRPB
113    *
114    c      N= 1
115    CC ml 20/04/05:
116    cc      Z= 0.
117    cc      Y= 0.
118    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
119    C #       Z=-TRPB(3)+0.11+TRCP(3)
120    C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
121    CC end ml.
122          Z=-0.005
123          y=-0.225-TPAS(2)+TRCP(2)
124    C # end CAF mod
125          DO I=1,4
126    CC ml 20/04/05:        
127    CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
128    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
129    C #          X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
130             X=-TPAS(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
131    CC end ml
132    c         CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')
133    C #          CALL GSPOS('TRCP',I,'TRPB',X,Y,Z,0,'ONLY')
134             CALL GSPOS('TRCP',I,'TPAS',X,Y,Z,0,'ONLY')
135    C # end CAF mod
136          ENDDO  
137    c      X= -TRSL(1) + TRCP(1)
138    c      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
139    C      N= 2
140    C      X= -X
141    C      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
142    *****
143    cc ml 20/04/05:
144    *
145    *** Positioning volume TRSL     into the mother TRPL
146    * Positioning volume TRSL     into the mother TRPB
147    *
148    c ml: 26/11/04:
149    C      Y= 0.
150    cc      Y=THBP(2)
151    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
152    C #       Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)
153    C      Y=2.75
154    CC      Z= 0.
155    C #       Z=-TRPB(3)+0.11+TRCP(3)
156          Z=-0.005
157        DO I=1, NTRSL        DO I=1, NTRSL
158           X= -TRPL(1) + TRSL(1) +2.*(I-1)*TRSL(1)           X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
159           CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')       +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
160             IF(I.LE.3) THEN
161                Y=TPAS(2)-TRSL(2)
162             ELSE
163                Y=TPAS(2)-3*TRSL(2)
164             ENDIF
165             CALL GSPOS('TRSL',I,'TPAS',X,Y,Z,0,'ONLY')
166        ENDDO        ENDDO
167    C # end CAF mod
168    **********
169    cc ml: 20/04/05
170    *
171    ********* Positioning volume THBP     into the mother TRPL
172    * Positioning volume THBP     into the mother TRPB
173    *
174    cc      Y=-TRSL(2)
175    C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
176    C #       Y=-TRPB(2)+1.5+THBP(2)
177    C # cc      Z= 0.
178    C #       Z=-TRPB(3)+0.11+TRCP(3)
179          Y=-TPAS(2)+THBP(2)
180          Z=-0.005
181          DO I=1, NTHBP
182    CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)
183    C #          X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)
184             X=-TPAS(1)+2*I*TRCP(1)+(2*I-1)*THBP(1)
185    cc         CALL GSPOS('THBP',I,'TRPL',X,Y,Z,0,'ONLY')
186    C #          CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')
187             CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
188          ENDDO
189    c end ml.
190    cc ml: 20/04/05
191    *
192    * Positioning volume TPAI into the mother TRPB
193    *
194    C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
195          X=0
196          Y=-TPAS(2)+TPAI(2)
197          Z=-TPAS(3)+TPAI(3)
198    C #       X=-TRPB(1)+2.29+TPAI(1)
199    C #       Y=-TRPB(2)+1.5+2.5+TPAI(2)
200    C #       Z=-TRPB(3)+TPAI(3)
201          CALL GSPOS('TPAI',I,'TPAS',X,Y,Z,0,'ONLY')
202    C # end CAF mod
203    *
204    * Positioning volume TPAS into the mother TRPB
205    *
206    C # Mod By Caf on 15th Jul 2005
207          X=0
208          Y=-TRPB(2)+1.5+TPAS(2)
209          Z=0
210    C #       X=-TRPB(1)+2.29+TPAS(1)
211    C #       Y=-TRPB(2)+1.5+TPAS(2)
212    C #       Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)
213          CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
214    c end ml.
215    
216  *  *
217  * Positioning volume MGPI     into the mother MGPL  * Positioning volume MGPI     into the mother MGPL
218  *  *
# Line 101  Line 222 
222        Z= 0.        Z= 0.
223        CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')        CALL GSPOS('MGPI',N,'MGPL',X,Y,Z,0,'ONLY')
224  *  *
225  * Positioning volume TRPL     into the mother TRPB  * Positioning volume MGPL     into the mother MGFR
226  *  *
227        N= 1        CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
228        Z= 0.  C Ml: 20/04/05:
229        CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')  c*
230    c* Positioning volume TRPL     into the mother TRPB
231    c*
232    c      N= 1
233    c      Z= 0.
234    c      CALL GSPOS('TRPL',N,'TRPB',X,Y,Z,0,'ONLY')
235    c end ml.
236  *  *
237  * Positioning volumes TRPB & MGPL into the mother SPEB  * Positioning volumes TRPB & MGPL into the mother SPEB
238  *  *
239        N= 1        N= 1
240        X= 0.        X= 0.
241        Y= 0.        Y= SPEB(2) - TRPB(2)
242        Z= SPEB(3) - TRPB(3)        Z= SPEB(3) - TRPB(3)
243        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')        CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
244        DO I=1, 5        DO I=1, 5
245           Z= Z - TRPB(3) - MGPL(3)           X= 0.
246             Y= 0.
247             Z= Z - TRPB(3) - MGFR(3)
248           N= I           N= I
249           CALL GSPOS('MGPL',N,'SPEB',X,Y,Z,0,'ONLY')           CALL GSPOS('MGFR',N,'SPEB',X,Y,Z,0,'ONLY')
250           Z= Z - MGPL(3) - TRPB(3)           Y= SPEB(2) - TRPB(2)
251             Z= Z - MGFR(3) - TRPB(3)
252           N= (I+1)           N= (I+1)
253           CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')           IF(I.EQ.5) THEN
254                CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,10,'ONLY')
255             ELSE
256                CALL GSPOS('TRPB',N,'SPEB',X,Y,Z,0,'ONLY')
257             ENDIF
258        ENDDO        ENDDO
259  *  *
260        RETURN        RETURN
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23