Name: Yan-Qing Ma Date: 01/08/09-08:05:25 AM Z
Instead of using GSD, maybe DiracGamma works.
e.g. DiracGamma[Momentum[a]]/.Momentum[a]:>s1*Momentum[p]+s2*Momentum[q]
You can change them to dimension D if you like. Good luck!