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!
This archive was generated by hypermail 2b29 : 09/04/20-12:55:05 AM Z CEST