subroutine pjetfilt_event

C=======================================================================
C
C===Implicit None Declaration:
C
$$IMPLICIT
C
C===Global Declarations:

$$INCLUDE 'YBOS$LIBRARY:ERRCOD.INC'    ! YBOS error codes
$$INCLUDE 'YBOS$LIBRARY:BNKTYP.INC'    ! bank type codes
$$INCLUDE 'UIPACK$LIBRARY:UINAMES.INC' ! UIPACK declarations
$$INCLUDE 'UIPACK$LIBRARY:UIERROR.INC'
$$INCLUDE 'UIPACK$LIBRARY:UIPARAM.INC'
$$INCLUDE 'C$INC:ERLEVL.INC'           ! error severity levels for error logger
$$INCLUDE 'C$INC:JOBSTA.INC'           ! job status parameters
$$INCLUDE 'C$INC:BCS.INC'              ! primary YBOS array
$$INCLUDE 'C$INC:PIVALU.INC'           ! Values of pi
$$INCLUDE 'C$INC:EVTINFO.INC'
$$INCLUDE 'C$INC:JETINFO.INC'
$$INCLUDE 'stntuple$source:my_jetinfo.inc'
$$INCLUDE 'stntuple$source:exo_phoinfo.inc'
C
         integer i,j,k
         integer nphotight,phoindex(10)
         real coriso
C My analysis variables...

          LOGICAL PASS
          INTEGER ANPTRG,status

          PASS=.FALSE.

               
         nphotight = 0
         do i = 1,10
         phoindex(i) = 0
         enddo
         cor_version = 1

          CALL TOPFND(.TRUE.,.TRUE.,.TRUE.)
          call exo_pho

C         find a good photon
C         photon quantities are filled by exo_phoval.cdf into the
C         common block exo_phoinfo.inc
   
         do i = 1,npho
           if (phodetect(i).eq.0) then
           if ((phoet(i)/phocor(i).gt.20.0)
     &         .or.(phoet(i).gt.20.0)) then
           if (abs(phocesz(i)).gt.14.0.and.
     &         abs(phocesz(i)).lt.217.0) then
           if (abs(phocesx(i)).lt.21.5) then
           if (phon3d(i).le.1.and.abs(phopt(i)).lt.1.0) then
           if (phohadem(i).lt.(0.055+0.00045*phoe(i))) then
            coriso = Phoco4(i)-Phovcor(i)-
     &               Pholcor(i)*Phoet(i)/Phocor(i)
            if (PHOBOXISO(i).LT.4.0.or.
     &      min(coriso,phoco4(i)).lt.
     &      max(0.1*phoet(i),0.1*phoet(i)/phocor(i)) ) then
             nphotight = nphotight + 1
             phoindex(nphotight) = i
           endif
           endif
           endif
           endif
           endif
           endif
           endif
         enddo
             
C        final event cuts
         if (good_run) then
         if (nphotight.gt.0) then
           pass = .true.
         endif
         endif

          IF (PASS) THEN
            STATUS=ANPTRG(.TRUE.)
          ELSE
            STATUS=ANPTRG(.FALSE.)
         ENDIF

          RETURN
          END