subroutine jetsort(et,eta,maxjet)
c
c Sort jet by energy
c input:  et,  transverse enegry of the jets (real)
c         eta, pseudorapidity of the jets (real)
c         maxjet, number of jets
c
c output: et1, transverse enegry of the jets after ordering(real)
c         eta1, pseudorapidity of the jets after ordering(real)

      IMPLICIT NONE
      integer i,j,k,l,ierr,maxjet
      integer max1
      real et0,eta0,et(*),eta(*)

      max1=maxjet+1
      do i=1,maxjet
         et0 = et(i)
         l=i+1
         k=0
c
c     check highest jet in the remaining sample
c
         do j=l,max1
            if (et(j).gt.et0) then
               et0 = et(j)
               eta0 = eta(j)
	       k=j
            endif
         enddo
c
c     classify jets
c         
         if (k.eq.0) goto 100
         et(k) = et(i)
         eta(k) = eta(i)
         et(i) = et0
         eta(i) = eta0
 100  enddo
         
      return
      end