LHE文件中错误报告重离子初始化

你能检查3.5.3吗？
您确定这不是Athena后处理问题吗（如果您使用该框架？）

我在文件中有这个：

399<初始化>
400 2212 2212 5.574400e+05 5.574400 e+05 0 0 0-4 1
401 1.461900e+05 3.975100e+02 1.461900e+05 1
402<生成器名称='MadGraph5_美国海军陆战队@NLO'版本='3.5.3'>请引用1405.0301
403</init>

对于以下脚本：

生成a>jj
输出
发射
设置lpp 2
设置ebeam 557440.0
设置pdlabel edff
设置nb_proton1 82
设置nb_neutron1 126
设置nb_proton2 82
设置nb_neutron2 126

干杯，

奥利维尔

谢谢奥利维尔！我确实看到，当独立运行时（在3.5.1和3.5.3中），质量能量中心得到了尊重，这表明我们的发射配置中存在一些可疑之处。话虽如此，但即使在你的例子中，粒子也肯定是错的（光束不是质子），对吗？

谢谢，
扎克

事实上，我从来没有尝试过按照重离子态可以计算的10位数代码设置这些数字。我可以修复它，但这将增加该部分中的数字数量，并且某些代码可能无法解析它（或拒绝刷新它/…）。

我会看看补丁。

奥利维尔

>2024年1月16日19:54，Zachary Marshall<隐藏的电子邮件地址>写道：
>
>谢谢奥利维尔！我在单机运行时看到了（在3.5.1和
>3.5.3）质量能量中心得到尊重，这表明
>我们的发射配置有些可疑。说到这里，因为
>当然，即使在你的例子中，粒子也是错误的（光束不是
>质子），对吗？
>
>谢谢，
>扎克
>
> --
>您收到此错误通知是因为您订阅了
> MadGraph5_aMC@NLO。
>https（https）：//错误.启动板.net/bugs/2049278
>
>标题：
>LHE文件中错误报告重离子初始化
>
>中的状态MadGraph5_aMC@NLO:
>新建
>
>错误描述：
>尊敬的作者：，
>
>我们已经看到MadGraph的LHE输出中明显存在的问题
>尝试以重离子初始状态运行时（如中所述
>https（https）：//阿西夫。组织/pdf/2207.03012.pdf格式). 该过程是光诱导的：
>
>生成a>y2>a
>
>并且卡被相应地设置为使用离子：
>
>2=lpp1！梁1型
>2=lpp2！梁2型
>557440.0=电子束1！束流1总能量（GeV）
>557440.0=ebeam2！光束2总能量（GeV）
>edff=pdlabel！PDF集合
>230000=lhaid！如果pdlabel=lhapdf，这是lhapdf编号
>82=nb_质子1
>126=nb_中子1
>82=nb_质子2
>126=nb_中子2
>
>LHE初始化行如下所示：
>
><初始化>
>2212 2212 6.500000e+03 6.500000e/03 0 0 247000 247000-4 1
>+1.5507000e-07+4.7623000e-10+1.5507000 e-07 1
><生成器名称='MadGraph5_美国海军陆战队@NLO'版本='3.5.1'>请引用1405.0301\
>运算器>
></init>
>
>这似乎是错误的-我们回到了13TeV质心质子质子
>碰撞。你能帮我们理解可能发生的事情吗
>错了吗？
>
>谢谢，
>扎克
>
>要管理有关此错误的通知，请访问：
>https（https）：//错误.启动板.net/mg5amcnlo公司/+错误/2049278/+订阅
>

只是一个简短的跟进：是的，看起来更复杂的跑步导致了重心的改变：

设置group_subprocesses自动
设置忽略_六夸克_进程错误
设置low_mem_多芯_nlo_generation错误
设置complex_mass_scheme为False
集合包含_轻子_已启动_进程错误
设定轨距单位
设置loop_optimized_输出True
设置loop_color_flows为False
设置max_npoint_对于通道0（_C）
设置默认值_非联轴器99
设置通道99的最大值
设置zerowidth_tchannel为True
将nlo_mixed_expansion设置为True
导入型号DMspin2已满
定义p=gucdsu~c~d~s~
定义j=gucdsu~c~d~s~
定义l+=e+mu+
定义l-=e-mu-
定义vl=ve vm vt
定义vl~=ve~vm~vt~
定义p=21 2 4 1 3-2-4-1-3 5-5
定义j=p
生成a>y2>a
输出aay2aa-nojpeg
启动aay2aa
更新ion_pdf
启动aay2aa
更新ion_pdf

最后，重心重置为6500 GeV。

干杯，
扎克

谢谢奥利维尔！我认为可能值得给Pythia8的人发一封电子邮件，看看他们打算做什么。我不知道他们是如何处理这些事情的。

最佳，
扎克

好的，让我随时更新。如果需要，我有可用的修补程序：

diff—git a/模板/LO/源/设置运行。f b/模板/LO/源/设置运行。（f）
指数ba99f895a。.876d325ed 100644
---a/模板/LO/源/设置运行。（f）
+++b/模板/LO/源/设置运行。（f）
@@-145,9+145,21@@c banner.py有望纠正此类错误的run_card。
C为Les Houches init信息填充公共块
do i=1,2
       if（lpp(i） 等式1。或.lpp(i） .eq.2），则
-idbmup（i）=2212
+if（nb_proton(i） 等式1。和nb_中子(i） .eq.0），然后
+idbmup（i）=2212
+埃尔塞夫（nb_proton）(i） .等式0。和nb_中子(i） .等式1），然后
+idbmup（i）=2112
+其他
+idbmup（i）=1000000000+（nb_proton(i） +nb_中子(i） ）*10
+$+nb_质子（i）*10000
+结尾
       埃尔塞夫(lpp（i）。等式-1或。lpp（i）。等式-2）则
-idbmup（i）=-2212
+if（nb_proton(i） 等式1。和nb_中子(i） .eq.0），然后
+idbmup（i）=-2212
+其他
+idbmup（i）=-1*（1000000000+（nb_proton(i） +nb_中子(i） ）*10
+$+nb_质子（i）*10000）
+结尾
       埃尔塞夫(lpp（i）。等式3）然后
       idbmup公司(i） =11
       埃尔塞夫(lpp（i）。等式-3）则
diff—git a/madgraph/各种各样的/横幅。py b/madgraph格式/各种各样的/横幅。第页
索引ece8d0a33。.bee9a0061 100755
---a/madgraph软件/各种各样的/横幅。第页
+++b/madgraph格式/各种各样的/横幅。第页
@@-2766,11+2766,22@@class RunCard（运行卡）(配置文件）:
                     fsock.writeline（锁定写入线）s（行）
       f锁定。关闭（）

-@staticmethod（静态方法）
-定义get_idbmup（lpp）：
+定义get_idbmup（self，lpp，beam=1）：
“”“返回粒子碰撞pdg代码”“”
如果lpp在（1,2，-1，-2）中：
-返回math.copysign（2212，lpp）
+目标=2212
+如果'nb_proton1'在self中：
+nbp=自['nb_proton%s'%beam]
+nbn=自['nb_中子%s'束]
+如果nbp==1且nbn==0：
+目标=2212
+elif nbp==0且nbn==1：
+目标=2112
+其他：
+目标=1000000000
+目标+=10*（nbp+nbn）
+目标+=10000*nbp
+返回math.copysign(目标，lpp）
（3，-3）中的elif lpp：
返回math.copysign（11，lpp）
（4，-4）中的elif lpp：
@@-2786,8+2797,8@@class RunCard(配置文件）:
lhe文件的<init>块的第一行。"""

输出={}
-output[“idbmup1”]=self.get_idbmup公司(自我['lpp1'])
-output[“idbmup2”]=self.get_idbmup公司(自我['lpp2'])
+output[“idbmup1”]=self.get_idbmup公司(自我['pp1']，梁=1）
+output[“idbmup2”]=self.get_idbmup公司(自我['lpp2']，梁=2）
       输出[“ebmup1”]=自己[“ebeam1”]
       输出[“ebmup2”]=自己[“ebeam2”]
       输出[“pdfgup1”] = 0

嗨，扎克，

我计划今天发布3.5.4，因此检查任何未解决的问题并面对这个问题？
你联系他们了吗？我应该按3.5.4中的补丁来包含来自光束的10位数信息吗？

现在我也意识到我们还没有解决能源问题。
在脚本中使用命令
更新_ pdf
它只在run_card.dat中添加一个块（默认值为p p）。
像这样，这样的命令毫无意义，因为您没有指定任何内容（至少从提供的脚本中）。
所以实际上，从你提供的脚本来看，我希望p p@13TeV

我不知道我们是否收到了皮提亚作者的回复；我刚刚检查了我们这边的联系人。我认为最好暂时不打补丁，以防这打破了他们的一些假设。

我相信删除启动命令中的多余垃圾会有所帮助，但我不能100%确定我们是否完全做到了。

感谢您的跟进！

最佳，
扎克