羊八井ARGO试验簇射芯位的重建-McGillphysics.ppt

羊八井ARGO实验簇射芯位的重建 刘 传 磊 物理与微电子学院 Centre Tree Mtree Dtree Wavelets 事例数随方位角的分布 质子事例数随方位角的分布 事例数随天顶角的分布 考虑到阵列的边缘效应 对于质子簇射事例 更加理想的重建方法 计算出中心地毯区的每个pad当作 r0时的 f(r。) Max(f(r0)) 簇射的芯位 把6×5个pad作为一组，用每一组pad中最左下角的pad作为r0 来计算。找出最大的pad组后，再在这一组中寻找 的最大值，则对应最大值的pad的中心为簇射的芯位。 --Wavelet algorithm ?芯位（假定芯区粒子的密度最高） 第四部分 对芯位重建方法的检验 对芯位重建方法的检验 大气簇射模拟 原初粒子 Gamma/Proton 能量(TeV) 0.1-10 模拟事例 100000 天顶角 (度) 7~45 能谱指数 －1.49 强相互作用模型 QGSJET 方位角（度） 0~360 强子能量cutoff(GeV) 0.3 Muon能量cutoff(GeV) 0.3 电子能量cutoff(GeV) 0.015 光子能量cutoff(GeV) 0.015 观测平面高度（m） 4300 对芯位重建方法的检验 ARGO探测器模拟 对芯位重建方法的检验 模拟后得到了27713和14073个 事例和质子事例 用8000个 事例和14073个质子事例来进行检验 对芯位重建方法的检验 事例的能量随击中PAD数的变化（ 事例） 对芯位重建方法的检验 事例的能量随击中PAD数的变化（质子事例） 对芯位重建方法的检验 各种重建方法所需CPU时间( 事例重建 ) 28’ Wavelet (a=5m) 28’ Wavelet (a=25m) 31’02’’ Dtree 27’42’’ Wavelet (a=20m) 31’04’’ Mtree 28’ Wavelet (a=15m) 32‘ Tree 28’ Wavelet (a=10m) 8’00“ centre 对芯位重建方法的检验 各种重建方法所需CPU时间(质子事例重建 ) 39’40“ Wavelet (a=5m) 42’ Wavelet (a=25m) 40‘ Dtree 43‘’’ Wavelet (a=20m) 40‘ Mtree 39’51“ Wavelet (a=15m) 38‘ Tree 39’53“ Wavelet (a=10m) 12‘ centre 0~30 ~60 ~90 ~120 ~150 ~180 ~210 ~240 ~270 ~300 ~330 ~360 Centre 705 657 645 661 667 649 646 659 648 675 671 706 Tree 708 653 650 656 670 646 647 660 647 679 671 702 Mtree 714 654 646 653 672 646 648 657 645 683 668 703 Dtree 707 653 648 656 674 647 647 658 643 682 670 704 Wavelet(a=10m) 718 653 649 644 675 647 649 664 640 686 664 700 对芯位重建方法的检验 对芯位重建方法的检验 0~30 ~60 ~90 ~120 ~150 ~180 ~210 ~240 ~270 ~300 ~330 ~360 Centre 1232 1172 1111 1246 1139 1157 1166 1170 1140 1186 1155 1191 Tree 1225 1174 1117 1242 1142 1158 1161 1171 1135 1189 1156 1195 Mtree 1212 1181 1110 1251 1136 1161 1165 1165 1132 1187 1164 1201 Dtree 1230 1177 1110 1249 1139 1157 1160 1173 1141 1187 1151 1191 Wavelet(a=20m) 1217 1169 1109 1255 1138 1152 1172 1163 1136 1182 1171 1201 对芯位重建方法的检验 事例的芯位重建精度随方位