应用Kalman filter建立时间投影室的径迹重建算法研究.pdfVIP

第 30卷 第 1期 计 算 物 理 VoI．3O．No．1 2013年 1月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS Jan．，2013 文章编号：1001．246X(2013)O1 127—07 应用 Kalmanfilter建立时间投影室的径迹重建算法 李 波 ， 高原宁 ， 杨振伟 ， 李玉兰 (1．粒子技术与辐射成像教育部重点实验室(清华大学)，北京 100084；2．清华大学工程物理系，北京 100084) 摘 要：建立直线径迹模型，详细验证将其用于Kalmanfilter小型时间投影室径迹重建算法的正确性．结果表明 径迹参数在一定的击中点下能够收敛到真实值．实际的径迹重建程序在MarlinTPC框架下实现．通过合理地选择 径迹重建参数 ，程序能够有效地对噪声进行甄别并重建出径迹．本算法对大型时间投影室的径迹重建有借鉴 意义． 关键词 ：时间投影室；径迹重建 ；Kalmanfilter；径迹模型 中图分类号：0572．21 文献标识码 ：A 0 引言 时间投影室 (TPC，timeprojectionchamber)是一种带 电粒子径迹探测器…，一般用气体作为探测介质． 当带电粒子经过 TPC时，使气体电离而产生电子一离子对 ，电子在 TPC内部均匀电场作用下，向读 出端盖漂 移并产生雪崩放大．电子信号由按照一定规则排列的若干读 出盘(pad)读出．TPC通过读出盘记录的信号确 定径迹 方 向的位置 ，方向(即电子漂移的方向)位置由 = 。(t一：)确定．其 中，。为电子在电场 中的漂 移速度，时间t，为带电粒子穿过TPC的时间，由外部时间探测器的触发信号给出；t为电子到达读出面的时 间．除了根据带电粒子的三维径迹得到其动量之外 ，TPC还能通过读出盘上收集到的电荷量得到带电粒子的 能量损失 dE／dx并以此进行粒子鉴别． 未来的国际直线对撞机 (ILC，internationallinearcollider)的物理 目标对中心径迹探测器的动量分辨率 的要求达到6(1／p)=9×10 (GeV ·C )～．如果采用微单元气体探测器 (MPGD，micropatterngaseous detector)作为电子倍增器件，TPC的位置分辨率有望达到 100p,m．并且由于具有较长的辐射长度、对径迹测 量的采样点多，以及有很好的探测效率等优势，TPC成为ILC中心径迹探测器的重要候选者 ．在 目前的研 究中，最重要的任务是用 TPC原型机 (prototype)实际验证位置分辨率是否能满足 ILC的要求．而在通过 TPC 的实验数据分析位置分辨率时，前提条件是需要重建程序能够从原始数据正确而有效地重建出径迹。本文针 对 目前的TPC原型机建立以Kalmanfilter为基础的径迹重建算法． 1 TPC的径迹重建算法 在 TPC的径迹重建中，首先需要将读出盘上的脉冲波形重建为击 中(hit)．击中可 以认为是某一读 出行 所记录到的粒子在 TPC工作气体中产生电离的平均位置．一个击中由在某一读 出行 内相邻读出盘上的脉冲 组成，同时各脉冲也要求在时间上相邻地到达 TPC读出面．设脉冲 i的电荷量 Q=∑q(其 中q为脉冲的波 形采样值)，读出盘的位置为 ，击 中 方 向坐标由重心法 = ∑(Q；)／zQ 确定．击中的Y方 向坐标为该方向读出盘的中心位置．和 方 向一样 ，击中的Z-坐标 由脉冲到达TPC的时间 (通过电子漂移速度得到对应的漂移距离)和电荷量由重心法计算得到． 收稿 日期 ：2012—03—19；修回 日期 ：2012—08—03 基金项 目：国家 自然科学基金 和清华大学 自主科研计划 (2009THZ02231)资助项 目 作者简介：李波 (1981一)，男，博士，从事高能物理探测器数据重建算法和相关数据分析等研究，E-mail：libei07@mails．tsinghua．edu．cn 128 计 算 物 理 第 3O卷 在重建出各读出行上所记录的击 中之后 ，接下来根据这些击中进行径迹重建