环路中的费米子流问题

你用check函数检查了born/loop/real矩阵元素了吗？
如果这是一个与费米子流有关的问题，它很可能会发现这个问题。

干杯，

奥利维尔

你好，奥利维尔，

我直接跑了
检查d d~>s s~NP=2 QCD=0 QED=0[QCD]
检查d~d>s s~NP=2 QCD=0 QED=0[QCD]
仪表测试失败。它适用于其他一些子流程。我不明白的是，它如何作为完全相同的顶点工作/反义词到处都在使用。有什么想法吗？

干杯，

本吉

>2020年10月13日22:33，Olivier Mattelaer<隐藏电子邮件地址>写道：
>
>你用check函数检查了born/loop/real矩阵元素了吗？
>如果这是一个与费米子流有关的问题，它很可能会发现这个问题。
>
>干杯，
>
>奥利维尔
>
> --
>您收到此错误通知是因为您订阅了该错误
>报告。
>https（https）：//错误.启动板.net/bugs/1899608
>
>标题：
>环路中的费米子流问题
>
>中的状态MadGraph5_aMC@NLO:
>新建
>
>错误描述：
>大家好，
>
>我们正在一个模型中执行一些新的物理计算，其中
>SM被一个喇叭（和一个无关的不可见状态）所扩展
>为我们的目的）。UFO在附件中提供。以下过程
>兴趣在于dijet生产。我使用MG5_aMC 2.6.8作为
>进程定义，我键入
>
>导入模型RPVMono_UFO_RG公司-纸张
>设置complex_mass_scheme为True
>生成p p>j j NP=2 QCD=0 QED=0[QCD]
>输出
>
>为了得到新物理分量的QCD修正，我修改了madgraph/循环_图表_生成。py文件。我已将筛选设置为true（大约l388）
>编辑过滤器_manual=真
>然后包括在内
>loop_pdgs=[abs（x）表示diag.get中的x_循环线（_L）_pdgs（）]
>is_loop_gluon=（在loop_pdgs中为21）
>is_loop_sigma=（9000002在loop_pdgs中）
>如果是loop_sigma而不是loop_glion：
>valid_diag=假
>我得到了所需的图表（这样就可以工作了）。
>
>然而，当启动横截面计算时，我遇到了极点不可取消的问题。奇怪的是，极点对于子流程来说很好
>d dbar到s sbar或s sbar到d dbar
>但不适用于子流程
>dbar d到s sbar
>这很奇怪，因为图表完全相同。我怀疑这里的费米子流有问题（因为我们有相互冲突的箭头图）。有什么想法吗？
>
>干杯，
>
>本吉
>
>要管理有关此bug的通知，请转到：
>https（https）：//错误.启动板.net/mg5amcnlo公司/+错误/18899608/+订阅

你好，奥利维尔，

这里有消息吗？我能帮我做点什么吗？

干杯，

本吉

在我的例子中，当尝试生成这些循环时，它们无法生成：

命令“check d d~>s s~NP=2 QCD=0 QED=0[QCD]”中断，出现错误：
ValueError:-5不在列表中
请报告此错误https（https）：//错误.启动板.net/mg5amcnlo
有关更多信息，请参阅“MG5_debug”。
请将此文件附加到您的报告中。

干杯，

奥利维尔

>2020年10月16日08:09，Benjamin Fuks<隐藏电子邮件地址>写道：
>
>你好，奥利维尔，
>
>这里有消息吗？我能帮我做点什么吗？
>
>干杯，
>
>本吉
>
> --
>您收到此错误通知是因为您订阅了
> MadGraph5_aMC@NLO。
>https（https）：//错误.启动板.net/bugs/1899608
>
>标题：
>环路中的费米子流问题
>
>中的状态MadGraph5_aMC@NLO:
>新建
>
>错误描述：
>大家好，
>
>我们正在一个模型中执行一些新的物理计算，其中
>SM由平方（和不相关的不可见状态）扩展
>为我们的目的）。UFO在附件中提供。以下过程
>兴趣在于dijet生产。我使用MG5_aMC 2.6.8作为
>进程定义，我键入
>
>导入模型RPVMono_UFO_RG公司-纸张
>设置complex_mass_scheme为True
>生成p p>j j NP=2 QCD=0 QED=0[QCD]
>输出
>
>为了得到新物理分量的QCD修正，我修改了madgraph/循环_图表_生成。py文件。我已将筛选设置为true（大约l388）
>编辑过滤器_manual=真
>然后包括在内
>loop_pdgs=[diag.get中x的abs（x）_循环线（_L）_pdgs（）]
>is_loop_gluon=（在loop_pdgs中为21）
>is_loop_sigma=（9000002在loop_pdgs中）
>如果是is_loop_sigma而不是is_loop_gluon：
>valid_diag=错误
>我得到了所需的图表（这样就可以工作了）。
>
>然而，当启动横截面计算时，我遇到了极点不可取消的问题。奇怪的是，极点对于子流程来说很好
>d dbar到s sbar或s sbar到d dbar
>但不适用于子流程
>dbar d到s sbar
>这很奇怪，因为图表完全相同。我怀疑费米子流有问题（因为我们有冲突的箭头图）。有什么想法吗？
>
>干杯，
>
>本吉
>
>要管理有关此bug的通知，请转到：
>https（https）：//错误.启动板.net/mg5amcnlo公司/+错误/1899608/+订阅

这太奇怪了。你和我使用的是同一版本的mg5（2.6.8版）吗？

干杯，

B>（B）

>2020年10月16日09:33，Olivier Mattelaer<隐藏电子邮件地址>写道：
>
>在我的例子中，当试图生成这些循环时，它们无法生成
>具有：
>
>命令“check d d~>s s~NP=2 QCD=0 QED=0[QCD]”中断，出现错误：
>ValueError:-5不在列表中
>请在上报告此错误https（https）：//错误.启动板.net/mg5amcnlo
>有关更多信息，请参阅“MG5_debug”。
>请将此文件附加到您的报告中。
>
>干杯，
>
>奥利维尔
>
>>2020年10月16日08:09，Benjamin Fuks<隐藏电子邮件地址>写道：
>>
>>你好，奥利维尔，
>>
>>这里有消息吗？我能帮我做点什么吗？
>>
>>干杯，
>>
>>本吉
>>
>> --
>>您收到此错误通知是因为您订阅了
>> MadGraph5_aMC@NLO。
>>https（https）：//错误.启动板.net/bugs/1899608
>>
>>标题：
>>环路中的费米子流问题
>>
>>中的状态MadGraph5_aMC@NLO:
>>新建
>>
>>错误描述：
>>大家好，
>>
>>我们正在一个模型中执行一些新的物理计算，其中
>>SM被一个喇叭（和一个无关的不可见状态）所扩展
>>为我们的目的）。UFO在附件中提供。以下过程
>>兴趣在于dijet生产。我使用MG5_aMC 2.6.8作为
>>进程定义，我键入
>>
>>导入模型RPVMono_UFO_RG公司-纸张
>>将complex_mass_scheme设置为True
>>生成p p>j j NP=2 QCD=0 QED=0[QCD]
>>输出
>>
>>为了只对新的物理组件进行QCD校正，我修改了madgraph/循环_图表_生成。py文件。我已将筛选设置为true（大约l388）
>>编辑过滤器_manual=真
>>然后包括在内
>>loop_pdgs=[diag.get中x的abs（x）_循环线（_L）_pdgs（）]
>>is_loop_gluon=（在loop_pdgs中为21）
>>is_loop_sigma=（9000002在loop_pdgs中）
>>如果是loop_sigma而不是loop_glion：
>>valid_diag=假
>>我得到了所需的图表（这样就可以工作了）。
>>
>>然而，当启动横截面计算时，我遇到了极点不可取消的问题。奇怪的是，极点对于子流程来说很好
>>d dbar到s sbar或s sbar到d dbar
>>但不适用于子流程
>>dbar d到s sbar
>>这很奇怪，因为图表完全一样。我怀疑费米子流有问题（因为我们有冲突的箭头图）。有什么想法吗？
>>
>>干杯，
>>
>>本吉
>>
>>要管理有关此bug的通知，请转到：
>>https（https）：//错误.启动板.net/mg5amcnlo公司/+错误/1899608/+订阅
>
> --
>您收到此错误通知是因为您订阅了该错误
>报告。
>https（https）：//错误.启动板.net/bugs/1899608
>
>标题：
>环路中的费米子流问题
>
>中的状态MadGraph5_aMC@NLO:
>新建
>
>错误描述：
>大家好，
>
>我们正在一个模型中执行一些新的物理计算，其中
>SM被一个喇叭（和一个无关的不可见状态）所扩展
>为我们的目的）。UFO在附件中提供。以下过程
>兴趣在于dijet生产。我使用MG5_aMC 2.6.8作为
>过程定义，我键入
>
>导入模型RPVMono_UFO_RG公司-纸张
>设置complex_。。。

这个问题似乎不是来自循环。但Benj将重新生成模型，并检查问题是否仍然存在

你好，
如果出现问题，请从失败的文件夹中附加checkpoles.log文件。

干杯，

马尔科

>2020年10月16日11:58，Olivier Mattelaer<隐藏电子邮件地址>写道：
>
>这个问题似乎不是来自循环。但Benj会重新-
>生成模型并检查问题是否仍然存在
>
>**更改时间：mg5amcnlo
>状态：新建=>未完成
>
> --
>您收到此错误通知是因为您是错误受让人。
>https（https）：//错误.启动板.net/bugs/1899608
>
>标题：
>环路中的费米子流问题
>
>中的状态MadGraph5_aMC@NLO:
>不完整
>
>错误描述：
>大家好，
>
>我们正在一个模型中执行一些新的物理计算，其中
>SM被一个喇叭（和一个无关的不可见状态）所扩展
>为我们的目的）。UFO在附件中提供。以下过程
>兴趣在于dijet生产。我使用MG5_aMC 2.6.8作为
>进程定义，我键入
>
>进口型号RPVMono_UFO_RG公司-纸张
>设置complex_mass_scheme为True
>生成p p>j j NP=2 QCD=0 QED=0[QCD]
>输出
>
>为了得到新物理分量的QCD修正，我修改了madgraph/循环_图表_生成。py文件。我已将筛选设置为true（大约l388）
>编辑过滤器_manual=真
>然后包括在内
>loop_pdgs=[diag.get中x的abs（x）_循环（_L）_pdgs（）]
>is_loop_gluon=（在loop_pdgs中为21）
>is_loop_sigma=（在loop_pdgs中为9000002）
>如果是loop_sigma而不是loop_glion：
>valid_diag=假
>我得到了所需的图表（这样就可以工作了）。
>
>然而，当启动横截面计算时，我遇到了极点不可取消的问题。奇怪的是，极点对于子流程来说很好
>d dbar到s sbar或s sbar到d dbar
>但不适用于子流程
>dbar d到s sbar
>这很奇怪，因为图表完全相同。我怀疑费米子流有问题（因为我们有冲突的箭头图）。有什么想法吗？
>
>干杯，
>
>本吉
>
>要管理有关此bug的通知，请转到：
>https（https）：//错误.启动板.net/mg5amcnlo公司/+错误/1899608/+订阅

你好，Marco，

我试过一个新的UFO，但没有什么变化。因此，UFO似乎正常。检查如下：
信息：check_poles的结果：
信息：电杆成功取消20分以上的20分（容差=1.0e-05）
信息：P0_ddx_ssx
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
信息：电杆成功取消20分以上的20分（容差=1.0e-05）
信息：P0_dsx_dsx
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
信息：极点成功取消20点以上的20点（容差=1.0e-05）
信息：P0_sd_sd
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
信息：电杆成功取消20分以上的20分（容差=1.0e-05）
信息：P0_sdx_sdx
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
信息：电杆成功取消20分以上的20分（容差=1.0e-05）
信息：P0_ssx_ddx
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
信息：电杆成功取消20分以上的20分（容差=1.0e-05）
信息：P0_dxd_ssx
信息：test_ME的结果：
信息：通过。
信息：check_poles的结果：
在“calculate_xsect”中检测到错误
写入调试文件/用户/fuks/工作/工具/madgraph/2.6.8/PROCNLO公司_RPV蒙诺_UFO_RG公司-纸张_7/运行_03_tag_1_debug。日志
如果您需要此问题的帮助，请联系我们https（https）：//答案。启动板。净值/mg5amcnlo
aMCatNLO错误：磁极未取消，运行无法继续

我会在下面附上最后一份日志，因为其他的都很好。

你好，瓦伦丁，

在与Marco在skype上进行离线讨论后，发现这个问题与回路中费米子流的处理有关。因此，我将此问题分配给您：-）

如果你想更快地治疗，请随时通过skype与我联系。

干杯，

本吉

PS：我已经用一个新鲜的和旧的UFO（来自我们常见的2015年论文）进行了测试，并且用2.6.8（python2）和2.8.1（python 3）都进行了测试。问题是一样的。

你好，瓦伦丁，
本杰明检查了当对虚拟对象进行独立输出时，它们返回相同但不同的极点/有限部分
d天~>秒~
和
d~d>s s~。

我恐怕我无能为力。。。。

谢谢！

干杯，

马尔科

大家好，

有什么消息吗？

干杯，

本吉

大家好，

这最终解决了吗？

干杯，

本吉

大家好，

在与Valentin进行了现场讨论后，以下修复（我无耻地从Valentin的电子邮件中复制了粘贴内容）解决了这个问题。请随意查看并将其合并到主分支。

干杯，

本吉

<><><><><><><><><><>

我们需要修改helas_objects.py（check_majorana中的第1349行_和_flip_流量）并更换管路

              new_wf=波函数[波函数。索引（新_wf）]

通过
if（非新wf。获取（'is_循环'）或（new_wf。获取（'pdg_代码'）>0):
              索引_wf=波函数。索引（新_世界粮食计划署）
              其他：
              对于i_wf，枚举中的wf(波函数）:
                     如果new_wf==wf和wf.get('pdg_code')==新_wf.get公司('pdg_code'):
                            index_wf=i_wf
                            打破
              其他：
                     提升值错误
              new_wf=波函数[索引wf]

这是必要的，因为对于循环波函数，外部分支的回收需要注意PDG代码，而不仅仅是粒子的属性/电荷：

#在优化输出模式下，启动回路波函数
#如果具有相同的pdg，则可以回收，因为无论其pdg是什么
#它具有相同的系数和环路动量零，
#因此，原则上无需将其添加到
#图表_波函数。然而，这对于功能来说是必要的
#检查和_修复终止_流以正确更新的依赖项
#外部L-切majorana波函数的先前图表
#需要翻转。
如果没有自我优化_输出：
              波函数Nnumber=波函数离子编号+1
              外部_循环wf。集合（“数字”，波函数Nu编号）
              图表_波函数。追加(外部_loop_wf）
其他：
              尝试：
              外部的，外部的_循环wf=\
                     外部_循环wfs_dict[外部_循环wf。获取（'pdg_代码'）]