Transcript BESIII MDC 事例重建与刻度

BESIII MDC 事例重建
• 径迹寻找
• 径迹拟合
• 时间刻度
• dE/dx粒子鉴别
• 总结
毛泽普
2010-4-19
于南昌
一
径迹寻找
BESIII MDC 结构
• 43丝层(24S +19A)
• 6796 信号丝
• XY = 130 m,
• Pt/Pt = 0.5 %(@1 GeV)
• dE/dx 分辨 0.6-0.7%
MDC 事例重建系统
Hit wire ID
MDC
TDC(时间)
Event Rec.
B(磁场)
Event –Test
Tracking
Kalman fit
径迹属性:
空间位置
动量
电荷
粒子种类
…….
dE/dx PID
事例起始时间计算
我们开发了五种独立的方法共同完
成事例起始时间计算:
1 MDC径迹TOF匹配法
2 EMC & TOF匹配法
3 MDC径迹段直线拟合法
4 MDC径迹法
5 MDC径迹拟合方法
• 效率：
– Bhabha，dimu：>99.8%
– Hadron: 99.6%
– Cosmic：99.9%
– 误判率~1%
• 时间分辨: 0.3ns~0.4ns
MDC 径迹快重建
Fast, simple tracking for Test
xy100um
z4mm
xy 100um
p/p56MeV
z  4mm
p  56MeV
Eff=99%
MDC 径迹寻找
• 在均匀z向磁场中(BX = BY = 0, Bz=C)
带电粒子运动轨迹(圆柱螺旋线Z)描述：
X(s) = x0 + R[cos(0+hscos/R) -cos0]
Y(s) = y0 + R[sin(0+hscos/R)-sin0]
Z(s) = z0 + ssin
x-y 平面投影: 园 (d2X/ds2=cos2 /Rh )
• 在一个特定的参考系中径迹参数为:
= (d , , , dz, tan )T , d : signed distance of helix from pivot in x-y plane,  :
The azimuthal angle to the helix center, : 1/pt, dz: signed z distance of the helix
from pinot in the z direction, tan: the slope of the track,  : dip angle
• 拟合方法 (Least Square Method)
2 =  (i/i)2, i=1, nhits ; i =d(xi(i) – di
E ={(1/2) (∂ 2 2 / (∂T∂ ) }-1
•
径迹属性:
Pt = 1/||,
Py = (1/||) cos(0 + ),
Px = (1/|| ) (-sin(0 + ),
Pz =( 1/||) tan 
(i = dfi - dmi )
MDC 径迹寻找
我们开发了两套独立的径迹寻找程序: MdcPatRec 和 MdcTsfRec
两套程序均已用于MDC的数据处理中, 其原理如下:
MDC 径迹寻找程序基本性能
•
当前基本性能: (from J/ψ Run 9947 Bhabha events )
Tracking efficiency : Barrel: 98.6%
CPU for MC data: about 10 ms/track
Efficiency vs P
pt/pt vs P
Efficiency vs  angle
二 MDC径迹拟合(Kalman-Filter方法)
• Kalman-Filter 基本原理：
对离散数据，用当前状态矢量预测下一个状态矢量的LSM方法
•
径迹拟合的目的：
径迹精细修正 NUMF, Multiple scattering, Energy loss …
• 径迹拟合的基本过程和原理：
1.预测：用径迹当前的状态矢量预测第k时候的状态矢量
2.过滤：加权组合第k时的预测信息和测量信息,估计k状态矢量信息
3.平滑：用全部时刻n(n>k)的测量信息回推, 估计k时刻的状态
MDC 径迹拟合
• 我们开发了五个不同用处的程序块，现已经用于物理分析中
1.分析流程：filter单向由外往里, 输出IP点参数，用于物理分析
2.刻度流程：双向filter迭加, 输出每点径迹参数
3.平滑流程：输出每个击中层以及最外点径迹参数、每小段的飞行时间总和，用
于dE/dx刻度，外推, TOF刻度等
4.宇宙线校准流程：用于校准流程，将宇宙线进行连接并进行拟合
5.次级顶点重建工具：ExtToSecondVertexTool按照用户指定的位置拟合.
• MDC 径迹拟合基本性能 by 
0.3GeV
cos=0.83
三 MDC 时间刻度
• 时间刻度的任务与方法(反复叠带，逐渐逼
近)
* 要达到空间分辨130μm, 动量分辨 0.5%@1GeV/c
* MDC 信号道数(信号数): 6796
* X-T关系刻度
* T0刻度
* 时幅关系刻度
•
* 几何位置校
• 数据样本：通常用Bhabha或dimu事例
• 刻度理论模型:
N hit
 
2
i 1
径迹残差法
(i )
(i )
2
(d meas
 dtrac
k)
 i2
Dmeas：径迹与信号丝间的测量距离
Dtrack: 拟合径迹与信号丝的距离(拟合距离）
σ i： 该测量点的权重(空间分辨）
刻度流程
MDC 时间刻度
• X-T关系刻度原因
小单元漂移室中电子漂移特性
– 单元内电场分布不均匀导致了电子漂移速度的
非均匀
– 信号丝附近，电场较强，漂移速度较大
– 随着漂移距离的增大，电场逐渐减弱，漂移速
度也逐渐减小
– 单元边界，漂移线严重弯曲，因而随着漂移距
离的增大漂移时间迅速增大
磁场下电子漂移线等时线图
• X-T关系刻度方法
– 时间谱积分法
– Δd-T迭代法
•
X-T 关系函数
采用5 阶多项
式
＋1阶多项式
漂移时间分布
漂移时间与漂移距离的关系
MDC 时间刻度
• T0刻度方法和原理
*
对时间分布前沿进行拟合，得到T0的初始值
* 径迹残差法对各丝层、单丝修正T0，反复
叠带，逐渐逼近
•
T0 的初步确定
时幅关系刻度
如果用二次函数描述信号上升沿：
t2  t1 
Vth  T  (
1

VT 1
1
)
VT 2
得出时间游动： t  p0 
p1
Q
实际情况中只能测量信号的电荷量Q，只能假
设 Vt ∞Q
，使用如下函数来修正时间
c1
游动： t  c0 
VT
信号脉冲幅度引起时间游动
MDC 几何位置校准
•
•
单丝位置刻度: 残差分布方法
漂移室端面板校准参数
MDC 时间测量主要性能现状
• 动量分辨:11.4 MeV/c
•
空间分辨: 134 µm
•
正负电贺对称性<5MeV/c
•
径迹空间分布正常
•
径迹质量稳定性好
from J/psi Bhabha
from J/psi Bhabha
P = 11.4MeV/c
from J/psi Bhabha
Spatial resolution vs run (Bhabha)
P vs phi
e+/e- momentum vs run
四 dE/dx 粒子鉴别
• dE/dx 粒子鉴别原理和方法
• dE/dx 刻度
修正dE/dx值的不一致性: 带电粒子的电离特性;信号丝气体放大的不均匀性; MDC磁场的
不 均匀性;取数过程环境温度、压强等条件的变化等
dE/dx 刻度项目
• 单丝级别刻度（只用电子样本刻度）
* 径迹长度修正
* run by run 修正外界环境影响，如气压、温度等
* 单丝增益：修正单元电场非均性、电子学增益差异
* 漂移距离和入射角联合修正
• 径迹级别刻度 (用各种粒子样本做刻度):
* 空间电荷效应的修正
* dE/dx 能损曲线刻度
* σdE/dx 的刻度
拟合函数σdE/dx=f(βγ)*g(sinθ)*h(Nhit)*I
使用bhabha数据样本得到拟合函数: g, h 和i :
利用强子样本得到 ( 与dE/dx能损曲线相似) ：f
MDC dE/dx测量主要性能现状
σdE/dx~6%(for π)
3σ π/ĸ 分离~760MeV/c

K
P
总结
• 在3年左右时间我们完成了MDC数据处理程序的（事例
起始时间计算，径迹寻找，径迹拟合，离线刻度，dE/dx
粒子鉴别）五个系统的设计、编程和调试.
• 经过BESIII 实验数据的调试和运行，证明了五个
系统设计达到了实际指标，基本满足物理分析的要求
。
• 更加精细的调整和参数的优化工作还将继续进行
谢谢各位
BESIII 探测器
MDC 时间测量数值—TDC
不是漂移时间
– 通过TDC测量到的原始时间并不等同于漂移时间，其中包含：
• TES: 事例起始时间 (EsTimeAlg)
• Tflight：粒子的飞行时间（重建中计算并修正）
• Tdrift：漂移时间
• Tprop：信号在丝上的传播时间（重建中计算并修正）
• Twalk：信号幅度差异引起的时间游动 → 时幅关系刻度
• Telec：信号在电子学通道上的传输时间 → T0刻度
TES
Tflight
Tdrift
TTDC
Tprop Twalk
Telec
北京谱仪
数
据在线事例选择
获
取
原始数据
BESIII数据处理流程
重建数据(DST)
离线刻度
calibration
事例重建
reconstruction
mc
MC Data 数
generator据
simulation产
生
重建数据(Rec)
数
据
处
理
物
数据分类
理
分
析
BESIII事例重建过程顺序
Event Data Flow
MC 数据
Generator
HepEvt
BES III
RawData
实验数据
Histograms
Ntuples
物理分析
Event
Converters
G4Event
Simulation
RawData
CnvSvc
Digits
Digitization
Calibration
Algorithms
RecHits
Reconstruction
Algorithms
DstTracks
Rec2DstAlg
Analysis
Tools
DstData
RootDstCnvSvc
McTruth
Hits
RecTracks
Calibration Framework
• Framework is based on GLAST ‘s scheme

The calibration constants for each sub-detector are
produced by the associated calibration algorithm

The framework provides reconstruction algorithms a
standard way to obtain the calibration data objects
MDC Tracking Module(1) R-  Tracking & S-Z finding
1) Conformal transformation
2) Segment Finding (TSF)
3) S-Z calculation in Z finding
MDC Tracking Module(2): R-  Tracking & S-Z finding
1）Segment finding: Search segments in each
super-layer using a pattern look-up table
7
6
4
a) Create Hit pattern by MC:
* 4 hit pattern
* 3 hit pattern
add stereo segments
3
2
1
b) Fill real hit pattern table and matching
Set one word for a group of 8 wire, each
bit for a wire. Set “1” for a hit wire,
others “0”
c) Link segments to 2D tracks
2）Track link (S-Z finding)
5
0
1 0 0 1 0 1 0 1
7 6 5 4 3 2 1 0
Wires No. 0-7
• R - 平面径迹园拟合(非迭代拟合),
S-Z平面直线拟合(非迭代拟合)
• 空间螺旋线拟合(迭代拟合)
dE/dx 能损曲线刻度
• 数据样本:
– 1.electron:
• bhabha:1.81.95GeV
• radee:0.5-1.8GeV
– 2.muon:
• dimu:1.8-1.94GeV
• cosmic:0.6-10GeV
– 3.pion:gama4pi:0.51.7GeV
– 4:kaon:gama4k:0.51.0GeV
– 5.pronton:0.3-0.9Ge