c--- This subroutine is a basic implementation of hep-ph/9801442 valid for 
c---- f(p1)+f(p2) --> gamma(p3)+gamma(p4)+f(p5)+f(p6)
      subroutine basic_2jet_frix(pin,passed,isub) 
      implicit none 
      include 'constants.f' 
      include 'frag.f'
      double precision pin(mxpart,4)
      integer isub
      logical passed,in_cone
      double precision ret_ET,R
      double precision R45,R35,R36,R46,ET3,ET4,ET5,ET6,
     & alpha_3,alpha_4

      passed=.true.
   
!------ THIS ROUTINE IS FOR LORD ONLY---------------------------------
   
c--- This routine also called for 'gmgmjt' -> always need second branch
!      if ( ((part .eq. 'lord') .or. (part .eq. 'virt')
!     &.or. (isub .eq. 1) ) .and.
!     &  ((case .eq. 'gamgam') .or. (case .eq. 'Higaga')) ) then
c---- LO, virtual and dipole pieces for gamgam and Higaga have no jets --> passed cut
!         return
 
c---  Real radiation (isub=0) and 'gmgmjt' should be checked
!      elseif (isub .eq. 0) then 
         
c---- Calculate angles and energies
         ET3=ret_ET(pin,3)
         ET4=ret_ET(pin,4)
         ET5=ret_ET(pin,5)
         ET6=ret_ET(pin,6) 
         
         R35=R(pin,3,5)
         R36=R(pin,3,6)
         R45=R(pin,4,5)
         R46=R(pin,4,6)
         
         
         alpha_3=ET3/(1d0-dcos(cone_ang)) 
         alpha_4=ET4/(1d0-dcos(cone_ang))
         
         
         if((R35 .lt. cone_ang) .and. (R45 .lt. cone_ang)) then  
            passed = in_cone(R35,ET5,alpha_3) 
            if (passed) then 
               passed = in_cone(R45,ET5,alpha_4)
            endif
         elseif((R36 .lt. cone_ang) .and. (R46 .lt. cone_ang)) then  
            passed = in_cone(R36,ET6,alpha_3) 
            if (passed) then 
               passed = in_cone(R46,ET6,alpha_4)
            endif
         elseif(R35.lt.cone_ang) then 
            passed=in_cone(R35,ET5,alpha_3)
         elseif(R36.lt.cone_ang) then 
            passed=in_cone(R36,ET6,alpha_3)
         elseif(R45.lt.cone_ang) then 
            passed=in_cone(R45,ET5,alpha_4)
         elseif(R46.lt.cone_ang) then 
            passed=in_cone(R46,ET6,alpha_4)
         else
            return 
         endif
            

!         if((R35 .lt. cone_ang) .and. (R45 .lt. cone_ang)) then  
!            passed = in_cone(R35,ET5,alpha_3) 
!            if (passed) then 
!               passed = in_cone(R45,ET5,alpha_4)
!            endif
!            return 
!         elseif (R35 .lt. cone_ang) then 
!            passed=in_cone(R35,ET5,alpha_3) 
!            return 
!         elseif (R45 .lt. cone_ang) then 
!            passed=in_cone(R45,ET5,alpha_4) 
!            return 
!         else
!            passed = .true. 
!            return 
!         endif                       

!      endif

      return 
      end