添加交互顺序时错误丢弃R2

席琳提供的其他信息：

亲爱的各位，

我们一直在测试与非直瞄SMEFT轴向异常相关的处理，尤其是g g>z。
当Z在壳层上时，这个过程取消了，所以我用重τ（177GeV）测试g g>ta+ta-，因为这个过程对于无质量费米子消失了。我使用处于最终状态而非顶部的τ来分离产生Z的耦合（与SM相比，顶部的额外贡献）和衰减到最终状态的耦合。我从一个更简单的EFT模型（SMEFTNLOaxial attached）开始，其中我只有修改夸克z顶点的算符（对左右耦合有额外贡献，但没有新的洛伦兹结构）和顶胶子相互作用（色磁算符）。我们已经手动检查了R2，它们是正确的

MG结果似乎存在几个问题：

brs检验只适用于处于最终状态的τ，而不适用于顶部，似乎存在Lambda^-4剩余项，但振幅本身不仅应为振幅的平方。

平方安培与预期值不匹配。循环是已知的（见所附的Cen注释，方程.1.2，它已经由Andres独立检查过），我添加了来自z传播子和τ衰变的额外因子，并取平方。平方安培与ga_ta^2 mta^2/mz^4 aS^2 ga_top^2成正比，其中ga_ta和ga_top是τ和顶部与z的轴向耦合，我没有发现eta和mz标度，而另一个很好。aS很明显，但不依赖于矢量耦合。我用一个具有相同洛伦兹结构的操作符替换了这个循环（1.2中的levi-civita位于（）之间），然后我得到了这个树级结果的完美一致。ME-go-MG比我自己的计算要小得多（5阶），所以我不认为我们在大幅度抵消后会被数值不确定性所支配。

我创建了一个新的NLO模型，其中SM的QCD部分加上4个独立的轴向联轴器到z，带有不同耦合顺序的向上（仅第1代）/向下（仅第一代）/tau/top（NPu/NPd/NPta/NPT，该模型作为轴测），以试图找出问题，但我当时无法进行检查。然而，对于向上或向下或两者一起，这两个矩阵元素的相隔系数为2.25（我的比MG大）。该系数为常数，因此在这种情况下，振幅具有正确的mta和mz相关性。另一方面，当我只做顶部时，这不起作用，我没有得到正确的mta和mz依赖性，并且在大mta限制下，我没有恢复无质量的情况（两个MG结果是一个系数4，而地雷有10%的差异）。

很抱歉收到这么长的邮件，但希望这能帮助您找到问题。

最佳，

赛琳

瓦伦丁，对此有何评论？

你好，Celine/Olivier，

Celine附上的模型似乎提前更新了您较长的电子邮件，这是需要测试的最新版本吗？

无论如何，我尝试在MG v3.1.1中运行它（为此我必须将模型转换为v3），我可以确认我得到了：

```
MG5_aMC>检查brs g g>g z QCD=99 QED=1 NP=2[virt=QCD]
[...]
测量结果：
过程[virt=QCD]矩阵BRS比率结果
g g>g z 1.9378359556e-03 4.2873465380e+00 2.2124403904e+03失败
摘要：0/1通过，1/1失败
失败的进程：g g>g z
```

当然，虽然删除了NP贡献，即返回到纯SM情况，但它仍然有效：

```
MG5_aMC>检查brs g>g z[virt=QCD]
[...]
测量结果：
过程[virt=QCD]矩阵BRS比率结果
g g>g z 9.0039116013e-02 1.6706986246e-19 1.8555253523e-18通过
摘要：1/1通过，0/1失败
```

然而，当我这样做时：
```
MG5_aMC>导入模型SMEFTNLOWZtaugsctGonly公司
MG5_aMC>集合已确认_v3.1_语法正确
MG5_aMC>生成g g>g z QCD=99 QED=1 NP=2[virt=QCD]
MG5_aMC>输出CelineTEST
MG5_aMC>发布-f
```

然后我检查生成的代码，即：
```
CelineTEST公司/子流程/P0_gg_gz号/helas呼叫_ampb_ 1。（f）
```
我发现以下情况：
```
[...]
C 58号回路图的计数器振幅
致电VVV21_0（宽（1,1），W（1,2），W（1,3），W（1,4），R2GC_127_18安培(1,34))
[...]
```

其中，`R2GC_127_18`是g g g Z的R2耦合，其耦合顺序为`NP=2`。
因此，它似乎被正确地包括在这个案例中。

那么@Celine，是什么让你认为它被从MadLoop生成的过程输出中删除/缺失了？