1 / 54

羊八井 ARGO 实验簇射芯位的重建

羊八井 ARGO 实验簇射芯位的重建. 刘 传 磊 物理与微电子学院. 羊八井 ARGO 实验簇射芯位的重建. Contents. 羊八井 ARGO 实验 离线数据分析软件介绍 不同的芯位重建方法 对芯位重建方法的检验 结论及展望. 第一部分 羊八井 ARGO 实验. 羊八井 ARGO 实验. ARGO 物理目的 γ 天文 弥散 γ 观测 γ 暴物理 反质子丰度测定 宇宙线质子谱 大气簇射结构的详细研究 太阳和日地物理. 羊八井 ARGO 实验.

kadeem-bell
Download Presentation

羊八井 ARGO 实验簇射芯位的重建

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. 羊八井ARGO实验簇射芯位的重建 刘 传 磊 物理与微电子学院

  2. 羊八井ARGO实验簇射芯位的重建 Contents • 羊八井ARGO实验 • 离线数据分析软件介绍 • 不同的芯位重建方法 • 对芯位重建方法的检验 • 结论及展望

  3. 第一部分 羊八井ARGO实验

  4. 羊八井ARGO实验 • ARGO 物理目的 • γ天文 • 弥散γ观测 • γ暴物理 • 反质子丰度测定 • 宇宙线质子谱 • 大气簇射结构的详细研究 • 太阳和日地物理

  5. 羊八井ARGO实验 111.01m 77.87m 98.84m 76.36m

  6. 羊八井ARGO实验

  7. 羊八井ARGO实验 1.25m 2.8m

  8. 羊八井ARGO实验 事例率高、数据量大 实验数据的处 理和存储 medea++

  9. 第二部分 离线数据分析软件介绍

  10. 离线数据分析软件介绍 • 分析软件的整体结构 • 探测器阵列的几何描述 • 事例信息的描述和重建 • 事例的存储 • 物理分析工具

  11. 离线数据分析软件介绍 • 探测器阵列的几何描述 • 采用了“复合-部件”的模式（Composite-component）： • ArgoComp类来描述一个部件标识号码，在阵列中的位置，在一个复合部件中的位置等信息。 • 模板类ArgocompT用于描述一个复合部件和复合部件的子部件的几何信息。

  12. 离线数据分析软件介绍 事例信息的描述和重建

  13. 离线数据分析软件介绍 重建过程事例信息流程 EvDAQ RecLevel0 EvRaw RecLevel1 EvRec1 RecLevel2 EvRec2

  14. 离线数据分析软件介绍 事例信息的存储 物理分析工具

  15. 第三部分 不同的芯位重建方法

  16. 不同的芯位重建方法 • Centre • Tree • Mtree • Dtree • Wavelets

  17. 不同的芯位重建方法 • N= 给定事例的PAD击中数 • 对于某一击中 PAD i (i=1,…,N): • xi , yi≡PAD 的中心坐标 • ti ≡粒子到达时间 • ri(2)≡(xi , yi) • ri(3)≡(xi, yi, zi) （with zi=ti·clight） • ni ≡击中Strip数

  18. 不同的芯位重建方法 1、重心法 rc = 芯位： 击中PAD位置的加权平均

  19. 不同的芯位重建方法 在 MTree、Tree和DTree 三种方法 中，都会用到Tree 长度这一概念。

  20. 不同的芯位重建方法 • 一个二维r-tree示例 r-tree 长度定义如下： l(r) =

  21. 不同的芯位重建方法 2、 Mtree 方法 芯位： N个Tree长度中最小那个 所对应的参照PAD的中心位置

  22. 不同的芯位重建方法 n1=1 l(r1)=3x2+5x5=31 5 3 n3=5 4 n2=2 l(r3)=5x1+4x2=13 l(r2)=3x1+4x5=23

  23. 不同的芯位重建方法 3、Tree 方法 rc= rc= where wi= 芯位： 考虑到Tree长度时击中PAD位置的加权平均

  24. ) 不同的芯位重建方法 4、Dtree 方法 rc= where wi= (1- 芯位： 考虑到Tree长度时击中PAD位置的加权平均

  25. 不同的芯位重建方法 5、平滑算法-----小波法 Fourier 变换: F(s)= f(x)=

  26. 不同的芯位重建方法 小波变换: Given a function , if it can meet the requirements as following，we call it mother wavelet. （x） ……….(a) C = < if （0） = 0 = 0 ……….(b)

  27. --Wavelet algorithm 对母小波 作伸缩、平移，得到一小波序列 ｛ a,b(x)｝ a,b(x) = ( ) 系数a (scale)表明小波函数的平滑程度 系数b表明小波函数的位置.

  28. --Wavelet algorithm A standard choice for the mother wavelet is the so called mexican hat function Where D means D-dimensional area

  29. --Wavelet algorithm

  30. --Wavelet algorithm 小波变换 • r = (x1, x2, …, xD) • r分离的、有限的

  31. --Wavelet algorithm 芯位（假定芯区粒子的密度最高） • 计算出中心地毯区的每个pad当作 r0时的 f(r。) Max(f(r0)) 簇射的芯位 • 把6×5个pad作为一组，用每一组pad中最左下角的pad作为r0 来计算。找出最大的pad组后，再在这一组中寻找 的最大值，则对应最大值的pad的中心为簇射的芯位。

  32. 第四部分 对芯位重建方法的检验

  33. 对芯位重建方法的检验 大气簇射模拟 原初粒子 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

  34. 对芯位重建方法的检验 ARGO探测器模拟

  35. 对芯位重建方法的检验 模拟后得到了27713和14073个 事例和质子事例 用8000个 事例和14073个质子事例来进行检验

  36. 对芯位重建方法的检验 事例的能量随击中PAD数的变化（ 事例）

  37. 对芯位重建方法的检验 事例的能量随击中PAD数的变化（质子事例）

  38. 对芯位重建方法的检验 各种重建方法所需CPU时间( 事例重建)

  39. 对芯位重建方法的检验 各种重建方法所需CPU时间(质子事例重建)

  40. 事例数随方位角的分布 对芯位重建方法的检验 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

  41. 质子事例数随方位角的分布 对芯位重建方法的检验 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

  42. 对芯位重建方法的检验 事例的芯位重建精度随方位角的变化

  43. 对芯位重建方法的检验 质子事例的芯位重建精度随方位角的变化

  44. 对芯位重建方法的检验 （a）簇射 （b）质子簇射 1 TeV的 -AS和P-AS的次级粒子 在羊八井观测平面处的横向分布

  45. 事例数随天顶角的分布 对芯位重建方法的检验 7～10 10～20 20～30 30～40 40～45 Centre 450 2111 2574 2221 619 Tree 453 2119 2570 2224 610 Mtree 455 2110 2574 2227 608 Dtree 454 2110 2578 2221 612 Wavelet(a=5m ) 455 2085 2589 2238 609 Wavelet(a=10m) 451 2098 2584 2233 608 Wavelet(a=15m) 451 2111 2573 2233 609 Wavelet(a=20m) 454 2113 2576 2227 606 Wavelet(a=25m) 456 2114 2573 2226 605

  46. 7～10 10～20 20～30 30～40 40～50 >50 Centre 706 3659 4547 3910 1231 12 Tree 701 3661 4561 3902 1228 12 Mtree 695 3668 4558 3894 1236 14 Dtree 704 3662 4549 3908 1230 12 Wavelet (a=10m) 690 3645 4583 3882 1247 18 对芯位重建方法的检验 质子事例数随天顶角的分布

  47. 对芯位重建方法的检验 事例的芯位重建精度随天顶角的变化

  48. 对芯位重建方法的检验 质子事例的芯位重建精度随天顶角的变化

  49. 对芯位重建方法的检验 事例的芯位重建精度随击中PAD数的变化 gamma proton

  50. 对芯位重建方法的检验 Scale不同时小波法的芯位重建情况 gamma proton

More Related