TID got wrong result

[xysxxs]

Hello,

When I calculate a hzz loop with TID, it does not get right result. As the figure shows, a and b should be the same, but somehow b is wrong. The only difference is replace k1->k2+k3 before or after TID. Anyone knows why?

Thanks in advance

SPD[k1, k1] = s; SPD[k2, k2] = mz^2; SPD[k3, k3] = mz^2; SPD[k1, k2] = (SPD[k1, k1] + SPD[k2, k2] - SPD[k3, k3])/2; SPD[k2, k3] = -((SPD[k3, k3] + SPD[k2, k2] - SPD[k1, k1])/2);

-I/Pi^2 FAD[{q, mh}, {q + k1, mh}, {q + k2,
mz}] (MTD[[Alpha], [Beta]] -
FVD[q + k2, [Alpha]] FVD[q + k2, [Beta]]/mz^2)
a = FCE[ExpandScalarProduct[
Simplify[TID[%, q, ToPaVe -> True, UsePaVeBasis -> True]] /.
k1 -> k2 + k3]] /. FVD[k2, [Alpha]] -> 0 /.
FVD[k3, [Beta]] -> 0;

-I/Pi^2 FAD[{q, mh}, {q + k1, mh}, {q + k2,
mz}] (MTD[[Alpha], [Beta]] -
FVD[q + k2, [Alpha]] FVD[q + k2, [Beta]]/mz^2) /. k1 -> k2 + k3
b = FCE[ExpandScalarProduct[
Simplify[TID[%, q, ToPaVe -> True, UsePaVeBasis -> True]]]] /.
FVD[k2, [Alpha]] -> 0 /. FVD[k3, [Beta]] -> 0;

[vsht]

Hi,

did you check that the difference between a and b doesn't vanish upon applying PaXEvaluate?