Name: john Date: 02/28/19-11:05:49 AM Z
I have just pasted what I have got in the output after running
amps = FCFAConvert[CreateFeynAmp[diags, Truncated ->
True,GaugeRules->{},PreFactor->1/((2^D)*(Pi)^(D/2))],
IncomingMomenta->{p},OutgoingMomenta->{p},LoopMomenta->{q},UndoChiralSplittings->True,ChangeDimension->D,List->False,SMP->True];
and then
ampsEval=(ChangeDimension[amps,D]//Contract//DiracSimplify)/.{DiracTrace->Tr}