BesIII软件系统强化课程

1、1BESIII EMC 事例重建与刻度刘春秀 （高能所实验物理中心）2010-04-19中国高能物理大会-南昌2主要内容主要内容 BESIII CsIBESIII CsI晶体量能器晶体量能器(EMC)(EMC)简介简介 EMCEMC事例重建事例重建簇射能量和位置重建簇射能量和位置重建簇射位置和能量修正簇射位置和能量修正TOFTOF能量重建能量重建, ,它是它是EMCEMC重建算法的补充重建算法的补充利用利用EMCEMC时间信息去除束流本底时间信息去除束流本底 EMC CsIEMC CsI晶体单元绝对能量刻度晶体单元绝对能量刻度 总结总结3BESCsI(Tl)BESCsI(Tl)晶体量能器的功能

2、与指标晶体量能器的功能与指标1. 测量光子的能量覆盖范围为：20 MeV2 GeV，重点能区低于500 MeV。2. 能量分辨: （BESII 201GeV)3. 位置分辨： （BESII 3cm1GeV） 4. 提供中性（）能量触发。5. 在能量大于200MeV的区域具有良好的e/分辨，和具有良好的双角分辩能力。6晶体探测单元的电子学噪声水平小于220KeV。E(GeV)2.5%/E/E)(/6,GeVEmmYX 电磁量能器电磁量能器(EMC)(EMC)在在BESBES中占有十分重要的地位，基本功中占有十分重要的地位，基本功能是测量能是测量光子和电子的能量和位置信息及穿过量能器的其光子和电子

3、的能量和位置信息及穿过量能器的其它带电粒子的沉积能量。它带电粒子的沉积能量。4BESCsI(Tl)BESCsI(Tl)晶体量能器的几何结构晶体量能器的几何结构Barrel: 120 x44=5280 w:21.6t Endcaps: 2x(96,96,80,80,64,64) =960 w: 4tTotal: 6240 w: 25.6 tBESIII EMC特有的支撑结构：桶部和端盖晶体均以吊挂方式固定在支撑架上 5CsICsI晶体量能器事例重建晶体量能器事例重建6CsI(Tl)晶体量能器事例重建晶体量能器事例重建重建原理: 电磁簇射在晶体中发展，把能量损失在一系列相邻的晶体中 寻找相邻的晶体

4、，进而计算出入射粒子的总能量和击中位置7EMC 事例重建流程事例重建流程Particles Energy(E)Raw dataADC (Elec. Calib. Applied)conversionabsolute constant Energy(Ei)clusteringBhabha Calib.DataBaseShower Energy(E)Correction 0 Calib.Rad. Bha. Calib.TOF EnergyShower energy correctionore+e- Calib.at high energy c ( J/ )at low energy初始刻度常数来自

5、相初始刻度常数来自相对光产额测量和无束对光产额测量和无束流宇宙线刻度流宇宙线刻度8重建算法流程重建算法流程晶体能量刻度能量簇团寻找能量簇团劈裂digithitclustershower9簇射寻找簇射寻找 能量簇团(Cluster)寻找: 一系列能量沉积大于某个阈值的相邻晶体的集合。 寻找种子：沉积能量的极大值种子数目=1, cluster=1个shower，对应一个入射粒子种子数目=m, m1, cluster=m个shower，对应多个入射粒子通常出现在高动量0形成的簇团或强相互作用形成的簇团中 能量簇团劈裂多个簇射1GeV光子形成的能量簇团晶体1.5GeV0形成的能量簇团10簇射能量簇射能

6、量 把簇射中包含的所有晶体能量求和，就得到了簇射的能量EAll 。不过参与求和的晶体数越多，包含进来的噪声也越多，这样会降低能量分辨。 由于电磁簇射能量沉积比较集中，通常表示为种子周围9块晶体能量求和E33 ；或25块晶体能量求和E55。 具体用E33还是E55要根据噪声的水平来选取。 如果该簇射是由一个能量簇团劈裂而来的，参与求和的晶体能量还需要乘以对应的权重。NiishowerEE111簇射位置簇射位置NNcjjjjjjjxwExwE()linearjjjwEE()00,lnlnlogjjjtotwEMaxaEE 加权平均 线性权重 对数权重 位置修正修正粒子打到晶体表面不同位置带来的偏差

7、用bhabha事例，根据漂移室的径迹外推信息进行修正修正公式：01234arctan()xppxpppx 线性权重对数权重12（a）线性权重，修正前；（b）线性权重，修正后；（c）对数权重，修正前；（d）对数权重，修正后。 13簇射能量修正簇射能量修正 原因: 量能器前物质能损晶体侧边、尾端的能量泄漏 簇射能量 入射粒子能量 目的：沉积能量修正到粒子的真实能量 方法：选用已知光子能量的物理过程MC: 单光子数据: 能量范围 40MeV-1.0GeV 0 （简称0刻度）数据: 能量范围 0.5GeV-1.5GeV e+e- e+e-（简称辐射Bhabha刻度）高能区e+e- 双光子事例14具体操

8、作：具体操作：首先用首先用MCMC单光子刻度簇射能量单光子刻度簇射能量( ( , ,E)E)依赖依赖, , 然后用然后用高能区高能区e+e-数据数据刻度簇射能量刻度簇射能量( )依赖，依赖， 最后用最后用 0 数据刻度不同能区簇射能量。数据刻度不同能区簇射能量。 利用利用cJ/ （）和）和 e+e- e+e- 过程检查光子重建性能过程检查光子重建性能能量分辨达到能量分辨达到2.4%1GeV （对于光子）。（对于光子）。 0 signal in inclusive DATA/MCDATA/MCEMC桶部的能量分辨桶部的能量分辨簇射能量修正簇射能量修正15电磁量能器前物质对其性能的影响电磁量能器前

9、物质对其性能的影响 EMC前有30% 辐射长度的物质，其中TOF占的比重最大 TOF上的能损延伸到几百个MeV，这将使能量分辨和探测效率变差Energy loss in TOF for 1GeV photons16方法：用方法：用e e+ +e e- - + + - - 事例刻度单根闪烁体能量和击中位置，事例刻度单根闪烁体能量和击中位置，将重建将重建TOFTOF沉积能量与量能器簇射匹配，并将该能量与沉积能量与量能器簇射匹配，并将该能量与EMCEMC簇簇射能量相加，簇射能量修正后射能量相加，簇射能量修正后，得到重建光子能量。，得到重建光子能量。结果：结果：对高能双光子事例能量分辨提高对高能双光子

10、事例能量分辨提高8%8%，低能量低能量(0.8GeV)(0.8GeV)光子探测效率提高光子探测效率提高6%, 6%, 0 0重建效率提高重建效率提高12%12%。 MCMC和数据符合的很好和数据符合的很好TOFTOF能量重建是能量重建是EMCEMC重建算法的补充，提高了低能光子的探重建算法的补充，提高了低能光子的探测效率，改善了高能光子的能量分辨。测效率，改善了高能光子的能量分辨。TOFTOF沉积能量重建并与沉积能量重建并与EMCEMC匹配（匹配（BESIIIBESIII首创）首创）17TOF能量重建17EMCTOFouterinnerIPCalibrate TOF energy using

11、di-mu events; add the TOF energy into the EMC shower energy to reconstruct photon energy18e+e- e+e- 和和 cJ/ （）过程中过程中光子探测效率Solid-Without TOF, circle-With TOFPhoton efficiency increased significantly when E0.8GeV19 0 0 0 0J/ 中中 0 0的探测效率的探测效率 circle: without TOF energydot: with TOF energycircle: without

12、 TOF energydot: with TOF energyMC efficiencyDATA efficiencyMC efficiency improvementDATA efficiency improvement12%12%M 0.12-0.145GeV)0 efficiency increase about 12% in low energy range 20BESIIIBESIII与其它国际上与其它国际上CSICSI晶体探测器不同，它除了有幅度晶体探测器不同，它除了有幅度信息信息ADCADC输出，还有时间信息输出，还有时间信息TDCTDC输出。利用输出。利用EMCEMC时间信息时

13、间信息可以部分去除随机产生的束流本底。可以部分去除随机产生的束流本底。方法方法: : 用用 0 0 J/ ， J/ l l+ +l l- - 过程选择低能光子样本。过程选择低能光子样本。比较好光子与假光子（来自随机触发本底）时间信息，比较好光子与假光子（来自随机触发本底）时间信息，即比较真、假光子的即比较真、假光子的EMC时间减去事例起始时间。时间减去事例起始时间。结果：结果： 低能下低能下(0.3GeV)束流相关本底排除率达到束流相关本底排除率达到75%。利用利用EMC时间信息去除束流本底（时间信息去除束流本底（BESIII首创）首创）21 物理光子的时间信息分布物理光子的时间信息分布 (

14、0 0J/ ，J/ ee, ) eEMC时间减事例起始时间的分布时间减事例起始时间的分布束流本底引起的假光子的时间束流本底引起的假光子的时间EMC时间减事例起始时间随能量的变化时间减事例起始时间随能量的变化事例起始时间事例起始时间 EMC时间伴随不同的触发有不同分布时间伴随不同的触发有不同分布 22利用利用EMC时间减事例起始时间排除束流本底前后的假光子时间减事例起始时间排除束流本底前后的假光子( 0 0J/ )BarrelEndcap利用利用EMC时间减事例起始时间可以提高低能量下假光子的排除率时间减事例起始时间可以提高低能量下假光子的排除率光子选择条件可设为光子选择条件可设为10几几MeV

15、True photonFake photon before Time cutFake photon after Time cutDATADATA2324将ADC道数转换成能量 其中 ADCi 经过了减台阶和电子学刻度 ，ci 是CsI晶体的绝对能量刻度常数 晶体探测单元安装之前，没有做绝对能量刻度 实验室只给出了晶体相对光产额的测量 利用无束流的宇宙线数据做首次的刻度iiicADCE25Cosmic ray宇宙线数据刻度原理宇宙线数据刻度原理l 事例筛选事例筛选 EMC背对背 能量大于10MeV的CsI晶体只有两块l 由于由于BESIII探测器获取的宇宙线主要探测器获取的宇宙线主要是是 子，它

16、的动量基本上大于子，它的动量基本上大于1GeV， 在在CsI晶体中的能量沉积可以认为与它晶体中的能量沉积可以认为与它在晶体中穿过的长度成正比在晶体中穿过的长度成正比l 对每块晶体检查沉积能量分布对每块晶体检查沉积能量分布 首先用相对光产额数据刻度首先用相对光产额数据刻度刻度常数刻度常数 然后与蒙特卡罗数据比较然后与蒙特卡罗数据比较 CiMC=能量平均值能量平均值(MC)/能量平均值（能量平均值（data）l 对每块晶体宇宙线刻度常数为对每块晶体宇宙线刻度常数为 Cicosmic= Ci光产额光产额 x CiMC晶晶体体光光产产额额光光产产额额平平均均值值光光产产额额iCi dE/dx of c

17、osmic ray in CsI crystal26Bhabha事例刻度事例刻度 无束流宇宙线刻度只能给出一个比较粗的初始刻度，更精确的刻度必须用对撞的物理事例完成。 晶体有不同的光产额和不同的光输出不均匀性，由于辐照损伤等原因，光产额和光输出不均匀性随时间发生变化，因此必须经常对单块晶体做能量刻度。采用采用Bhabha事例做刻度，得到更精确的晶体单元绝对事例做刻度，得到更精确的晶体单元绝对能量刻度常数。能量刻度常数。27 NkikiikEgE12exp2),( NkkkikiNkkjkikijijijjEEREEQRQg12exp12 , ,其中求2 最小，可以给出如下的矩阵方程：Q 是62

18、40阶稀疏矩阵如果积累足够统计量的Bhabha事例，我们可以通过求矩阵反演得到所有晶体的刻度常数用稀疏矩阵软件包 (SLAP )求解矩阵方程.Ee(,) ：电子或正电子的运动学能量 f (Ee, ) ：沉积在EMC中的能量比分Ekexp ： 预期能量 , (,) ：能量分辨i 晶体编号 k shower编号 gii块晶体的刻度常数目前f & 由模拟数据给出构造 2 ：Bhabha 刻度原理刻度原理),(),(exp eekEfEE 其中e+e-11mrad281.843GeV1.843GeV位置分辨位置分辨(boss651)(boss651)BarrelEast EndcapWest EndcapMC4.6mm5.6mm5.6mmData4.5mm6.0mm5.9mm MC