***********************************************************
*
      subroutine prthit97(prtdata,prt1ene,prt2ene,ierr)
*
*        returns: mean 'energy' values for the 7 PRT counter
*                 pairs which have a coincidence signal and
*                 good timing for at least on counter.
*
*        comment: in 1997 all PRT counters were operational;
*                 for now the same cuts as in 1996 will be
*                 used (needs further studies).
*
*        author: K.Desler
*        last change: 12/08/97
*
***********************************************************

      implicit none
      integer ierr,i
      real prt1ene,prt2ene
      integer energy

      integer ishft

      real prtdata(14,2)
      real e(14),t(2)

      data e /
c pulse height cuts PRT2, channels 1 - 10
     +     10., 10., 7., 8., 5., 10., 10., 11., 13., 12.,
c pulse height cuts PRT1, channels 11 - 14
     +     1.5, 2.0, 2.0, 2.5/

c timing cut
      data t / -20.0, 10.0/

      ierr = 0
      prt1ene = 0
      prt2ene = 0

c  PRT2: channels 1, 2
      if (prtdata(1,1).gt.e(1).and.prtdata(2,1).gt.e(2).and.
     +  ((prtdata(1,2).gt.t(1).and.prtdata(1,2).lt.t(2)).or.
     +   (prtdata(2,2).gt.t(1).and.prtdata(2,2).lt.t(2)))) then
          prt2ene = prt2ene + (prtdata(1,1)+prtdata(2,1))/2.
      endif

c  PRT2: channels 3, 4
      if (prtdata(3,1).gt.e(3).and.prtdata(4,1).gt.e(4).and.
     +  ((prtdata(3,2).gt.t(1).and.prtdata(3,2).lt.t(2)).or.
     +   (prtdata(4,2).gt.t(1).and.prtdata(4,2).lt.t(2)))) then
          prt2ene = prt2ene + (prtdata(3,1)+prtdata(4,1))/2.
      endif


c  PRT2: channels 7, 8
      if (prtdata(7,1).gt.e(7).and.prtdata(8,1).gt.e(8).and.
     +  ((prtdata(7,2).gt.t(1).and.prtdata(7,2).lt.t(2)).or.
     +   (prtdata(8,2).gt.t(1).and.prtdata(8,2).lt.t(2)))) then
          prt2ene = prt2ene + (prtdata(7,1)+prtdata(8,1))/2.
      endif


c  PRT1: channel 11, 12
      if (prtdata(11,1).gt.e(11).and.prtdata(12,1).gt.e(12).and.
     +  ((prtdata(11,2).gt.t(1).and.prtdata(11,2).lt.t(2)).or.
     +   (prtdata(12,2).gt.t(1).and.prtdata(12,2).lt.t(2)))) then
          prt1ene = prt1ene + (prtdata(11,1)+prtdata(12,1))/2.
      endif

c  PRT1: channel 13, 14
      if (prtdata(13,1).gt.e(13).and.prtdata(13,1).gt.e(13).and.
     +  ((prtdata(13,2).gt.t(1).and.prtdata(13,2).lt.t(2)).or.
     +   (prtdata(14,2).gt.t(1).and.prtdata(14,2).lt.t(2)))) then
          prt1ene = prt1ene + (prtdata(13,1)+prtdata(14,1))/2.
      endif

      return
      end