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}