double complex function amp_qqagg_ql(i1,qh,i2,h2,i3,h3,i4,h4,
     . i5,lh,j6,j7)
      implicit none
      include 'constants.f'
      include 'sprods_com.f'
      include 'zprods_com.f'
      integer i1,i2,i3,i4,i5,i6,i7,j,k,qh,h2,h3,h4,lh,g1,g2,g3,g4,lg,
     . j6,j7
      double complex xa(mxpart,mxpart),xb(mxpart,mxpart)
      double precision t,iqcd
      double complex A7h1,A7h2,A7h3,A7h4,A7h5,A7h6,A7h7,A7h8
c-----
      iqcd=one
C---A(qh,h2,h3,h4,lh)
C---h=1 LH
C---h=2 RH
      g1=qh
      g2=h2
      g3=h3
      g4=h4
      lg=lh
      i6=j6
      i7=j7
      if (g1*lg.eq.2)then
      i7=j6
      i6=j7
      lg=3-lg
      endif

      if (g1.eq.2) then
      do j=1,mxpart
      do k=1,mxpart
      xa(j,k)=za(j,k)
      xb(j,k)=zb(j,k)
      enddo
      enddo
      elseif (g1.eq.1) then
      g1=2
      g2=3-g2
      g3=3-g3
      g4=3-g4
      lg=3-lg
      do j=1,mxpart
      do k=1,mxpart
      xa(j,k)=zb(k,j)
      xb(j,k)=za(k,j)
      enddo
      enddo

      endif

c-----A7h1
c-----A(2,2,2,2,2)
      if (
     . (g1.eq.2).and.(g2.eq.2).and.(g3.eq.2).and.(g4.eq.2).and.(lg.eq.2)
     . ) then
c-----A(2,2,2,2,2)
      A7h1=czip
      A7h1=xa(i5,i6)**2
     ./(xa(i1,i3)*xa(i3,i4)*xa(i4,i5)*xa(i2,i6)*xa(i2,i7))
      amp_qqagg_ql = A7h1
c-----A7h2
c-----A(2,2,2,1,2) =
      elseif
     .((g1.eq.2).and.(g2.eq.2).and.(g3.eq.2).and.(g4.eq.1).and.(lg.eq.2)
     . ) then
c-----A(2,2,2,1,2)
      A7h2=czip
      A7h2=1/(s(i3,i4)*t(i2,i6,i7)*xa(i2,i6)*xa(i2,i7))*(
     . xa(i4,i1)*xb(i1,i3)*( xa(i4,i1)*xb(i1,i2)*xa(i2,i6)
     . +xa(i4,i1)*xb(i1,i7)*xa(i7,i6)+xa(i4,i3)*xb(i3,i2)*xa(i2,i6)
     . +xa(i4,i3)*xb(i3,i7)*xa(i7,i6) )*xa(i5,i6)/xa(i1,i3)/t(i1,i3,i4)
     .-xa(i4,i5)*xb(i5,i3)*( xa(i6,i4)*xb(i4,i3)+xa(i6,i5)*xb(i5,i3) )*
     . ( xa(i6,i2)*xb(i2,i1)+xa(i6,i7)*xb(i7,i1) )/xb(i4,i5)/t(i3,i4,i5)
     .-( xa(i4,i1)*xb(i1,i2)*xa(i2,i6)+xa(i4,i1)*xb(i1,i7)*xa(i7,i6)
     .  +xa(i4,i3)*xb(i3,i2)*xa(i2,i6)+xa(i4,i3)*xb(i3,i7)*xa(i7,i6) )*
     . ( xa(i6,i4)*xb(i4,i3)+xa(i6,i5)*xb(i5,i3) )/xa(i1,i3)/xb(i4,i5)
     .)
      amp_qqagg_ql = A7h2
c-----A7h3
c-----A(2,2,1,2,2) =
      elseif
     .((g1.eq.2).and.(g2.eq.2).and.(g3.eq.1).and.(g4.eq.2).and.(lg.eq.2)
     . ) then
c-----A(2,2,1,2,2)
      A7h3=czip
      A7h3=1/(s(i3,i4)*t(i2,i6,i7)*xa(i2,i6)*xa(i2,i7))*(
     .-xb(i1,i4)**2*( xa(i3,i1)*xb(i1,i2)*xa(i2,i6)
     . +xa(i3,i1)*xb(i1,i7)*xa(i7,i6)+xa(i3,i4)*xb(i4,i2)*xa(i2,i6)
     . +xa(i3,i4)*xb(i4,i7)*xa(i7,i6) )*xa(i5,i6)/xb(i1,i3)/t(i1,i3,i4)
     .+xa(i3,i5)**2*( xa(i6,i3)*xb(i3,i4)+xa(i6,i5)*xb(i5,i4) )*
     . ( xa(i6,i2)*xb(i2,i1)+xa(i6,i7)*xb(i7,i1) )/xa(i4,i5)/t(i3,i4,i5)
     .+xb(i1,i4)*xa(i3,i5)*xa(i5,i6)*( xa(i6,i2)*xb(i2,i1)
     . +xa(i6,i7)*xb(i7,i1) )/xb(i1,i3)/xa(i4,i5) )
      amp_qqagg_ql = A7h3
c-----A7h4
c-----A(2,1,2,2,2)
      elseif
     .((g1.eq.2).and.(g2.eq.1).and.(g3.eq.2).and.(g4.eq.2).and.(lg.eq.2)
     . ) then
c-----A(2,1,2,2,2)
      A7h4=czip
      A7h4=(xa(i5,i2)*xb(i2,i7)+xa(i5,i6)*xb(i6,i7))**2
     ./(t(i2,i6,i7)*xa(i1,i3)*xa(i3,i4)*xa(i4,i5)*xb(i2,i6)*xb(i2,i7))
      amp_qqagg_ql = A7h4
c-----A7h5
c-----A(2,2,1,1,2) =
      elseif
     .((g1.eq.2).and.(g2.eq.2).and.(g3.eq.1).and.(g4.eq.1).and.(lg.eq.2)
     . ) then
c-----A(2,2,1,1,2)
      A7h5=czip
      A7h5=-(xa(i6,i2)*xb(i2,i1)+xa(i6,i7)*xb(i7,i1))**2
     ./(t(i2,i6,i7)*xb(i1,i3)*xb(i3,i4)*xb(i4,i5)*xa(i2,i6)*xa(i2,i7))
      amp_qqagg_ql = A7h5
c-----A7h6
c-----A(2,1,2,1,2) =
      elseif
     .((g1.eq.2).and.(g2.eq.1).and.(g3.eq.2).and.(g4.eq.1).and.(lg.eq.2)
     . ) then
c-----A(2,1,2,1,2)
      A7h6=czip
      A7h6=1/(s(i3,i4)*t(i2,i6,i7)*xb(i2,i6)*xb(i2,i7))*(
     . xa(i4,i1)*xb(i1,i3)*( xa(i4,i1)*xb(i1,i7)+xa(i4,i3)*xb(i3,i7) )*
     . ( xa(i5,i2)*xb(i2,i7)+xa(i5,i6)*xb(i6,i7) )/xa(i1,i3)/t(i1,i3,i4)
     .-xa(i4,i5)*xb(i5,i3)*( xb(i3,i4)*xa(i4,i2)*xb(i2,i7) 
     . + xb(i3,i4)*xa(i4,i6)*xb(i6,i7) + xb(i3,i5)*xa(i5,i2)*xb(i2,i7) 
     . + xb(i3,i5)*xa(i5,i6)*xb(i6,i7) )*xb(i1,i7)/xb(i4,i5)/t(i3,i4,i5)
     .-( xb(i3,i4)*xa(i4,i2)*xb(i2,i7) + xb(i3,i4)*xa(i4,i6)*xb(i6,i7) 
     .  +xb(i3,i5)*xa(i5,i2)*xb(i2,i7) + xb(i3,i5)*xa(i5,i6)*xb(i6,i7))*
     . ( xa(i4,i1)*xb(i1,i7)+xa(i4,i3)*xb(i3,i7) )/xa(i1,i3)/xb(i4,i5))
      amp_qqagg_ql = A7h6
c-----A7h7
c-----A(2,1,1,2,2) =
      elseif
     .((g1.eq.2).and.(g2.eq.1).and.(g3.eq.1).and.(g4.eq.2).and.(lg.eq.2)
     . ) then
c-----A(2,1,1,2,2)
      A7h7=czip
      A7h7=1/(s(i3,i4)*t(i2,i6,i7)*xb(i2,i6)*xb(i2,i7))*(
     .-xb(i1,i4)**2*( xa(i3,i1)*xb(i1,i7)+xa(i3,i4)*xb(i4,i7) )*
     . ( xa(i5,i2)*xb(i2,i7)+xa(i5,i6)*xb(i6,i7) )/xb(i1,i3)/t(i1,i3,i4)
     .+xa(i3,i5)**2*xb(i1,i7)*( xb(i4,i3)*xa(i3,i2)*xb(i2,i7)
     . +xb(i4,i3)*xa(i3,i6)*xb(i6,i7)+xb(i4,i5)*xa(i5,i2)*xb(i2,i7)
     . +xb(i4,i5)*xa(i5,i6)*xb(i6,i7) )/xa(i4,i5)/t(i3,i4,i5)
     .-xb(i1,i4)*xa(i3,i5)*( xa(i5,i2)*xb(i2,i7)+xa(i5,i6)*xb(i6,i7) )*
     . xb(i7,i1)/xb(i1,i3)/xa(i4,i5) )
      amp_qqagg_ql = A7h7
c-----A7h8
c-----A(2,1,1,1,2)
      elseif
     .((g1.eq.2).and.(g2.eq.1).and.(g3.eq.1).and.(g4.eq.1).and.(lg.eq.2)
     . ) then
c-----A(2,1,1,1,2)
      A7h8=czip
      A7h8=-xb(i1,i7)**2
     ./(xb(i1,i3)*xb(i3,i4)*xb(i4,i5)*xb(i2,i6)*xb(i2,i7))
      amp_qqagg_ql = A7h8
c----
      else
      amp_qqagg_ql = czip
      endif
c      write(*,*) qh,h2,h3,h4,lh
c      write(*,*) amp_qqagg_ql
c-----done
      return
      end