BESIII EMC 事例重建与刻度知识培训.ppt

*求?2最小，可以给出如下的矩阵方程：Q是6240阶稀疏矩阵如果积累足够统计量的Bhabha事例，我们可以通过求矩阵反演得到所有晶体的刻度常数用稀疏矩阵软件包(SLAP)求解矩阵方程.Ee(?,?)：电子或正电子的运动学能量f(Ee,?,?)：沉积在EMC中的能量比分Ekexp：预期能量,?(?,?)：能量分辨i—晶体编号k－shower编号gi－i块晶体的刻度常数目前f?由模拟数据给出构造?2：Bhabha刻度原理其中e+e-11mrad*1.843GeV位置分辨(boss651)BarrelEastEndcapWestEndcap?MC4.6mm5.6mm5.6mmData4.5mm6.0mm5.9mm?MC4.7mm6.8mm6.8mmData5.1mm7.6mm7.2mm1.843GeV桶部e5X5能量分辨2.3%boss6511.843GeV端盖e5X5能量分辨4.2%boss6513.7-3.27Psip4.2-4.14Psip5.25-6.23.65GeV6.7-7.2JpsiMachinestudy晶体刻度常数平均值（相对于3月7日至11日的刻度常数）随时间的变化。Bhabha刻度结果：能量分辨：桶部2.3%@1.843GeV位置分辨：桶部～5mm@1.843GeV模拟和数据符合的很好晶体刻度常数随时间的变化反映了晶体辐照损伤情况对撞取数阶段辐照较小加速器机器研究时辐照程度大一些*总结BESIII晶体量能器运行状态很好，性能稳定。BESIIIEMC已经达到了设计/验收指标蒙特卡罗模拟和数据符合很好能量分辨2.5%@1GeV位置分辨6mm@1GeVEMC重建刻度软件已经成功的处理了BESIII采集的所有数据*谢谢各位*接下来，我给大家讲BESIIICsI晶体量能器的刻度与重建*主要内容包括，BESIII电磁量能器简介，电磁量能器重建电磁量能器能量刻度 单块晶体能量刻度 簇射能量修正****Emc刻度流程图从数据库中读取单块晶体的能量刻度常数，将ADC转换成沉积能量。单块晶体的初始刻度常数来自于实验室的相对光产额测量和无束流宇宙线刻度。在簇团寻找过程中，重建簇射能量和位置。利用对撞物理事例，Bhabha事例对单块晶体能量做精确刻度。然后利用已知光子能量的物理过程对簇射能量做修正。可以用pi0事例刻度、单光子辐射Bhabha事例、，在高能区还可以用e+e-???，在低能区，用?????c(?J/?)。最后重建出粒子能量由于辐照损伤等原因，晶体的光产额和光输出不均匀性随时间发生变化，因此需要经常利用bhabha事例做晶体能量刻度。*量能器重建算法的流程，首先通过刻度服务，将晶体的ADC信息转换为真实的沉积能量；再根据晶体编号，用迭代法寻找能量簇团；然后在能量簇团中寻找种子，并且根据种子的数目和簇团的形状将簇团劈裂为一个或多个簇射；最后计算簇射的能量和位置并进行修正。*在能量簇团寻找的算法中，定义相邻的晶体是非常重要的，并且跟几何结构关系密切。对于桶部量能器，每一块晶体有八个邻居；而对于端盖量能器，由于每层晶体数目不完全相同，晶体的邻居数也不一样，这是端盖一个扇区的晶体排列，其中带○的晶体就有七个邻居。*把簇射中包含的所有晶体能量求和，就得到了簇射的能量。不过参与求和的晶体数越多，包含进来的噪声也越多，这样会降低能量分辨。由于电磁簇射能量沉积比较集中，通常表示为种子周围9块晶体能量求和；或25块晶体能量求和。具体用还是要根据噪声的水平来选取。如果该簇射是由一个能量簇团劈裂而来的话，参与求和的晶体能量还需要乘以对应的权重。***Becauseofleakageandfrontmaterial,theshowerenergydon’tequaltotheincidentenergyofphoton.Soacorrectionisneeded.InMC,thisprocedureisstraightforwardbecauseweknowtheincidentenergyofeachsinglephoton.InData,below1GeVpi0calibrationisusedandradiativebhabhacancalibrateenergygreaterthan0.