--- gpamela/gpcard/gpcardv.F 2006/04/09 23:28:47 1.1 +++ gpamela/gpcard/gpcardv.F 2006/11/30 12:29:04 1.2 @@ -4,260 +4,254 @@ * Volumes definition for the side anticoincidence counters * * Called by: GPGEO * * * -* Author: Francesco Cafagna, 08/02/96 16.58.49 * +* Authors: Francesco Cafagna, 08/02/96 16.58.49 * +* Alessandro Bruno, 20/11/06 20.30.00 * * * ************************************************************************ #include "gconst.inc" #include "gpgeo.inc" #include "gpmed.inc" #include "gprot.inc" - INTEGER IROT,IVOLU,N,NMED,NUM + INTEGER IROT,IVOLU,N,NMED,NUM,I REAL X,Y,Z REAL CAR1(11),CR1P(11) - REAL C1D1(11),C1N1(11),C1N2(11) + REAL C1D1(11),C2D1(11) REAL CAR2(11),CR2P(11) - REAL C2D1(11),C2N1(11),C2N2(11) - REAL CAR1H,CAR1D,CAR1U,CAR1T - REAL CAR2H,CAR2D,CAR2U,CAR2T - REAL BETA, ALPHA, ANGLCAR1, ANGLCAR2 + REAL C1N1(11),C1N2(11) + REAL ANGLCAR1,ANGLCAR2 + REAL H1,H2 + REAL A1,B1,C1,D1,F1,ANG1C + REAL A2,B2,C2,D2,F2,ANG2C + REAL AA, BB,ZH +* +* Scint sheets +* + H1=CAR(2)+CAR(3) + H2=CAR(1)-H1 +* +* Greater Scintillator CAR1 (type B) +* + B1=CARB(2)+CAR(4)+CAR(5)+H1/TAN(CAR(7)*DEGRAD) + C1=H2/TAN(CAR(7)*DEGRAD) + A1=CARB(1)-C1 + D1=CARB(1)-CARB(2)-CAR(4)-CAR(5) + F1=C1+H1/TAN(CAR(7)*DEGRAD) + ANG1C=ATAN((B1/2.+D1-A1/2.)/CAR(1)) + + C1D1(1) = CAR(1)/2. + C1D1(2) = -ANG1C*RADDEG + C1D1(3) = 0. + C1D1(4) = CAR(6) + C1D1(5) = A1/2. + C1D1(6) = A1/2. + C1D1(7) = 0. + C1D1(8) = CAR(6) + C1D1(9) = B1/2. + C1D1(10) = B1/2. + C1D1(11) = 0. + +* +* Smaller Scintillator CAR2 (type A) +* + B2=CARA(2)+CAR(4)+CAR(5)+H1/TAN(CAR(7)*DEGRAD) + C2=H2/TAN(CAR(7)*DEGRAD) + A2=CARA(1)-C2 + D2=CARA(1)-CARA(2)-CAR(4)-CAR(5) + F2=C2+H1/TAN(CAR(7)*DEGRAD) + ANG2C=ATAN((B2/2.+D2-A2/2.)/CAR(1)) + + C2D1(1) = CAR(1)/2. + C2D1(2) = -ANG2C*RADDEG + C2D1(3) = 0. + C2D1(4) = CAR(6) + C2D1(5) = A2/2. + C2D1(6) = A2/2. + C2D1(7) = 0. + C2D1(8) = CAR(6) + C2D1(9) = B2/2. + C2D1(10) = B2/2. + C2D1(11) = 0. * * CAR1 Al box sheet * NMED=MAL - CAR1CF=3.012 - CAR1ANGL=16.83 - CAR1H=6.7*CAR1CF - CAR1D=12*CAR1CF - CAR1U=11.3*CAR1CF - CAR1T=1.2 - CAR1(1) = 0.5*CAR1H - CAR1(2) = CAR1ANGL + + CAR1(1) = C1D1(1)+ATZ + CAR1(2) = C1D1(2) CAR1(3) = 0. - CAR1(4) = 0.5*CAR1T - CAR1(5) = 0.5*CAR1D + CAR1(4) = C1D1(4)+ATY + CAR1(5) = C1D1(5)+ATZ*(1./COS(ATAN(D1/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CAR1(6) = CAR1(5) CAR1(7) = 0. - CAR1(8) = 0.5*CAR1T - CAR1(9) = 0.5*CAR1U + CAR1(8) = C1D1(8)+ATY + CAR1(9) = C1D1(9)+ATZ*(1./COS(ATAN(D1/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CAR1(10) = CAR1(9) CAR1(11) = 0. + CALL GSVOLU('CAR1','TRAP',NMED,CAR1,11,IVOLU) * * CR1P PLASTIC BOX TO CREATE RIM * NMED=MPLAS - CR1P(1) = 0.5*(CAR1H-.7) - CR1P(2) = CAR1ANGL + + CR1P(1) = C1D1(1)+PT + CR1P(2) = C1D1(2) CR1P(3) = 0. - CR1P(4) = 0.5*(CAR1T-.4) - CR1P(5) = 0.5*(CAR1D-0.7) + CR1P(4) = C1D1(4) + CR1P(5) = C1D1(5)+PT*(1./COS(ATAN(D1/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CR1P(6) = CR1P(5) CR1P(7) = 0. - CR1P(8) = 0.5*(CAR1T-.4) - CR1P(9) = 0.5*(CAR1U-.7) + CR1P(8) = C1D1(8) + CR1P(9) = C1D1(9)+PT*(1./COS(ATAN(D1/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CR1P(10) = CR1P(9) CR1P(11) = 0. + CALL GSVOLU('CR1P','TRAP',NMED,CR1P,11,IVOLU) * -* Scint sheet -* NMED=MSCIN - C1D1(1) = 0.5*(CAR1H-1.3) - C1D1(2) = CAR1ANGL - C1D1(3) = 0. - C1D1(4) = 0.5*(CAR1T-.4) - C1D1(5) = 0.5*(CAR1D-1.3) - C1D1(6) = C1D1(5) - C1D1(7) = 0. - C1D1(8) = 0.5*(CAR1T-.4) - C1D1(9) = 0.5*(CAR1U-1.3) - C1D1(10) = C1D1(9) - C1D1(11) = 0. +* CALL GSVOLU('C1D1','TRAP',NMED,C1D1,11,IVOLU) * -* Empty volume to be subtracted from C1D1 sheet +* Empty volume to be subtracted from C1D1 and C2D1 sheets * NMED=MN2 - BETA=ATAN( - + 2.*C1D1(1)/ - + (C1D1(9)-C1D1(5)+2.*C1D1(1)*TAN(C1D1(2)*3.14/180.)) ) - A=3.*3.02 - B=3.*3.02 - C1N1(1) = 0.5*A*SIN(BETA) - C1N1(9) = 0.5*B - C1N1(2) = -RADDEG*ATAN( - + (C1N1(9)-2.*C1N1(1)/TAN(BETA))*1./(2.*C1N1(1)) ) + + AA=(ATZ/SIN(DEGRAD*CAR(7)))/(1./TAN(CAR(7)*DEGRAD)+CAR(4)/CAR(2)) + AA=CAR(2)+CAR(3)+ATZ+AA + BB=AA/TAN(CAR(7)*DEGRAD)+AA*CAR(4)/CAR(2) + + C1N1(1) = 0.5*AA + C1N1(2) = RADDEG*ATAN((BB/2.-AA/TAN(CAR(7)*DEGRAD))/AA) C1N1(3) = 0. - C1N1(4) = 0.5*(CAR1T) + C1N1(4) = CAR1(4) C1N1(5) = 0. C1N1(6) = C1N1(5) C1N1(7) = 0. - C1N1(8) = 0.5*(CAR1T) + C1N1(8) = CAR1(8) + C1N1(9) = 0.5*BB C1N1(10) = C1N1(9) C1N1(11) = 0. + CALL GSVOLU('C1N1','TRAP',NMED,C1N1,11,IVOLU) - ALPHA=ATAN( - + ABS(A*SIN(BETA)/(B - A*COS(BETA))) ) - BETA=3.14-ALPHA - A=.6*3.02 - B=1.1*3.02 - C1N2(1) = ABS(0.5*A*SIN(BETA)) - C1N2(9) = 0.5*B - C1N2(2) = -(90. - RADDEG* - + ATAN( 2.*C1N2(1) /( C1N2(9) + (2.*C1N2(1)/ABS(TAN(BETA)))) )) + AA=CAR(3) + ATZ + BB=(ATZ+CAR(3))*(CAR(5)/CAR(3)-CAR(4)/CAR(2)) + + C1N2(1) = 0.5*AA + C1N2(2) = RADDEG*ATAN((BB/2.+(CAR(3)+ATZ)*CAR(4)/CAR(2))/AA) C1N2(3) = 0. - C1N2(4) = 0.5*CAR1T + C1N2(4) = CAR1(4) C1N2(5) = 0. C1N2(6) = C1N2(5) C1N2(7) = 0. - C1N2(8) = 0.5*CAR1T + C1N2(8) = CAR1(8) + C1N2(9) = 0.5*BB C1N2(10) = C1N2(9) C1N2(11) = 0. + CALL GSVOLU('C1N2','TRAP',NMED,C1N2,11,IVOLU) * * CAR2 Al box sheet * NMED=MAL - CAR2CF=3.012 - CAR2ANGL=16.24 - CAR2H=6.7*CAR2CF - CAR2D=10.35*CAR2CF - CAR2U=10.2*CAR2CF - CAR2T=1.2 - CAR2(1) = 0.5*CAR2H - CAR2(2) = CAR2ANGL + + CAR2(1) = C2D1(1)+ATZ + CAR2(2) = C2D1(2) CAR2(3) = 0. - CAR2(4) = 0.5*CAR2T - CAR2(5) = 0.5*CAR2D + CAR2(4) = C2D1(4)+ATY + CAR2(5) = C2D1(5)+ATZ*(1./COS(ATAN(D2/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CAR2(6) = CAR2(5) CAR2(7) = 0. - CAR2(8) = 0.5*CAR2T - CAR2(9) = 0.5*CAR2U + CAR2(8) = C2D1(8)+ATY + CAR2(9) = C2D1(9)+ATZ*(1./COS(ATAN(D2/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CAR2(10) = CAR2(9) CAR2(11) = 0. + CALL GSVOLU('CAR2','TRAP',NMED,CAR2,11,IVOLU) * * CR1P PLASTIC BOX TO CREATE RIM * NMED=MPLAS - CR2P(1) = 0.5*(CAR2H-.7) - CR2P(2) = CAR2ANGL + + CR2P(1) = C2D1(1)+PT + CR2P(2) = C2D1(2) CR2P(3) = 0. - CR2P(4) = 0.5*(CAR2T-.4) - CR2P(5) = 0.5*(CAR2D-0.7) + CR2P(4) = C2D1(4) + CR2P(5) = C2D1(5)+PT*(1./COS(ATAN(D2/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CR2P(6) = CR2P(5) CR2P(7) = 0. - CR2P(8) = 0.5*(CAR2T-.4) - CR2P(9) = 0.5*(CAR2U-.7) + CR2P(8) = C2D1(8) + CR2P(9) = C2D1(9)+PT*(1./COS(ATAN(D2/CAR(1)))+ + + 1./SIN(DEGRAD*CAR(7)))/2. CR2P(10) = CR2P(9) CR2P(11) = 0. + CALL GSVOLU('CR2P','TRAP',NMED,CR2P,11,IVOLU) * -* Scint sheet -* NMED=MSCIN - C2D1(1) = 0.5*(CAR2H-1.3) - C2D1(2) = CAR2ANGL - C2D1(3) = 0. - C2D1(4) = 0.5*(CAR2T-.4) - C2D1(5) = 0.5*(CAR2D-1.3) - C2D1(6) = C2D1(5) - C2D1(7) = 0. - C2D1(8) = 0.5*(CAR2T-.4) - C2D1(9) = 0.5*(CAR2U-1.3) - C2D1(10) = C2D1(9) - C2D1(11) = 0. +* CALL GSVOLU('C2D1','TRAP',NMED,C2D1,11,IVOLU) * -* Empty volume to be subtracted from C2D1 sheet * - NMED=MN2 - BETA=ATAN( - + 2.*C2D1(1)/ - + (C2D1(9)-C2D1(5)+2.*C2D1(1)*TAN(C2D1(2)*3.14/180.)) ) - A=3.*3.02 - B=3.*3.02 - C2N1(1) = 0.5*A*SIN(BETA) - C2N1(9) = 0.5*B - C2N1(2) = -RADDEG*ATAN( - + (C2N1(9)-2.*C2N1(1)/TAN(BETA))*1./(2.*C2N1(1)) ) - C2N1(3) = 0. - C2N1(4) = 0.5*(CAR2T) - C2N1(5) = 0. - C2N1(6) = C2N1(5) - C2N1(7) = 0. - C2N1(8) = 0.5*(CAR2T) - C2N1(10) = C2N1(9) - C2N1(11) = 0. - CALL GSVOLU('C2N1','TRAP',NMED,C2N1,11,IVOLU) - ALPHA=ATAN( - + ABS(A*SIN(BETA)/(B - A*COS(BETA))) ) - BETA=3.14-ALPHA - - A=.6*3.02 - B=1.1*3.02 - C2N2(1) = ABS(0.5*A*SIN(BETA)) - C2N2(9) = 0.5*B - C2N2(2) = -(90. - RADDEG* - + ATAN( 2.*C2N2(1) /( C2N2(9) + (2.*C2N2(1)/ABS(TAN(BETA)))) )) - C2N2(3) = 0. - C2N2(4) = 0.5*CAR2T - C2N2(5) = 0. - C2N2(6) = C2N2(5) - C2N2(7) = 0. - C2N2(8) = 0.5*CAR2T - C2N2(10) = C2N2(9) - C2N2(11) = 0. - CALL GSVOLU('C2N2','TRAP',NMED,C2N2,11,IVOLU) +* The zenihtal angle of the paddles * + ANGLCAR2=CARA(3)*DEGRAD !CAR2 inclination + ANGLCAR1=CARB(3)*DEGRAD !CAR1 inclination + NMED=MN2 - CARDB(1)=1.5*CAR2(9) - CARDB(2)=1.5*CAR1(9) - CARDB(3)=CAR1(1) CALL GSVOLU('CARD','BOX ',NMED,CARDB, 3,IVOLU) * -* The zenihtal angle of the paddles -* - ANGLCAR2=ASIN( 0.5*(CAR1D-(CAR1U-2.*(C1N2(9)+C1N1(9)))) / CAR1H ) - ANGLCAR1=ASIN( 0.5*(CAR2D-(CAR2U-2.*(C2N2(9)+C2N1(9)))) / CAR2H ) -* * Define ANGLCAR rotation around Z axis * - IRZCAR2=11 - CALL GSROTM(IRZCAR2,90.,0.,(90.-ANGLCAR2*RADDEG), - + 90.,(180.-ANGLCAR2*RADDEG),90.) - IRZCAR22=12 - CALL GSROTM(IRZCAR22,90.,180.,(270.-ANGLCAR2*RADDEG), - + 270.,(180.-ANGLCAR2*RADDEG),270.) - IRZCAR1=13 - CALL GSROTM(IRZCAR1,90.,270.,(90.-ANGLCAR1*RADDEG), - + 0.,(180.-ANGLCAR1*RADDEG),0.) - IRZCAR12=14 - CALL GSROTM(IRZCAR12,90.,90.,(90.-ANGLCAR1*RADDEG), - + 180.,(180.-ANGLCAR1*RADDEG),180.) +* CAR1A + IRZCAR1B=11 + CALL GSROTM(IRZCAR1B,90.,0.,(90.-ANGLCAR1*RADDEG), + + 90.,(180.-ANGLCAR1*RADDEG),90.) +* CAR1B + IRZCAR1A=12 + CALL GSROTM(IRZCAR1A,90.,180.,(270.-ANGLCAR1*RADDEG), + + 270.,(180.-ANGLCAR1*RADDEG),270.) +* CAR2A + IRZCAR2B=13 + CALL GSROTM(IRZCAR2B,90.,270.,(90.-ANGLCAR2*RADDEG), + + 0.,(180.-ANGLCAR2*RADDEG),0.) +* CAR2B + IRZCAR2A=14 + CALL GSROTM(IRZCAR2A,90.,90.,(90.-ANGLCAR2*RADDEG), + + 180.,(180.-ANGLCAR2*RADDEG),180.) + * * Positioning in CAR1 * N=1 - X=0 + X=(ATZ-PT)*(1./COS(ATAN(D1/CAR(1)))- + + 1./SIN(DEGRAD*CAR(7)))/2. Y=0 Z=0 CALL GSPOS('CR1P',N,'CAR1',X,Y,Z,0,'MANY') N=1 - X=( CAR1(9)+CAR1(1)*TAN(CAR1(2)*DEGRAD) ) - - + (C1N1(9)+C1N1(1)*TAN(C1N1(2)*DEGRAD)) Y=0 Z=CAR1(1)-C1N1(1) + X=CAR1(9)+CAR1(1)*TAN(ANG1C)- + + (C1N1(9)-C1N1(1)*TAN(C1N1(2)*DEGRAD)) + X=-X !ALEX MOD CALL GSPOS('C1N1',N,'CAR1',X,Y,Z,0,'ONLY') N=1 - X=( CAR1(9)+CAR1(1)*TAN(CAR1(2)*DEGRAD) ) - 2.*C1N1(9) - - + (C1N2(9)+C1N2(1)*TAN(C1N2(2)*DEGRAD)) Y=0 Z=CAR1(1)-C1N2(1) + X=CAR1(9)+CAR1(1)*TAN(ANG1C)- + + (2.*C1N1(9)+C1N2(9)-C1N2(1)*TAN(C1N2(2)*DEGRAD)) + X=-X !ALEX MOD CALL GSPOS('C1N2',N,'CAR1',X,Y,Z,0,'ONLY') N=1 - X=0 + X=PT*(1./COS(ATAN(D1/CAR(1)))-1./SIN(DEGRAD*CAR(7)))/2. Y=0 Z=0 CALL GSPOS('C1D1',N,'CR1P',X,Y,Z,0,'ONLY') @@ -265,53 +259,336 @@ * Positioning in CAR2 * N=1 - X=0 + X=(ATZ-PT)*(1./COS(ATAN(D2/CAR(1)))- + + 1./SIN(DEGRAD*CAR(7)))/2. Y=0 Z=0 CALL GSPOS('CR2P',N,'CAR2',X,Y,Z,0,'MANY') - N=1 - X=( CAR2(9)+CAR2(1)*TAN(CAR2(2)*DEGRAD) ) - - + (C2N1(9)+C2N1(1)*TAN(C2N1(2)*DEGRAD)) + N=2 Y=0 - Z=CAR2(1)-C2N1(1) - CALL GSPOS('C2N1',N,'CAR2',X,Y,Z,0,'ONLY') - N=1 - X=( CAR2(9)+CAR2(1)*TAN(CAR2(2)*DEGRAD) ) - 2.*C2N1(9) - - + (C2N2(9)+C2N2(1)*TAN(C2N2(2)*DEGRAD)) + Z=CAR2(1)-C1N1(1) + X=CAR2(9)+CAR2(1)*TAN(ANG2C)- + + (C1N1(9)-C1N1(1)*TAN(C1N1(2)*DEGRAD)) + X=-X !ALEX MOD + CALL GSPOS('C1N1',N,'CAR2',X,Y,Z,0,'ONLY') Y=0 - Z=CAR2(1)-C2N2(1) - CALL GSPOS('C2N2',N,'CAR2',X,Y,Z,0,'ONLY') + Z=CAR2(1)-C1N2(1) + X=CAR2(9)+CAR2(1)*TAN(ANG2C)- + + (2.*C1N1(9)+C1N2(9)-C1N2(1)*TAN(C1N2(2)*DEGRAD)) + X=-X !ALEX MOD + CALL GSPOS('C1N2',N,'CAR2',X,Y,Z,0,'ONLY') N=1 - X=0 + X=PT*(1./COS(ATAN(D2/CAR(1)))-1./SIN(DEGRAD*CAR(7)))/2. Y=0 Z=0 CALL GSPOS('C2D1',N,'CR2P',X,Y,Z,0,'ONLY') * * Positioning in card * - +* +* CAR1A +* CAR2B CAR2A +* CAR1B +* + ZH=49.229-ZTPLA-TPLA(3) +*--- CAR1B N=1 - X=-(C2N1(9) + C2N2(9)) + CAR2(1)*ABS(TAN(CAR2(2)*DEGRAD)) - Y=(CAR1(9) - (C1N1(9)+C1N2(9))) - + + CAR2(1)*ABS(COS(ANGLCAR2*RADDEG)) - 2.*CAR2(4) - Z=0 - CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR22,'ONLY') + Y=-CARB(5)+CAR(8)/2.*SIN(ANGLCAR1) + Z=CARB(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR1)/2. + X=((CARB(7)/2.)-TAN(ANG1C)*CAR(8)/2.)-CARB(4) + CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR1B,'ONLY') +*--- CAR1A + N=2 + X=-X + Y=-Y + CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR1A,'ONLY') +*--- CAR2A + N=1 + X=CARA(4)-CAR(8)/2.*SIN(ANGLCAR2) + Z=CARA(6)+ZH-CARDB(3)-CAR(8)*COS(ANGLCAR2)/2. + Y=(CARA(7)/2.-TAN(ANG2C)*CAR(8)/2.)-CARA(5) + CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR2A,'ONLY') +*--- CAR2B + N=2 + X=-X + Y=-Y + CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR2B ,'ONLY') +c--------------------------------------------------------------- +* +* LEGS +* + NMED=MAL + CALL GSVOLU('LEGB','BOX ',NMED,LEGB, 3,IVOLU) +* +* empty volumes to be subtracted from LEGB + NMED=MN2 +* + CALL GSVOLU('LEG1','TRAP',NMED,LEG1, 11,IVOLU) + X=-LEGB(1)+LEG1(1) + Y=-LEGB(2)+LEG1(4) + Z=0. + CALL GSPOS('LEG1',1,'LEGB',X,Y,Z,RY90,'ONLY') + Y=-LEGB(2)+2.*1.95-LEG1(4) + CALL GSPOS('LEG1',2,'LEGB',X,Y,Z,RY90,'ONLY') +* + CALL GSVOLU('LEG2','BOX ',NMED,LEG2, 3,IVOLU) + X=-LEGB(1)+LEG2(1) + Y=-LEGB(2)+1.95 + Z=0. + CALL GSPOS('LEG2',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG3','BOX ',NMED,LEG3, 3,IVOLU) + X=-LEGB(1)+0.65*2.+LEG3(1)+2./10 + Y=LEGB(2)-LEG3(2) + Z=LEGB(3)-LEG3(3) + CALL GSPOS('LEG3',1,'LEGB',X,Y,Z,0,'ONLY') + Z=-LEGB(3)+LEG3(3) + CALL GSPOS('LEG3',2,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG4','BOX ',NMED,LEG4, 3,IVOLU) + X=-LEGB(1)+LEG4(1) + Y=LEGB(2)-LEG4(2) + Z=0. + CALL GSPOS('LEG4',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG5','BOX ',NMED,LEG5, 3,IVOLU) + X=LEGB(1)-LEG5(1) + Y=LEGB(2)-LEG5(2) + Z=0. + CALL GSPOS('LEG5',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG6','BOX ',NMED,LEG6, 3,IVOLU) + X=-LEGB(1)+0.65*2.+LEG3(1) + Y=LEGB(2)-LEG6(2) + Z=0. + CALL GSPOS('LEG6',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG7','BOX ',NMED,LEG7, 3,IVOLU) + X=LEGB(1)-1.125-0.375 + Y=-LEGB(2)+LEG7(2) + Z=11.85-.6-LEG7(3) + CALL GSPOS('LEG7',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG8','BOX ',NMED,LEG8, 3,IVOLU) + Z=0. + CALL GSPOS('LEG8',1,'LEGB',X,Y,Z,0,'ONLY') +* + CALL GSVOLU('LEG9','BOX ',NMED,LEG9, 3,IVOLU) + Z=-(11.85-.6-LEG7(3)) + CALL GSPOS('LEG9',1,'LEGB',X,Y,Z,0,'ONLY') +* +* Positioning LEG BOXES in CARD +* + N=1 + X=-CARDB(1)+LEGB(1)+3.8 + Y=-CARDB(2)+LEGB(2)+0.7 + Z=-0.5 + CALL GSPOS('LEGB',N,'CARD',X,Y,Z,0,'MANY') N=2 - X=-X + 2.*CAR1(4) - Y=-Y + 2.*CAR1(4) - CALL GSPOS('CAR2',N,'CARD',X,Y,Z,IRZCAR2 ,'ONLY') + X=X + Y=-Y + CALL GSPOS('LEGB',N,'CARD',X,Y,Z,RX180,'MANY') + N=3 + X=-X + Y=Y + CALL GSROTM(RZALEX180,90.,180.,90., + + 270.,180.,0.) + CALL GSPOS('LEGB',N,'CARD',X,Y,Z,RZ180,'MANY') + N=4 + X=X + Y=-Y + CALL GSPOS('LEGB',N,'CARD',X,Y,Z,RY180,'MANY') +* +c--------------------------------------------------------------- +* +* CARD TOP PLATE +* +* TPLT +* + NMED=MAL + CALL GSVOLU('TPLT','BOX ',NMED,TPLT, 3,IVOLU) +* +* Empty volumes to be subtracted from TPLT +* + NMED=MN2 + CALL GSVOLU('TOPC','BOX ',NMED,TOPC, 3,IVOLU) +* + CALL GSVOLU('TH11','BOX ',NMED,TH11, 3,IVOLU) + CALL GSVOLU('TH12','BOX ',NMED,TH12, 3,IVOLU) + CALL GSVOLU('TH21','BOX ',NMED,TH21, 3,IVOLU) + CALL GSVOLU('TH22','BOX ',NMED,TH22, 3,IVOLU) +* +* XY plane + CALL GSVOLU('TOPT','TRAP',NMED,TOPT,11,IVOLU) +* + CALL GSVOLU('TOHO','TUBE',NMED,TOHO,3,IVOLU) +* +* Positioning empty volumes... +* +* Positioning TOPC in TPLT +* N=1 - Y=(C1N1(9) + C1N2(9)) - CAR1(1)*ABS(TAN(CAR1(2)*DEGRAD)) - + + 2.* CAR2(4) - X=(CAR2(9) - (C2N1(9)+C2N2(9))) - + + CAR1(1)*ABS(SIN(ANGLCAR1*RADDEG)) + 2.*CAR1(4) - Z=0 - CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR12,'ONLY') + X=0. + Y=0. + Z=0. + CALL GSPOS('TOPC',N,'TPLT',X,Y,Z,0,'ONLY') +* +* +* Positioning H11 in TPLT +* + N=1 + X=-23.6 + Y=14.4 + Z=0. + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + N=2 + X=-23.6 + Y=-14.4 + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + N=3 + X=23.6 + Y=14.4 + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + N=4 + X=23.6 + Y=-14.4 + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') +* +* Positioning H12 in TPLT +* + N=1 + X=-23.6 + Y=5.5 + Z=0. + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + N=2 + X=-23.6 + Y=-5.5 + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + N=3 + X=23.6 + Y=5.5 + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + N=4 + X=23.6 + Y=-5.5 + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') +* +* Positioning H21 in TPLT +* + N=1 + X=-17.2 + Y=19.4 + Z=0. + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + N=2 + X=-17.2 + Y=-19.4 + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + N=3 + X=17.2 + Y=19.4 + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + N=4 + X=17.2 + Y=-19.4 + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') +* +* Positioning H22 in TPLT +* + N=1 + X=-10.2 + Y=19.4 + Z=0. + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + N=2 + X=-10.2 + Y=-19.4 + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + N=3 + X=10.2 + Y=19.4 + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + N=4 + X=10.2 + Y=-19.4 + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + N=5 + X=0. + Y=19.4 + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + N=6 + X=0. + Y=-19.4 + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') +* +* Positioning TOPT in TPLT +* + N=1 + X=-CARDB(1)+4.5/4. + Y=-CARDB(2)+4.5/2. + Z=0. + CALL GSPOS('TOPT',N,'TPLT',X,Y,Z,0,'ONLY') + N=2 + X=X + Y=-Y + CALL GSPOS('TOPT',N,'TPLT',X,Y,Z,RX180,'ONLY') + N=3 + X=-X + Y=Y + CALL GSPOS('TOPT',N,'TPLT',X,Y,Z,RZ180,'ONLY') + N=4 + X=X + Y=-Y + CALL GSPOS('TOPT',N,'TPLT',X,Y,Z,RY180,'ONLY') +* +* Positioning TOHO in TPLT +* + N=1 + X=TPLT(1)-6.2-TOHO(1) + Y=TPLT(2)-4.1-TOHO(2) + CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY') N=2 - Y=-Y + 2.* CAR2(4) - X=-X + 2.* CAR2(4) - CALL GSPOS('CAR1',N,'CARD',X,Y,Z,IRZCAR1,'ONLY') + X=-X + Y=Y + CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY') + N=3 + X=X + Y=-Y + CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY') + N=4 + X=-X + Y=Y + CALL GSPOS('TOHO',N,'TPLT',X,Y,Z,0,'ONLY') +* +* Positioning TPLT in CARD volume +* + N=1 + X=0. + Y=0. + Z=CARDB(3)-TPLT(3) + CALL GSPOS('TPLT',N,'CARD',X,Y,Z,0,'ONLY') +C--------------------------------------------------------------- +* +* MAGNETIC SCREEN +* + NMED=10 + CALL GSVOLU('MGSC','BOX ',NMED,MGSC, 3,IVOLU) +* Empty volume to be subtracted from MGSC + NMED=MN2 + CALL GSVOLU('MGSH','BOX ',NMED,MGSH, 3,IVOLU) +* +* Positioning in CARD... +* + N=1 + X=0. + Y=0. + Z=0. + CALL GSPOS('MGSH',N,'MGSC',X,Y,Z,0,'ONLY') + Z=-CARDB(3)+MGSC(3) + CALL GSPOS('MGSC',N,'CARD',X,Y,Z,0,'ONLY') +c--------------------------------------------------------------- RETURN END