subroutine gen6_rap(r,p,wt6,*)
      implicit none
      include 'constants.f'
      include 'masses.f'
      include 'mxdim.f'
      include 'breit.f'

      integer nu
      double precision r(mxdim)
      double precision wt6,p(mxpart,4),tp(4),tm(4),bp(4),bm(4)
      double precision wtepnn,wtnbem,ep(4),em(4),nn(4),nb(4),wp(4),wm(4)
      double precision wtttb,wtwp,wtwm,s3min,wt0
      parameter(wt0=1d0/twopi**4)
      data wp,wm,ep,em,nn,nb,bp,bm/32*0d0/


*     q(-p1) +qbar(-p2)=t(nu(p3)+e^+(p4)+b(p5))                        *
*                       +t~(b~(p6)+e^-(p7)+nu(p8))+g(p9)               *
*                                                                      * 


      wt6=0d0
      mass2=mt
      mass3=mt
      call gen2m(r,p,wtttb,*999)
      wtttb=(mt*twidth*pi)**2*wtttb
      do nu=1,4
      tp(nu)=p(3,nu)
      tm(nu)=p(4,nu)
      enddo
      s3min=0d0
      n3=1
      mass3=wmass      
      width3=wwidth      

      call phi1_2m(mb,r(5),r(6),r(7),s3min,tp,bm,wp,wtwp,*999)
      call phi1_2m(mb,r(8),r(9),r(10),s3min,tm,bp,wm,wtwm,*999)
      call phi3m0(r(11),r(12),wp,nn,ep,wtepnn,*999)
      call phi3m0(r(13),r(14),wm,em,nb,wtnbem,*999)
      wt6=wt0*wtepnn*wtnbem*wtwp*wtwm*wtttb
*     q(-p1) +qbar(-p2)=t(nu(p3)+e^+(p4)+b(p5))                        *
*                       +t~(b~(p6)+e^-(p7)+nu(p8))               *
*                                                                      * 
      do nu=1,4
      p(3,nu)=nn(nu)
      p(4,nu)=ep(nu)
      p(5,nu)=bm(nu)

      p(6,nu)=bp(nu)
      p(7,nu)=em(nu)
      p(8,nu)=nb(nu)

      enddo

      return
 999  wt6=0d0
      return 1
      
      end