! functions for counting particles 

      integer function count_photo()
      implicit none
      include 'constants.f'
      integer j
      logical is_photon

      count_photo=0
      do j=1,mxpart
         if (is_photon(j)) then 
            count_photo=count_photo+1
         endif
      enddo 

      return
      end

      
      integer function count_jets()
      implicit none
      include 'constants.f'
      include 'npart.f'
      integer j
      logical is_hadronic

c---- Count final state jets
      count_jets=0

      do j=3,npart+2
         if (is_hadronic(j)) then 
            count_jets=count_jets+1
         endif
      enddo 

      return 
      end