SUBROUTINE GPCARDV
************************************************************************
*                                                                      *
* Volumes definition for the side anticoincidence counters             *
* Called by: GPGEO                                                     *
*                                                                      *
* 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,I
      REAL X,Y,Z
      REAL CAR1(11),CR1P(11)
      REAL C1D1(11),C2D1(11)
      REAL CAR2(11),CR2P(11)
      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
      
      CAR1(1) = C1D1(1)+ATZ
      CAR1(2) = C1D1(2)
      CAR1(3) = 0.
      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) = 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) = C1D1(1)+PT
      CR1P(2) = C1D1(2)
      CR1P(3) = 0.
      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) = 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)
*
      NMED=MSCIN
*
      CALL GSVOLU('C1D1','TRAP',NMED,C1D1,11,IVOLU)
*
* Empty volume to be subtracted from C1D1 and C2D1 sheets
*
      NMED=MN2
      
      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) = CAR1(4)
      C1N1(5) = 0.
      C1N1(6) = C1N1(5)
      C1N1(7) = 0.
      C1N1(8) = CAR1(8)
      C1N1(9) = 0.5*BB
      C1N1(10) = C1N1(9)
      C1N1(11) = 0.
      
      CALL GSVOLU('C1N1','TRAP',NMED,C1N1,11,IVOLU)

      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) = CAR1(4)
      C1N2(5) = 0.
      C1N2(6) = C1N2(5)
      C1N2(7) = 0.
      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
      
      CAR2(1) = C2D1(1)+ATZ
      CAR2(2) = C2D1(2)
      CAR2(3) = 0.
      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) = 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) = C2D1(1)+PT
      CR2P(2) = C2D1(2)
      CR2P(3) = 0.
      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) = 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)
*
      NMED=MSCIN
*     
      CALL GSVOLU('C2D1','TRAP',NMED,C2D1,11,IVOLU)
*
*
* The zenihtal angle of the paddles
*
      ANGLCAR2=CARA(3)*DEGRAD   !CAR2 inclination
      ANGLCAR1=CARB(3)*DEGRAD   !CAR1 inclination

      NMED=MN2
      CALL GSVOLU('CARD','BOX ',NMED,CARDB, 3,IVOLU)
*
*  Define ANGLCAR rotation around Z axis
*
*    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=(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
      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
      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=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')
*
* Positioning in CAR2
*
      N=1
      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=2
      Y=0
      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)-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=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
      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
      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
      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
      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