subroutine detkt(vtx)
      IMPLICIT NONE
c
c  selects jet using the kT algorithm. The routine which selects jets is
c  ktjzes(vtx,sira_flag,sira_prob,ebeam,pbeam). The input parameters
c  can be entered by input cards
c  The jets are then sorted by Et (subroutine sort) and a fla for forward
c  jets is set (bfkl_flag).
c
#include "ktjjets.inc"
#include "ktcut.inc"
#include "jdiszes.inc"
#include "zecpbm.inc"

      real vtx(3),sira_prob
      real ebeam,pbeam,x_jet(20)
      real et_jet(20),eta_jet(20)
      integer sira_flag,ierr,i,flag
c
c Call kt package (put together by S.Eisenhart)
c
c T. Haas: Use proper beam energies
c     ebeam = -27.52
c     pbeam = 820.
      ebeam = -BEAMOM(2)
      pbeam =  BEAMOM(1)
      sira_flag = 1
      sira_prob = 0.9
c
      call ktjzes(vtx,sira_flag,sira_prob,ebeam,pbeam)
c
c     The result of this routine is store in the common black
c     ktjjets.
c
c
c     First sort jets by transverse energy with ycut = 0.1
c
      call vfill(et_jet,20,-999.99)
      call vfill(eta_jet,20,-999.99)
      call vfill(x_jet,20,-999.99)
      njet_b=njkt(1)
      do i=1,njkt(1)
         et_jet(i) = pjkt(1,i,1)
         eta_jet(i) = pjkt(1,i,2)	 
	 x_jet(i) = pjkt(1,i,4)	 
      enddo
c
      call jetsort(et_jet,eta_jet,njkt(1))
c
c     store the 2 highest Et jets in the Zes common block
c
      et1_b = et_jet(1)
      et2_b = et_jet(2)
      eta1_b = eta_jet(1)
      eta2_b = eta_jet(2)
c
c     set the bfkl flag
c     The cuts are defined in forjets.inc
c 
      flag = 0
      call bfkl_flag(et_jet,eta_jet,x_jet,njkt(1),flag)
      Fjetb = flag
c
c     Then sort jets with ycut = 0.5
c
      call vfill(et_jet,20,-999.99)
      call vfill(eta_jet,20,-999.99)      
      njet_c=njkt(2)
      do i=1,njkt(2)
         et_jet(i) = pjkt(2,i,1)
         eta_jet(i) = pjkt(2,i,2)
	 x_jet(i) = pjkt(2,i,4)	 
      enddo
c
      call jetsort(et_jet,eta_jet,njkt(2))
c
c     store the 2 highest Et jets in the Zes common block
c
      et1_c = et_jet(1)
      et2_c = et_jet(2)
      eta1_c = eta_jet(1)
      eta2_c = eta_jet(2)
c
c     set the bfkl flag
c     The cuts are defined in forjets.inc
c 
      flag = 0
      call bfkl_flag(et_jet,eta_jet,x_jet,njkt(2),flag)
      Fjetc = flag
c
c
      return
      end