real function qsumall(id)
      integer id
      include 'nt.inc'
      logical hexist
      external hexist
      integer i
      character*80 title
      integer nbins
      real xmax,rmax
      data nbins,xmax/400,5./
      rmax=xmax*sqrt(2.)
      if (.not.hexist(id)) then
        do i=1,8
          write(title,*) 'qsum for sector',i
          call hbook1(id+i-1,title,nbins,-rmax,rmax,0)
          write(title,*) 'qsumb for sector',i
          call hbook1(id+i+9,title,nbins,-rmax,rmax,0)
          write(title,*) 'qsump for sector',i
          call hbook1(id+i+19,title,nbins,-rmax,rmax,0)
        enddo
      endif
      xr=(+offset(1)+offset(2))/sqrt(2.)
      yr=(-offset(1)+offset(2))/sqrt(2.)
      anorm=(2*rmax/nbins)*0.1/100
      pA=1.6e-7/2.3
      if (abs(yr).lt.0.05) then
        do i=1,8
          call hfill(id+i-1,xr,yr,qsum(i)*pA/anorm)
          call hfill(id+i+9,xr,yr,qsumb(i)*pA/anorm)
          call hfill(id+i+19,xr,yr,qsump(i)*pA/anorm)
        enddo
      endif
      qsumall=qsum(1)+qsum(2)+qsum(3)+qsum(4)
     +       +qsum(5)+qsum(6)+qsum(7)+qsum(8)
      end