PROGRAM TST
C
C      test program
C      list standard particle definitions
C
#include "stdlun.inc"
      character*1000 filrd


C initialize HEP logical units
      lnhwrt=0
      lnhrd=23
      lnhrd=6
      lnhdcy=0
      lnhout=6
      lnhpdf=0
      lnhdmp=0

      nread=0
      maxread=100
C...initialize MCFio once only
      call mcfio_init()
C...now open and read events

          filrd='sample.stdhep'
          print *,'Enter input file name'
          read(5,*) filrd
          print *,'Input file is ',trim(filrd)
          print *,'Enter number of read records'
          read(5,*) maxread
          print *,'Max read records is ',maxrecord

          call stdxropen(filrd,ntries,istr,lok)
          if(lok.ne.0) write(lnhout,1001) filrd
          if(lok.ne.0) go to 200
          write(lnhout,1002) filrd
  10      call stdxrd(ifl,istr,lok)

          nread=nread+1
          if( ifl .ne. 4 .and.ifl.ne.1) then
             print *,' IFL .ne. 4 : nread=',nread,' ifl=',ifl
          endif

          print *,' nread=',nread,' ... ifl=',ifl

          if(lok.ne.0) go to 100
          if( nread .gt. maxread ) then 
           print *,'Number of read records exceeds maxread(',maxread,')'
           go to 100
          end if
C... histogram cross section if it exists

          if(ifl.eq.100) write(lnhout,1003) filrd, generatorname
C... keep looping if this is not a STDHEP record
          if(ifl.ne.4.and.ifl.ne.1) go to 10
C          call tstflem
          call heplst(1)
          call prtev4()
          if ( nread .lt. 100 ) go to 10
 100      continue
          call stdxend(istr)
 200      continue

C...print histograms
C      call hrput(0,'sample.hbook','N')
1001  format(' could not open file ',a20)
1002  format(' file ',a20,' has been opened')
1003  format(/' file ',a20,' was generated by ',a20)

      STOP
      END

      subroutine prtev4
#include "stdhep.inc"
#include "hepev4.inc"

c      common/HEPEVT/NEVHEP,NHEP,ISTHEP(NMXHEP),IDHEP(NMXHEP),
c     &JMOHEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP)
c      common/hepev4/eventweightlh, alphaqedlh, alphaqcdlh, scalelh(10),
c     1              spinlh(3,NMXHEP), icolorflowlh(2,NMXHEP), idruplh

       print *,' ********* hepev4 contents ************************'
       print *,' nevhep=',nevhep,' nhep=',nhep
       print *,' eventweightlh=',eventweightlh,' idruplh=',idruplh
       print *,' alphaqedlh=',alphaqedlh,' alphaqcdlh=',alphaqcdlh
       print *,' scalelh()=',(scalelh(k),k=1,5)
       print *,'           ',(scalelh(k),k=6,10)
       do 100 i=1, nhep
         print *,'i=',i,' spin=',spinlh(1,i),spinlh(2,i),spinlh(3,i),
     1   ' colorflow=',icolorflowlh(1,i),icolorflowlh(2,i)
100    continue
   
       return
       end