subroutine digi_cer(iclust)
      implicit none
      include 'beta_geom.inc'
      include 'sane_cwn.inc'
      integer iclust
      integer i
      real predicted_x,predicted_y
      integer mirror(8)
      real CER_SANE_GEOM_CUT_LOW(8) 
      real CER_SANE_GEOM_CUT_HI(8) 
      integer CER_SANE_GEOM_CUT_X(8)
      data CER_SANE_GEOM_CUT_LOW /0,   0,  10,  10,  
     ,     22, 22,  35,  35/     
      data CER_SANE_GEOM_CUT_HI /22,  22,  42, 42, 
     ,     54,  54,  56, 56/
      data CER_SANE_GEOM_CUT_X /1,  -1,  1, -1, 
     ,	1,  -1,  1, -1/
      integer cer_n
      mirror = 0
      cer_h(iclust)=0
      predicted_x = cwn_x_ur(iclust)
      predicted_y = cwn_y_ur(iclust)
      do i=1,cwn_cergood
         do cer_n=1,8
         if(predicted_y.gt.(CER_SANE_GEOM_CUT_LOW(cer_n)*4-1)-120..and.
     ,        predicted_y.lt.(CER_SANE_GEOM_CUT_HI(cer_n)*4+1)-120.and.
     ,              predicted_x*CER_SANE_GEOM_CUT_X(cer_n).gt.-20)then
            if(cwn_cg_y(i).gt.
     ,           ((CER_SANE_GEOM_CUT_LOW(cer_n)*4-1)-120)*205/335..and.
     ,           cwn_cg_y(i).lt.
     ,           ((CER_SANE_GEOM_CUT_HI(cer_n)*4-1)-120)*205/335..and.
     , cwn_cg_x(i)*CER_SANE_GEOM_CUT_X(cer_n).gt.-20*205/335.
     ,       )then
c            write(*,*)predicted_y,predicted_x,
c     ,           cwn_cg_y(i),cwn_cg_x(i),
c     ,           cwn_cergood,i,cer_n
               cer_h(iclust)=cer_h(iclust)+1
            endif
         endif
      enddo
      enddo
      

      

      end
      subroutine clear_cer()
      include 'sane_cwn.inc'
      cwn_cergood=0

      do i=1,1000
         cwn_cg_y(i)=0
         cwn_cg_x(i)=0
         cwn_cg_xp(i)=0
         cwn_cg_yp(i)=0
      enddo
      end