PROGRAM TSTJET
C
C      test program
C      compare event list in jetset and standard structure
C
#include "pydat1.inc"
#include "stdlun.inc"
      integer istr

      INTEGER NEVT,I,NPRNT,LOK
C
      NEVT=100
      NPRNT=1
C initialize HEP logical units
      LNHWRT=23
      LNHRD=0
      LNHDCY=0
      LNHOUT=22
      MSTU(11)=LNHOUT
      OPEN(UNIT=LNHOUT,FILE='stdtstjx.lpt',STATUS='NEW')
      call stdxwinit('stdtstjx.io','STDHEP/Jetset example',
     1               NEVT,istr,LOK)
      call mcfio_FPrintDictionary(LNHOUT)
C
      CALL TSTBKEM(20)
C...
      CALL HPTRLSTJ
C          Write begin-run record   
      call stdxwrt(100,istr,LOK)
      IF(LOK.NE.0) WRITE(LNHOUT,*) ' Problem writing begin run record'
C
      DO 100 I=1,NEVT
	CALL PY2ENT(0.D0,2,-2,20.D0)
        CALL LUNHEP(1)
        IF(I.LE.NPRNT)THEN
          CALL PYLIST(2)
          CALL HEPLST(1)
          CALL HEPTREE
        ENDIF
        CALL TSTFLEM(20)
        call stdxwrt(1,istr,LOK)
 100  CONTINUE
C
C          Write end-of-run record  
      call stdxwrt(200,istr,LOK)
      IF(LOK.NE.0) WRITE(LNHOUT,*) ' Problem writing end run record'
c...close event file
      call stdxend(istr)
C...print histograms
      CALL HOUTPU(LNHOUT)
      CALL HISTDO
c...close output file
      CLOSE(UNIT=LNHOUT)
      STOP
      END