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

Presentation on theme: "羊八井ARGO实验簇射芯位的重建 刘 传 磊 物理与微电子学院."— Presentation transcript:

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

2 羊八井ARGO实验簇射芯位的重建 Contents 羊八井ARGO实验 离线数据分析软件介绍 不同的芯位重建方法 对芯位重建方法的检验
结论及展望 ARGO实验合作组开发研制的离线数据分析软件-MEDEA++

3 第一部分 羊八井ARGO实验

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

5 羊八井ARGO实验 111.01m 77.87m 98.84m 76.36m 一万平米的RPC地毯式阵列，中央五千米的地毯

6 羊八井ARGO实验

7 羊八井ARGO实验 1.25m 2.8m

8 羊八井ARGO实验 事例率高、数据量大 实验数据的处 理和存储 medea++
事例率高、数据量大 实验数据的处 理和存储 medea++ 使得对实验数据 面向对象的编程语言C++开发了离线数据分析软件MEDEA++

9 第二部分 离线数据分析软件介绍 下面介绍一下medea++，这个软件可以对ARGO所探测到的EAS簇射事例进行重建并可以对数据进行存储和管理，还可以做物理分析。

10 离线数据分析软件介绍 分析软件的整体结构 探测器阵列的几何描述 事例信息的描述和重建 事例的存储 物理分析工具
软件的整体结构包括下面4部分。 探测器阵列的几何描述：给出探测器的组成和几何尺寸以及这些组成在阵列中的位置等信息。事例信息的描述和重建部分是整个程序的核心部分。他主要是对实验上探测器获得的或者MC数据以及重建后的事例信息的描述。同时它控制事例的整个重建过程。存储是采用程序提供的ROOT数据结构来存储事例信息。分析工具，他提供了某些物理量的直方图或NTUPLE文件，所以用户可以根据需要产生一些文件来进行物理分析工作。

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

12 离线数据分析软件介绍 事例信息的描述和重建 1、 下面4个类完成对事例信息的描述，EVDAQ类描述来自阵列数据获取系统的事例信息。
EVRAW类描述可直接用于重建的事例的信息。EVREC1类包括出去躁声的击中 EVREC2重建后的事例的信息 RECLEVEL0完成由EVDAQ到EVRAW信息的传输 RECLEVEL1去掉非活动区的击中和躁声击中 RECLEVEL2完成重建工作，平面拟合然后芯位重建然后圆锥拟合

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

14 离线数据分析软件介绍 事例信息的存储 物理分析工具
存储是采用程序提供的ROOT数据结构来存储事例信息。分析工具，他提供了某些物理量的直方图或NTUPLE文件，所以用户可以根据需要产生一些文件来进行物理分析工作。

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 长度这一概念。
Mtree tree and dtree 三种方法都要用到tree 长度这个定义。 ???????

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

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

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

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

27 --Wavelet algorithm 对母小波 作伸缩、平移，得到一小波序列 系数a (scale)表明小波函数的平滑程度
对母小波 作伸缩、平移，得到一小波序列 ｛ 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分离的、有限的 r = (x1, x2, …, xD)
小波变换如式子所示，对于探测器阵列来说，变量r是不连续和有限的。小波变换是一个密度和小波函数的求和。其中小波函数的是这样定义的：SCALE

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

32 第四部分 对芯位重建方法的检验 检验是通过MC方法得到一些事例数据，然后MC给出的芯位与重建的芯位差值，这一个差值分别在不同的添丁角范围、方位角范围

33 对芯位重建方法的检验 大气簇射模拟 原初粒子 Gamma/Proton 能量(TeV) 0.1-10
模拟事例 天顶角 (度) ~45 能谱指数 － 强相互作用模型 QGSJET 方位角（度） 0~ 强子能量cutoff(GeV) 0.3 Muon能量cutoff(GeV) 电子能量cutoff(GeV) 光子能量cutoff(GeV) 观测平面高度（m）

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

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

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

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

38 对芯位重建方法的检验 各种重建方法所需CPU时间( 事例重建 ) centre 8’00“ Wavelet (a=10m) 28’ Tree
32‘ (a=15m) Mtree 31’04’’ (a=20m) 27’42’’ Dtree 31’02’’ (a=25m) (a=5m)

39 对芯位重建方法的检验 各种重建方法所需CPU时间(质子事例重建 ) centre 12‘ Wavelet (a=10m) 39’53“
Tree 38‘ (a=15m) 39’51“ Mtree 40‘ (a=20m) 43‘’’ Dtree (a=25m) 42’ (a=5m) 39’40“

40 对芯位重建方法的检验 事例数随方位角的分布
0~30 ~60 ~90 ~120 ~150 ~180 ~210 ~240 ~270 ~300 ~330 ~360 Centre Tree Mtree Dtree Wavelet(a=10m)

41 对芯位重建方法的检验 质子事例数随方位角的分布
0~30 ~60 ~90 ~120 ~150 ~180 ~210 ~240 ~270 ~300 ~330 ~360 Centre Tree Mtree Dtree Wavelet(a=20m)

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

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

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

45 对芯位重建方法的检验 事例数随天顶角的分布 7～10 10～20 20～30 30～40 40～45
7～ ～ ～ ～ ～45 Centre Tree Mtree Dtree Wavelet(a=5m ) Wavelet(a=10m) Wavelet(a=15m) Wavelet(a=20m) Wavelet(a=25m)

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

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

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

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

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

51 第五部分 结论及展望

52 结论及展望 结论 对于 簇射事例 小波法的重建结果最好，精度最高。其次是Mtree方法。 对于质子簇射事例
对于 簇射事例 小波法的重建结果最好，精度最高。其次是Mtree方法。 对于质子簇射事例 小波法的重建结果不是令人满意，这是由于这种平滑算法不能很好地描述质子簇射的横向分布特征。而其它几种方法的重建效果没有很明显的区别。

53 结论及展望 下一步需要做的工作 考虑到阵列的边缘效应 对于质子簇射事例 更加理想的重建方法 ???????

54 谢 大 家 感谢张学尧教授的指导！ 感谢各位教授的帮助！ 感谢同学们的支持！