《粒子物理与核物理实验方法（二）》数据处理与物理分析.pptVIP

数据处理与物理分析 何康林 hekl@ihep.ac.cn 内容提要 数据处理与物理分析流程的简单介绍 探测器刻度 事例重建:径迹重建和事例顶点重建 粒子鉴别技术 衰变顶点重建和运动学拟合 物理分析 包括一个分析实例 误差分析 粒子物理实验:有哪些类型的数据 事例数据 记带 探测器获取的事例数据 通常意义上的”数据处理”,指处理事例数据 探测器状态数据 数据库 探测器运行环境数据:温度,气压,等等 slow control 探测器刻度数据,几何定位校准数据 calibration 描述性数据 XML文件 运行配置文件 取数条件,电子学读出配置等 configuration 探测器物质,几何描述等 service 数据处理流程 数据过滤,标记,存储 探测器刻度,探测器几何定位的精确校准 事例的重建 事例模拟 物理分析 数据简并 内容进行压缩 丢弃不需要的事例 分步处理后的数据,内容越来越”物理” 有些步骤中,数据要反复处理,直到满足物理要求 BES3数据处理的挑战 非常高的物理事例率 ~3000Hz@J/y 数据量大,尽管单一事例长度不长 B-工厂,LHC,事例率~100Hz 本底高 高亮度,束流相关本底 物理事例,但是非人们感兴趣的”信号”事例 分布式数据处理与分析 国际合作组,物理学家来自不同国家,不同大学和研究机构 同一物理学家参加不同物理课题,同一物理课题有不同物理学家加入 BES3上主要的物理过程 内容十分丰富 数据过滤 初识碰撞率:~125M Hz 事例记录率:~3000 Hz 实现方法 多级触发: YES/NO 非常快的一级触发:硬件 可编程技术 利用各个子探测器的模拟信号,快速处理 慢的二级触发:软件方法 快CPU 利用各个子探测器的数字信号,完整事例分析 算法可靠性:丢失的事例永远找不回来!!! High Level Trigger 在线事例筛选 Hardware trigger L1 125M Hz?~4000 2000 Hz Software trigger HLT ,进一步压缩本底 4000 2000 Hz?3000 1000 Hz 事例标记与分类: 亮度测量: Bhabha ee , di-muon mm , di-photon gg 物理事例: 强子, t对,… 本底事例: 宇宙线,beam gas, … “噪声”事例: random beam gate 算法复杂度介于一级触发判选与事例重建算法之间,在速度上进行了优化 判选依据:各种事例的”形态” shape 分布 共面角分布 飞行时间差 CPU时间 每个事例平均耗时~10ms Intel? Xeon? CPU 2.8GHz 对4000Hz的事例率,在线计算机机群CPU节点数 40 探测器刻度与事例重建 通过HLT的原始数据将记录到磁带或磁盘上，并很快进行事例重建 原始数据记录的是ADC和TDC计数，事例重建就是把这些计数“变成”粒子的能量，动量，质量，位置等物理量 事例重建之前，这些计数首先要转换成能量，位置和时间等物理量，转换方式与探测器相关 每个探测器都需要进行刻度，刻度常数通常存储在数据库，并随时可能更新 各个探测器的刻度完成之后，才可以进行完整事例的重建 数据产生与存储 事例重建生成物理对象,内容包括 径迹:带电径迹和中性径迹 顶点:初级顶点和次级顶点 粒子鉴别信息 事例标记信息 物理对象依赖于 探测器刻度,几何定位 重建程序版本,以及算法所使用的参数 Monte Carlo数据也需要进行事例重建 以上数据都需要永久保存 记盘,记带,或者记入数据库 ,供物理分析人员使用 物理分析 事例选择与分析:多维变量的分析与选择 信号提取,研究本底,确定信号的显著性 必要的修正 利用模拟数据,计算接收度,重建效率,… 数据与模拟不一致时,改进分析,并做必要的效率修正 物理量的计算 质量,宽度,寿命,物理截面,分支比,自旋宇称,… 相应的误差分析 统计误差:样本的有限性造成的统计涨落 系统误差:由于对物理规律和探测器性能认识的局限性,所做的一些关键的假设.由此造成的不确定性,称为”系统误差”.包括:探测器模型,理论模型,本底研究,探测器分辨能力,等等 常见的分析方法 简单分析手段 多维变量的选择与分析 事例选择,交互式处理平台 直方图累积与显示 曲线与模型的拟合 实验结果的估计,置信水平与区间,上下限的确定 实验误差处理与合并 常用分析工具 粒子鉴别 衰变顶点重建,运动学拟合 等等 较为复杂的分析手段 Dalitz图分析 分波分析方法 等等 公用物理量的系统误差分析 径迹重建效率 光子探测效率 粒子鉴别效率 等等 探测器刻度 量能器刻度 漂移室脉冲幅度刻度 飞行时间计数器的刻度 漂移室的时间刻度 探测器几何位置的精确校