/[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.5 by pam-ba, Thu May 18 10:52:32 2006 UTC revision 3.7 by pamela, Tue Jan 29 18:25:18 2008 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gptofv.F,v 3.4 2006/05/02 10:40:13 pam-ba Exp $  * $Id: gptofv.F,v 3.6 2006/11/16 18:45:29 pam-ba Exp $
3  *  *
4  * $Log: gptofv.F,v $  * $Log: gptofv.F,v $
5    * Revision 3.6  2006/11/16 18:45:29  pam-ba
6    * Simulated an aluminum container for S4
7    *
8    * Revision 3.5  2006/05/18 10:52:32  pam-ba
9    * TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
10    *
11  * Revision 3.4  2006/05/02 10:40:13  pam-ba  * Revision 3.4  2006/05/02 10:40:13  pam-ba
12  * TOF geometry and position updated and a new material, the mylar, added  * TOF geometry and position updated and a new material, the mylar, added
13  *  *
# Line 98  Line 104 
104        CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU)        CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU)
105  *  *
106  *  *
107    *  S4 container volume.
108    *
109          NMED=MAL
110          CALL GSVOLU('S4AL','BOX ',NMED,S4AL, 3,IVOLU)
111    *
112    *
113  *  Change medium just for S4.  *  Change medium just for S4.
114  *  *
115        NMED=MSCIN        NMED=MSCIN
# Line 173  c$$$>>>>>>> 3.3 Line 185  c$$$>>>>>>> 3.3
185        Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)        Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
186        DO I=1,NS12X        DO I=1,NS12X
187           N= I           N= I
188           Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)  c ab mod         Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)
189             Y= -S1A(2) + (2*I-1)*S12M(2)
190           CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
191        ENDDO        ENDDO
192  *  *
# Line 181  c$$$>>>>>>> 3.3 Line 194  c$$$>>>>>>> 3.3
194        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
195        DO I=1,NS11Y        DO I=1,NS11Y
196           N= I           N= I
197           X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)  c ab mod         X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)
198             X= -S1A(1) + (2*I-1)*S11M(1)
199           CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
200        ENDDO        ENDDO
201  *  *
# Line 191  c$$$>>>>>>> 3.3 Line 205  c$$$>>>>>>> 3.3
205        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+2.*S11M(3)+HPVCMAX/2.        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+2.*S11M(3)+HPVCMAX/2.
206        CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY')        CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY')
207  C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)  C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
208  C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
209  C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)  C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
210  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
211  *        *
212        X=0.        X=0.
213        Y=0.        Y=0.
214        Z=0.        Z=0.
215        CALL GSPOS('S2A ',N,'S2  ',X,Y,Z,0,'ONLY')            CALL GSPOS('S2A ',N,'S2  ',X,Y,Z,0,'ONLY')
216  *  *
217        Y= 0.        Y= 0.
218        Z=-S2A(3)+HGLUTOF+S22M(3)        Z=-S2A(3)+HGLUTOF+S22M(3)
219        DO I=1,NS22Y        DO I=1,NS22Y
220           N= I           N= I
221           X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1)          c ab mod         X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1)
222             X= -S2A(1) + (2*I-1)*S22M(1)
223           CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
224        ENDDO        ENDDO
225  *  *
# Line 212  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z, Line 227  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,
227        Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)        Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
228        DO I=1,NS21X        DO I=1,NS21X
229           N= I           N= I
230           Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)  c ab mod         Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)
231             Y= -S2A(2) + (2*I-1)*S21M(2)
232           CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
233        ENDDO        ENDDO
234  *  *
# Line 222  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z, Line 238  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,
238        Z=-S2A(3)+HGLUTOF+2*S22M(3)+2.*S21M(3)+HPVCMIN/2.        Z=-S2A(3)+HGLUTOF+2*S22M(3)+2.*S21M(3)+HPVCMIN/2.
239        CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY')        CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY')
240  C      Z=-S2A(3)+HGLUTOF+S22M(3)  C      Z=-S2A(3)+HGLUTOF+S22M(3)
241  C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')
242  C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)  C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
243  C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
244  *  *
245        X= 0.        X= 0.
246  ccC      Z=-S3(3)+HAIRS3+S32M(3)  ccC      Z=-S3(3)+HAIRS3+S32M(3)
# Line 246  ccC      Z= -S3(3)+HAIRS3+2*S32M(3)+HAIR Line 262  ccC      Z= -S3(3)+HAIRS3+2*S32M(3)+HAIR
262    
263  *  *
264  C      Z=-S3(3)+HAIRS3+S32M(3)  C      Z=-S3(3)+HAIRS3+S32M(3)
265  C      CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')
266  C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)  C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
267  C      CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')
268  ****************************************  ****************************************
269  c ml: 21/12/04:  c ml: 21/12/04:
270        X=0.        X=0.
# Line 291  C      ENDDO Line 307  C      ENDDO
307        N=1        N=1
308  C      DO I=1,NS22Y  C      DO I=1,NS22Y
309  C         N= I  C         N= I
310  C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)          C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)
311        CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')        CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
312  C      ENDDO  C      ENDDO
313        Y= 0.        Y= 0.
# Line 313  C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S3 Line 329  C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S3
329        Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)        Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
330        CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')        CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
331  C      ENDDO  C      ENDDO
332          X=0.
333          Y=0.
334          Z=0.
335          CALL GSPOS('S4  ',N,'S4AL',X,Y,Z,0,'ONLY')
336  *  *
337        RETURN        RETURN
338        END        END
Legend:
Removed from v.3.5  
changed lines
  Added in v.3.7
  ViewVC Help
Powered by ViewVC 1.1.23