• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    暗物質(zhì)粒子探測器徑跡重建方法的研究?

    2016-06-27 08:14:07魯同所雷仕俊藏京京常進伍
    天文學(xué)報 2016年3期
    關(guān)鍵詞:徑跡暗物質(zhì)卡爾曼濾波

    魯同所雷仕俊藏京京常 進伍 健

    (1中國科學(xué)院紫金山天文臺南京210008)

    (2中國科學(xué)院暗物質(zhì)與空間天文重點實驗室南京210008)

    (3中國科學(xué)院大學(xué)北京100049)

    暗物質(zhì)粒子探測器徑跡重建方法的研究?

    魯同所1,2,3?雷仕俊1,2?藏京京1,2常 進1,2伍 健1,2

    (1中國科學(xué)院紫金山天文臺南京210008)

    (2中國科學(xué)院暗物質(zhì)與空間天文重點實驗室南京210008)

    (3中國科學(xué)院大學(xué)北京100049)

    暗物質(zhì)粒子探測器(DAMPE,Dark Matter Particle Explorer)旨在通過對空間高能粒子的大能量范圍、高能量分辨和高空間分辨的觀測間接探究暗物質(zhì)的存在和分布.徑跡重建是利用探測器不同位置探測通道的觀測結(jié)果重建高能粒子入射位置和方向,其結(jié)果的準(zhǔn)確性決定了探測器的空間分辨本領(lǐng).暗物質(zhì)粒子探測器徑跡重建主要依靠其中的硅陣列探測器(STK,Silicon Tracker)和BGO(鍺酸鉍,Bi4Ge3O12)量能器兩個子探測器的觀測.結(jié)合兩者的設(shè)計原理和結(jié)構(gòu),利用多次束流試驗和地面宇宙線實驗數(shù)據(jù)細(xì)致探討了徑跡重建的基本方法,主要包含3個基本步驟:徑跡點的選取,徑跡點的擬合以及最佳徑跡的判選.粒子在每層STK和BGO晶體中往往留下多個擊中點,首先給出了一種利用BGO重心法得到粒子粗略徑跡并以此限定STK徑跡點的方法.接著對于選定的一組可能的STK徑跡點嘗試了Kalman濾波和直線最小二乘擬合兩種擬合法,并發(fā)現(xiàn)兩種方法分別獨立得到比較一致的結(jié)果,從而驗證了結(jié)果的有效性.最后,對于重建得到的多條可能徑跡,討論了一些利用BGO量能器和STK結(jié)果進行聯(lián)合判選的方法.利用提出的徑跡重建方法對多次束流試驗數(shù)據(jù)的處理表明,暗物質(zhì)粒子探測器的空間分辨完全達(dá)到了設(shè)計要求.

    儀器:探測器,行星和衛(wèi)星:探測,伽馬射線:星際介質(zhì),宇宙學(xué):暗物質(zhì),方法:數(shù)據(jù)分析

    1 引言

    有關(guān)暗物質(zhì)的問題是困擾天文學(xué)和現(xiàn)代物理學(xué)發(fā)展的一大難題,目前對暗物質(zhì)的探測正處在高速發(fā)展階段.在中國,以暗物質(zhì)為目標(biāo)的暗物質(zhì)粒子探測器于2015年12月17日發(fā)射升空,它采用的是大縱深、三維全吸收型的量能器,具有測量能量范圍大、能量分辨率高和本底抑制能力強等優(yōu)勢.根據(jù)物理設(shè)計的要求和模擬的結(jié)果,暗物質(zhì)粒子探測器由4個子探測器組成,從頂層到底層分別為:塑料閃爍體探測器(PSD,Plastic Scintillator Detector)、硅陣列探測器、BGO量能器和中子探測器(NUD,Neutron Detector).如圖1所示.

    塑閃陣列探測器分系統(tǒng)有效探測面積為820 mm×820 mm,整個分系統(tǒng)由兩層共82個X、Y方向互相垂直放置的塑閃單元模塊組成,其中78個塑閃單元模塊的尺寸為884 mm×28 mm×10 mm,另外4個塑閃單元模塊的尺寸為884 mm×25 mm×10 mm.塑閃陣列探測器探測單元的兩端采用光電倍增管將光信號轉(zhuǎn)換為電信號再進一步處理.

    硅陣列探測器有6大層,每個大層由相互垂直的分別沿X、Y方向的兩層硅微條探測器構(gòu)成,主要用于測量入射粒子的方向和徑跡.在第二、第三和第四大層上方各設(shè)置了一層1.0 mm厚的鎢板,用來將伽馬光子轉(zhuǎn)換為正負(fù)電子.考慮到光子轉(zhuǎn)換后,需測量產(chǎn)生的次級粒子徑跡,為實現(xiàn)粒子電荷的區(qū)分,特別是區(qū)分伽馬光子和電子[1],鎢板下方應(yīng)有多層硅微條探測器[2].

    BGO探測器主體由308根BGO晶體構(gòu)成,如圖2所示,其探測面積為60 cm×60 cm. BGO探測器共分為7大層,每大層由X、Y兩小層構(gòu)成.每小層22個探測單元,由22根晶體排列組成,每個探測單元(BGO晶體)的尺寸為2.5 cm×2.5 cm×60 cm.BGO探測器探測單元的兩端采用1個R5610A-01光電倍增管(PMT)將光信號轉(zhuǎn)換為電信號[3?6].

    圖1 子探測器概念圖Fig.1 A sketch of 4 sub-detectors

    圖2 BGO量能器位置示意圖Fig.2 The detailed geometry of BGO calorimeter

    中子探測器位于暗物質(zhì)粒子探測器底層,采用厚度為1 cm、摻硼(B)的塑料閃爍體探測器(Saint Gobian公司生產(chǎn)的BC454).探測器由外形輪廓693 mm×693 mm的塑料閃爍體分割成4個獨立的正方形,并各自去除一個角以耦合一個PMT讀出構(gòu)成一個完整的探測平面.它是為測量宇宙線中的強子(主要為質(zhì)子)與中子探測器上層的物質(zhì)發(fā)生作用產(chǎn)生的次級粒子而設(shè)計的,根據(jù)這些中子在探測器內(nèi)的能量沉積,可以判斷入射粒子的類型,配合BGO量能器進一步區(qū)分質(zhì)子和電子[7?8].

    2 徑跡重建的流程

    考慮探測器的幾何構(gòu)型和精度,為了近似描述粒子的數(shù)學(xué)模型,我們在徑跡擬合中采用的方法如圖3所示.

    圖3 徑跡重建擬合的流程圖Fig.3 A fitting fl ow diagram for reconstructing the incident particle trajectory

    3 BGO量能器中粒子的徑跡重建

    BGO量能器中的徑跡重建思路如下:當(dāng)入射粒子通過量能器時,由于粒子與介質(zhì)原子的相互作用(主要是電磁相互作用和強相互作用),入射粒子把能量沉積在量能器中,根據(jù)簇射形狀,我們大致可以判斷出徑跡方向.具體實現(xiàn)如下:首先是BGO層的選取,設(shè)定一個能量閾值以排除背景噪聲,根據(jù)多次統(tǒng)計結(jié)果,初步設(shè)定2%,也就是當(dāng)每層BGO的能量占總能量超過2%時,該層才作為參考層,參與BGO的徑跡重建.其次,選取每一層中參與徑跡重建的BGO晶體,我們選取每層中能量沉積最高的晶體及左右兩側(cè)各一根共3根晶體計算該層能量沉積的重心,詳見(1)式.如果其中一層中,能量沉積最高的晶體恰好位于邊緣,也就是左側(cè)或者右側(cè)沒有晶體,為了減小誤差,我們只選取這根能量最高的晶體參與重建;如果其中一層中,能量沉積最高的晶體不在邊緣,但是其左側(cè)和右側(cè)的晶體能量相差很大,也就是一側(cè)的能量相對于另一側(cè)基本忽略不計,我們依然選取能量較高的2根晶體參與重建,因為當(dāng)粒子同時穿過相鄰的2根晶體,接著從下方穿出的這種情況是可能存在的.然后做直線擬合:采用的是能量加權(quán)的線性擬合,用于加權(quán)的能量為參與重建的晶體能量之和.(1)式中,xc、xmax、Emax分別表示BGO量能器某層的能量重心位置、該層能量最大的那根晶體的位置和能量.xmax?1、xmax+1、Emax?1、Emax+1分別表示xmax左、右兩側(cè)晶體的位置和能量.

    圖4為垂直入射400 GeV的質(zhì)子徑跡重建效果.該圖為XZ平面的投影圖,圖形右側(cè)幾排點為量能器中實際擊中的BGO晶體的能量重心位置,每個點橫向誤差棒長度表示沉積能量大小,直線為擬合的粒子徑跡.圖中(X<0)6個點是STK的團簇(cluster).

    圖4 從暗物質(zhì)粒子探測器頂部垂直入射的400 GeV的質(zhì)子徑跡重建效果Fig.4A reconstructed particle event display.The particle is a proton with a momentum of 400 GeV vertically impacting on the top surface of DAMPE.

    4 硅陣列探測器(STK)的徑跡重建

    4.1 Cluster的選取

    BGO晶體尺寸比較大,為2.5 cm×2.5 cm×60 cm,所以在重建過程難免會存在一定的誤差,這是系統(tǒng)本身造成的,我們可以利用這個誤差范圍來限制硅探測器團簇的選取.

    圖5為利用BGO量能器信息重建出的粒子徑跡方向的分布圖,在STK徑跡重建過程中選取團簇時可以給出參考.左上圖表示BGO中粒子徑跡在XZ平面內(nèi)的投影與Z軸的夾角θx隨事例累計的分布情況;右下圖表示BGO中粒子徑跡在Y Z平面內(nèi)的投影與Z軸的夾角θy隨事例累計的分布情況;左下圖橫坐標(biāo)表示粒子徑跡在XZ平面內(nèi)的投影與Z軸的夾角θx,縱坐標(biāo)表示徑跡在Y Z平面內(nèi)的投影與Z軸的夾角θy.右上圖表示粒子徑跡與Z軸方向的夾角θ隨事例累計的分布情況.

    根據(jù)測算,為減小偏差,硅探測器中重建的徑跡斜率(傾角)取偏離BGO中重建的徑跡斜率(傾角)±5?,從圖5我們可以清楚地看出,5?以內(nèi)基本包含了絕大多數(shù)的擊中點,截距取±25 mm(一根BGO晶體的厚度),也就是限定這樣一個范圍在硅探測器中尋找徑跡.假定BGO重建的徑跡在XZ平面的投影方程為lc:z=ax+b,如圖6所示.則限定團簇選取范圍的直線方程為:

    其中a、b分別表示投影方程的斜率和截距.

    圖5 BGO中徑跡方向分布圖Fig.5 The track reconstruction performance of BGO calorimeter

    4.2 尋跡

    尋跡是一種模式識別過程,屬于人工智能的范疇[9?10],它的任務(wù)是把一組實驗測量點分成兩類:第1類可分為許多子類,每個子類中的點都是由同一個粒子引起,即對應(yīng)于同一徑跡;第2類則是在充分可信的程度下無法劃分給任何粒子的那些點,即噪聲.

    目前徑跡識別的方法有兩類,即所謂的局域方法和全局方法.局域方法中一次只選擇一條候選徑跡,典型的做法是從少數(shù)幾個點開始(候選徑跡初始化),然后對屬于該候選徑跡的點進行預(yù)測,例如基于已經(jīng)找出的候選徑跡按照當(dāng)前使用的徑跡模型進行內(nèi)插或者外延來預(yù)測徑跡的其他擊中點.如果找到了額外的點,就將它們加入到候選徑跡中去,否則,經(jīng)過一定次數(shù)的嘗試以后(嘗試次數(shù)取決于算法所容許的探測器的探測失效程度),該候選徑跡被丟棄.為了尋找候選徑跡,局域方法總是必須進行無結(jié)果的嘗試,于是在不同的組合中會用到同一個點,因此隨著點數(shù)的增加,計算時間的增加速率遠(yuǎn)高于線性增加的速率.如果所有的對象(擊中點)以相同的方式出現(xiàn)在算法之中,則稱為全局方法.這一算法產(chǎn)生一個徑跡表,利用它能夠更容易地從原始數(shù)據(jù)中找出徑跡,該算法可以被視為一種事例坐標(biāo)或者空間點的全部數(shù)據(jù)集的一般的變換.全局方法的計算用時原則上應(yīng)當(dāng)正比于一個事例中的擊中點數(shù).

    具體到該探測器,當(dāng)高能電子進入硅微條探測器會發(fā)生電離損失,損失很小的能量,但由于硅陣列探測器中共有3 mm厚度的鎢板,所以很大一部分電子在鎢板中會開始發(fā)生電磁簇射[11?12],產(chǎn)生更多的電子,轉(zhuǎn)換為電信號后被探測到.或者在接下來幾層硅微條中會產(chǎn)生很多團簇,當(dāng)然也可能存在反散射粒子與上一層硅微條作用產(chǎn)生的團簇,但相比電磁簇射要小得多.我們選擇團簇的重心位置作為我們處理數(shù)據(jù)過程中的一個點.高能質(zhì)子進入硅微條探測器中,主要發(fā)生電離損失,由于鎢板厚度雖小,但密度較大,也有可能會發(fā)生少量的強子簇射,但是概率比較小[13].

    根據(jù)以上物理原理,結(jié)合限定條件,遍歷尋找所有可能的徑跡.其中一種方法為:遍歷過程中,在限定的范圍內(nèi),我們?nèi)我膺x取2個點去擬合一條直線,然后判斷其他層上是否有靠近該直線的點(依據(jù)是:如圖6,如果點到直線ls的距離小于1 mm,我們就認(rèn)為該點靠近直線),找出靠近直線的點最多,并且方差σ2最小的那條直線,然后利用該直線上的點和靠近直線的點,重新擬合直線,那么重新擬合后的直線就是我們要尋找的硅微條探測器中重建的徑跡.為了減少計算量、提高準(zhǔn)確度,最終我們選取的候選徑跡至少包含3個點,因為理想情況下應(yīng)該有6個點,即便有一個壞點,還有5個點可以參與擬合,如果選取3個以下的點,盡管有2個點便可以構(gòu)成一條直線,但是這種情況出現(xiàn)的概率比較小,幾乎可以忽略不計.

    圖6 團簇的尋找范圍Fig.6 The geometric range for searching the cluster

    在尋跡過程中主要考慮以下兩點:一是徑跡要盡可能長、盡可能是一條直線,通俗地講就是選取的點盡可能多,方差盡可能小.二是要與BGO徑跡的方向最匹配,盡可能一致,也就是斜率和截距盡可能接近,但是實際處理過程中,比對斜率和截距的方式是比較困難的.我們采取一種更加合理、方便的方法,如圖7所示,圖中兩條直線分別表示利用BGO信息重建擬合的徑跡和利用硅微條信息并且遍歷所有可能點后得到的一條可能的徑跡,它們與硅微條的交點分別為X1,X2,···,X6和X′1,X′2,···,X′6.我們很容易得到:

    越小時,這兩條直線越靠近.也就是當(dāng)它最小時,基本上認(rèn)為兩條直線最匹配.另外一種方法:因為任意兩點即可確定一條直線,我們在限定的范圍內(nèi)任意選取2點擬合一條直線,排除偏離利用BGO重建的徑跡5?以上的直線,如圖5所示,因為5?以內(nèi)基本包含了絕大多數(shù)的點(可能的徑跡).然后判斷其他層上是否有靠近該直線的點,找出靠近直線的點最多并且χ2最小的直線,再利用該條直線上的點和靠近直線的點重新擬合直線,那么,新擬合出的直線就是我們要尋找的硅微條探測器中重建的徑跡.

    圖7 硅陣列探測器重建的一條可能的徑跡在XZ平面內(nèi)的投影Fig.7 A projected line of track candidate in XZ plane reconstructed with STK

    4.3 徑跡擬合

    4.3.1 帶有能量加權(quán)的最小二乘法直線擬合

    已知數(shù)據(jù)測量點(xi,yi)(i=1,2,···,n),能量加權(quán)就是令xi=Exixi,yi=Eyiyi,擬合y=b′′+a′′x,求出a′′、b′′:

    其中Exi、Eyi分別表示第xi和第yi個團簇或者晶體對應(yīng)的能量值,a′′和b′′分別表示直線的斜率和截距.

    重建的結(jié)果如圖8和圖9.圖8表示30?方向入射的100 GeV的電子徑跡重建結(jié)果,其中實線是BGO重心法擬合的結(jié)果,虛線為BGO給出的限制,點劃線為硅陣列探測器中粒子徑跡的最佳擬合結(jié)果,這是在XZ平面的投影圖,橫坐標(biāo)表示沿Z軸方向各個擊中的重心位置,縱坐標(biāo)表示沿X方向各個擊中點的重心位置.圖9顯示的是0?方向入射的400 GeV的質(zhì)子徑跡重建效果.

    圖10能很好地檢驗重建效果,左上圖表示BGO中粒子徑跡在XZ平面內(nèi)的投影與Z軸的夾角θx隨事例累計的分布情況;右下圖表示BGO中粒子徑跡在Y Z平面內(nèi)的投影與Z軸的夾角θy隨事例累計的分布情況;左下圖橫坐標(biāo)表示粒子徑跡在XZ平面內(nèi)的投影與Z軸的夾角θx,縱坐標(biāo)表示徑跡在Y Z平面內(nèi)的投影與Z軸的夾角θy.右上圖表示粒子徑跡與Z軸方向的夾角θ隨事例累計的分布情況.圖中對于垂直入射的粒子,理想的重建結(jié)果角度分布的中心位置應(yīng)該在0?,但是對于圖9 0?入射400 GeV質(zhì)子的重建結(jié)果, Y方向基本在零點位置,而X方向偏離中心零點0.35?左右,這可能有兩方面的原因:(1)由于探測器安裝偏差及儀器測量等系統(tǒng)性問題,造成粒子可能不是精確的0?入射.(2)可能是我們的算法還有待完善,沒能找出真實的徑跡.從左上角和右下角的圖形中我們可以看出,在偏離中心位置的還有尾巴,這些基本是重建出的假徑跡,我們還沒能將其全部排除.左下角的圖形也能清楚地看出偏離中心位置還存在一些散點.

    圖8 硅微條探測器與BGO量能器聯(lián)合重建的結(jié)果(100 GeV電子,30?方向入射)Fig.8A trajectory of an electron with a momentum of 100 GeV and an incident angle of 30?.The track is the reconstructed result jointly using BGO and STK.

    圖9 硅微條探測器與BGO量能器聯(lián)合重建的結(jié)果(400 GeV質(zhì)子,0?方向入射)Fig.9 A trajectory of a proton with a momentum of 400 GeV and an incident angle of 0?.The track is the reconstructed result jointly using BGO and STK.

    4.3.2 基于卡爾曼濾波的直線擬合

    卡爾曼濾波是匈牙利數(shù)學(xué)家Rudof Email Kalman在1960年首先提出的一種線形濾波與預(yù)測問題的方法[14].在直線擬合過程中,經(jīng)常用到最小二乘法[15?16],對于有偽點的線性擬合,由于部分偽點與擬合直線偏差大,出現(xiàn)了擬合徑跡偏離真實徑跡的現(xiàn)象,有些偽點的出現(xiàn)是由人為造成的,這種情況可將偽點剔除再進行擬合.但是根據(jù)經(jīng)驗或者直覺將偽點剔除,此種做法缺乏說服力,所以要通過有效方法去判斷.卡爾曼濾波是一種統(tǒng)計估算方法,由于其具有最小無偏方差性,能夠去除實驗中的隨機誤差來獲取更接近真實值的信息,目前被廣泛應(yīng)用于各領(lǐng)域.對于解決大部分的問題來說,它是一種最有用、效率最高的算法.

    圖10 硅微條探測器中的徑跡方向分布圖Fig.10 The track reconstruction performance of the Silicon Tracker Detector

    圖11和圖12為宇宙射線μ子擊中探測器后,利用帶有能量加權(quán)的最小二乘法和卡爾曼濾波分別處理后得到徑跡的斜率和截距的散點分布圖,其中圖11的橫坐標(biāo)表示利用帶有能量加權(quán)的最小二乘法直線擬合后得到的徑跡的斜率(傾角),單位為度;縱坐標(biāo)表示通過卡爾曼濾波處理后得到的徑跡的斜率(傾角),單位為度;圖12的橫坐標(biāo)表示利用帶有能量加權(quán)的最小二乘法直線擬合后得到的徑跡的截距,單位為mm;縱坐標(biāo)表示通過卡爾曼濾波處理后得到的徑跡的截距,單位為mm.從圖11~12中,我們可以清楚地看出,兩者斜率和截距的分布比較集中,大致都在斜率為1的直線型區(qū)域,說明:基于卡爾曼濾波的直線擬合和帶有能量加權(quán)的最小二乘法重建效果大體相當(dāng),基本差別不算太大.

    圖11 兩種方法擬合得到的斜率(傾角)的散點圖Fig.11 A scattered plot of the slopes(inclinations) fitted with two methods

    圖12 兩種方法擬合得到的截距散點圖Fig.12 A scattered plot of the intercepts fitted with two methods

    在0?入射的100 GeV電子和400 GeV的質(zhì)子擊中探測器后產(chǎn)生的數(shù)據(jù)中隨機選取兩個事例,經(jīng)過帶有能量加權(quán)的最小二乘法與卡爾曼濾波方法處理后的徑跡擬合參數(shù)比較如表1.

    表1 徑跡擬合參數(shù)Table 1 The parameters of track fitting

    重建結(jié)果如圖13和圖14所示.圖13和圖14中,實線是利用卡爾曼濾波法重建的徑跡,點線是帶有能量加權(quán)的最小二乘法直線擬合的結(jié)果,虛線是利用BGO量能器簇射情況重建的結(jié)果,黑色點為擊中探測器的位置.我們統(tǒng)計了多個能量段的徑跡重建結(jié)果后發(fā)現(xiàn):對于低能區(qū),基于卡爾曼濾波的擬合結(jié)果要比最小二乘法有明顯的優(yōu)勢,如圖13,實線與點線并不重合,且實線更加接近于真實的徑跡.對于高能區(qū),帶有能量加權(quán)的最小二乘法與卡爾曼濾波結(jié)果差別不大,如圖14,實線與點線基本重合,經(jīng)過探測器位置標(biāo)定后,基本都接近于真實的徑跡.

    5 誤差分析

    造成重建結(jié)果誤差的因素有很多,我們對最主要的幾個可能做一下介紹:

    (1)束流試驗中,探測器位置本身就是粗略的,并不是很精確,或者說探測器某個部件可能會偏離中心位置,而且偏離方向還可能不太一致,這可能會影響我們對坐標(biāo)系的選取,造成探測器某一個部件到坐標(biāo)軸或者坐標(biāo)原點的計算出現(xiàn)誤差.這要求我們處理數(shù)據(jù)之前進行位置校正,確保在同一個坐標(biāo)系下進行.

    (2)由于探測器本身材料的問題,當(dāng)粒子擊中探測器后,可能會產(chǎn)生很多次級粒子,比如多重散射,包含帶電粒子,我們無法判斷團簇是由哪些粒子造成的.當(dāng)然也可能有反沖粒子與探測器作用產(chǎn)生的團簇.

    (3)我們處理數(shù)據(jù)時選取的擊中點是團簇的重心位置,但是粒子實際擊中探測器的位置不一定就是重心位置,可能有偏差.

    圖13 0?入射的100 GeV電子徑跡重建的結(jié)果Fig.13A reconstructed trajectory of electron with a momentum of 100 GeV and an incident angle of 0?

    圖14 0?入射的400 GeV質(zhì)子徑跡重建的結(jié)果Fig.14 A reconstructed trajectory of proton with a momentum of 400 GeV and an incident angle of 0?

    6 總結(jié)

    徑跡重建在暗物質(zhì)粒子探測方面有著很重要的意義,再現(xiàn)入射粒子徑跡,有助于得到入射粒子的種類、能量、強度、入射方向等相關(guān)信息,進而判斷分析入射粒子的產(chǎn)生以及源的方向.如果這些粒子是由暗物質(zhì)衰變或者湮滅產(chǎn)生的,還可以了解暗物質(zhì)的分布情況及與其他物質(zhì)分布的相關(guān)性.徑跡重建最重要的幾點是:偽點的剔除、尋跡、徑跡擬合和相關(guān)的修正.過程中要考慮幾個關(guān)鍵問題,比如:探測器性能(主要是幾何構(gòu)形和精度(分辨)),還有能量損失、多重散射等等,然后以足夠好的精度構(gòu)建一個近似描述粒子軌跡的數(shù)學(xué)模型.為了驗證模型和算法的準(zhǔn)確性,在對真實的實驗數(shù)據(jù)進行分析之前,整個重建鏈過程應(yīng)當(dāng)利用模擬數(shù)據(jù)進行測試,模擬數(shù)據(jù)中所估計的參數(shù)的真值應(yīng)當(dāng)是已知的,以便能夠?qū)w一化殘差進行檢查.擬合參數(shù)的精確誤差應(yīng)當(dāng)加以仔細(xì)研究,它們可能受到偽點或者偽徑跡及其去除方法的影響,以及探測器位置安裝偏差的影響.

    本文主要采用能量加權(quán)的最小二乘法和基于卡爾曼濾波的直線擬合方法重建入射粒子徑跡,通過對比,重建徑跡的效果(誤差)跟入射粒子能量大小存在一定關(guān)系.大體上說,入射能量越高,重建結(jié)果誤差可能會越小.另外我們發(fā)現(xiàn):在高能區(qū)利用加權(quán)的最小二乘法直線擬合與卡爾曼濾波經(jīng)過預(yù)測、過濾、平滑得到的結(jié)果基本沒有太大差別,但是在低能區(qū),卡爾曼濾波還是擁有比較明顯的優(yōu)勢.卡爾曼濾波是一種強有力、靈活和有效的工具,它不僅僅適用于徑跡擬合,還適用于最優(yōu)外延和內(nèi)插的計算、偽點的搜索和排除、利用真實的數(shù)據(jù)進行誤差調(diào)整和探測器校準(zhǔn)以及徑跡段的合并.但是卡爾曼濾波也有自身的缺陷,有時候很難找到卡爾曼濾波的一個適當(dāng)?shù)某跏贾?如果初始探測器的角度分辨和動量分辨很差,并且與下一個探測器不是很靠近,第1次的預(yù)測步驟會存在以下問題:

    (1)很難準(zhǔn)確計算粒子真實軌跡穿過的物質(zhì)量,因為預(yù)測的路徑可能與實際路徑完全不符.

    (2)如果初始值偏離真實徑跡過遠(yuǎn),徑跡模型的線性近似將會無效.

    (3)當(dāng)預(yù)測的路徑與下一個探測器表面不相交,預(yù)測步驟將無法進行.

    對于這種問題的一個解決方法是將模式識別方法求得的參考徑跡取為線性近似展開點.首先計算參考徑跡與所有表面的交點(稱為參考狀態(tài)),然后利用狀態(tài)和參考狀態(tài)的差值進行預(yù)測、濾波和平滑化運算.

    此外,儀器在放置時位置測量并不是完全精確的,儀器的標(biāo)定對徑跡重建的結(jié)果也會產(chǎn)生一定的影響,因此,為提高徑跡重建的效果,對儀器進行校準(zhǔn)也是有必要的.

    [1]Chang J.ChJSS,2014,34:550

    [2]董亦凡.暗物質(zhì)粒子探測衛(wèi)星硅探測器數(shù)據(jù)獲取系統(tǒng)的研究.北京:中國科學(xué)院高能物理研究所,2015

    [3]張云龍.空間暗物質(zhì)探測電磁量能器的研究.合肥:中國科學(xué)技術(shù)大學(xué),2011

    [4]謝明剛,郭建華,伍健,等.天文學(xué)報,2014,55:170

    [5]Xie M G,Guo J H,Wu J,et al.ChA&A,2015,39:129

    [6]郭建華,蔡明生,胡一明,等.天文學(xué)報,2012,53:72

    [7]張磊,郭建華,張永強.天文學(xué)報,2014,55:522

    [8]Zhang L,Guo J H,Zhang Y Q.ChA&A,2015,39:380

    [9]Fruhwirth R,Regle M.高能物理數(shù)據(jù)分析.朱永生,譯.合肥:中國科學(xué)技術(shù)大學(xué)出版社,2011

    [10]劉靖,毛澤普,李衛(wèi)國,等.高能物理與核物理,1998,22:587

    [11]汪曉蓮.粒子探測技術(shù).合肥:中國科學(xué)技術(shù)大學(xué),2009

    [12]Mao Z P,Wang T J,Wang S Q.HEPNP,1993,17:193

    [13]Grupen C,Shwartz B,Spieler H.Particle Detectors.Cambridge:Cambridge University Press,2011

    [14]Kalman R E.Transactions of the ASME-Journal of Basic Engineering,1960:35

    [15]王江榮.自動化與儀器儀表,2013:140

    [16]高飛,童恒慶.系統(tǒng)工程與電子技術(shù),2006,28:775

    Researches on Track Reconstruction for DAMPE

    LU Tong-suo1,2,3LEI Shi-jun1,2ZANG Jing-jing1,2Chang Jin1,2WU Jian1,2
    (1 Purple Mountain Observatory,Chinese Academy of Sciences,Nanjing 210008)
    (2 Key Laboratory of Dark Matter and Space Astronomy,Chinese Academy of Sciences,Nanjing 210008)
    (3 University of Chinese Academy of Sciences,Beijing 100049)

    The Dark Matter Particle Explorer(DAMPE)is aimed to study the existence and distribution of dark matter via observation of high energy particles in space with unprecedented large energy bandwidth,high energy resolution,and high space resolution.The track reconstruction is to restore the positions and angles of the incident particles using the multiple observations of di ff erent channels at di ff erent positions, and its accuracy determines the angular resolution of the detector.The track reconstruction is mainly based on the observations of two sub-detectors,namely,the Silicon Tracker(STK)detector and the BGO(Bi4Ge3O12)calorimeter.In accordance with the design and structure of the two sub-detectors and using the data collected during the beam tests and ground tests,we provide a detailed introduction of the track reconstruction of DAMPE data,including three basic steps,the selection of track hits, the fitting of track hits,and the judgement of the best track among(most probably) many of them.Since a high energy particle most probably leaves more than one hit in each level of the STK and BGO,we first provide a method to constrain the STK clusters for the track reconstruction using the rough result of the BGO reconstruction. We apply two di ff erent algorithms,the Kalman filter and the least square linear fitting, to fit the track hits.The consistency of the results obtained independently via the two algorithms con firms the validity of our track reconstruction results,and we discuss the advantages/disadvantages of each method.Several criteria combining the BGO and STK detection are discussed for picking out the most possible track among all the tracks found in the track reconstruction.Using the track reconstruction methods mentioned in this article and the beam test data,we con firm that the angular resolution of DAMPE satis fies the requirement in design.

    instrumentation:detectors,planets and satellites:detection,gamma rays: ISM,cosmology:dark matter,methods:data analysis

    P111

    :A

    10.15940/j.cnki.0001-5245.2016.03.011

    2015-10-16收到原稿,2015-11-16收到修改稿

    ?國家自然科學(xué)基金項目(11303105、11303106、11303107)資助

    ?tslu@pmo.ac.cn

    ?sjlei@pmo.ac.cn

    猜你喜歡
    徑跡暗物質(zhì)卡爾曼濾波
    基于蒙特卡羅模擬方法的圓筒形固體核徑跡氡探測器探測效率的研究
    裂變徑跡LA-ICP-MS/FT法原理、實驗流程和應(yīng)用
    暗能量VS暗物質(zhì)(下)
    什么是暗物質(zhì)
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    什么是暗物質(zhì)
    太空探索(2016年2期)2016-07-12 09:57:45
    基于模糊卡爾曼濾波算法的動力電池SOC估計
    基于擴展卡爾曼濾波的PMSM無位置傳感器控制
    撲朔迷離的暗物質(zhì)
    太空探索(2014年5期)2014-07-12 09:53:42
    核乳膠中質(zhì)子徑跡計算機輔助識別方法研究
    亚洲精品日韩av片在线观看| av国产免费在线观看| 久久精品国产亚洲av香蕉五月| 免费电影在线观看免费观看| 人人妻,人人澡人人爽秒播| 免费观看在线日韩| 成人av一区二区三区在线看| 看十八女毛片水多多多| 99久久精品一区二区三区| 人妻丰满熟妇av一区二区三区| 欧美国产日韩亚洲一区| 级片在线观看| 中文字幕人妻熟人妻熟丝袜美| 日本三级黄在线观看| 麻豆av噜噜一区二区三区| 97超级碰碰碰精品色视频在线观看| 夜夜夜夜夜久久久久| 精品日产1卡2卡| 亚洲男人的天堂狠狠| 亚洲三级黄色毛片| 搞女人的毛片| 中国美女看黄片| 嫩草影院精品99| 国产欧美日韩精品亚洲av| 深夜精品福利| 免费不卡的大黄色大毛片视频在线观看 | 久久中文看片网| 长腿黑丝高跟| 久久人人爽人人爽人人片va| 国产男人的电影天堂91| 丰满乱子伦码专区| 欧美zozozo另类| 在线免费观看不下载黄p国产 | 丰满的人妻完整版| 欧美黑人欧美精品刺激| 日韩欧美免费精品| 99久久精品国产国产毛片| 1024手机看黄色片| 如何舔出高潮| 国产成人a区在线观看| 美女被艹到高潮喷水动态| 成人精品一区二区免费| 国产精品国产高清国产av| 在线a可以看的网站| 999久久久精品免费观看国产| 一进一出抽搐gif免费好疼| 神马国产精品三级电影在线观看| 一a级毛片在线观看| 国产精品av视频在线免费观看| 一进一出好大好爽视频| or卡值多少钱| 热99在线观看视频| 看免费成人av毛片| 亚洲avbb在线观看| 午夜免费激情av| a级毛片免费高清观看在线播放| 成人毛片a级毛片在线播放| 精品国内亚洲2022精品成人| 国产69精品久久久久777片| 干丝袜人妻中文字幕| 亚洲人成网站在线播| 尤物成人国产欧美一区二区三区| 日本一本二区三区精品| 成人综合一区亚洲| 搡老熟女国产l中国老女人| 熟妇人妻久久中文字幕3abv| 一级黄色大片毛片| eeuss影院久久| 午夜免费成人在线视频| 亚洲一区二区三区色噜噜| 麻豆久久精品国产亚洲av| 成年女人看的毛片在线观看| 99久久久亚洲精品蜜臀av| 一区二区三区高清视频在线| 色在线成人网| 五月玫瑰六月丁香| 女人十人毛片免费观看3o分钟| 露出奶头的视频| 全区人妻精品视频| 日本免费一区二区三区高清不卡| 日韩欧美精品v在线| 久久久久久伊人网av| 久久精品国产清高在天天线| 国产成年人精品一区二区| АⅤ资源中文在线天堂| 亚洲av免费在线观看| 日韩欧美三级三区| 亚洲一级一片aⅴ在线观看| 天堂动漫精品| 亚洲av一区综合| 少妇猛男粗大的猛烈进出视频 | 99在线视频只有这里精品首页| 热99re8久久精品国产| 很黄的视频免费| 日本黄色片子视频| 国产精品永久免费网站| 99久久精品一区二区三区| 国产精品国产高清国产av| 国产亚洲精品久久久久久毛片| 国产精品久久久久久久电影| 国内少妇人妻偷人精品xxx网站| 精品人妻1区二区| 热99在线观看视频| 变态另类丝袜制服| 天美传媒精品一区二区| 最后的刺客免费高清国语| 91精品国产九色| 国产伦一二天堂av在线观看| 国产不卡一卡二| 欧美成人a在线观看| 亚洲精品乱码久久久v下载方式| 日本成人三级电影网站| 亚洲一级一片aⅴ在线观看| 国内精品宾馆在线| 精品久久久噜噜| av福利片在线观看| 一区二区三区激情视频| 一区二区三区免费毛片| 亚洲国产欧美人成| 又粗又爽又猛毛片免费看| 久久久午夜欧美精品| 亚洲av中文av极速乱 | 欧美日本亚洲视频在线播放| 一区二区三区高清视频在线| 午夜福利在线观看免费完整高清在 | 国产亚洲欧美98| 欧美成人性av电影在线观看| 日本撒尿小便嘘嘘汇集6| 国产大屁股一区二区在线视频| 亚洲欧美日韩东京热| 精品人妻1区二区| 欧美黑人欧美精品刺激| 国产成人影院久久av| 中文字幕免费在线视频6| 成人欧美大片| 在线a可以看的网站| 美女大奶头视频| 欧美性感艳星| 高清毛片免费观看视频网站| 十八禁国产超污无遮挡网站| 人妻制服诱惑在线中文字幕| 成人国产一区最新在线观看| 成人欧美大片| 赤兔流量卡办理| 能在线免费观看的黄片| 日本三级黄在线观看| 又黄又爽又刺激的免费视频.| 国产亚洲精品av在线| 日韩中文字幕欧美一区二区| 五月伊人婷婷丁香| 色哟哟·www| 毛片女人毛片| 一卡2卡三卡四卡精品乱码亚洲| 午夜精品一区二区三区免费看| 黄色配什么色好看| 国产乱人视频| 亚洲美女视频黄频| 中文字幕熟女人妻在线| 国模一区二区三区四区视频| 91精品国产九色| 亚洲无线观看免费| 精品一区二区免费观看| 99热这里只有是精品在线观看| 精品福利观看| 亚洲最大成人av| 老司机福利观看| 中文字幕人妻熟人妻熟丝袜美| 在线a可以看的网站| 亚洲熟妇熟女久久| 精品午夜福利视频在线观看一区| 日本-黄色视频高清免费观看| 少妇的逼好多水| 日本欧美国产在线视频| 午夜精品在线福利| 91久久精品国产一区二区三区| 国产麻豆成人av免费视频| 成人美女网站在线观看视频| 老熟妇乱子伦视频在线观看| eeuss影院久久| 日日干狠狠操夜夜爽| 99热这里只有精品一区| 又紧又爽又黄一区二区| 哪里可以看免费的av片| 俺也久久电影网| 人妻丰满熟妇av一区二区三区| 国产一级毛片七仙女欲春2| 色吧在线观看| 亚洲精品国产成人久久av| 日韩高清综合在线| 亚洲人成网站在线播| 欧美日韩瑟瑟在线播放| 欧美日本亚洲视频在线播放| 舔av片在线| 国产中年淑女户外野战色| 小蜜桃在线观看免费完整版高清| 欧美区成人在线视频| 成人国产麻豆网| 我要搜黄色片| 国产v大片淫在线免费观看| 午夜日韩欧美国产| 午夜福利18| 久久亚洲真实| 中文字幕久久专区| 村上凉子中文字幕在线| 精品久久久久久久久av| 夜夜夜夜夜久久久久| 日韩精品中文字幕看吧| 美女免费视频网站| 国产黄色小视频在线观看| 如何舔出高潮| 黄色一级大片看看| 一进一出好大好爽视频| 成人性生交大片免费视频hd| 午夜精品一区二区三区免费看| 日韩精品青青久久久久久| 日本免费a在线| 亚洲av电影不卡..在线观看| 夜夜看夜夜爽夜夜摸| 九色成人免费人妻av| 1024手机看黄色片| 国产伦在线观看视频一区| 午夜福利欧美成人| 夜夜爽天天搞| 午夜久久久久精精品| 日韩欧美 国产精品| 国内毛片毛片毛片毛片毛片| 内地一区二区视频在线| 白带黄色成豆腐渣| 香蕉av资源在线| 97热精品久久久久久| 国产在视频线在精品| 色哟哟哟哟哟哟| 长腿黑丝高跟| 99精品久久久久人妻精品| 日韩欧美三级三区| 91麻豆精品激情在线观看国产| 午夜免费激情av| 美女被艹到高潮喷水动态| 一区二区三区高清视频在线| 少妇的逼好多水| 国产人妻一区二区三区在| 精品久久久噜噜| 免费在线观看影片大全网站| 国产av不卡久久| 成人性生交大片免费视频hd| 精品久久久久久久人妻蜜臀av| 久久欧美精品欧美久久欧美| 在线观看美女被高潮喷水网站| 十八禁国产超污无遮挡网站| 国产蜜桃级精品一区二区三区| 久久99热6这里只有精品| 精品一区二区三区视频在线观看免费| 成人午夜高清在线视频| 国产亚洲精品av在线| 精品午夜福利在线看| 成人亚洲精品av一区二区| 亚洲 国产 在线| 国产成人影院久久av| 欧美人与善性xxx| 亚洲18禁久久av| 精品99又大又爽又粗少妇毛片 | 国产精品久久久久久av不卡| 免费av观看视频| 在线免费观看的www视频| 男人狂女人下面高潮的视频| 日本五十路高清| 日韩亚洲欧美综合| 熟女电影av网| 性插视频无遮挡在线免费观看| 中文资源天堂在线| a级毛片a级免费在线| 国产精品一区二区免费欧美| 搡老岳熟女国产| 国产成人aa在线观看| 欧美区成人在线视频| 亚洲av第一区精品v没综合| 国产伦精品一区二区三区四那| 日日撸夜夜添| 欧洲精品卡2卡3卡4卡5卡区| eeuss影院久久| 韩国av在线不卡| 色哟哟哟哟哟哟| 色精品久久人妻99蜜桃| 婷婷亚洲欧美| 男女之事视频高清在线观看| 欧美人与善性xxx| 国产精品日韩av在线免费观看| 又爽又黄a免费视频| 亚洲精品成人久久久久久| 午夜精品久久久久久毛片777| 欧美区成人在线视频| 岛国在线免费视频观看| 欧美激情国产日韩精品一区| 欧美一区二区亚洲| 亚洲国产日韩欧美精品在线观看| 琪琪午夜伦伦电影理论片6080| 国产精品一区二区免费欧美| 国产精品国产三级国产av玫瑰| 女人被狂操c到高潮| 一区二区三区高清视频在线| 国产淫片久久久久久久久| 婷婷亚洲欧美| 一区二区三区四区激情视频 | 日韩欧美在线二视频| 亚洲精品在线观看二区| 热99re8久久精品国产| av.在线天堂| 欧美成人a在线观看| 亚洲性久久影院| 天堂动漫精品| 国产高清三级在线| 老司机福利观看| 一级黄色大片毛片| 国产一区二区激情短视频| 夜夜夜夜夜久久久久| 精品福利观看| 美女xxoo啪啪120秒动态图| 99精品在免费线老司机午夜| 九色成人免费人妻av| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 欧美中文日本在线观看视频| 日韩欧美 国产精品| 给我免费播放毛片高清在线观看| 日本黄色视频三级网站网址| 美女cb高潮喷水在线观看| 成人国产综合亚洲| 俄罗斯特黄特色一大片| 欧美最黄视频在线播放免费| 九色成人免费人妻av| 精品一区二区三区av网在线观看| 精品午夜福利在线看| 老女人水多毛片| eeuss影院久久| 啦啦啦啦在线视频资源| 久久国产乱子免费精品| 精品日产1卡2卡| 人人妻,人人澡人人爽秒播| 午夜精品久久久久久毛片777| 国产一级毛片七仙女欲春2| 国产成人影院久久av| 老熟妇乱子伦视频在线观看| 亚洲av第一区精品v没综合| 亚洲国产日韩欧美精品在线观看| 美女xxoo啪啪120秒动态图| 午夜日韩欧美国产| 久久久久久国产a免费观看| 亚洲欧美日韩高清在线视频| 国产成年人精品一区二区| 国产 一区 欧美 日韩| 国产精品久久久久久精品电影| 免费黄网站久久成人精品| 精品久久久久久久久亚洲 | 他把我摸到了高潮在线观看| 99久久成人亚洲精品观看| 久久这里只有精品中国| 中文字幕人妻熟人妻熟丝袜美| 午夜福利欧美成人| 亚洲男人的天堂狠狠| 少妇被粗大猛烈的视频| 日本三级黄在线观看| 亚洲美女视频黄频| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 九九爱精品视频在线观看| 日韩精品有码人妻一区| 亚洲国产精品久久男人天堂| 国产精品一及| 天堂影院成人在线观看| 他把我摸到了高潮在线观看| 久久久久久久亚洲中文字幕| 床上黄色一级片| 久久久久久久亚洲中文字幕| 国产男人的电影天堂91| a级一级毛片免费在线观看| 国国产精品蜜臀av免费| 国产亚洲精品av在线| av在线天堂中文字幕| 国产伦在线观看视频一区| 在线观看美女被高潮喷水网站| 欧美人与善性xxx| 噜噜噜噜噜久久久久久91| 小蜜桃在线观看免费完整版高清| 99九九线精品视频在线观看视频| 欧美+日韩+精品| 免费看美女性在线毛片视频| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 国产精品1区2区在线观看.| 最近最新免费中文字幕在线| 夜夜看夜夜爽夜夜摸| 久久九九热精品免费| 麻豆久久精品国产亚洲av| 99精品久久久久人妻精品| 中文字幕av成人在线电影| 无人区码免费观看不卡| 精品久久久久久,| 特大巨黑吊av在线直播| 黄色女人牲交| 精品久久久噜噜| 亚洲精品粉嫩美女一区| 韩国av一区二区三区四区| 99九九线精品视频在线观看视频| 亚洲国产精品合色在线| 亚洲精品久久国产高清桃花| 我的老师免费观看完整版| 可以在线观看的亚洲视频| 又粗又爽又猛毛片免费看| xxxwww97欧美| 在线观看一区二区三区| 一本精品99久久精品77| 国产黄a三级三级三级人| 天堂av国产一区二区熟女人妻| 色哟哟·www| 日本一本二区三区精品| 一区福利在线观看| 国产精品无大码| 最近中文字幕高清免费大全6 | 日韩在线高清观看一区二区三区 | 色精品久久人妻99蜜桃| 亚洲欧美日韩高清在线视频| 99久久久亚洲精品蜜臀av| 欧美在线一区亚洲| 中文在线观看免费www的网站| 国产视频一区二区在线看| а√天堂www在线а√下载| 99久久精品一区二区三区| 99热网站在线观看| 一卡2卡三卡四卡精品乱码亚洲| 日日啪夜夜撸| 久久久国产成人免费| 亚洲人成网站在线播放欧美日韩| 成人午夜高清在线视频| 国产伦在线观看视频一区| 亚洲国产欧美人成| 特级一级黄色大片| 欧美成人一区二区免费高清观看| 999久久久精品免费观看国产| 乱系列少妇在线播放| 国产男靠女视频免费网站| 久久精品国产99精品国产亚洲性色| 少妇被粗大猛烈的视频| av在线老鸭窝| 国产精品日韩av在线免费观看| 色5月婷婷丁香| 日韩av在线大香蕉| 亚洲中文字幕一区二区三区有码在线看| 乱系列少妇在线播放| 2021天堂中文幕一二区在线观| 人妻久久中文字幕网| 黄色女人牲交| 美女被艹到高潮喷水动态| 欧美性感艳星| av福利片在线观看| 亚洲成人久久爱视频| 12—13女人毛片做爰片一| 一本一本综合久久| 2021天堂中文幕一二区在线观| 一卡2卡三卡四卡精品乱码亚洲| 色吧在线观看| 中文字幕av在线有码专区| 一个人免费在线观看电影| av福利片在线观看| 乱系列少妇在线播放| 精品久久久久久,| 国内久久婷婷六月综合欲色啪| 亚州av有码| 久久久久久久久中文| 狠狠狠狠99中文字幕| 中文亚洲av片在线观看爽| 成人性生交大片免费视频hd| 国产 一区精品| 日韩欧美国产在线观看| 亚洲无线在线观看| 婷婷亚洲欧美| 久久久国产成人免费| 久久人人精品亚洲av| 内地一区二区视频在线| 少妇的逼好多水| 亚洲,欧美,日韩| 97热精品久久久久久| 12—13女人毛片做爰片一| 欧美3d第一页| 美女黄网站色视频| 亚洲国产日韩欧美精品在线观看| 国产精品久久视频播放| 大型黄色视频在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 88av欧美| 91狼人影院| 久久天躁狠狠躁夜夜2o2o| 国产高清视频在线播放一区| 又爽又黄无遮挡网站| 亚洲av.av天堂| 赤兔流量卡办理| 欧美在线一区亚洲| 12—13女人毛片做爰片一| 成人二区视频| 亚洲av美国av| 欧美日韩国产亚洲二区| 一本久久中文字幕| 日本五十路高清| 午夜福利成人在线免费观看| 国产精品一区二区三区四区久久| 国产精品99久久久久久久久| 禁无遮挡网站| 国产男靠女视频免费网站| 国产av一区在线观看免费| 国产精品三级大全| 黄色一级大片看看| 99精品久久久久人妻精品| 91av网一区二区| 国产 一区 欧美 日韩| 嫩草影院精品99| 国产乱人伦免费视频| 国产淫片久久久久久久久| 97人妻精品一区二区三区麻豆| 国产aⅴ精品一区二区三区波| 久久久久久大精品| 国产精品98久久久久久宅男小说| 亚洲国产色片| 男人舔女人下体高潮全视频| 国产欧美日韩精品一区二区| 岛国在线免费视频观看| 12—13女人毛片做爰片一| 亚洲欧美精品综合久久99| 亚洲国产欧洲综合997久久,| 久久精品国产亚洲av涩爱 | 免费搜索国产男女视频| 国内少妇人妻偷人精品xxx网站| 成人亚洲精品av一区二区| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 成人性生交大片免费视频hd| 久久精品久久久久久噜噜老黄 | av视频在线观看入口| 永久网站在线| 如何舔出高潮| 99在线人妻在线中文字幕| 亚洲国产精品合色在线| 国产人妻一区二区三区在| 一a级毛片在线观看| 成人欧美大片| 亚洲精品影视一区二区三区av| 精品人妻偷拍中文字幕| 别揉我奶头 嗯啊视频| 亚洲精品国产成人久久av| av天堂中文字幕网| 欧美色欧美亚洲另类二区| 一进一出好大好爽视频| 色5月婷婷丁香| 免费观看的影片在线观看| 亚洲乱码一区二区免费版| 又黄又爽又刺激的免费视频.| 欧美日韩中文字幕国产精品一区二区三区| 少妇高潮的动态图| 热99re8久久精品国产| 丰满的人妻完整版| 色5月婷婷丁香| 亚洲四区av| 国产精品福利在线免费观看| 亚洲在线自拍视频| 女人十人毛片免费观看3o分钟| 精品人妻偷拍中文字幕| 国产av一区在线观看免费| 乱人视频在线观看| av.在线天堂| 极品教师在线视频| 亚洲国产欧美人成| 麻豆成人av在线观看| 久久99热6这里只有精品| 亚洲自偷自拍三级| 日韩精品中文字幕看吧| 欧美xxxx性猛交bbbb| 日日撸夜夜添| 日韩国内少妇激情av| av女优亚洲男人天堂| 男女做爰动态图高潮gif福利片| 在线国产一区二区在线| 日本一二三区视频观看| www日本黄色视频网| 99精品在免费线老司机午夜| 亚洲五月天丁香| 欧美日本亚洲视频在线播放| 国产色婷婷99| 嫩草影院入口| 亚洲乱码一区二区免费版| 啦啦啦观看免费观看视频高清| 久久亚洲真实| 国产三级在线视频| 精品国内亚洲2022精品成人| 亚洲精华国产精华精| 日本与韩国留学比较| 成人综合一区亚洲| 男女视频在线观看网站免费| av中文乱码字幕在线| 又黄又爽又刺激的免费视频.| 亚洲av熟女| 嫩草影视91久久| 国产高清三级在线| 丝袜美腿在线中文| av中文乱码字幕在线| 成人综合一区亚洲| 久久天躁狠狠躁夜夜2o2o| 日韩欧美精品v在线| 国产伦精品一区二区三区四那| 亚洲va在线va天堂va国产| 五月伊人婷婷丁香| 午夜日韩欧美国产| 在现免费观看毛片| 久久精品国产亚洲av香蕉五月| 国产成人一区二区在线| 精品一区二区三区视频在线| 免费av毛片视频| 少妇熟女aⅴ在线视频|