Name: Yan-Qing Ma Date: 02/10/09-02:40:08 PM Z
Maybe you can calculate them like what I did.
For example, when you get the following expression,
In[5]:= %/.DiracTrace->Tr
Out[5]= -2 (4 (m2 m4 - p2.p4) 2 p3.ep[p1])
p3.ep(*)[p1]
replace ep[p1] by LorentzIndex[mu] and ep(*)[p1] by LorentzIndex[nu], and then multiply the expression by a tensor which is the sum of polarizaion (maybe it is MT[mu,nu] or something else).