Thanks for the reply, though the problem persists. I actually tried your solution before but I ran into the same problem.
res=((M12 // FermionSpinSum) /. DiracTrace -> Tr) // Contract
I still get a product of traces, because there are still dummy indexes inside them:
tr(...GA[alpha]...GA[alphalin]..) tr (... GA[alpha]...GA[alphalin]...)
Everything is multiplying, so there are no "free indexes" M12 is a Lorentz scalar.
This archive was generated by hypermail 2b29 : 09/20/18-12:40:01 AM Z CEST