program read_03
       dimension xp(3)
       dimension pp(3)
       dimension pol(3)
       character cdum*72
      OPEN(UNIT=12,FILE='for003.dat',FORM='FORMATTED',STATUS='OLD')
       nmupos = 0
       npipos = 0
       nkpos = 0
       nmuneg = 0
       npineg = 0
       nkneg = 0
   10  continue
       read(12,'(1a72)') cdum
       print 02, cdum
02    format(1a72)
       read(12,'(1a72)') cdum
       print 02, cdum
       do 20 j=1,9999999
       read(12,*,end=16)ievt,ipnum,iptyp,ipflg,tp,wt,xp,pp,pol
      if(iptyp.eq.2) nmupos = nmupos + 1
      if(iptyp.eq.3) npipos = npipos + 1
      if(iptyp.eq.4) nkpos = nkpos + 1
      if(iptyp.eq.-2) nmuneg = nmuneg + 1
      if(iptyp.eq.-3) npineg = npineg + 1
      if(iptyp.eq.-4) nkneg = nkneg + 1
20     continue
16     continue
       print 333, nmupos,npipos,nkpos,nmuneg,npineg,nkneg
333    format(' mu+  pi+ k+  mu- pi- k-', 6i8)
       close(12)
       stop
       end