--- gpamela/gpcard/gpcardv.F 2006/12/17 14:53:05 1.3 +++ gpamela/gpcard/gpcardv.F 2008/01/29 18:25:16 1.4 @@ -24,14 +24,14 @@ H2=CAR(1)-H1 * * Greater Scintillator CAR1 (or type B) -* +* B1=CARB(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7)) C1=H2/TAN(DEGRAD*CAR(7)) A1=CARB(1)-C1 D1=CARB(1)-CARB(2)-CAR(4)-CAR(5) F1=C1+H1/TAN(DEGRAD*CAR(7)) ANG1C=ATAN((B1/2.+D1-A1/2.)/CAR(1)) - + C1D1(1) = CAR(1)/2. C1D1(2) = -ANG1C*RADDEG C1D1(3) = 0. @@ -46,9 +46,9 @@ NMED=MSCIN CALL GSVOLU('C1D1','TRAP',NMED,C1D1,11,IVOLU) -* +* * Smaller Scintillator CAR2 (or type A) -* +* B2=CARA(2)+CAR(4)+CAR(5)+H1/TAN(DEGRAD*CAR(7)) C2=H2/TAN(DEGRAD*CAR(7)) A2=CARA(1)-C2 @@ -95,7 +95,7 @@ NMED=MAL CALL GSVOLU('CAR1','TRAP',NMED,CAR1,11,IVOLU) * -* CR1P PLASTIC BOX TO CREATE RIM +* CR1P PLASTIC BOX TO CREATE RIM * BB=C1D1(9)+PT*(1./COS(ATAN(D1/CAR(1)))- + D1/CAR(1)+1./TAN(DEGRAD*CAR(7))+ @@ -139,11 +139,11 @@ CAR2(9) = BB CAR2(10) = CAR2(9) CAR2(11) = 0. - + NMED=MAL CALL GSVOLU('CAR2','TRAP',NMED,CAR2,11,IVOLU) * -* CR1P PLASTIC BOX TO CREATE RIM +* CR1P PLASTIC BOX TO CREATE RIM * BB=C2D1(9)+PT*(1./COS(ATAN(D2/CAR(1)))- + D2/CAR(1)+1./TAN(DEGRAD*CAR(7))+ @@ -164,14 +164,14 @@ CR2P(10) = CR2P(9) CR2P(11) = 0. - NMED=MPLAS + NMED=MPLAS CALL GSVOLU('CR2P','TRAP',NMED,CR2P,11,IVOLU) * * Volumes to be subtracted from C1D1 and C2D1 * AA=CAR(2)+CAR(3) BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2) - + VSN1(1) = 0.5*AA VSN1(2) = RADDEG*ATAN((BB/2.-AA/TAN(DEGRAD*CAR(7)))/AA) VSN1(3) = 0. @@ -184,9 +184,9 @@ VSN1(10) = VSN1(9) VSN1(11) = 0. - AA=CAR(3) + AA=CAR(3) BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2)) - + VSN2(1) = 0.5*AA VSN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA) VSN2(3) = 0. @@ -209,7 +209,7 @@ AA=(PT/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+ + CAR(4)/CAR(2))+CAR(2)+CAR(3)+PT BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2) - + VPN1(1) = 0.5*AA VPN1(2) = VSN1(2) VPN1(3) = 0. @@ -221,12 +221,12 @@ VPN1(9) = 0.5*BB VPN1(10) = VPN1(9) VPN1(11) = 0. - - AA=CAR(3) + PT - + + AA=CAR(3) + PT - + COS(ATAN(CAR(4)/CAR(2)))*PT/SIN(ATAN(CAR(5)/CAR(3))- + ATAN(CAR(4)/CAR(2))) BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2)) - + VPN2(1) = 0.5*AA VPN2(2) = RADDEG*ATAN((BB/2.+AA*CAR(4)/CAR(2))/AA) VPN2(3) = 0. @@ -249,7 +249,7 @@ AA=(ATZ/SIN(DEGRAD*CAR(7)))/(1./TAN(DEGRAD*CAR(7))+ + CAR(4)/CAR(2))+CAR(2)+CAR(3)+ATZ BB=AA/TAN(DEGRAD*CAR(7))+AA*CAR(4)/CAR(2) - + VAN1(1) = 0.5*AA VAN1(2) = VSN1(2) VAN1(3) = 0. @@ -262,11 +262,11 @@ VAN1(10) = VAN1(9) VAN1(11) = 0. - AA=CAR(3) + ATZ - + AA=CAR(3) + ATZ - + COS(ATAN(CAR(4)/CAR(2)))*ATZ/SIN(ATAN(CAR(5)/CAR(3))- + ATAN(CAR(4)/CAR(2))) BB=AA*(CAR(5)/CAR(3)-CAR(4)/CAR(2)) - + VAN2(1) = 0.5*AA VAN2(2) = RADDEG*ATAN((BB/2.+(AA)*CAR(4)/CAR(2))/AA) VAN2(3) = 0. @@ -279,7 +279,7 @@ VAN2(10) = VAN2(9) VAN2(11) = 0. - NMED=MN2 + NMED=MN2 CALL GSVOLU('VAN1','TRAP',NMED,VAN1,11,IVOLU) CALL GSVOLU('VAN2','TRAP',NMED,VAN2,11,IVOLU) * @@ -354,7 +354,7 @@ X=-(CAR1(9)-CAR1(1)*TAN(CAR1(2)*DEGRAD)- + (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD))) CALL GSPOS('VAN2',N,'CAR1',X,Y,Z,0,'ONLY') -* +* N=1 X=PT*(1./COS(ATAN(D1/CAR(1)))- + (1./SIN(DEGRAD*CAR(7))))/2. @@ -412,7 +412,7 @@ X=-(CAR2(9)-CAR2(1)*TAN(CAR2(2)*DEGRAD)- + (2.*VAN1(9)+VAN2(9)-VAN2(1)*TAN(VAN2(2)*DEGRAD))) CALL GSPOS('VAN2',N,'CAR2',X,Y,Z,0,'ONLY') -* +* N=1 X=PT*(1./COS(ATAN(D2/CAR(1)))-1./SIN(DEGRAD*CAR(7)))/2. Y=0 @@ -434,13 +434,13 @@ * CAR1B * ZH=49.229-ZTPLA-TPLA(3) -*--- CAR1B +*--- CAR1B N=1 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 +*--- CAR1A N=2 X=-X Y=-Y @@ -451,7 +451,7 @@ 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 +*--- CAR2B N=2 X=-X Y=-Y @@ -463,62 +463,62 @@ NMED=MAL CALL GSVOLU('LEGB','BOX ',NMED,LEGB, 3,IVOLU) * -* empty volumes to be subtracted from LEGB +* empty volumes to be subtracted from LEGB * NMED=MN2 - CALL GSVOLU('LEG1','TRAP',NMED,LEG1, 11,IVOLU) + 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') + CALL GSPOS('LEG1',1,'LEGB',X,Y,Z,RY90,'ONLY') Y=-LEGB(2)+2.*LEGP(1)-LEG1(4) - CALL GSPOS('LEG1',2,'LEGB',X,Y,Z,RY90,'ONLY') -* - CALL GSVOLU('LEG2','BOX ',NMED,LEG2, 3,IVOLU) + 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)+LEGP(1) Z=0. CALL GSPOS('LEG2',1,'LEGB',X,Y,Z,0,'ONLY') * - CALL GSVOLU('LEG3','BOX ',NMED,LEG3, 3,IVOLU) + CALL GSVOLU('LEG3','BOX ',NMED,LEG3, 3,IVOLU) X=-LEGB(1)+LEGP(2)*2.+LEG3(1)+0.2 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 GSPOS('LEG3',2,'LEGB',X,Y,Z,0,'ONLY') * - CALL GSVOLU('LEG4','BOX ',NMED,LEG4, 3,IVOLU) + 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) +* + 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) + CALL GSVOLU('LEG6','BOX ',NMED,LEG6, 3,IVOLU) X=-LEGB(1)+LEGP(2)*2.+LEG3(1) Y=LEGB(2)-LEG6(2) Z=0. - CALL GSPOS('LEG6',1,'LEGB',X,Y,Z,0,'ONLY') + CALL GSPOS('LEG6',1,'LEGB',X,Y,Z,0,'ONLY') * - CALL GSVOLU('LEG7','BOX ',NMED,LEG7, 3,IVOLU) + CALL GSVOLU('LEG7','BOX ',NMED,LEG7, 3,IVOLU) X=LEGB(1)-LEGP(3)-LEGP(4) Y=-LEGB(2)+LEG7(2) Z=LEGP(5)-LEGP(8)-LEG7(3) - CALL GSPOS('LEG7',1,'LEGB',X,Y,Z,0,'ONLY') + CALL GSPOS('LEG7',1,'LEGB',X,Y,Z,0,'ONLY') * - CALL GSVOLU('LEG8','BOX ',NMED,LEG8, 3,IVOLU) + CALL GSVOLU('LEG8','BOX ',NMED,LEG8, 3,IVOLU) Z=0. - CALL GSPOS('LEG8',1,'LEGB',X,Y,Z,0,'ONLY') + CALL GSPOS('LEG8',1,'LEGB',X,Y,Z,0,'ONLY') * - CALL GSVOLU('LEG9','BOX ',NMED,LEG9, 3,IVOLU) + CALL GSVOLU('LEG9','BOX ',NMED,LEG9, 3,IVOLU) Z=-(LEGP(5)-LEGP(8)-LEG7(3)) - CALL GSPOS('LEG9',1,'LEGB',X,Y,Z,0,'ONLY') + CALL GSPOS('LEG9',1,'LEGB',X,Y,Z,0,'ONLY') * * Positioning LEG BOXES in CARD * @@ -558,8 +558,8 @@ * 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) + CALL GSVOLU('TH21','BOX ',NMED,TH21, 3,IVOLU) + CALL GSVOLU('TH22','BOX ',NMED,TH22, 3,IVOLU) * CALL GSVOLU('TOPT','TRAP',NMED,TOPT,11,IVOLU) * @@ -573,8 +573,8 @@ N=1 X=0. Y=0. - Z=0. - CALL GSPOS('TOPC',N,'TPLT',X,Y,Z,0,'ONLY') + Z=0. + CALL GSPOS('TOPC',N,'TPLT',X,Y,Z,0,'ONLY') * * * Positioning H11 in TPLT @@ -583,19 +583,19 @@ X=-TOPP(1) Y=TOPP(2) Z=0. - CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') N=2 X=-TOPP(1) Y=-TOPP(2) - CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') N=3 X=TOPP(1) Y=TOPP(2) - CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') N=4 X=TOPP(1) Y=-TOPP(2) - CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH11',N,'TPLT',X,Y,Z,0,'ONLY') * * Positioning H12 in TPLT * @@ -603,19 +603,19 @@ X=-TOPP(1) Y=TOPP(3) Z=0. - CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') N=2 X=-TOPP(1) Y=-TOPP(3) - CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') N=3 X=TOPP(1) Y=TOPP(3) - CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') N=4 X=TOPP(1) Y=-TOPP(3) - CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH12',N,'TPLT',X,Y,Z,0,'ONLY') * * Positioning H21 in TPLT * @@ -623,19 +623,19 @@ X=-TOPP(4) Y=TOPP(5) Z=0. - CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') N=2 X=-TOPP(4) Y=-TOPP(5) - CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') N=3 X=TOPP(4) Y=TOPP(5) - CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') N=4 X=TOPP(4) Y=-TOPP(5) - CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH21',N,'TPLT',X,Y,Z,0,'ONLY') * * Positioning H22 in TPLT * @@ -643,27 +643,27 @@ X=-TOPP(6) Y=TOPP(5) Z=0. - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') N=2 X=-TOPP(6) Y=-TOPP(5) - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') N=3 X=TOPP(6) Y=TOPP(5) - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') N=4 X=TOPP(6) Y=-TOPP(5) - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') N=5 X=0. Y=TOPP(5) - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') N=6 X=0. Y=-TOPP(5) - CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') + CALL GSPOS('TH22',N,'TPLT',X,Y,Z,0,'ONLY') * * Positioning TOPT in TPLT * @@ -705,34 +705,34 @@ 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--------------------------------------------------------------- + CALL GSPOS('TPLT',N,'CARD',X,Y,Z,0,'ONLY') +C--------------------------------------------------------------- * * MAGNETIC SCREEN * - NMED=10 - CALL GSVOLU('MGSC','BOX ',NMED,MGSC, 3,IVOLU) + NMED=MIRON + 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) + 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--------------------------------------------------------------- + CALL GSPOS('MGSC',N,'CARD',X,Y,Z,0,'ONLY') +c--------------------------------------------------------------- RETURN END