$LoadFeynArts=True;
<<FeynCalc`
$FAVerbose=False;

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

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

Check with Pascual and Tarrach

amps2=(amps/.{p2->-p2,p3->-p3});

ampsPT=(TID[#,q,ToPaVe->True,UsePaVeBasis->True]&/@amps2);

uvParts={PaVe[1,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[0,0,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->(1/(64 EpsilonUV (π^4))),PaVe[1,1,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[1,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[2,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[0,0,1,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->-(1/(192 EpsilonUV (π^4))),PaVe[0,0,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->-(1/(192 EpsilonUV (π^4))),PaVe[1,1,1,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[1,1,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[1,2,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0,PaVe[2,2,2,{SPD[p2,p2],SPD[p2,p2]+2 SPD[p2,p3]+SPD[p3,p3],SPD[p3,p3]},{0,0,0},PaVeAutoOrder->True,PaVeAutoReduce->True]->0};

ampsPT2=ampsPT/.FCI[uvParts]

ampsPTUVpart1=Total[ampsPT2]//FCE//Collect2[#,{Epsilon,MTD}]&

(* p2+p3 = -p1 *)
ampsPTUVpart2=ampsPTUVpart1/.{FVD[p2,i_]+2FVD[p3,i_]:>FVD[p3,i]-FVD[p1,i],
FVD[p3,i_]+2FVD[p2,i_]:>FVD[p2,i]-FVD[p1,i]}

PTVertexFuLorentzStruct[{p_,q_,k_},{mu_,nu_,si_},{a_,b_,c_}]:=
-I SUNF[a,b,c](MTD[mu,nu]FVD[p-q,si]+MTD[nu,si]FVD[q-k,mu]+MTD[si,mu]FVD[k-p,nu])

Cancel[FCI[ExpandScalarProduct[ampsPTUVpart2]]/ExpandScalarProduct[FCI[PTVertexFuLorentzStruct[{p1,p2,p3},{λ,μ,ν},{a,b,c}]]]]//Simplify
res=%/(I SMP["g_s"])

PTGhostTriangleResult=SMP["g_s"]^2/(4Pi)^2  CA/8 1/(3EpsilonUV)

res-PTGhostTriangleResult