使用面向对象技术实现的BESIII主漂移径迹重建程序.docVIP

BES III主漂移径迹重建软件MdcPatRec的设计与实现( 张 瑶1;1)，李卫东2;2)，张学尧1;3)，毛泽普2;4)，马秋梅2，，，，，，，，，，，，，，，，，，，，，，，，，，，，，，(1.山东大学物理系济南 250100中国科学院高能物理研究所实验物理中心北京 100049北京大学物理学院技术物理系北京 100871湖南大学应用物理系长沙 410082中国科学院研究生院北京 100049引言 北京谱仪BES III[1]是北京正负电子对撞机BEPC II[2]上的新一代高精度探测器。BES III的主漂移室(MDC)主要用来精确测量从相互作用点产生的带电粒子的动量，并测量粒子能量损失以进行粒子鉴别。MDC径迹重建软件MdcPatRec是基于BaBar漂移室径迹重建软件开发的[3]，负责MDC中带电径迹的动量和方向的精确重建，是BES III离线数据处理的最重要环节之一。本文首先描述BES III主漂移室几何结构和离线软件框架的设计；然后，介绍径迹重建算法数据流和事例模型的设计，以及重建流程中的主要步骤；最后，使用MdcPatRec算法对蒙特卡罗模拟数据进行重建，给出算法的主要性能指标。 BES III主漂移室与BES III离线软件系统 主漂移室内径为126cm，外径为1620cm包含43个同心丝层的圆柱型精密丝室外室的端面板阶梯状，每四个信号读出层为一个超层漂移单元基于小单元设计每个单元分布单元按90°对称排列外室其它部分的单元按22.5°对称排列氦基混合气作为工作气体π×93% 。 BES III离线软件[4](BES III Offline Software System，BOSS)是根据BES III实验的需开发的离线数据处理。基于框架Gaudi []开发的，以数据为中心瞬态数据各种服务模块(Service)转换器Converters)等模块之间通过标准的接口数据数据处理相对独立瞬态数据宿存数据独立管理φ平面投影为圆弧线，s-z平面上为直线，其中s为r-φ平面上径迹的飞行距离。带电粒子在穿过的路径上留下的电离粒子在电场的作用下漂移到信号丝上，留下一系列击中信号，提供粒子的位置信息。MdcPatRec算法依据原始击中信号，使用模式匹配的方法进行径迹的寻找，并使用最小二乘法拟合，得到径迹的物理信息。重建算法MdcPatRec继承服务模块BEPC II采用多束团对撞机制，在线数据获取系统无法给出每一个事例的准确时间起点(EvTime)，所以事例的起始时间由算法EsTimeAlg进行计算。MdcPatRec软件包计算出带电径迹的径迹参数、动量等物理量，并输出为重建后的径迹MdcTrack。  图 1 MdcPatRec数据流 MdcPatRec的事例数据模型 主漂移室的数据重建利用寻迹算法，从大量探测单元的原始击中信息中识别出穿越主漂移室的圆或螺旋线径迹，然后对这些径迹进行拟合，得出相应带电粒子的动量和运动方向。从数据的角度看，主漂移室重建是一个由原始数据MdcDigi组合成带电径迹MdcTrack的过程。因此，建立一个统一的事例数据模型很重要，它可以被应用于数据准备、径迹寻找与拟合等数据重建的不同阶段。MdcPatRec事例数据的设计模型UML图如图2所示。  图 2 MdcPatRec事例数据的模型UML图 在重建前的数据准备过程中，每个MdcDigi对象被转换成主漂移室击中MdcHit，转换后的MdcHit将保留它与MdcDigi的关联。原始探测器击中MdcDigi所包含的Identifier唯一标识了主漂移室的一根信号丝，MdcGeomSvc根据Identifier定位探测单元，用以提供该信号丝的几何位置信息。 径迹重建以MdcHit为输入，通过径迹寻找与拟合算法形成径迹。寻迹算法首先在r-φ平面进行圆径迹的寻迹和拟合，然后在三维空间进行螺旋线径迹的寻找和拟合。为统一不同空间的寻迹和拟合，引入径迹表示TrkRep管理不同拟合过程中的径迹。TrkRep具有一个径迹上击中的列表，还包含一个描述路径属性的Trajectory和进行径迹拟合的接口。当径迹分别进行二维和三维拟合时，TrkRep会根据径迹路径的性质包含不同类型的径迹参数和Trajectory对象：圆弧线表示TrkCircleRep和螺旋线表示TrkHelixRep分别包含圆径迹TrkCircleTraj对象和螺旋线径迹TrkHelixTraj对象。 Trajectory表示三维空间中的路径，用于计算在路径上给定路径长度处的位置、方向、曲率矢量等信息。Trajectory可以描述粒子的飞行路径，还可以描述击中单元的几何信息，所有与几何相关的计算由基类Trajectory的接口实现。所以，径迹重建模块与主漂移室几何间的通信能够通过T