清苇大净工程物理系 Department of Engineering Physics,Tsinghua Universlty 实验三.宇宙线缪子寿命测量 前言 本实验的目的是：学习一种测量微观粒子寿命的方法 1.从真实的实验测量中理解微观粒子的寿命。 2.学习如何在数据拟合中扣除偶然符合本底的贡献。 一.实验原理 本实验是宇宙线缪子的寿命测量。首先讲一下实验原理，宇宙线粒子与大气中的原子核碰撞大量次 级粒子，次级粒子产生缪子。缪子到达我们的闪烁体探测器，会在闪烁体中沉积能量，部分低能的 缪子会在闪烁体中沉积完所有能量并停止在闪烁体中，然后发生衰变， u+→et+e+yu h→e十e+u 通过测量初始缪子进入闪烁体的信号和衰变产物的电子在闪烁体内发出的第二个信号之间的时间 差，可以测量缪子的寿命。 在粒子物理标准模型理论中，缪子的衰变几率为 An=1 Gimic2 T4192m3i7 其中T4是缪子的寿命，也就是平均衰变时间，缪子的静止寿命为2.24s，静止衰变时放出的电子的 动能分布为从0MeV至53MeV(平均值为37MeV),衰变时间成指数分布。平均衰变时间目前 最好的测量结果如下所示： T4=(2.19703±0.00004)×106s, 我们采用双端符合的方法来压低热噪声，同时缪子的能量沉积导致的光电子数也要大于热噪声的光 电子的单个计数，通过提高触发阈值也可以大量抑制热噪声。这样，缪子打到闪烁体上能观察到一 个符合信号，衰变出来的电子可以观察到一个延迟符合信号，两个信号的时间差预计在微秒量级。 对于实验室中使用的大约5cm厚度的闪烁体，有1%的缪子能量较低（动能小于10MeV),大约 为1个/分钟或更低，因此应采用体积更大或更多的闪烁体。另外，经估计，如果缪子有1MV的能 量沉积，可以采用较高的阈值，例如100毫伏，进行实验。 二.实验主要内容 1,按照示意图连接好设备，检查信号情况，查看是否有延迟的衰变电子的信号，事例率是否与预 估结果相近。 2.利用示波器测量缪子发信号和延迟电子信号之间的时间差，计算平均时间，做寿命测量的粗 略估计。 3.采集约30分钟数据，利用示波器测量时间，做统计直方图。 4.利用最小二乘法拟合实验数据。 三.需要观察或测量的问题 1.如何初始检查信号的质量？ 2.如何利用最小二乘法去寻找参数的最优值？ 3.为什么可以用常数项作为偶然符合的本底项？
实验三. 宇宙线缪子寿命测量 前言 本实验的目的是：学习一种测量微观粒子寿命的方法 1. 从真实的实验测量中理解微观粒子的寿命。 2. 学习如何在数据拟合中扣除偶然符合本底的贡献。 一. 实验原理 本实验是宇宙线缪子的寿命测量。首先讲一下实验原理，宇宙线粒子与大气中的原子核碰撞大量次 级粒子，次级粒子产生缪子。缪子到达我们的闪烁体探测器，会在闪烁体中沉积能量，部分低能的 缪子会在闪烁体中沉积完所有能量并停止在闪烁体中，然后发生衰变， 通过测量初始缪子进入闪烁体的信号和衰变产物的电子在闪烁体内发出的第二个信号之间的时间 差，可以测量缪子的寿命。 在粒子物理标准模型理论中，缪子的衰变几率为 其中 是缪子的寿命，也就是平均衰变时间，缪子的静止寿命为 ，静止衰变时放出的电子的 动能分布为从 至 （平均值为 ），衰变时间成指数分布。平均衰变时间目前 最好的测量结果如下所示： 我们采用双端符合的方法来压低热噪声，同时缪子的能量沉积导致的光电子数也要大于热噪声的光 电子的单个计数，通过提高触发阈值也可以大量抑制热噪声。 这样，缪子打到闪烁体上能观察到一 个符合信号，衰变出来的电子可以观察到一个延迟符合信号，两个信号的时间差预计在微秒量级。 对于实验室中使用的大约 厚度的闪烁体，有1%的缪子能量较低（动能小于 ），大约 为1个/分钟或更低，因此应采用体积更大或更多的闪烁体。 另外，经估计，如果缪子有1 MeV的能 量沉积，可以采用较高的阈值，例如100毫伏，进行实验。 二. 实验主要内容 1. 按照示意图连接好设备，检查信号情况，查看是否有延迟的衰变电子的信号，事例率是否与预 估结果相近。 2. 利用示波器测量缪子触发信号和延迟电子信号之间的时间差，计算平均时间，做寿命测量的粗 略估计。 3. 采集约30分钟数据，利用示波器测量时间，做统计直方图。 4. 利用最小二乘法拟合实验数据。 三. 需要观察或测量的问题 1. 如何初始检查信号的质量？ 2. 如何利用最小二乘法去寻找参数的最优值？ 3. 为什么可以用常数项作为偶然符合的本底项？ μ + → e+ + ν¯e + νμ， μ − → e− + νe + ν¯μ。 λμ = = , 1 τμ G2 F m2 μ c 2 192π 3ℏ 7 τμ 2.2 μs 0 MeV 53 MeV 37 MeV τμ = (2.19703 ± 0.00004) × 10 s, −6 5 cm 10 MeV
4.如何给出测量结果的误差？ 5.讨论如何改进我们的实验以克服上面的误差 6.考察世界上还有那些实验在测量缪子寿命？ 7.为什么还要继续测量缪子寿命？ 四.实验步骤 本实验的实验装置连接图如图1所示， 闪烁体探测器 闪烁体探测器 甄别器 高压电源 甄别器 扇出单元 逻辑单元 扇出单元 示波器 图1：实验装置连接示意图 经过甄别器的双端信号经过扇出单元分别扇出为两路，一路做符合，一路接到示波器上观察。符合 后的信号也要接到示波器上。 实验步骤如下，连接好设备后，把每格的时间宽度调节为一微秒，把示波器波形显示的余辉时间调 成无穷大。用示波器光标测量缪子触发信号和延迟电子信号之间的时间差。采集30分钟数据，利用 示波器的光标测量时间，做统计直方图。用含常数本底项的指数分布来拟合直方图。 ni=A×exp(-t:/T)+o 其中每个时间bin中观察到的事例数为m,t为该bin的平均时间，A,0,T为待拟合常数。其中 A为信号部分指数分布前的归一化因子，为一常数本底项，描述偶然符合本底，τ为待测缪子寿 命。 五.实验注意事项 这个实验有几处需要注意的地方。同学们可以通过适当调节甄别器阈值来压低余波，做实验的时候 可以把高压加到1700伏，甄别器设置为100毫伏。本次实验使用光标来测量时间，光标测量的误差为 0.02微秒。因为两侧的光标都会贡献误差，所以总的误差为0.03微秒。 V01+号≈0.03s 连好装置之后，在示波器上可以采样到一个缪子触发的符合信号和微秒量级之后紧跟着的一个符合 信号，如图2所示。证明我们观察到了一个缪子衰变的事例
4. 如何给出测量结果的误差？ 5. 讨论如何改进我们的实验以克服上面的误差 6. 考察世界上还有那些实验在测量缪子寿命？ 7. 为什么还要继续测量缪子寿命？ 四. 实验步骤 本实验的实验装置连接图如图1所示， 图1：实验装置连接示意图 经过甄别器的双端信号经过扇出单元分别扇出为两路，一路做符合，一路接到示波器上观察。符合 后的信号也要接到示波器上。 实验步骤如下，连接好设备后，把每格的时间宽度调节为一微秒，把示波器波形显示的余辉时间调 成无穷大。用示波器光标测量缪子触发信号和延迟电子信号之间的时间差。采集30分钟数据，利用 示波器的光标测量时间，做统计直方图。用含常数本底项的指数分布来拟合直方图。 其中每个时间bin中观察到的事例数为 ， 为该bin的平均时间， ， ， 为待拟合常数。其中 为信号部分指数分布前的归一化因子， 为一常数本底项，描述偶然符合本底， 为待测缪子寿 命。 五. 实验注意事项 这个实验有几处需要注意的地方。同学们可以通过适当调节甄别器阈值来压低余波，做实验的时候 可以把高压加到1700伏，甄别器设置为100毫伏。本次实验使用光标来测量时间，光标测量的误差为 0.02微秒。因为两侧的光标都会贡献误差，所以总的误差为0.03微秒。 连好装置之后，在示波器上可以采样到一个缪子触发的符合信号和微秒量级之后紧跟着的一个符合 信号，如图2所示。证明我们观察到了一个缪子衰变的事例。 ni = A × exp(−ti/τ) + y0 mi ti A y0 τ A y0 τ σ + ≈ 0.03 μs 2 1 σ 2 2 −−−−−− √
ek停h M10.0us缩放系数：10Y 噪声滤波器送闭 1■ 1500mw 2500mW 3500rmW 21.00us0.00000s3/-740mW <10H2 找到的搜索事件：2 ●幅值860mY 1●Pise3.892s☐165735 图2：缪子符合信号和缪子衰变信号 下面是用光标测量缪子信号和电子信号时间差的演示视频，通过调节余辉时间为无穷大，显示亮度 为100%，我们可以使得缪子衰变的信号停留在示波器上，看到缪子衰变出的电子信号后马上移动光 标测量时间差。 0:00/0:10 六.实验数据分析 在做实验数据的指数拟合时，有其他的本底存在，例如偶然符合的本底，这些本底的时间结构常数 远大于缪子寿命，可以用常数项代替。 拟合时是通过最小化最小二乘法的X来得到最优的拟合参数。最小二乘法的X可以按如下公式来 构造： X2=∑m:-m)2/m 其中m:对应每个时间区间(bi)的统计量。下图为一典型的时间差分布图及拟合结果(1TDC=1 ns):
图2：缪子符合信号和缪子衰变信号 下面是用光标测量缪子信号和电子信号时间差的演示视频，通过调节余辉时间为无穷大，显示亮度 为100%，我们可以使得缪子衰变的信号停留在示波器上，看到缪子衰变出的电子信号后马上移动光 标测量时间差。 六. 实验数据分析 在做实验数据的指数拟合时，有其他的本底存在，例如偶然符合的本底，这些本底的时间结构常数 远大于缪子寿命，可以用常数项代替。 拟合时是通过最小化最小二乘法的 来得到最优的拟合参数。最小二乘法的 可以按如下公式来 构造： 其中 对应每个时间区间（bin）的统计量。 下图为一典型的时间差分布图及拟合结果（1 TDC=1 ns）： 0:00 / 0:10 χ 2 χ 2 χ = / 2 ∑ i (ni − mi) 2 ni mi
500 400 □4c ExpDecayFit +A"ep(t 100 2doF-1395 R=0.957 A-5243+1878 【=2,10236+ 100 中中中9 5,000 1e+04 15e+04 2e+04 25e+04 3e04 TDC 图3：实验数据拟合结果 从图中我们可以看到，拟合结果中包含指数分布的三个参数以及各自的不确定度，还有X比上自由 度(degrees of freedom)，用来衡量拟合结果的好坏。对于固定自由度而言，X2相比于自由度的差 异越小，说明拟合结果越好。 清第火
图3：实验数据拟合结果 从图中我们可以看到，拟合结果中包含指数分布的三个参数以及各自的不确定度，还有 比上自由 度（degrees of freedom），用来衡量拟合结果的好坏。对于固定自由度而言， 相比于自由度的差 异越小，说明拟合结果越好。 χ 2 χ 2