double complex function BSYALslggpppm(e1,p2,e4,p3,za,zb,zab,zba)
      implicit none
C-----Authors: John Campbell and Keith Ellis, November 2011
C---- arXiv:1101.5947 [hep-ph], Eq. (97)
C---- (These are twiddle functions, c.f.arXiv:1101.5947[hep-ph],Eq.(91))
      include 'constants.f'
      include 'masses.f'
      include 'zprods_decl.f'
      include 'sprods_com.f'
      include 'zabprods_decl.f'
      include 'qdef.f'
      include 'massiveintegrals.f'
      double precision s12,s23,mt2,mt3
      integer e1,p2,p3,e4,j
      j=p2-1
      mt2=mt**2
      mt3=mt**3
      s12=mt2+s(1,p2)
      s23=s(p2,p3)
      BSYALslggpppm =  + I3m12x3x4(j) * (  - 1.D0/2.D0/(zab(p2,q1,p3))*
     &    za(e1,p3)*za(e4,p3)*zb(p3,p2)*zab(p2,q1,p2)*mt + 1.D0/2.D0/(
     &    zab(p2,q1,p3))*za(e4,e1)*zab(p2,q1,p2)*zab(p3,q1,p2)*mt - 1.D0
     &    /2.D0/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p3)*za(e4,p3)*zab(p2,
     &    q1,p2)*mt3 + 1/(zab(p2,q1,p3))/(zab(p2,q1,p3))*za(e4,p3)*zab(
     &    p2,q1,p2)*zab(e1,q1,p3)*mt3 + 1/(zab(p2,q1,p3))/(zab(p2,q1,p3
     &    ))*za(e4,e1)*zab(p2,q1,p2)*s12*mt3 - 1/(zab(p2,q1,p3))/(zab(
     &    p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4,p3)*zab(p2,q1,p2)**2*
     &    mt3 + 1.D0/2.D0/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(zab(p2,q1,p3
     &    ))/(za(p3,p2))*za(e1,p2)*za(e4,p2)*zab(p2,q1,p2)**2*zab(p3,q1
     &    ,p3)*mt3 )
      BSYALslggpppm = BSYALslggpppm + I312x3x4(j) * ( 1.D0/2.D0/(zab(p2
     &    ,q1,p3))*za(e1,p3)*za(e4,p3)*zb(p3,p2)*zab(p2,q1,p2)*mt - 1.D0
     &    /2.D0/(zab(p2,q1,p3))*za(e4,e1)*zab(p2,q1,p2)*zab(p3,q1,p2)*
     &    mt + 1.D0/2.D0/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p3)*za(e4,p3
     &    )*zab(p2,q1,p2)*mt3 - 1/(zab(p2,q1,p3))/(zab(p2,q1,p3))*za(e4
     &    ,p3)*zab(p2,q1,p2)*zab(e1,q1,p3)*mt3 - 1/(zab(p2,q1,p3))/(
     &    zab(p2,q1,p3))*za(e4,e1)*zab(p2,q1,p2)*s12*mt3 + 1/(zab(p2,q1
     &    ,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4,p3)*zab(p2,
     &    q1,p2)**2*mt3 - 1.D0/2.D0/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(
     &    zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4,p2)*zab(p2,q1,p2)
     &    **2*zab(p3,q1,p3)*mt3 )
      BSYALslggpppm = BSYALslggpppm + I3m13x2x4(j) * (  - 1.D0/2.D0/(
     &    zab(p2,q1,p3))*za(e1,p3)*za(e4,p3)*zb(p3,p2)*zab(p3,q1,p3)*mt
     &     + 1.D0/2.D0/(zab(p2,q1,p3))*za(e4,e1)*zab(p3,q1,p2)*zab(p3,
     &    q1,p3)*mt - 3.D0/2.D0/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p3)*
     &    za(e4,p3)*zab(p3,q1,p3)*mt3 + 1/(zab(p2,q1,p3))/(zab(p2,q1,p3
     &    ))*za(e1,p2)*za(e4,p3)*zb(p3,p2)*zab(p3,q1,p3)*mt3 + 1/(zab(
     &    p2,q1,p3))/(zab(p2,q1,p3))*za(e4,e1)*zab(p3,q1,p3)*s12*mt3 - 
     &    2.D0/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*
     &    za(e4,p3)*zab(p2,q1,p2)*zab(p3,q1,p3)*mt3 + 1.D0/2.D0/(zab(p2
     &    ,q1,p3))/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2
     &    )*za(e4,p2)*zab(p2,q1,p2)*zab(p3,q1,p3)**2*mt3 )
      BSYALslggpppm = BSYALslggpppm + I313x2x4(j) * ( 1.D0/2.D0/(zab(p2
     &    ,q1,p3))*za(e1,p3)*za(e4,p3)*zb(p3,p2)*zab(p3,q1,p3)*mt - 1.D0
     &    /2.D0/(zab(p2,q1,p3))*za(e4,e1)*zab(p3,q1,p2)*zab(p3,q1,p3)*
     &    mt + 3.D0/2.D0/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p3)*za(e4,p3
     &    )*zab(p3,q1,p3)*mt3 - 1/(zab(p2,q1,p3))/(zab(p2,q1,p3))*za(e1
     &    ,p2)*za(e4,p3)*zb(p3,p2)*zab(p3,q1,p3)*mt3 - 1/(zab(p2,q1,p3)
     &    )/(zab(p2,q1,p3))*za(e4,e1)*zab(p3,q1,p3)*s12*mt3 + 2.D0/(
     &    zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4,p3
     &    )*zab(p2,q1,p2)*zab(p3,q1,p3)*mt3 - 1.D0/2.D0/(zab(p2,q1,p3))
     &    /(zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4,
     &    p2)*zab(p2,q1,p2)*zab(p3,q1,p3)**2*mt3 )
      BSYALslggpppm = BSYALslggpppm + I41x2x4x3 * (  - 1.D0/2.D0/(zab(
     &    p2,q1,p3))*za(e1,p3)*za(e4,p3)*zb(p3,p2)*s23*mt3 + 1.D0/2.D0
     &    /(zab(p2,q1,p3))*za(e1,p3)*za(e4,p3)*zb(p3,p2)*zab(p2,q1,p3)*
     &    zab(p3,q1,p2)*mt - 1.D0/2.D0/(zab(p2,q1,p3))*za(e4,e1)*zab(p2
     &    ,q1,p3)*zab(p3,q1,p2)**2*mt + 1.D0/2.D0/(zab(p2,q1,p3))*za(e4
     &    ,e1)*zab(p3,q1,p2)*s23*mt3 - 3.D0/2.D0/(zab(p2,q1,p3))/(za(p3
     &    ,p2))*za(e1,p3)*za(e4,p3)*zab(p2,q1,p2)*zab(p3,q1,p3)*mt3 + 
     &    1/(zab(p2,q1,p3))/(zab(p2,q1,p3))*za(e1,p2)*za(e4,p3)*zb(p3,
     &    p2)*zab(p2,q1,p2)*zab(p3,q1,p3)*mt3 + 1/(zab(p2,q1,p3))/(zab(
     &    p2,q1,p3))*za(e4,e1)*zab(p2,q1,p2)*zab(p3,q1,p3)*s12*mt3 - 2.D
     &    0/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)*za(e4
     &    ,p3)*zab(p2,q1,p2)**2*zab(p3,q1,p3)*mt3 + 1.D0/2.D0/(zab(p2,
     &    q1,p3))/(zab(p2,q1,p3))/(zab(p2,q1,p3))/(za(p3,p2))*za(e1,p2)
     &    *za(e4,p2)*zab(p2,q1,p2)**2*zab(p3,q1,p3)**2*mt3 )

c---- SB overall sign change to put in agreement with eq. (97)
      BSYALslggpppm = -BSYALslggpppm

      return
      end