/[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.4 by pam-ba, Tue May 2 10:40:13 2006 UTC revision 3.6 by pam-ba, Thu Nov 16 18:45:29 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gptofv.F,v 3.3 2006/04/09 23:29:22 cafagna Exp $  * $Id: gptofv.F,v 3.5 2006/05/18 10:52:32 pam-ba Exp $
3  *  *
4  * $Log: gptofv.F,v $  * $Log: gptofv.F,v $
5    * Revision 3.5  2006/05/18 10:52:32  pam-ba
6    * TOF geometry completed and a new material, the polystyrene (density 35 g/l), added
7    *
8    * Revision 3.4  2006/05/02 10:40:13  pam-ba
9    * TOF geometry and position updated and a new material, the mylar, added
10    *
11  * Revision 3.3  2006/04/09 23:29:22  cafagna  * Revision 3.3  2006/04/09 23:29:22  cafagna
12  * Several new things, among this: ND and CARD  * Several new things, among this: ND and CARD
13  *  *
# Line 51  Line 57 
57  *  *
58        CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)        CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU)
59  *  *
60  *  Define the S3A     volume  *  Define the S3     volume
61  *  *
62        CALL GSVOLU('S3A ','BOX ',NMED,S3A, 3,IVOLU)        CALL GSVOLU('S3  ','BOX ',NMED,S3, 3,IVOLU)
63  *  *
64        NMED=MMYL        NMED=MMYL
65  *  *
# Line 81  Line 87 
87  *  *
88        CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)        CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU)
89  *  *
90    *  Define the POL1     volume
91    *
92          NMED=MPOL
93          CALL GSVOLU('POL1','BOX ',NMED,POL1, 3,IVOLU)
94    *
95    *  Define the POL2     volume
96    *
97          CALL GSVOLU('POL2','BOX ',NMED,POL2, 3,IVOLU)
98    *
99    *  Define the POLY     volume
100    *
101          CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU)
102    *
103    *
104    *  S4 container volume.
105    *
106          NMED=MAL
107          CALL GSVOLU('S4AL','BOX ',NMED,S4AL, 3,IVOLU)
108    *
109  *  *
110  *  Change medium just for S4.  *  Change medium just for S4.
111  *  *
# Line 139  c$$$* Line 164  c$$$*
164  c$$$      CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)  c$$$      CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU)
165  c$$$  c$$$
166  c$$$>>>>>>> 3.3  c$$$>>>>>>> 3.3
 c*  
 c*  Define the STOF volume  
 c*        
 c      NMED=MAL  
 c      CALL GSVOLU('STOF','BOX ',NMED,STOF, 3,IVOLU)  
 c*  
 c*  Define the SGLU volume  
 c*  
 cC COLLA ARALDITICA (DA VERIFICARNE LA COMPOSIZIONE)  
 c      NMED=MSIC  
 c      CALL GSVOLU('SGLU','BOX ',NMED,SGLU, 3,IVOLU)  
 c*  
 c*  Define the SPV1 volume  
 c*  
 c      NMED=? (PVC espanso)  
 c      CALL GSVOLU('SPV1','BOX ',NMED,SPV1, 3,IVOLU)  
 c*  
 c*  Define the SMYL volume  
 c*  
 c      NMED=? (MYLAR)  
 c      CALL GSVOLU('SMYL','BOX ',NMED,SMYL, 3,IVOLU)  
 c*  
 c*  Define the SPV2 volume  
 c*  
 c      NMED=? (PVC espanso)  
 c      CALL GSVOLU('SPV2','BOX ',NMED,SPV2, 3,IVOLU)  
 cc end ml.  
167  *  *
168  * Positioning volumes  * Positioning volumes
169  *  *
# Line 174  cc end ml. Line 172  cc end ml.
172        Z=0.        Z=0.
173        CALL GSPOS('S1A ',N,'S1  ',X,Y,Z,0,'ONLY')        CALL GSPOS('S1A ',N,'S1  ',X,Y,Z,0,'ONLY')
174  *  *
175          N=1
176          X=0.
177          Y=0.
178          Z=-S1A(3)+HGLUTOF+HPVCMIN/2.
179          CALL GSPOS('POL2',N,'S1A ',X,Y,Z,0,'ONLY')
180    *
181        X=0.        X=0.
182        Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)        Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3)
183        DO I=1,NS12X        DO I=1,NS12X
# Line 181  cc end ml. Line 185  cc end ml.
185           Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)           Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2)
186           CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')
187        ENDDO        ENDDO
188    *
189        Y= 0.        Y= 0.
190        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)        Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3)
191        DO I=1,NS11Y        DO I=1,NS11Y
# Line 188  cc end ml. Line 193  cc end ml.
193           X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)           X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1)
194           CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY')
195        ENDDO        ENDDO
196    *
197          N=1
198          X=0.
199          Y=0.
200          Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+2.*S11M(3)+HPVCMAX/2.
201          CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY')
202  C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)  C      Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3)
203  C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY')    
204  C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)  C      Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3)
# Line 213  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z, Line 224  C      CALL GSPOS('S11M',N,'S1A ',X,Y,Z,
224           Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)           Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2)
225           CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')
226        ENDDO        ENDDO
227    *
228          N=1
229          X=0.
230          Y=0.
231          Z=-S2A(3)+HGLUTOF+2*S22M(3)+2.*S21M(3)+HPVCMIN/2.
232          CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY')
233  C      Z=-S2A(3)+HGLUTOF+S22M(3)  C      Z=-S2A(3)+HGLUTOF+S22M(3)
234  C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY')    
235  C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)  C      Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3)
236  C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY')    
237  *  *
238        X= 0.        X= 0.
239        Z=-S3A(3)+HAIRS3+S32M(3)  ccC      Z=-S3(3)+HAIRS3+S32M(3)
240          Z=-S3(3)+S32M(3)
241        DO I=1,NS32X        DO I=1,NS32X
242           N= I           N= I
243           Y= -S3A(2)+ (2*I-1)*S32M(2)           Y= -S3(2)+ (2*I-1)*S32M(2)
244           CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')
245        ENDDO        ENDDO
246  *  *
247        Y= 0.        Y= 0.
248        Z= -S3A(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)  ccC      Z= -S3(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
249          Z=-S3(3)+2*S32M(3)+HAIRS3+S31M(3)
250        DO I=1,NS31Y        DO I=1,NS31Y
251           N= I           N= I
252           X= -S3A(1)+(2*I-1)*S31M(1)           X= -S3(1)+(2*I-1)*S31M(1)
253           CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')           CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')
254        ENDDO        ENDDO
255    
256  *  *
257  C      Z=-S3A(3)+HAIRS3+S32M(3)  C      Z=-S3(3)+HAIRS3+S32M(3)
258  C      CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S32M',N,'S3  ',X,Y,Z,0,'ONLY')    
259  C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)  C      Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3)
260  C      CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY')      C      CALL GSPOS('S31M',N,'S3  ',X,Y,Z,0,'ONLY')    
261  ****************************************  ****************************************
262  c ml: 21/12/04:  c ml: 21/12/04:
263        X=0.        X=0.
# Line 283  C         N= I Line 303  C         N= I
303  C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)          C      X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1)        
304        CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')        CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY')
305  C      ENDDO  C      ENDDO
 c*  
 c*  Positioning volumes into S11:  
 c+  
 c      X=0.  
 c      Y=0.  
 c      Z=-S11(3)+SMYL(3)  
 c      N=1  
 c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')  
 c      Z=Z+SMYL(3)+S11Y(3)  
 c      DO I=1,NS11Y  
 c         N= I  
 c         X= -S11(1) + (2*I-1)*S11Y(1)  
 c         CALL GSPOS('S11Y',N,'S11 ',X,Y,Z,0,'ONLY')  
 c      ENDDO  
 c      N=2  
 c      Z=Z+S11Y(3)+SMYL(3)  
 c      CALL GSPOS('SMYL',N,'SMYL',X,Y,Z,0,'ONLY')  
 c      Z=Z+SMYL(3)+SPV1(3)  
 c      N=1  
 c      CALL GSPOS('SPV1',N,'SPV1',X,Y,Z,0,'ONLY')  
 c      Z=Z+SPV1(3)+SGLU(3)  
 c      N=1  
 c      CALL GSPOS('SGLU',N,'SGLU',X,,Y,Z,0,'ONLY')  
 c      Z=Z+SGLU(3)+STOF(3)  
 c      N=1  
 c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')  
 c*  
 c* Positioning volumes into S12:  
 c*  
 c      Z=-S12(3)+STOF(3)  
 c      N=2  
 c      CALL GSPOS('STOF',N,'STOF',X,Y,Z,0,'ONLY')  
 c      Z=Z+STOF(3)+SGLU(3)  
 c      CALL GSPOS('SGLU',N,'SGLU',X,Y,Z,0,'ONLY')  
 c      Z=Z+SGLU(3)+SPV2  
306        Y= 0.        Y= 0.
307        Z= -S31M(3)+HMYLTOF+S31Y(3)        Z= -S31M(3)+HMYLTOF+S31Y(3)
308        N=1        N=1
# Line 337  C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S3 Line 322  C      Y= -S32M(2) + HMYLTOF+ (2*I-1)*S3
322        Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)        Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2)
323        CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')        CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY')
324  C      ENDDO  C      ENDDO
325          X=0.
326          Y=0.
327          Z=0.
328          CALL GSPOS('S4  ',N,'S4AL',X,Y,Z,0,'ONLY')      
329  *  *
330        RETURN        RETURN
331        END        END
Legend:
Removed from v.3.4  
changed lines
  Added in v.3.6
  ViewVC Help
Powered by ViewVC 1.1.23