BESIII实验的顶点及共振态粒子重建研究
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
北京正负电子对撞机(BEPC)及其大型通用探测器——北京谱仪(BES)改造工程BEPCⅡ/BESⅢ已经成功联调运行,并于2008年获取首批对撞数据。事例顶点重建是BESⅢ数据分析系统的重要组成部分,是物理分析的前提条件。BESⅢ的整体设计中未包含顶点探测器,粒子在探测器上的最内层击中信息位于主漂移室内,共振态粒子的衰变顶点大部分不能被直接探测。事例顶点重建过程中需要解决的主要问题是,如何利用共振态粒子衰变顶点的条件优化末态粒子径迹参数,如何高精度地重建共振态粒子的质量和寿命,如何有效去除前期径迹重建过程中产生的伪径迹,获得高精度的事例初级顶点。高效和精确的顶点重建能够提高径迹的质量分辨,有效地压低本底事例,以便观察信号事例特征。高精度的事例初级顶点重建能够准确地定位对撞束团的中心位置和分辨。
     本论文针对上述问题,围绕BESⅢ上的顶点和共振态粒子重建展开以下研究工作:
     (1)相对长寿命共振态粒子K_s~0,∧/(?)和Ξ重建研究
     共振态粒子K_s~0和∧/(?)是BESⅢ上大量产生的次级衰变粒子,它们的重建关键是衰变顶点的重建和粒子在衰变之前飞行过程中的运动学拟合。本论文详细讨论了全局顶点拟合方法的原理,研究了顶点限定性条件、粒子在飞行过程中满足的运动学限定性条件,以及共振态粒子参数的估算方程;描述了共振态粒子重建软件的设计和实现。利用蒙特卡罗数据进行重建,给出共振态粒子重建的质量,并通过对衰变长度的研究给出重建K_s~0,人的本底压制条件。详细分析了共振态粒子K_s~0寿命测量过程中探测器分辨、探测效率函数,设计开发了不借助模拟信息而直接运用数据估算探测效率的算法。蒙特卡罗数据研究表明拟合的寿命结果在误差范围内是合理的,从而为真实数据的共振态粒子寿命测量构建了基础框架。还讨论了级联衰变粒子Ξ的重建,通过对多次衰变粒子的重建研究,进一步验证了算法设计的合理性与可行性。
     (2)基于Kalman滤波方法的顶点重建算法及顶点寻找算法研究
     基于Kalman滤波方法,在BESⅢ离线软件框架BOSS下开发了顶点拟合算法。事例初级顶点重建的精确度直接影响对撞束团的离线重建位置和分辨,并影响产生自初级顶点的次级粒子重建。提出基于Kalman滤波方法的BESⅢ顶点拟合算法,通过对径迹的过滤、预测和平滑三个步骤,有效去除主漂移室径迹重建过程中的伪径迹信息,直接改善事例初级顶点拟合结果。详细讨论了Kalman顶点拟合软件设计与实现,描述了算法在BESⅢ上的具体运用,并结合蒙特卡罗样本优化了算法中相关的参数。提出了BESⅢ上的几何顶点寻找算法,其目的是将事例的末态粒子径迹进行归类。详细阐述了该算法的原理并描述了相应的软件设计和实现,蒙特卡罗数据研究表明经过顶点寻找之后的顶点拟合,其顶点分辨有明显改善。
     (3)BESⅢ实验数据的事例初级顶点和次级顶点重建及数据分析
     BESⅢ实验数据分析是对本论文中顶点重建算法,包括顶点寻找、全局顶点拟合以及基于Kalman滤波的顶点拟合算法的可靠检验。详细描述了针对BESⅢ取数初期,顶点重建的具体方案,通过实验数据中事例初级顶点重建分析了束团对撞中心位置的偏移结果,并根据实际情况改进了物理分析中径迹选择的条件。同时给出了实验数据中共振态粒子K_s~0,人和Ξ的初步重建结果。此外,为了满足BESⅢ数据质量监测系统(DQM)对事例顶点重建运行时间性能的要求,讨论了快速顶点拟合算法,设计并实现了相应的软件包。蒙特卡罗样本测试的结果表明该算法完全可以满足DQM对顶点重建的要求。
The Beijing Electron Positron Collider(BEPC) and Beijing Spectrometer (BES) have been upgraded to BEPCⅡ/BESⅢ.The raw collision data have been collected at a high luminosity since July,2008.Event reconstruction,and especially vertex reconstruction,is an indispensable and important part of BESⅢdata analysis system,and also a precondition for subsequent physics analysis.BESⅢdoes not include a vertex detector in its overall design.The hit information of tracks that closest to the origin locates at the Main Drift Chamber(MDC) of BESⅢ,thus the decay vertices of resonance particles cannot be detected directly. The key problems of vertex reconstruction are:how to optimize track parameters of daughter particles based on the fact that they derive from a common decay point;how to precisely reconstruct the mass and lifetime of relatively long lifetime resonance particles;how to determine strict primary vertex positions and effectively remove the ghost tracks generated from track reconstruction,to suppress backgrounds so as to analysis signal events.High-precision primary vertex reconstruction will be able to determine the average of beam profile.
     In order to solve the above problems on vertex and resonance particle reconstruction at BESⅢ,The research works of this dissertation are:
     (1) Reconstruction of relatively long lifetime resonance particles
     The relatively long lifetime resonance particles,such as K_s~0 and∧/(?),etc,are abundantly produced atτ-charm energy region at BESⅢ.The reconstruction of resonance particles is actually reconstruction of decay vertices and kinematic fitting between the production and decay vertices.In this thesis,a global vertex fitting algorithm based on the Least Squares Method is presented,in which we also study the effects of vertex constraints,kinematic equations and the calculation of track parameters of resonance particles.A global vertex reconstruction program is designed and implemented using C++ language under the framework of BESⅢ offline Software System,BOSS.The results of resonance particle reconstruction based on Monte Carlo simulation are given.The cut on decay length is studied in order to suppress the backgrounds in physics analysis,and the reconstruction efficiency for their lifetime measurements is discussed.In addition,the reconstruction of cascading decay particleΞ,is also studied.It further verifies the applicability and feasibility of the algorithm.
     (2) Vertex reconstruction based on Kalman filter technique and vertex finding algorithm
     The purpose of primary vertex reconstruction is to determine the interaction points of events,and then further determine the average of beam profile in high energy physics experiments.The information of a general interaction point is a critical parameter for reconstruction of secondary particles decayed from that point.In the thesis,a BESⅢvertex fitting algorithm based on Kalman filter technique is developed for handling the problems of ghost tracks and badly fitted tracks.The corresponding program is discussed,and its algorithm parameters are optimized using simulation data.A vertex finding algorithm which is to sort a set of tracks into subsets that share a common point is also proposed.
     (3) Analysis of vertex reconstruction using raw collision data at BE-SⅢ
     For checking the performance of vertex reconstruction algorithms,the beam bunch positions reconstructed from raw collision data collected at BESⅢare discussed.The procedure of raw collision data processing is described in detail, and for physics analysis,an optimized condition of selecting fitted tracks based on the results of primary vertex reconstruction is proposed.The preliminary results of resonance particle reconstruction from raw data is also given and analyzed.In order to achieve the time requirements in BESⅢData Quality Monitoring(DQM) system,a fast vertex fitting algorithm is developed and its results are compared with those of the global vertex fitting using simulation data.
引文
[1]章乃森编著.1987.粒子物理学:上册[M].北京:科学出版社.
    [2]Perkins D H.2000.Introduction to high energy physics[M].Cambridge University Press.
    [3]Barger V,Phillips R.1987.Collider Physics[M].Addison-Wesley Publishing Company.
    [4]Fritzsch H,Gell-Mann M,and Leutwyler H.1973.Phys.Lett.B 47,365.
    [5]李政道.1987.粒子物理和场论简引,科学出版社.
    [6]赖得L.1983.基本粒子与对称性,科学出版社.
    [7]Gross D J and Wilczek F.1973.Phys.Rev.Lett.30,1343.
    [8]Politzer H D.1973.Phys.Rev.Lett.30,1346.
    [9]Eichten T,et al.1973.Phys.Lett.B 46,274-280.
    [10]Arnison G,et al.1983.Phys.Lett.B 126,398.
    [11]Bagnaia P,et al.1983.Phys.Lett.B 126,130.
    [12]郑志鹏,等.1998.北京谱仪正负电子物理[M].广西科学技术出版社.
    [13]PEPII-An Asymmetric B Factory.1993.Conceptual Design Report[G].SLAC-418,LBL-5379.
    [14]BaBar Collaboration.1995.BaBar Technical Design Report[G].SLAC-R-457.
    [15]KEKB B-Factory Design Report.1995.KEK Report[G]95-7.
    [16]Belle Collaboration.1995.Belle Technical Design Report[G].KEK Report 95-1.
    [17]唐孝威,等.1982.粒子物理试验方法[M].人民教育出版社.
    [18]Kisel I.2006.Event reconstruction in the CBM experiment[J].Nuclear Instruments and Methods in Physics Research Section A,Volume 566,Issue 1,Pages 85-88.
    [19]Naur P,Randell B(editors).1968.Software engineering:Report of a conference sponsored by the NATO Science Committee[G],Garmisch,Germany:Scientific Affairs Division,NATO,7-11 October.
    [20]Backus J W,et al.1957.The FORTRAN Automatic Coding System for the IBM 704 EDPM:Preliminary Operator's Manual[M].International Business Machines.
    [21]Dijkstra E.1970.Notes on Structured Programming[G].
    [22]Meyer B.1999.Object-Oriented Software Construction[M].Prentice-Hall,New York,NY,2/ed.
    [23]Stroustrup B.2000.The C++ Programming Language[M].Special Edition(3rd Edition),Addison-Wesley Professional.
    [24]Gaudi project homepage[EB/OL]:http://proj-gaudi.web.cern.ch/proj-gaudi/.
    [25]Root-an object-oriented data analysis framework[EB/OL],http://root.cern.ch/.
    [26]Li W G.BEPCⅡ/BESⅢ Upgrade and the Prospective Physics.
    [27]IHEP.2004.Preliminary Design Report:The BESⅢ Detector[G].
    [28]IHEP.2004.Preliminary Design Report:The BESⅢ Detector[G].Technical Report IHEP-BEPCⅡ-SB-13.
    [29]Allison W.1980.Relativistic Charged Particle Identification by Energy Loss[J].Ann.Rev.,Part Sci.,30:253.
    [30]Geant4[EB/OL],http://geant4.web.cern.ch/geant4.
    [31]Li W D,Liu H M,et al.2006.The Offline Software for the BESⅢ Experiment[G].Proceeding of CHEP06,Mumbai.
    [32]邓子艳,等.2006.面向对象的BESⅢ探测器模拟系统[J].高能物理核物理,30(05):371-377.
    [33]eXtended Markup Language(XML)[EB/OL],http://www.xml.com/.
    [34]Geometry Description Markup Language(GDML)[EB/OL],http://gdml.web.cern.ch/GDML/.
    [35]臧石磊.2005.B~-→J/ψΛ(?)的寻找和BESⅢ漂移室重建软件的研究[Ph.D.].中科院高能物理研究所博士论文.
    [36]张瑶.2007.BESⅢ主漂移室径迹重建软件的设计与开发[Ph.D.].山东大学博士论文.
    [37]张瑶,等.2007.BESⅢ主漂移径迹重建软件MdcPatRec的设计与实现[J].核电子学与探测技术,27(04):665-669.
    [38]张瑶,等.2007.BESⅢ主漂移室基于模式匹配的径迹重建[J].高能物理与核物理,31(06):570-575.
    [39]Billoir P.1984.Track fitting with multiple scatteing:A new method[J].Nucl.Instrum.Method,225(2):352-366.
    [40]Avery P.1992.The charged particle energy correction in CLEO[OL].CLEO Note CBX92-40.
    [41]Zeini F and Avery P.1995.UFTRAK:new energy correction for charged tracks[G].CLEO note CBX 95-34,CLEO.
    [42]He K L.2007.Status and Plan for Analysis and Tools software[R].BESIII Collaboration Meeting.
    [43]Yan L.2008.Status of Kinematics Fitting Software[R].BESIII Collaboration Meeting.
    [44]Qin G,et al.2008.Particle identification using artificial neural networks at BESIII.Chinese Physics C,32(1).
    [45]莫晓虎,朱永生.2001.ψ(2S)扫描实验中衰变宽度误差与亮度的确定[J].高能物理核物理,25(12):1133-1139
    [46]Dobbs S,et al.2005.Measurement of Absolute Hadronic Branching Fractions of D Mesons and e+e—→D(?) Cross Sections at Ec.m.=3773MeV[J].Phys.Rev.Lett.95,121801.
    [47]Sharpe S R.1998.Progress in Lattice Gauge Theory[G].29th International Conference on High-Energy Physics(ICHEP 98),Vancouver,Canada,hep-lat/9811006.
    [48]Toernqvist N A.1981.Found.of Phys,11,171.
    [49]Metcalf M,Regler M,and Broll C.1973.A Split Field Magnet Geometry Fit Program[G].NICOLE,CERN,Geneva,CEAN 73-2.
    [50]Billoir P,Fr(u|¨)hwirth R and Regler M.1985.Track Element Merging Strategy and Vertex Fitting in Complex Modular Detectors[J].Nucl.Instrum.Methods in.Physics Research A241,115.
    [51]Regler M and Fr(u|¨)hwirth It.1990.Reconstruction of Charged Tracks[M].In:Concepts and Techniques in High Energy Physics V,Plenum Publishing Corporation,New York,Ed.T.Ferbel.
    [52]Wang J K.2008.Progress in MDC Kalman Track Fitting[R].BESIII Collaboration Meeting.
    [53]Ohnishi Y.1996.Track parameterization[G].Belle note 148.
    [54]Avery P.1998.Applied Fitting Theory Ⅵ:formulas for kinematic fitting[OL].CLEO Note CBX98-37.
    [55]Kendall M G and Stuart A.1967.The Advanced Theory of Statistics[M].Volume Ⅱ(Interference and Relationship,second edition).Charles Griffin and Company Limited,London.
    [56]Eadie W T,Drijard D,James F E,Roos M and Sadoulet B.1971.Statistical Methods in Experimental Physics[M].North-Holland Publishing Company,Amsterdam.
    [57]Strandlie A.2004.Track reconstruction from bubble chambers to the LHC[J].Nucl.Instrum.Methods A,535(1-2):57-64.
    [58]Avery P.1991.Applied Fitting Theory Ⅰ:general least squares theory[OL].CLEO Note CBX91-72.
    [59]Meyer S L.1975.Data Analysis for Scientists and engineers[M].John Wiley and Sons.
    [60]Aitken A C.1934.On least squares and linear combinations of observations[J1.Proceedings of the Royal Society Edinburgh A,55,42.
    [61]He M,et al.2008.Energy loss correction for a crystal calorimeter[J].Chinese Physics C,32(04):269-274
    [62]Zhu Z A.2008.Magnet:results of field mapping[R].BESIII Collaboration Meeting.
    [63]Gamma E,Helm R,Johnson R,and Vlissides J.1995.Design Patterns:Elements of Reusable Object-Oriented Software[M],Addison Wesley.
    [64]Ma X,et al.2008.Determination of event start time at BESIII[J],Chinese Physics C,32(09):744-749.
    [65]Wang D Y.2006.dE/dx offline software and the performance.BESIII technical note[G].Institute of High Energy Physics,CAS.
    [66]Wang D Y,et al.2005.Software developments and studies of BESIII dE/dx reconstruction [J].In Proceedings of the 5th National Conference of Nuclear Instrumentation and Methods.
    [67]王亮亮,等.2007.BESIII主漂移室径迹的外推和径迹匹配[J],高能物理与核物理,31(02):183-188.
    [68]Hu J F,et al.2007.Correlation Analysis in Time-of-Flight Calibration of BESIII [J].High Energy Physics and Nuclear Physics,31(10):893-899.
    [69]He M,et al.2008.Energy loss correction for a crystal calorimeter[J].Chinese Physics C,32(04):269-274.
    [70]Xie Y G,et al.2007.Performance Study of RPC Prototypes for the BESIII Muon Detector[J].High Energy Physics and Nuclear Physics,31(01):70-75.
    [71]Astier P,Cardini A,Robert D C,et al.2000.Kalman filter track fits and track breakpoint analysis[J].Nucl.Instrum.Methods A,450(1):138-154.
    [72]Billoir P and Qian S.1990.Simultaneous patten recognition and track fitting by the Kalman filtering method[J].Nucl.Instrum.Methods A,294(1-2):2190228.
    [73]Emeliyanov D,Gorbounov I and Kisel I.2001.OTIt/ITR-CATS:Tracking based on cellular automaton and Kalman filter[G],HERA-B Note HERA-B 01-137,DESY.
    [74]王喆,袁野,张长春.2003.K_S~0重建及其寿命测量在BESII上的研究[J].高能物理与核物理,Vol 27.
    [75]Ping R G,Pang C Y.2007.Monte Carlo Generators for Tau-Charm-Physics at BESIII[G].BAD 522 V6,Based on EvtGen V00-11-07.
    [76]Back J,et al.2001.A User's Guide to the RooFitTools Package for Unbinned Maximum Likelihood Fitting[G].Babar Analysis Document,Version 12.
    [77]Verkerke W,Kirkby D.2006.RooFit Users Manual[OL].v2.07.
    [78]Andersen E,et al.1998.Enhancement of central A,Ξ and Ω yields in Pb-Pb collisions at 158 A GeV/c[J].Physics Letters B,433:209-216.
    [79]Ushida N,et al.1980.Measurement of D~+,F~+,and Λ_c~+ Charmed-Particle Lifetimes [J].Phys.Rev.Lett.45:1053-1056.
    [80]Aubert B,et al.2005.Precision measurement of the Λ_c~+ baryon mass[J].Phys.Rev.D 72:052006.
    [81]Particle Data Group.2008.Particle physics boooklet[S].Journal of Physics G,33,1.http://pdg.lbl.gov/.
    [82]Fr(u|¨)hwirth R,Prokofiev K,Speer T,et al.2003.New developments in vertex reconstruction for CMS[J].Nuclear Instruments and Methods in Physics Research A,502:699-701.
    [83]Fi(u|¨)hwirth R,Speer T.2004.A Gaussian-sum filter for vertex reconstruction[J].Nuclear Instruments and Methods in Physics Research A,534:217-221.
    [84]Jackson D J.1997.A Topological vertex reconstruction algorithm for hadronic jets[J],Nuclear instruments and methods in physics research A,247-253.
    [85]Thomas L.1995.Vertex Reconstruction and Fitting[EB].ftp://x4u2.desy.de/pub/herab/vertex/vertex.ps.
    [86]Garcia E,et al.2007.Vertex reconstruction using a single layer silicon detector[J]. Nuclear Instruments and Methods in Physics Research A,570:536-542.
    [87]Frühwirth R,Kubinec P,Mitaroff W,Regler M.1996.Vertex reconstruction and track bundling at the LEP collider using robust algorithms [J].Computer Physics Communications,Volume 96,Issues 2-3,Pages 189-208.
    [88]Frühwirth R,et al.1996.Vertex reconstruction and track bundling at the LEP collider using robust algorithms [J].Computer Physics Communications 96:189-208.
    [89]D'Hondt J,et al.2003.Properties of robust vertex fitting algorithms at high luminosity HEP experiments [J],Nuclear Science Symposium Conference Record,Volume 2:858-862.
    [90]Narain M.2001.Vertex Reconstruction by means of the Kalman Filter at the DO Detector [J],American Physical Society,April Meeting,Washington,DC Bulletin of the American Physical Society,Vol.46,No.2.
    [91]Wouter D H.2005.Decay chain fitting with a Kalman filter [J].Nuclear Instruments and Methods in Physics Research A 552:566-575.
    [92]Frühwirth R,Regler M,Bock R K,et al.2000.Data Analysis Techniques for High-Energy Physics[M].2nd edition.Number ISBN 0-521-63219-6,Cambridge University Press.
    [93]Rudolph E K.1960.A new approach to linear filtering and prediction problems[J].Transactions of the ASME-Journal of Basic Engineering,82(Series D):35-45.
    [94]Robert G B.1983.Introduction to random signal analysis and Kalman filtering[M].Number ISBN 0-471-08732-7,John Wiley and Sons.
    [95]Wolin E J,Ho L L.1993.Covariance matrices for track fitting with the Kalman filter[J].Nuclear Instruments and Methods in Physics Research A,329(3):493-500.
    [96]Robert K K.1998.Example of a Kalman filter on a straight line[EB].from Kutschke's web homepage.
    [97]Bischof H,Frühwirth R.1998.Recent developments in pattern recognition with applications in high-energy physics[J].Nuclear Instruments and Methods in Physics Research A,419:259-269.
    [98]Frühwirth R.1987.Application of Kalman filtering to track and vertex fitting[J].Nuclear Instruments and Methods in Physics Research A,262:444-450.
    [99]Nelder J,et al.1965.A simplex method for function minimization[M].Comp.J.
    [100]Rosenbrock H.1962.An automatic method for finding the greatest or least value of a function[M].Comp.J.
    [101]Powell M.1964.An efficient method for finding the minimum of a function of several variables without calculating derivatives[M].Comp.J.
    [102]MySQL[EB/OL].Http://www.mysql.com/.
    [103]Huang B.2009.The VertexDbSvc in BOSS6.4.3[R].BESⅢ Physics and Software Meeting.
    [104]Luo T,Wu L H and Wu Z.2008.MDC calibration and Alignment[R]miniworkshop on BESⅢ detector performance.
    [105]Fu C D,et al.2008.Study of the online event filtering algorithm for BESⅢ[J].Chinese Physics C,32(05):329-337.
    [106]Wang Z Y.2009.ψ'scan.total number.and the luminosity[R].BESⅢ Collaboration Meeting.
    [107]Ma X.2009.Progress of event time determination[R].BESⅢ Collaboration Meeting.
    [108]Zhang Y.2009.Status of MDC Tracking[R].BESⅢ Collaboration Meeting.
    [109]Wu L H.2009.Status of MDC Calibration[R].BESⅢ Collaboration Meeting.
    [110]Cao X X.2009.Status of dE/dx Calibration[R].BESⅢ Collaboration Meeting.
    [111]Sun Shengsen.2009.Status of TOF calibration[R].BESⅢ Collaboration Meeting.
    [112]He M.2008.Hadronic events selection[R].BESⅢ Physics and Software Meeting.
    [113]Sun X D.2009.Status of Data Quality monitoring[R],BESⅢ Collaboration Meeting.
    [114]Billoir P,Qlan S.1992.Fast vertex fitting with a local parametrization of tracks[J].Nuclear Instruments and Methods in Physics Research A,311:139-150.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700