羊八井ARGO试验簇射芯位的重建.pptVIP

CentreTreeMtreeDtreeWavelets不同的芯位重建方法不同的芯位重建方法N=给定事例的PAD击中数对于某一击中PADi(i=1,…,N):xi,yi≡PAD的中心坐标ti≡粒子到达时间ri(2)≡(xi,yi)ri(3)≡(xi,yi,zi)（withzi=ti·clight）ni≡击中Strip数不同的芯位重建方法1、重心法rc=芯位：击中PAD位置的加权平均不同的芯位重建方法在MTree、Tree和DTree三种方法中，都会用到Tree长度这一概念。一个二维r-tree示例r-tree长度定义如下：不同的芯位重建方法l(r)=不同的芯位重建方法2、Mtree方法芯位：N个Tree长度中最小那个所对应的参照PAD的中心位置不同的芯位重建方法l(r1)=3x2+5x5=31345n1=1n2=2n3=5l(r2)=3x1+4x5=23l(r3)=5x1+4x2=133、Tree方法rc=wherewi=不同的芯位重建方法rc=芯位：考虑到Tree长度时击中PAD位置的加权平均不同的芯位重建方法4、Dtree方法rc=wherewi=(1-)芯位：考虑到Tree长度时击中PAD位置的加权平均不同的芯位重建方法5、平滑算法-----小波法Fourier变换:F(s)=f(x)=小波变换:Givenafunction,ifitcanmeettherequirementsasfollowing，wecallitmotherwavelet.（x）C=if（0）=0=0不同的芯位重建方法……….(a)……….(b)｛a,b(x)｝a,b(x)=()系数a(scale)表明小波函数的平滑程度系数b表明小波函数的位置.--Waveletalgorithm对母小波作伸缩、平移，得到一小波序列--WaveletalgorithmAstandardchoiceforthemotherwaveletisthesocalledmexicanhatfunctionWhereDmeansD-dimensionalarea--Waveletalgorithm--Waveletalgorithm小波变换r=(x1,x2,…,xD)r分离的、有限的计算出中心地毯区的每个pad当作r0时的f(r。)Max(f(r0))簇射的芯位把6×5个pad作为一组，用每一组pad中最左下角的pad作为r0来计算。找出最大的pad组后，再在这一组中寻找的最大值，则对应最大值的pad的中心为簇射的芯位。--Waveletalgorithm?芯位（假定芯区粒子的密度最高）第四部分对芯位重建方法的检验对芯位重建方法的检验大气簇射模拟原初粒子 Gamma/Proton能量(TeV) 0.1-10模拟事例 100000 天顶角(度) 7~45能谱指数 －1.49 强相互作用模型QGSJET方位角（度） 0~360 强子能量cutoff(GeV) 0.3Muon能量cutoff(GeV) 0.3 电子能量cutoff(GeV) 0.015光子能量cutoff(GeV) 0.015 观测平面高度（m） 4300对芯位重建方法的检验ARGO探测器模拟对芯位重建方法的检验模拟后得到了27713和14073个事例和质子事例用8000个事例和14073个质子事例来进行检验对芯位重建方法的检验事例的能量随击中PAD数的变化（事例）对芯位重建方法的检验事例的能量随击中P