--- gpamela/gptof/gptofv.F 2006/05/02 10:40:13 3.4 +++ gpamela/gptof/gptofv.F 2008/01/29 18:25:18 3.7 @@ -1,7 +1,16 @@ * -* $Id: gptofv.F,v 3.3 2006/04/09 23:29:22 cafagna Exp $ +* $Id: gptofv.F,v 3.6 2006/11/16 18:45:29 pam-ba Exp $ * * $Log: gptofv.F,v $ +* Revision 3.6 2006/11/16 18:45:29 pam-ba +* Simulated an aluminum container for S4 +* +* Revision 3.5 2006/05/18 10:52:32 pam-ba +* TOF geometry completed and a new material, the polystyrene (density 35 g/l), added +* +* Revision 3.4 2006/05/02 10:40:13 pam-ba +* TOF geometry and position updated and a new material, the mylar, added +* * Revision 3.3 2006/04/09 23:29:22 cafagna * Several new things, among this: ND and CARD * @@ -51,9 +60,9 @@ * CALL GSVOLU('S2A ','BOX ',NMED,S2A, 3,IVOLU) * -* Define the S3A volume +* Define the S3 volume * - CALL GSVOLU('S3A ','BOX ',NMED,S3A, 3,IVOLU) + CALL GSVOLU('S3 ','BOX ',NMED,S3, 3,IVOLU) * NMED=MMYL * @@ -81,6 +90,25 @@ * CALL GSVOLU('S32M','BOX ',NMED,S32M, 3,IVOLU) * +* Define the POL1 volume +* + NMED=MPOL + CALL GSVOLU('POL1','BOX ',NMED,POL1, 3,IVOLU) +* +* Define the POL2 volume +* + CALL GSVOLU('POL2','BOX ',NMED,POL2, 3,IVOLU) +* +* Define the POLY volume +* + CALL GSVOLU('POLY','BOX ',NMED,POLY, 3,IVOLU) +* +* +* S4 container volume. +* + NMED=MAL + CALL GSVOLU('S4AL','BOX ',NMED,S4AL, 3,IVOLU) +* * * Change medium just for S4. * @@ -139,33 +167,6 @@ c$$$ CALL GSVOLU('SC22','BOX ',NMED,SC22, 3,IVOLU) c$$$ 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. * * Positioning volumes * @@ -174,35 +175,51 @@ Z=0. CALL GSPOS('S1A ',N,'S1 ',X,Y,Z,0,'ONLY') * + N=1 + X=0. + Y=0. + Z=-S1A(3)+HGLUTOF+HPVCMIN/2. + CALL GSPOS('POL2',N,'S1A ',X,Y,Z,0,'ONLY') +* X=0. Z=-S1A(3)+HGLUTOF+HPVCMIN+S12M(3) DO I=1,NS12X N= I - Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2) +c ab mod Y= -S1A(2) + HMYLTOF + (2*I-1)*S12M(2) + Y= -S1A(2) + (2*I-1)*S12M(2) CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY') ENDDO +* Y= 0. Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+S11M(3) DO I=1,NS11Y N= I - X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1) +c ab mod X= -S1A(1) + HMYLTOF + (2*I-1)*S11M(1) + X= -S1A(1) + (2*I-1)*S11M(1) CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY') ENDDO +* + N=1 + X=0. + Y=0. + Z=-S1A(3)+HGLUTOF+HPVCMIN+2*S12M(3)+2.*S11M(3)+HPVCMAX/2. + CALL GSPOS('POL1',N,'S1A ',X,Y,Z,0,'ONLY') C Z=-S1A(3)+HGLUTOF+HPVCMAX+S12M(3) -C CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY') +C CALL GSPOS('S12M',N,'S1A ',X,Y,Z,0,'ONLY') C Z=-S1A(3)+HGLUTOF+0.4+2*S12M(3)+S11M(3) -C CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY') -* +C CALL GSPOS('S11M',N,'S1A ',X,Y,Z,0,'ONLY') +* X=0. Y=0. Z=0. - CALL GSPOS('S2A ',N,'S2 ',X,Y,Z,0,'ONLY') + CALL GSPOS('S2A ',N,'S2 ',X,Y,Z,0,'ONLY') * Y= 0. Z=-S2A(3)+HGLUTOF+S22M(3) DO I=1,NS22Y N= I - X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1) +c ab mod X= -S2A(1) + HMYLTOF + (2*I-1)*S22M(1) + X= -S2A(1) + (2*I-1)*S22M(1) CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY') ENDDO * @@ -210,34 +227,44 @@ Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3) DO I=1,NS21X N= I - Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2) +c ab mod Y= -S2A(2) + HMYLTOF+ (2*I-1)*S21M(2) + Y= -S2A(2) + (2*I-1)*S21M(2) CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY') ENDDO +* + N=1 + X=0. + Y=0. + Z=-S2A(3)+HGLUTOF+2*S22M(3)+2.*S21M(3)+HPVCMIN/2. + CALL GSPOS('POLY',N,'S2A ',X,Y,Z,0,'ONLY') C Z=-S2A(3)+HGLUTOF+S22M(3) -C CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY') +C CALL GSPOS('S22M',N,'S2A ',X,Y,Z,0,'ONLY') C Z=-S2A(3)+HGLUTOF+2*S22M(3)+S21M(3) -C CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY') +C CALL GSPOS('S21M',N,'S2A ',X,Y,Z,0,'ONLY') * X= 0. - Z=-S3A(3)+HAIRS3+S32M(3) +ccC Z=-S3(3)+HAIRS3+S32M(3) + Z=-S3(3)+S32M(3) DO I=1,NS32X N= I - Y= -S3A(2)+ (2*I-1)*S32M(2) - CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY') + Y= -S3(2)+ (2*I-1)*S32M(2) + CALL GSPOS('S32M',N,'S3 ',X,Y,Z,0,'ONLY') ENDDO * Y= 0. - Z= -S3A(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3) +ccC Z= -S3(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3) + Z=-S3(3)+2*S32M(3)+HAIRS3+S31M(3) DO I=1,NS31Y N= I - X= -S3A(1)+(2*I-1)*S31M(1) - CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY') + X= -S3(1)+(2*I-1)*S31M(1) + CALL GSPOS('S31M',N,'S3 ',X,Y,Z,0,'ONLY') ENDDO + * -C Z=-S3A(3)+HAIRS3+S32M(3) -C CALL GSPOS('S32M',N,'S3A ',X,Y,Z,0,'ONLY') +C Z=-S3(3)+HAIRS3+S32M(3) +C CALL GSPOS('S32M',N,'S3 ',X,Y,Z,0,'ONLY') C Z=-S32(3)+HAIRS3+2*S32M(3)+HAIRS3+S31M(3) -C CALL GSPOS('S31M',N,'S3A ',X,Y,Z,0,'ONLY') +C CALL GSPOS('S31M',N,'S3 ',X,Y,Z,0,'ONLY') **************************************** c ml: 21/12/04: X=0. @@ -280,44 +307,9 @@ N=1 C DO I=1,NS22Y C N= I -C X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1) +C X= -S22M(1) + HMYLTOF + (2*I-1)*S22Y(1) CALL GSPOS('S22Y',N,'S22M',X,Y,Z,0,'ONLY') 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 Y= 0. Z= -S31M(3)+HMYLTOF+S31Y(3) N=1 @@ -337,6 +329,10 @@ Y= -S32M(2) + 2*(HMYLTOF/2)+ S32X(2) CALL GSPOS('S32X',N,'S32M',X,Y,Z,0,'ONLY') C ENDDO + X=0. + Y=0. + Z=0. + CALL GSPOS('S4 ',N,'S4AL',X,Y,Z,0,'ONLY') * RETURN END