/[PAMELA software]/gpamela/gptof/gptofv.F
ViewVC logotype

Diff of /gpamela/gptof/gptofv.F

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

revision 3.1.1.1 by cafagna, Thu Jul 11 16:02:01 2002 UTC revision 3.3 by cafagna, Sun Apr 9 23:29:22 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gptofv.F,v 3.2 2005/06/21 02:42:29 cafagna Exp $
3    *
4    * $Log: gptofv.F,v $
5    * Revision 3.2  2005/06/21 02:42:29  cafagna
6    * Major modification to the geometry and to the random number chain
7    *
8    * Revision 3.1.1.1  2002/07/11 16:02:01  cafagna
9    * First GPAMELA release on CVS
10  *  *
 * $Log$  
11  *  *
12  *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
13  *CMZ :  2.00/00 03/03/2000  15.39.05  by  Francesco Cafagna  *CMZ :  2.00/00 03/03/2000  15.39.05  by  Francesco Cafagna
# Line 58  Line 64 
64  *  *
65  *  Start subvolume definition.  *  Start subvolume definition.
66  *  *
67    c ml: 27/10/04:
68  *  *
69  *  Define the S11X     volume  c*  Define the S11X     volume
70    *  Define the S11Y     volume
71  *  *
72        CALL GSVOLU('S11X','BOX ',NMED,S11X, 3,IVOLU)  c      CALL GSVOLU('S11X','BOX ',NMED,S11X, 3,IVOLU)
73          CALL GSVOLU('S11Y','BOX ',NMED,S11Y, 3,IVOLU)
74  *  *
75  *  Define the S12Y     volume  c*  Define the S12Y     volume
76    *  Define the S12X     volume
77  *  *
78        CALL GSVOLU('S12Y','BOX ',NMED,S12Y, 3,IVOLU)  c      CALL GSVOLU('S12Y','BOX ',NMED,S12Y, 3,IVOLU)
79          CALL GSVOLU('S12X','BOX ',NMED,S12X, 3,IVOLU)
80    c end ml.
81  *  *
82  *  Define the S21X     volume  *  Define the S21X     volume
83  *  *
# Line 83  Line 95 
95  *  *
96        CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)        CALL GSVOLU('S32X','BOX ',NMED,S32X, 3,IVOLU)
97    
98        NMED=MCF  c ml: 21/12/04:
99          NMED=MAL
100  *  *
101  *  Define the SC11     volume  *  Define the SC11     volume
102  *  *
# Line 101  Line 114 
114  *  *
115        CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)        CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
116    
117    c*
118    c*  Define the STOF volume
119    c*      
120    c      NMED=MAL
121    c      CALL GSVOLU('STOF','BOX ',NMED,STOF, 3,IVOLU)
122    c*
123    c*  Define the SGLU volume
124    c*
125    cC COLLA ARALDITICA (DA VERIFICARNE LA COMPOSIZIONE)
126    c      NMED=MSIC
127    c      CALL GSVOLU('SGLU','BOX ',NMED,SGLU, 3,IVOLU)
128    c*
129    c*  Define the SPV1 volume
130    c*
131    c      NMED=? (PVC espanso)
132    c      CALL GSVOLU('SPV1','BOX ',NMED,SPV1, 3,IVOLU)
133    c*
134    c*  Define the SMYL volume
135    c*
136    c      NMED=? (MYLAR)
137    c      CALL GSVOLU('SMYL','BOX ',NMED,SMYL, 3,IVOLU)
138    c*
139    c*  Define the SPV2 volume
140    c*
141    c      NMED=? (PVC espanso)
142    c      CALL GSVOLU('SPV2','BOX ',NMED,SPV2, 3,IVOLU)
143    cc end ml.
144  *  *
145  *  Define the SC31     volume  *  Define the SC31     volume
146  *  *
# Line 113  Line 153 
153  *  *
154  * Positioning volumes  * Positioning volumes
155  *  *
156    c ml: 21/12/04:
157        Y= 0        Y= 0
158        Z= 0        Z= 0
159        DO I=1,NS11X  c ml: 27/10/04:
160           N= I  c      DO I=1,NS11X
161           X= -S11(1) + (2*I-1)*S11X(1)        DO I=1,NS11Y
162           CALL GSPOS('S11X',N,'S11 ',X,Y,Z,0,'ONLY')           N= I
163    c         X= -S11(1) + (2*I-1)*S11X(1)
164             X= -S11(1) + (2*I-1)*S11Y(1)
165    c         CALL GSPOS('S11X',N,'S11 ',X,Y,Z,0,'ONLY')
166             CALL GSPOS('S11Y',N,'S11 ',X,Y,Z,0,'ONLY')
167        ENDDO        ENDDO
168    c end ml.
169        X= 0.        X= 0.
170        Y= 0.        Y= 0.
171        DO I=1,2        DO I=1,2
# Line 130  Line 176 
176  *  *
177        X= 0        X= 0
178        Z= 0        Z= 0
179        DO I=1,NS12Y  c ml: 27/10/04
180           N= I  c      DO I=1,NS12Y
181           Y= -S12(2) + (2*I-1)*S12Y(2)        DO I=1,NS12X
182           CALL GSPOS('S12Y',N,'S12 ',X,Y,Z,0,'ONLY')           N= I
183    c         Y= -S12(2) + (2*I-1)*S12Y(2)
184             Y= -S12(2) + (2*I-1)*S12X(2)
185    c         CALL GSPOS('S12Y',N,'S12 ',X,Y,Z,0,'ONLY')
186             CALL GSPOS('S12X',N,'S12 ',X,Y,Z,0,'ONLY')
187        ENDDO        ENDDO
188    c end ml.
189        X= 0.        X= 0.
190        Y= 0.        Y= 0.
191        DO I=1,2        DO I=1,2
# Line 147  Line 198 
198        Z= 0        Z= 0
199        DO I=1,NS21X        DO I=1,NS21X
200           N= I           N= I
201           X= -S21(1) + (2*I-1)*S21X(1)  c ml: 27/10/04         X= -S21(1) + (2*I-1)*S21X(1)
202             Y= -S21(2) + (2*I-1)*S21X(2)
203    c end ml.
204           CALL GSPOS('S21X',N,'S21 ',X,Y,Z,0,'ONLY')           CALL GSPOS('S21X',N,'S21 ',X,Y,Z,0,'ONLY')
205        ENDDO        ENDDO
206        X= 0.        X= 0.
# Line 162  Line 215 
215        Z= 0        Z= 0
216        DO I=1,NS22Y        DO I=1,NS22Y
217           N= I           N= I
218           Y= -S22(2) + (2*I-1)*S22Y(2)  c ml: 27/10/04         Y= -S22(2) + (2*I-1)*S22Y(2)
219             X= -S22(1) + (2*I-1)*S22Y(1)        
220    c end ml.  
221           CALL GSPOS('S22Y',N,'S22 ',X,Y,Z,0,'ONLY')           CALL GSPOS('S22Y',N,'S22 ',X,Y,Z,0,'ONLY')
222        ENDDO        ENDDO
223        X= 0.        X= 0.
# Line 172  Line 227 
227           Z= (-1.)**I*S22(3) + (-1.)**(I+1)*SC22(3)           Z= (-1.)**I*S22(3) + (-1.)**(I+1)*SC22(3)
228           CALL GSPOS('SC22',N,'S22 ',X,Y,Z,0,'ONLY')           CALL GSPOS('SC22',N,'S22 ',X,Y,Z,0,'ONLY')
229        ENDDO        ENDDO
230  *  c*
231    c*  Positioning volumes into S11:
232    c+
233    c      X=0.
234    c      Y=0.
235    c      Z=-S11(3)+SMYL(3)
236    c      N=1
237    c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
238    c      Z=Z+SMYL(3)+S11Y(3)
239    c      DO I=1,NS11Y
240    c         N= I
241    c         X= -S11(1) + (2*I-1)*S11Y(1)
242    c         CALL GSPOS('S11Y',N,'S11 ',X,Y,Z,0,'ONLY')
243    c      ENDDO
244    c      N=2
245    c      Z=Z+S11Y(3)+SMYL(3)
246    c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')
247    c      Z=Z+SMYL(3)+SPV1(3)
248    c      N=1
249    c      CALL GSPOS('SPV1',N,'SPV1',X,Y,Z,0,'ONLY')
250    c      Z=Z+SPV1(3)+SGLU(3)
251    c      N=1
252    c      CALL GSPOS('SGLU',N,'SGLU',X,,Y,Z,0,'ONLY')
253    c      Z=Z+SGLU(3)+STOF(3)
254    c      N=1
255    c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
256    c*
257    c* Positioning volumes into S12:
258    c*
259    c      Z=-S12(3)+STOF(3)
260    c      N=2
261    c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')
262    c      Z=Z+STOF(3)+SGLU(3)
263    c      CALL GSPOS('SGLU',N,'SGLU',X,Y,Z,0,'ONLY')
264    c      Z=Z+SGLU(3)+SPV2
265        X= 0        X= 0
266        Z= 0        Z= 0
267        DO I=1,NS31Y        DO I=1,NS31Y
268           N= I           N= I
269           Y= -S31(2) + (2*I-1)*S31Y(2)  c ml: 27/10/04:
270    c         Y= -S31(2) + (2*I-1)*S31Y(2)
271             X= -S31(1) + (2*I-1)*S31Y(1)
272    c end ml.
273           CALL GSPOS('S31Y',N,'S31 ',X,Y,Z,0,'ONLY')           CALL GSPOS('S31Y',N,'S31 ',X,Y,Z,0,'ONLY')
274        ENDDO        ENDDO
275        X= 0.        X= 0.
# Line 192  Line 284 
284        Z= 0        Z= 0
285        DO I=1,NS32X        DO I=1,NS32X
286           N= I           N= I
287           X= -S32(1) + (2*I-1)*S32X(1)  c ml: 27/10/04:
288    c         X= -S32(1) + (2*I-1)*S32X(1)
289             Y= -S32(2) + (2*I-1)*S32X(2)
290    c end ml.
291           CALL GSPOS('S32X',N,'S32 ',X,Y,Z,0,'ONLY')           CALL GSPOS('S32X',N,'S32 ',X,Y,Z,0,'ONLY')
292    
293        ENDDO        ENDDO
294        X= 0.        X= 0.
295        Y= 0.        Y= 0.
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23