PROGRAM ICD
**************************************************
* This calculates parameters for ICD geometry for
* d0gstar. 
*
* The geometry tree is
*
*   EC+ -- ICDN -- DICN(16) -- ICD1     (ieta 12)
*                              ICD2     (ieta 13)
*                              ICD3     (ieta 14)
*
*   EC- -- ICDS -- DICN(16) -- ICD1     (ieta 12)
*                              ICD2     (ieta 13)
*                              ICD3     (ieta 14)
*
* Shape of volumes are PGON.
*
* A. Stone & S.Kunori, July,2000.
**************************************************

*      IMPLICIT NONE

      DEG70=3.14159*70.0/180.
      S70=SIN(DEG70)
      C70=COS(DEG70)
      
*     -- center of ICD....
      R0=104.0
      Z0=162.0

*********************************************
*  container box.
*********************************************

*     -- half thickness of box...
      boxXW=2.2/2.0
*     -- half length of box...
      boxXL=50.0/2.0

*     -- coordinate at the center of two end planes....
      boxZA=Z0-boxXL*C70
      boxRA=R0-boxXL*S70

      boxZB=Z0+boxXL*C70
      boxRB=R0+boxXL*S70

*     -- displacement...
      boxDZ=boxXW*S70
      boxDR=boxXW*C70
      boxDD=2.0*boxXW/C70

*     -- tile boundary A.
      boxZ1=boxZA-boxDZ
      boxRI1=boxRA+boxDR
      boxRO1=boxRI1

      boxZ2=boxZA+boxDZ
      boxRI2=boxRA-boxDR
      boxRO2=boxRI2+boxDD

*     -- tile boundar D.
      boxZ3=boxZB-boxDZ
      boxRO3=boxRB+boxDR
      boxRI3=boxRO3-boxDD

      boxZ4=boxZB+boxDZ
      boxRI4=boxRB-boxDR
      boxRO4=boxRI4
*     -- now write...
      WRITE(6,*) 'ICD '
      WRITE(6,100) boxZ1,boxRI1,boxRO1, boxZ2,boxRI2,boxRO2, 
     &             boxZ3,boxRI3,boxRO3, boxZ4,boxRI4,boxRI4

*********************************************
*  scintillator volumes
*********************************************

*     -- half thickness of scintillator.  0.5 inches
      XW=(2.54*0.5)/2.0
*     -- half length of scintillator...
      XL=48.6/2.0

*     -- distance from lower edge of scintillator...
*        calculate tile boundaies using same fraction shown
*        in ICD TDR.
      A=5.031+6.6063+7.531
      XL1=2.0*XL*5.031/A
      XL2=2.0*XL*(5.031+6.6063)/A
      XL3=2.0*XL

*     -- coordinate at the center of boundary plane...
      ZA=Z0-XL*C70
      RA=R0-XL*S70

      ZB=ZA+XL1*C70
      RB=RA+XL1*S70

      ZC=ZA+XL2*C70
      RC=RA+XL2*S70

      ZD=ZA+XL3*C70
      RD=RA+XL3*S70

*     -- displacement...
      DZ=XW*S70
      DR=XW*C70
      DD=2.0*XW/C70

*     -- tile boundary A.
      Z1=ZA-DZ
      RI1=RA+DR
      RO1=RI1

      Z2=ZA+DZ
      RI2=RA-DR
      RO2=RI2+DD
*     -- tile boundary B.
      Z3=ZB-DZ
      RO3=RB+DR
      RI3=RO3-DD

      Z4=ZB+DZ
      RI4=RB-DR
      RO4=RI4+DD
*     -- tile boundary C.
      Z5=ZC-DZ
      RO5=RC+DR
      RI5=RO5-DD

      Z6=ZC+DZ
      RI6=RC-DR
      RO6=RI6+DD
*     -- tile boundar D.
      Z7=ZD-DZ
      RO7=RD+DR
      RI7=RO7-DD

      Z8=ZD+DZ
      RI8=RD-DR
      RO8=RI8
*     -- now write...
      WRITE(6,*) 'ICD1'
      WRITE(6,100) Z1,RI1,RO1, Z2,RI2,RO2, Z3,RI3,RO3, Z4,RI4,RI4

      WRITE(6,*) 'ICD2'
      WRITE(6,100) Z3,RO3,RO3, Z4,RI4,RO4, Z5,RI5,RO5, Z6,RI6,RI6

      WRITE(6,*) 'ICD3'
      WRITE(6,100) Z5,RO5,RO5, Z6,RI6,RO6, Z7,RI7,RO7, Z8,RI8,RO8

 100  FORMAT(20X,3F8.2,2X,3F8.2)

      END