Trace Calculation Using FORM

FORM by J.Vermaseren,version 3.0(Jan 28 2001) 
Run at: Sat Mar 10 23:14:35 2001
    *
    SYMBOL m;
    VECTORS pi,pf,k,kp,e,ep;
    OFF STATS;
    *
    * B&D T1 p130.
    LOCAL M11 = g_(1,pf,ep,e,k,pi,k,e,ep);
    LOCAL M12 = g_(1,pf,ep,e,k,   k,e,ep);
    LOCAL M13 = g_(1,   ep,e,k,pi,k,e,ep);
    LOCAL M14 = g_(1,   ep,e,k,   k,e,ep);
    LOCAL T1 = M11 + M12 + M13 + M14;
    * B&D T2 p131.
    LOCAL M21 = g_(1,pf,e,ep,kp,pi,kp,ep,e);
    LOCAL M22 = g_(1,pf,e,ep,kp,   kp,ep,e);
    LOCAL M23 = g_(1,   e,ep,kp,pi,kp,ep,e);
    LOCAL M24 = g_(1,   e,ep,kp,   kp,ep,e);
    LOCAL T2 = M21 + M22 + M23 + M24;
    * B&D T3 p131.
    LOCAL M31 = g_(1,pf,ep,e,k,pi,kp,ep,e);
    LOCAL M32 = g_(1,pf,ep,e,k,   kp,ep,e);
    LOCAL M33 = g_(1,   ep,e,k,pi,kp,ep,e);
    LOCAL M34 = g_(1,   ep,e,k,   kp,ep,e);
    LOCAL T3 = M31 + m*M32 + m*M33 + m^2*M34;
    * B&D T3 p131.
    LOCAL M41 = g_(1,pf,e,ep,kp,pi,k,e,ep);
    LOCAL M42 = g_(1,pf,e,ep,kp,   k,e,ep);
    LOCAL M43 = g_(1,   e,ep,kp,pi,k,e,ep);
    LOCAL M44 = g_(1,   e,ep,kp,   k,e,ep);
    LOCAL T4 = M41 + m*M42 + m*M43 + m^2*M44;
    *
    TRACE4,1;
    *
    ID pf    = pi + k - kp;
    ID k.k   =  0;
    ID kp.kp =  0;
    ID e.e   = -1;
    ID ep.ep = -1;
    ID k.e   =  0;
    ID kp.ep =  0;
    ID pi.e  =  0;
    ID pi.ep =  0;
    ID k.kp  = k.pi - k.pf;
    ID k.pf  = kp.pi;
    ID pi.pi = m^2;
    *
    *
    PRINT T1,T2,T3,T4;
    .SORT

   T1 =
      8*pi.k*pi.kp + 16*pi.k*k.ep^2;

   T2 =
      8*pi.k*pi.kp - 16*pi.kp*kp.e^2;

   T3 =
      16*pi.k*pi.kp*e.ep^2 - 8*pi.k*pi.kp + 8*pi.k*kp.e^2 - 8*pi.kp*k.ep^2;

   T4 =
      16*pi.k*pi.kp*e.ep^2 - 8*pi.k*pi.kp + 8*pi.k*kp.e^2 - 8*pi.kp*k.ep^2;

    .END