C
C count qq->tt or gg->tt
C
C     -lpdflib -lmathlib -lpacklib
C
      PROGRAM HWIGPR
C---COMMON BLOCKS ARE INCLUDED AS FILE HERWIG65.INC
      INCLUDE 'HERWIG65.INC'
      INTEGER N
      integer iii
      EXTERNAL HWUDAT
C---MAX NUMBER OF EVENTS THIS RUN
      MAXEV=100000
C      MAXEV=2000000
C---BEAM PARTICLES
      PART1='PBAR'
      PART2='P'
C---BEAM MOMENTA
      PBEAM1=980.
      PBEAM2=980.
C---PROCESS
      IPROC=1706
C---INITIALISE OTHER COMMON BLOCKS
      CALL HWIGIN
C---USER CAN RESET PARAMETERS AT
C   THIS POINT, OTHERWISE DEFAULT
C   VALUES IN HWIGIN WILL BE USED.
CCCCCCC  Spin correlation
      SYSPIN=.true.
      THREEB=.true.

C      AUTPDF(1)='MRS'
C      AUTPDF(2)='MRS'
C      MODPDF(1)=41
C      MODPDF(2)=41
      AUTPDF(1)='CTEQ'
      AUTPDF(2)='CTEQ'
      MODPDF(1)=46
      MODPDF(2)=46

C      PTMIN=100.
      PTMIN=0.
C---COMPUTE PARAMETER-DEPENDENT CONSTANTS
      CALL HWUINC
C---CALL HWUSTA TO MAKE ANY PARTICLE STABLE
C      CALL HWUSTA('PI0     ')
C---USER'S INITIAL CALCULATIONS
      CALL HWABEG
C---INITIALISE ELEMENTARY PROCESS
      CALL HWEINI
C---LOOP OVER EVENTS

C. force
      do iii=1,10
         BRFRAC(iii)=0.0
      enddo
      BRFRAC(1)=1.0
      BRFRAC(6)=1.0


      DO 100 N=1,MAXEV

C---INITIALISE EVENT
         CALL HWUINE
C---GENERATE HARD SUBPROCESS
         CALL HWEPRO
C---GENERATE PARTON CASCADES
         CALL HWBGEN
C---DO HEAVY OBJECT DECAYS
CCCCCCC  force W->e nu
         CALL HWDHOB
C---DO CLUSTER FORMATION
         CALL HWCFOR
C---DO CLUSTER DECAYS
         CALL HWCDEC
C---DO UNSTABLE PARTICLE DECAYS
         CALL HWDHAD
C---DO HEAVY FLAVOUR HADRON DECAYS
         CALL HWDHVY
C---ADD SOFT UNDERLYING EVENT IF NEEDED
         CALL HWMEVT
C---FINISH EVENT
         CALL HWUFNE
C---USER'S EVENT ANALYSIS
         CALL HWANAL
 100  CONTINUE
C---TERMINATE ELEMENTARY PROCESS
      CALL HWEFIN
C---USER'S TERMINAL CALCULATIONS
      CALL HWAEND
      STOP
      END
C----------------------------------------------------------------------
      SUBROUTINE HWABEG
C     USER'S ROUTINE FOR INITIALIZATION
C----------------------------------------------------------------------
      END
C----------------------------------------------------------------------
      SUBROUTINE HWAEND
C     USER'S ROUTINE FOR TERMINAL CALCULATIONS, HISTOGRAM OUTPUT, ETC
C----------------------------------------------------------------------
      END
C----------------------------------------------------------------------
      SUBROUTINE HWANAL
C     USER'S ROUTINE TO ANALYSE DATA FROM EVENT
      INCLUDE 'HERWIG65.INC'

      integer nu/0/, nd/0/, ng/0/, no/0/, nev/0/
C.. ------------------------------------------------------------
      nev=nev+1

      if ( IDHEP(4) .eq. -1 .and. IDHEP(5) .eq. 1 .or.
     1     IDHEP(4) .eq. 1 .and. IDHEP(5) .eq. -1 ) then
         print *, '@ddbar'
         nd=nd+1
      else if (IDHEP(4) .eq. -2 .and. IDHEP(5) .eq. 2 .or.
     1         IDHEP(4) .eq. 2 .and. IDHEP(5) .eq. -2 ) then
         print *, '@uubar'
         nu=nu+1
      else if (IDHEP(4) .eq. 21 .and. IDHEP(5) .eq. 21 ) then
         print *, '@gg'
         ng=ng+1
      else
         print *, '@', IDHEP(4), IDHEP(5)
         no=no+1
      endif
      if ( mod(nev, 100) .eq. 0) then
         print *, 'u=', nu, ' d=', nd, ' g=', ng, ' other=', no
      endif
      END

C...  EOF
C----------------------------------------------------------------------
C      END
C----------------------------------------------------------------------