Subroutine book_nu(infile1,infile2,infile3,infile4,infile5, 
     c                  outfile)

*     ntuple ID=200
*
*     Usage: PAW> for/call book_nu.for('data-1.txt',..., 'data.hbook')
*

      integer count, nf
      parameter (nf=3)
      character*28 infile1,infile2,infile3,infile4,infile5,outfile
      integer ievent,nch,trueid,ndirb,nfit
      real zenith,azimuth,delang,truezen,trueazi,trueen,
     c     smootallphit,ldirc,jkprob,jkchi,jkrchi,
     c     r14,intwght,atmswght

      common/ntupr/ievent,nch,zenith(nf),azimuth(nf),delang(nf),
     c     truezen(nf),trueazi(nf),trueen(nf),trueid(nf),
     c     ndirb(nf),smootallphit(nf),ldirc(nf),jkprob(nf),
     c     nfit,jkchi(nf),jkrchi(nf),
     c     r14,intwght,atmswght


      do i=1,nf
         zenith(i)=0.0
         azimuth(i)=0.0
         delang(i)=0.0
         truezen(i)=0.0
         trueazi(i)=0.0
         trueen(i)=0.0
         trueid(i)=0.0
         ndirb(i)=0
         smootallphit(i)=0.0
         ldirc(i)=0.0
         jkprob(i)=0.0
         jkchi(i)=0.0
         jkrchi(i)=0.0
      enddo
      
      lin1=40
      lin2=41
      lin3=42
      lin4=43
      lin5=44
      lout=45
      id=200
      
      open(unit=lin1,file=infile1,status='old')
      open(unit=lin2,file=infile2,status='old')
      open(unit=lin3,file=infile3,status='old')
      open(unit=lin4,file=infile4,status='old')
      open(unit=lin5,file=infile5,status='old')
      
      call hropen(lout,'NTUPLE',outfile,'N',1024,istat)
     
      call hbnt(id,'New Ntuple',' ')
      call hbname(id,'ntupr',ievent,
     c     'ievent:I,nch:I,zenith(3):R,azimuth(3):R,delang(3):R,
     c     truezen(3):R,trueazi(3):R,trueen(3):R,trueid(3):I,
     c     ndirb(3):I,smootallphit(3):R,ldirc(3):R,jkprob(3):R, 
     c     nfit:I,jkchi(3):R,jkrchi(3):R,
     c     r14:R,intwght:R,atmswght:R')
      
      count = 0
 10   read(lin1,*,end=20) 
     c     ievent,nch,zenith(2),azimuth(2),delang(2)
      count = count + 1
      
 20   read(lin2,*,end=30) 
     c     truezen(2),trueazi(2),trueen(2),trueid(2) 

 30   read(lin3,*,end=40)
     c     ndirb(2),smootallphit(2),ldirc(2),jkprob(2) 

 40   read(lin4,*,end=50)
     c     nfit,jkchi(2),jkchi(3),jkrchi(2)

 50   read(lin5,*,end=60)
     c     r14,intwght,atmswght

      call hfnt(id)
      go to 10
      
     
 60   call hrout(id,icycle,'D')
      call hrend('NTUPLE')
      close(lin1)
      close(lin2)
      close(lin3)
      close(lin4)
      close(lin5)
      close(lout)
      end