PROGRAM TBHTEST

C...Double precision and integer declarations.
      IMPLICIT DOUBLE PRECISION(A-H, O-Z)
      IMPLICIT INTEGER(I-N)
      INTEGER PYK,PYCHGE,PYCOMP
      EXTERNAL PYDATA

C...Commonblocks.
      COMMON/PYSUBS/MSEL,MSELPD,MSUB(500),KFIN(2,-40:40),CKIN(200)
      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
      COMMON/PYDAT1/MSTU(200),PARU(200),MSTJ(200),PARJ(200)
      COMMON/PYDAT2/KCHG(500,4),PMAS(500,4),PARF(2000),VCKM(4,4)

C...Set processes to run
      MSEL = 0
      MSUB(401) = 1
      MSUB(402) = 1

C     Number of tries to find highest Xsect with 2->3 variables (D=10)
      MSTP(129) = 1000

C     Set H+ mass and tan(beta)
      PMAS(37,1) = 250.
      PARU(141)  = 30.

      niter = 1000
      ECM   = 1.4d4

C...Initialize PYTHIA
      CALL PYINIT('CMS','p+','p+',ECM)

C...Start PYTHIA loop (generate events)
      DO 20 J=1,niter

C     Generate one event
         CALL PYEVNT

C     Print first 3 events
         IF(J.le.3) CALL PYLIST(2)

 20   CONTINUE

      CALL PYSTAT(1)

      END