$LoadAddOns={"FeynHelpers"};
$LoadFeynArts=True;
<<FeynCalc`
$FAVerbose=False;

$KeepLogDivergentScalelessIntegrals=True;

top=CreateTopologies[1,1->2,ExcludeTopologies->{WFCorrections,SelfEnergies,Tadpoles}];
diags=InsertFields[top,{V[5]}->{V[5],V[5]},Model->"SMQCD",InsertionLevel -> {Particles},ExcludeParticles->{F[__],S[_],V[_],U[1|2|3|4]}];
Paint[diags, ColumnsXRows -> {2, 1},SheetHeader -> False,SheetHeader->None,Numbering -> None,ImageSize->{512,256}];

amps=FCFAConvert[CreateFeynAmp[diags, Truncated -> True,GaugeRules->{},PreFactor->1],IncomingMomenta->{p1},
OutgoingMomenta->{p2,p3},LoopMomenta->{q},DropSumOver->True,UndoChiralSplittings->True,ChangeDimension->D,List->True,SMP->True,LorentzIndexNames->{λ,μ,ν},
SUNIndexNames->{a,b,c}]//Contract//SUNSimplify

amps2=amps[[2]]/.p3->-p2/.p2->p

amps1=amps[[1]]/.p3->-p2/.p2->p

xxx=TID[amps1+amps2,q,ToPaVe->True]

xxx[[1]]//PaXEvaluateUVIRSplit[#,q]&
xxx[[2]]//PaXEvaluateUVIRSplit[#,q]&
xxx[[3]]//PaXEvaluateUVIRSplit[#,q]&

res=PaXEvaluateUVIRSplit[xxx,q]

res//FCHideEpsilon