real function fill_ptot()
      include 'pass2nt.inc'
      if(abs(cos3d).lt.0.5.and.(tpi-tk).gt.2)then
         call hfill(1000,ptot,0.,1.0)
      end if
      return
      end