application comis quit
      subroutine newnt
      integer part1,part2
      real en1,px1,py1,pz1,en2,px2,py2,pz2
      character*100 chform
      call hropen(42,'NTUPLE','out.rz','N',1024,istat)
      open(unit=41,file='in.dat',status='old')
      call hbnt(1,'pairs',' ')
      chform = 'part1:I,part2:I,en1,px1,py1,pz1,en2,px2,py2,pz2'
      call hbname(1,'pairs',part1,chform)
10    read(41,fmt=30,end=20) part1,en1,px1,py1,
     +pz1,part2,en2,px2,py2,pz2
30    format(I1,4E14.6,I2,4E14.6)
      call hfnt(1)
      goto 10
20    close(41)
      call hrout(1,icycle,' ')
      call hrend('NTUPLE')
      close(42)
      end
      call newnt
      end
quit