logical function filterWbbmas()
c--- this routine is specific to the "Wbbmas" processes 401-408;
c--F also for processes 520-529
c--- it inspects the jets to check whether an event should be included 
c--- routine returns FALSE if event does not pass the process-specific cuts
      implicit none
      include 'constants.f'
      include 'clustering.f'
      include 'jetlabel.f'
      logical veto3jets
      integer nproc
      common/nproc/nproc

      if (jets.lt.1) then
         filterWbbmas=.false.
         return
      endif


      filterWbbmas=.true.

c--F  3-jet veto to match with ATLAS analysis. In general we do not want such a veto
      veto3jets = .false.

c---  20 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4)) +b(p5)+b~(p6) [massive]'
      if ((nproc .eq. 20) .or. (nproc .eq. 25)) then
        if (jets .lt. 2) then
          filterWbbmas=.false.
          return
c--- note: jets should already contain bq and ba because of bbproc
        endif
      endif

c--  New final state slicing: Wb, W(bb~), Wbb

c--  401 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+b(p5) [1, 2 or 3 jets, 4FNS]'
      if ((nproc .eq. 401) .or. (nproc .eq. 406)) then
         if ((inclusive .eqv. .false.) .and. (jets .ne. 1)) then
            filterWbbmas=.false.
            return
         endif
         if (jets .eq. 1) then
            if((jetlabel(1).ne.'bq').and.(jetlabel(1).ne.'ba')) then
               filterWbbmas=.false.
               return
            endif
         endif
         if (jets .eq. 2) then
            if ((jetlabel(1).eq.'bb').or.(jetlabel(2).eq.'bb')) then
               filterWbbmas=.false.
               return
            endif
         endif
         if ((jets .eq. 3) .and. (veto3jets .eqv. .true.)) then
            filterWbbmas=.false.
            return
         endif
      endif

c--  402 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+(b+b~)(p5) [1 or 2 jets, 4FNS]'

      if ((nproc .eq. 402) .or. (nproc .eq. 407)) then
         if ((inclusive .eqv. .false.) .and. (jets .ne. 1)) then
            filterWbbmas=.false.
            return
         endif
         if ((jets .eq. 1) .and. (jetlabel(1) .ne. 'bb')) then
            filterWbbmas=.false.
            return
         endif
         if (jets .eq. 2) then
            if ((jetlabel(1).ne.'bb').and.(jetlabel(2).ne.'bb')) then
               filterWbbmas=.false.
               return
            endif
         endif
         if (jets .eq. 3) then
            filterWbbmas=.false.
            return
         endif
      endif

c--  403 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+b(p5)+b~(p6) [2 or 3 jets, 4FNS]'
      if ((nproc .eq. 403) .or. (nproc .eq. 408)) then
         if ((jets .eq. 3) .and.
     &       ((inclusive.eqv..false.).or.(veto3jets.eqv..true.))) then
            filterWbbmas=.false.
            return
         endif
         if (jets .eq. 1) then
            filterWbbmas=.false.
            return
         endif
         if (jets .eq. 2) then
            if ((jetlabel(1).eq.'pp').or.(jetlabel(2).eq.'pp')) then
               filterWbbmas=.false.
               return
            endif
         endif
      endif

      
      return
      end
      

c--- This code is now deprecated as of v6.1.
c--- The functionality is (mostly) replicated in processes 520-527.

cc--- 420 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+b(p5) [1 or 2 jets only]'
c      if ((nproc .eq. 420) .or. (nproc .eq. 425)) then
c        if ((jets .gt. 2) .or. (jets .lt. 1)) then
c          filterWbbmas=.false.
c          return
c        endif
c        if ((inclusive .eqv. .false.) .and. (jets .eq. 2)) then
c          filterWbbmas=.false.
c          return
c        endif
c        if (jets .eq. 1) then
c         if ((jetlabel(1) .ne. 'bq') .and. (jetlabel(1) .ne. 'ba')) then
c            filterWbbmas=.false.
c            return
c          endif
c        endif
c        if (jets .eq. 2) then
c          if ((jetlabel(1) .eq. 'bb') .or. (jetlabel(2) .eq. 'bb')) then
c            filterWbbmas=.false.
c            return
c          endif
c        endif
c      endif
c      
cc--- 421 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+(b+b~)(p5) [1 or 2 jets only]'
c      if ((nproc .eq. 421) .or. (nproc .eq. 426)) then
c        if ((jets .gt. 2) .or. (jets .lt. 1)) then
c          filterWbbmas=.false.
c          return
c        endif
c        if ((inclusive .eqv. .false.) .and. (jets .eq. 2)) then
c          filterWbbmas=.false.
c          return
c        endif
c        if (jets .eq. 1) then
c          if ((jetlabel(1) .ne. 'bb')) then
c            filterWbbmas=.false.
c            return
c          endif
c        endif
c        if (jets .eq. 2) then
c         if ((jetlabel(1) .ne. 'bb') .and. (jetlabel(2) .ne. 'bb')) then
c            filterWbbmas=.false.
c            return
c         endif
c        endif
c      endif
c      
cc--- 422 '  f(p1)+f(p2) --> W^+(-->nu(p3)+e^+(p4))+b(p5)+jet(p6) [2 or 3 jets only]'
c      if ((nproc .eq. 422) .or. (nproc .eq. 427)) then
c        if (jets .lt. 2) then
c          filterWbbmas=.false.
c          return
c        endif
c        if ((inclusive .eqv. .false.) .and. (jets .eq. 3)) then
c          filterWbbmas=.false.
c          return
c        endif
c        if (jets .eq. 2) then
c          if ((jetlabel(1) .eq. 'bb') .or. (jetlabel(2) .eq. 'bb')) then
c            filterWbbmas=.false.
c            return
c          endif
c        endif
c      endif