跳转至

控制衰变------衰变文件(Decfiles)

学习目标

  • 了解衰变文件(decfile)的基本结构

衰变文件(DecFile)控制衰变本身(即EvtGen的行为),并提供任何事件类型特定的配置(例如生成器切割)。 它们位于Gen/DecFile包中,我们在开始时已经检出并构建了该包,以.dec文件的形式存在于Gen/DecFiles/dkfiles目录中。 例如,事件类型27175000的.dec文件是Dst_D0pi,KKmumu=DecProdCut.dec

# EventType: 27175000
#
# Descriptor:[D*(2010)+ -> (D0 -> {K+ K- mu+ mu-}) pi+]cc
#
# NickName: Dst_D0pi,KKmumu=DecProdCut
#
# Cuts: DaughtersInLHCb
#
# Documentation: Forces the D* decay in generic b-bbar / c-cbar events + Requires products to be in LHCb acceptance
# EndDocumentation
#
# PhysicsWG: Charm
# Tested: Yes
# Responsible: Luisa Arrabito
# Email: unknown@<nospam>cern.ch
# Date: 20091215
#

Alias MyD0 D0
Alias MyantiD0 anti-D0
ChargeConj MyD0 MyantiD0

Decay D*+sig
  1.000 MyD0 pi+    VSS;
Enddecay
CDecay D*-sig

Decay MyD0
  1.000 K+ K- mu+ mu- PHSP;
Enddecay
CDecay MyantiD0
#
End 
注释部分实际上包含关键信息,在编译过程中会被解析,以创建<event-type>.py文件,该文件作为选项之一提供给Gauss! EventType是一系列标志,控制生成过程。其规则在LHCb-2005-034中有详细描述,还有一个网站允许你构建/解析事件类型。 例如,第一位数字的含义:1=包含b夸克,2=包含c夸克,3=最小偏置事件...

同样,该文档规定了"NickName"的命名约定------它也必须与文件名一致。请注意,一旦从某个衰变文件生成了蒙特卡洛(MC)样本,该衰变文件就不允许被修改,因此当你试图理解现有样本时,无需担心正在查看的DecFiles版本。 "Cuts"字段指定使用预设切割工具集中的哪一个。额外的Python代码可以添加到生成的配置文件中:

# InsertPythonCode:
# code ...
# EndInsertPythonCode

衰变文件的底部部分指定了衰变本身: 该衰变文件定义了一个信号D*+粒子,它100%衰变为D0pi+,而D0又100%衰变为K+K-mu+mu-。重要的是" MyD0"的定义。如果衰变目标是"D0"而不是"MyD0",那么D0会通过DECAY.DEC中实现的所有衰变模式衰变。 每个衰变的最后部分是使用的实际物理模型------在这种情况下,"VSS"是矢量到标量-标量跃迁的通用矩阵元,"PHSP"是仅相空间模型(矩阵元=常数)。请注意,使用PHSP时,衰变产物是完全非极化的------对于除(自旋0)到(自旋0 自旋0)之外的任何衰变,这会导致角衰变分布不正确!