REAL FUNCTION hitslocal()
*********************************************************
*                                                       *
* This file was generated by HUWFUN.                    *
*                                                       *
*********************************************************
*
*     Ntuple Id:      3    
*     Ntuple Title:   blhits
*     Creation:       27/01/:0 00.11.32
*
*********************************************************
*
      LOGICAL         CHAIN
      CHARACTER*128   CFILE
      INTEGER         IDNEVT,NCHEVT,ICHEVT
      REAL            OBS(13)
*
      COMMON /PAWIDN/ IDNEVT,OBS
      COMMON /PAWCHN/ CHAIN, NCHEVT, ICHEVT
      COMMON /PAWCHC/ CFILE
*
*--   Ntuple Variable Declarations
*
      REAL barrel(2500),clusid(2500),detid(2500),hitphi(2500)
     + ,layer(2500),nmchits(2500),nmcprim(2500),nmctrks(2500)
     + ,nstrip(2500),pulse(2500),view(2500),xdiff(2500),xerrl(2500)
     + ,xhmc(2500),xlocal(2500),yhmc(2500),zhmc(2500)
      INTEGER nhits
*
      COMMON /PAWCR4/ nhits,barrel,clusid,detid,hitphi,layer,nmchits
     + ,nmcprim,nmctrks,nstrip,pulse,view,xdiff,xerrl,xhmc,xlocal,yhmc
     + ,zhmc
*

*
*--   Enter user code here
*
      logical first
      data first/.true./
      save first
*
*--   Enter user code here
*
      if (first) then
         first=.false.
         call hbook1(139,' Residuals ',100,-0.01,0.01,0.)
         call hbook1(140,' Residuals - match with a primary',100,
     &        -0.01,0.01,0.)
         call hbook1(150,' Residuals - Nprimary.ge.5'
     &        ,100,-0.01,0.01,0.)

         do ilay=1,4
         call hbook1(100+ilay,' Layer  Nstrips ',50,0.,50.,0.)
         call hbook2(110+ilay,' Layer  Nstrips vs NMc Hits',
     &       50,0.,50.,25,0.,25.,0.)
         call hbook2(120+ilay,' Layer  Nstrips vs NMc Primary Hits',
     &       50,0.,50.,25,0.,25.,0.)
         call hbook2(130+ilay,' Layer  Nstrips vs NMc Tracks',
     &       50,0.,50.,25,0.,25.,0.)
         call hbook1(140+ilay,' Layer Residuals - primary',
     &        100,-0.01,0.01,0.)
         call hbook1(150+ilay,'  Residuals - Nprimary=#',
     &        100,-0.01,0.01,0.)
         end do
      endif
      do i=1,8
         ilay = (i+1)/2
c         print*, i,ilay
      end do
      do  i=1,nhits
         ilay=(layer(i)+1)/2 
         if (view(i).eq.1) then
            call hfill(100+ilay, nstrip(i),0.,1.)
            call hfill(110+ilay, nstrip(i),nmchits(i),1.)
            call hfill(120+ilay, nstrip(i),nmcprim(i),1.)
            call hfill(130+ilay, nstrip(i),nmctrks(i),1.)
            call hfill(139, xdiff(i),0.,1.)
            if (nmcprim(i).ne.0) then
               call hfill(140, xdiff(i),0.,1.)
               call hfill(140+ilay, xdiff(i),0.,1.)
*               print*, nmcprim(i)
               if (nmcprim(i).le.4) then
                  nnn=nmcprim(i)
                  call hfill(150+nnn, xdiff(i),0.,1.)
               else
                  call hfill(150, xdiff(i),0.,1.)
               endif
            endif
         else 
            print*, ' wrong view ', view(i)
         endif
      enddo
      
      hitslocal = 1. 
*
      END