/[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.2 by cafagna, Tue Jun 21 02:42:27 2005 UTC revision 3.3 by cafagna, Mon Jul 25 11:53:21 2005 UTC
# Line 1  Line 1 
1    
2  *  *
3  * $Id: gpspev.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpspev.F,v 3.2 2005/06/21 02:42:27 cafagna Exp $
4  *  *
5  * $Log: gpspev.F,v $  * $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  * Revision 3.1.1.1  2002/07/11 16:02:14  cafagna
10  * First GPAMELA release on CVS  * First GPAMELA release on CVS
11  *  *
# Line 31  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 54  cc*  Define the TRPL     volume Line 62  cc*  Define the TRPL     volume
62        NMED= MN2        NMED= MN2
63  CC      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)        CALL GSVOLU('TPAS','BOX ',NMED,TPAS, 3,IVOLU)
65          NMED= MAL
66        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)        CALL GSVOLU('TPAI','BOX ',NMED,TPAI, 3,IVOLU)
67  CC end ml.  CC end ml.
68  ****************  ****************
# Line 67  c      NMED= MN2 Line 76  c      NMED= MN2
76  *  *
77  * Define the THBP volume  * Define the THBP volume
78  *  *
79        NMED=MAL        NMED=MG10C
80        CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)        CALL GSVOLU('THBP','BOX ',NMED,THBP, 3,IVOLU)
81  C end ml.  C end ml.
82  *  *
# Line 87  C end ml. Line 96  C end ml.
96        X= 0.        X= 0.
97  c ml: 26/11/66:  c ml: 26/11/66:
98  c      Y= TRSL(2) - TSPA(2)  c      Y= TRSL(2) - TSPA(2)
99        Y=+TRSL(2)/2  C #       Y=+TRSL(2)/2
100          Y= 0.
101        Z= 0.        Z= 0.
102        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
103        N= 2  C #       N= 2
104  C      Y= Y - 2.*TSPA(2)  C # C      Y= Y - 2.*TSPA(2)
105        Y=-TRSL(2)/2  C #       Y=-TRSL(2)/2
106        CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')  C #       CALL GSPOS('TSPA',N,'TRSL',X,Y,Z,0,'ONLY')
107  c ml: 26/11/04:  c ml: 26/11/04:
108  *  *
109  c* Positioning volume TRCP     into the mother TRSL  c* Positioning volume TRCP     into the mother TRSL
# Line 105  c      N= 1 Line 115  c      N= 1
115  CC ml 20/04/05:  CC ml 20/04/05:
116  cc      Z= 0.  cc      Z= 0.
117  cc      Y= 0.  cc      Y= 0.
118        Z=-TRPB(3)+0.11+TRCP(3)  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
119        Y=-TRPB(2)+1.5-0.225+TRCP(2)  C #       Z=-TRPB(3)+0.11+TRCP(3)
120    C #       Y=-TRPB(2)+1.5-0.225+TRCP(2)
121  CC end ml.  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        DO I=1,4
126  CC ml 20/04/05:          CC ml 20/04/05:        
127  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  CC         X=-TRPL(1)+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)
128           X=-TRPB(1)+2.29+(2*I-1)*TRCP(1)+2*(I-1)*TRSL(1)  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  CC end ml
132  c         CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')  c         CALL GSPOS('TRCP',I,'TRPL',X,Y,Z,0,'ONLY')
133           CALL GSPOS('TRCP',I,'TRPB',X,Y,Z,0,'ONLY')  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          ENDDO  
137  c      X= -TRSL(1) + TRCP(1)  c      X= -TRSL(1) + TRCP(1)
138  c      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')  c      CALL GSPOS('TRCP',N,'TRSL',X,Y,Z,0,'ONLY')
# Line 130  cc ml 20/04/05: Line 148  cc ml 20/04/05:
148  c ml: 26/11/04:  c ml: 26/11/04:
149  C      Y= 0.  C      Y= 0.
150  cc      Y=THBP(2)  cc      Y=THBP(2)
151        Y=-TRPB(2)+1.5+2*THBP(2)+TRSL(2)  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  C      Y=2.75
154  CC      Z= 0.  CC      Z= 0.
155        Z=-TRPB(3)+0.11+TRCP(3)  C #       Z=-TRPB(3)+0.11+TRCP(3)
156          Z=-0.005
157        DO I=1, NTRSL        DO I=1, NTRSL
158  c         X= -TRPL(1) + TRSL(1) +2.*(I-1)*TRSL(1)           X=-TPAS(1)+2*(MOD(I-1,3)+1)*TRCP(1)+
159  cc         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*TRSL(1)       +        (2*(MOD(I-1,3)+1)-1)*TRSL(1)
160           X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*TRSL(1)           IF(I.LE.3) THEN
161  cc         CALL GSPOS('TRSL',I,'TRPL',X,Y,Z,0,'ONLY')              Y=TPAS(2)-TRSL(2)
162           CALL GSPOS('TRSL',I,'TRPB',X,Y,Z,0,'ONLY')           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  cc ml: 20/04/05
170  *  *
# Line 148  cc ml: 20/04/05 Line 172  cc ml: 20/04/05
172  * Positioning volume THBP     into the mother TRPB  * Positioning volume THBP     into the mother TRPB
173  *  *
174  cc      Y=-TRSL(2)  cc      Y=-TRSL(2)
175        Y=-TRPB(2)+1.5+THBP(2)  C # Caf mod on 19 Jul 2005 all volumes are placed inside the TPAS one
176  cc      Z= 0.  C #       Y=-TRPB(2)+1.5+THBP(2)
177        Z=-TRPB(3)+0.11+TRCP(3)  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        DO I=1, NTHBP
182  CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)  CC         X= -TRPL(1) + 2*I*TRCP(1)+(2*I-1)*THBP(1)
183           X=-TRPB(1)+2.29+2*I*TRCP(1)+(2*I-1)*THBP(1)  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')  cc         CALL GSPOS('THBP',I,'TRPL',X,Y,Z,0,'ONLY')
186           CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')  C #          CALL GSPOS('THBP',I,'TRPB',X,Y,Z,0,'ONLY')
187             CALL GSPOS('THBP',I,'TPAS',X,Y,Z,0,'ONLY')
188        ENDDO        ENDDO
189  c end ml.  c end ml.
190  cc ml: 20/04/05  cc ml: 20/04/05
191  *  *
192  * Positioning volume TPAI into the mother TRPB  * Positioning volume TPAI into the mother TRPB
193  *  *
194        X=-TRPB(1)+2.29+TPAI(1)  C # Mod By Caf on 19th Jul 2005 all volumes inside the TPAS
195        Y=-TRPB(2)+1.5+2.5+TPAI(2)        X=0
196        Z=-TRPB(3)+TPAI(3)        Y=-TPAS(2)+TPAI(2)
197        CALL GSPOS('TPAI',I,'TRPB',X,Y,Z,0,'ONLY')        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  * Positioning volume TPAS into the mother TRPB
205  *  *
206        X=-TRPB(1)+2.29+TPAS(1)  C # Mod By Caf on 15th Jul 2005
207          X=0
208        Y=-TRPB(2)+1.5+TPAS(2)        Y=-TRPB(2)+1.5+TPAS(2)
209        Z=-TRPB(3)+0.11+2*TRCP(3)+TPAS(3)        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')        CALL GSPOS('TPAS',I,'TRPB',X,Y,Z,0,'ONLY')
214  c end ml.  c end ml.
215    
# Line 183  c end ml. Line 221  c end ml.
221        Y= 0.        Y= 0.
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 MGPL     into the mother MGFR
226    *
227          CALL GSPOS('MGPL',N,'MGFR',X,Y,Z,0,'ONLY')
228  C Ml: 20/04/05:  C Ml: 20/04/05:
229  c*  c*
230  c* Positioning volume TRPL     into the mother TRPB  c* Positioning volume TRPL     into the mother TRPB
# Line 196  c end ml. Line 238  c end ml.
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.2  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23