全天區(qū)掃描法尋找gamma源(二)

全天區(qū)掃描法尋找gamma源（二）山東大學(xué)高能物理組李愛鳳2005.10.18對前段時間我的工作的總結(jié)：利用AS實驗的數(shù)據(jù)樣本，基于全天區(qū)掃描法尋找gamma源，計算精度比以前精細(xì)。內(nèi)容：一：全天區(qū)掃描法尋找gamma源的方法簡介。二：對數(shù)據(jù)進(jìn)行的方位角修正和大尺度各向異性修正。三：介紹我做出的結(jié)果。四：下一步的工作。一：全天區(qū)掃描法尋找gamma源的方法簡介在地平坐標(biāo)系中,將天區(qū)按恒星時、天頂角、方位角劃分成細(xì)小單元（l,m,n）。統(tǒng)計每個bin的事例樣本，以bin的中心值（lst,,）作為事例的位置信息。在赤道坐標(biāo)系中，將羊八井天區(qū)在赤經(jīng)、赤緯方向分成細(xì)小的單元（i,i）。統(tǒng)計每個bin的事例樣本，以bin的中心值（i,i）作為源的位置信息。地平坐標(biāo)系：Time(l):localsiderealday,3600bins(24s/bin)Zenith(m):<40o,500bins(0.08operbin)Azimuth(n):0o~360o,360o/0.08o.sin(zenith)binsNorthPoleZenithNonNoffNoffNon(l,m,n)--------I(i,j)總共：3600×754261個binsRa.(j)赤道坐標(biāo)系：Dec.(i):0o~80onorthsky（-9.89o——70.11o）Ra.(j):0o~360oCells:0.1o(Dec,)×0.1o(Ra,),800×3600binsDec(i)Non(m,n,l)--------I(i,j)

地平坐標(biāo)系：某一時刻，天區(qū)在地平坐標(biāo)系天頂角方位角的剖面展開圖示天頂角方向方位角方向0360040NonNoff對off源區(qū)的選擇規(guī)則：1：天頂角與on區(qū)相等。2：分部在on區(qū)的左右兩側(cè)所有的窗口。3：與on區(qū)窗口大小相等。估計本底的方法：等天頂角法(統(tǒng)計背源窗口的事例數(shù)來估計向源窗口的本底)；利用最小2方法擬合數(shù)據(jù)樣本給出各單元（赤經(jīng)赤緯）的相對宇宙線強度（I(I,j)

）。offoffoffoffoffonoffoffoffoffoffAzim.(l)迭代法求其中k是迭代次數(shù)其中k為迭代次數(shù)，k=0時，＝1

對2取極小（20多次的迭代）獲得各天區(qū)的相對宇宙線強度及其誤差I(lǐng)i,j

省略從宇宙線相對強度的角度考慮，計算各個天區(qū)的顯著性歸一化Keep:IfI(i,j)isacompletesolution,thenconst*I(i,j)isalsoasolution（1）全天區(qū)歸一化（2）一個赤緯帶上的所有的格子進(jìn)行歸一化。平滑處理（考慮到實驗的角分辨率）以天區(qū)單元為中心，以角分辨率為角半徑的范圍內(nèi)天區(qū)總強度的變化情況實驗的角分辨率：0.9o，考慮事例在不同統(tǒng)計量重的權(quán)重，加權(quán)平均的平滑處理方法。Smooth后的顯著性：R三：方位角修正

事例隨方位角分布不均勻性(2%)（AS實驗有一定的系統(tǒng)誤差，目前這種系統(tǒng)誤差理解為：陣列所處的地勢有一定的偏斜，探測器沒有處在一定的水平面上，東南方向較低，西北方向較高，傾斜度大約為1.3o。）方位角的分布的不均勻性依賴天頂角（天頂角越大，分布越不均勻）。事例方位角分部不均勻隨不同天頂角的變化修正方法：按天頂角環(huán)帶修正。方位角方向天頂角方向r三：大尺度各向異性的修正從宇宙線知識上來說，宇宙線事例在大尺度上是各向同性的。但作出結(jié)果宇宙線強度在大尺度上的各向異性的。全天區(qū)顯著性分部圖這種大尺度的各向異性必須加以修正。具體做法是：對背景進(jìn)行修正赤經(jīng)（00-3600）赤緯00-800444修正之后的天區(qū)顯著性平滑前的平滑后的四：做出的結(jié)果：全天區(qū)顯著性各向異性修正前各向異性修正后Grab3.413.07Mrk4212.942.73各向異性修正前各向異性修正后GrabMrk421方位角修正前方位角修正后赤經(jīng)赤緯顯著性（）329.25－7.344.56304.35－1.744.5838.7513.965.02221.4513.464.6838.8514.164.52341.5559.265.07(沒做出來)得到了大約30個顯著性高于4.0的位置五：下一步