兰州重离子储存环外靶实验终端时间投影室的动量分辨率模拟

李 贺　张 松　卢 飞　钟 晨　马余刚（中国科学院上海应用物理研究所 嘉定园区　上海 20800）2（中国科学院大学　北京 00049）



兰州重离子储存环外靶实验终端时间投影室的动量分辨率模拟

李 贺1，2张 松1卢 飞1钟 晨1马余刚1
1（中国科学院上海应用物理研究所 嘉定园区上海 201800）2（中国科学院大学北京 100049）

时间投影室（Time Projection Chamber， TPC）是在多丝正比室（Multi Wire Proportional Chamber， MWPC）和多丝漂移室（Multi Wire Drift Chamber， MWDC）的基础上发展起来的一种高空间分辨粒子径迹探测器。因其具有很高的位置分辨能力，从而可以获得很好的动量分辨率。在开展兰州重离子储存环外靶实验终端 CEE （Cooler-Storage-Ring External-target Experiment）-TPC的建造之前，估算CEE-TPC的动量分辨能力是一项很有必要的工作。基于Geant4和Kalman Filter等模拟程序，对时间投影室探测器测量带电粒子的整个过程进行了蒙特卡罗模拟，给出了CEE-TPC探测不同种类的带电粒子的动量分辨率，并且对影响探测器动量分辨率的因素进行了分析。通过模拟计算，得到了CEE-TPC对π介子、质子和氘核的典型动量分辨率在5%左右，能够较好地鉴别这三种粒子，为未来CEE-TPC的建造提供了可供参考的技术参数。

时间投影室，动量分辨率，蒙特卡罗模拟

时间投影室（Time Projection Chamber， TPC）是一种应用于粒子物理与核物理大型实验中的大体积气体探测器［1］，它是在多丝正比室（Multi Wire Proportional Chamber， MWPC）和多丝漂移室（Multi Wire Drift Chamber， MWDC）［2-3］的基础上发展起来的一种高分辨粒子径迹探测器，能同时对多个带电粒子的三维空间径迹和能量损失进行测量，成为很多大型对撞机上高能物理实验的中心径迹探测器，如美国RHIC （Relativistic Heavy Ion Collider）高能重离子对撞机上的 STAR-TPC［4］和日本理化研究所大型放射性束装置上的SAMURAI-TPC［5］。

CEE（Cooler-Storage-RingExternal-target Experiment）合作组计划在兰州重离子加速器装置（Heavy Ion Research Facility in Lanzhou， HIRFL）的冷却储存环（Cooler-Storage-Ring， CSR）［6-7］上建造一台低温高密核物质测量谱仪［8-16］，实现对撞能量100-800MeV·u-1范围的重离子碰撞，对碎片产物进行接近全空间的测量，为致密天体性质、核团簇现象、重离子核反应动力学、同位旋非对称及奇异性核物质状态方程、核物质相变和相图等重要科学问题［17-28］提供基础实验数据。该谱仪包括大接受度超导磁铁、大型时间投影室 CEE-TPC、硅像素探测器［29］、高计数率飞行时间谱仪［30］等。CEE-TPC是其中非常重要的探测器。

CEE-TPC为立方体型的结构，上部为阳极，下部为阴极，两极之间即为电子漂移腔，漂移腔四壁有分压电阻和场腔壁，分立的场腔壁存在均匀的阶梯电压，作用是在漂移区产生匀强电场。CEE-TPC放置在大型二级磁铁中间，以保证工作时有与均匀电场平行的均匀磁场。CEE-TPC在阳极附近设有进行电子雪崩放大和读出的端盖探测器。当带电粒子入射到漂移场腔内时，由于电离作用在整条径迹上发生初级电离，产生电子-离子对，初级电离产生电子会在电场的作用下向顶部阳极匀速漂移。电子信号到达端盖探测器后被放大和收集，产生输出电信号。信号经过后续的处理可以对粒子三维径迹进行重建，并且根据能损、磁刚度和粒子在磁场中的偏转得到粒子类型与粒子动量。CEE-TPC漂移区设计长、宽、高分别为110 cm×90 cm×70 cm，整个设备位于超导二级磁铁的磁场区，竖直方向均匀磁场强度为0.5 T。探测器内部的工作气体使用90%的氩气和10%的甲烷。动量分辨率是衡量时间投影室探测器性能的一个重要标准，本文对 CEE-TPC进行了蒙特卡罗模拟并得到了动量分辨率等指标。模拟程序分为入射粒子源的产生、粒子与探测器介质相互作用、电子漂移和读出、三维径迹重建和后续数据处理等几个部分。

1　入射粒子源的产生

CEE-TPC主要用于测量带电粒子的三维径迹，在实际的重离子碰撞核物理实验中，进入探测器的粒子种类繁多，主要有π介子、质子、氘核等，实验室系下出射角度一般小于30°。用碰撞模拟程序［31］

模拟重核的对撞过程，按照能量均匀分布从碰撞产物里挑选出π介子、质子和氘核，将这些粒子的类型、动量、位置、出射方向等信息保存在root［32］文件中，作为下一步模拟粒子与探测器介质相互作用过程的初始粒子源。

2　粒子与探测器介质的相互作用

Geant4［33］由欧洲核子中心开发，可以模拟带电粒子与探测器相互作用的物理过程。它采用了C++面向对象的编程技术，其庞大的数据库包含了各类粒子与各类材料物质相互作用的物理过程和反应截面等数据。根据CEE-TPC的设计参数，在Geant4程序中定义探测器的结构。在Geant4程序中初始化粒子发生器（Particle gun）并注册了π介子、质子和氘核这三种入射粒子及其与气体原子发生相互作用的物理过程列表：如电磁相互作用、强相互作用、输运过程和各种能量的参数化相互作用。运行程序，完成探测器结构的初始化，读入上一步碰撞模型产生的粒子源数据并且作为粒子发生器的初始化参数。当带电粒子经过入射窗进入漂移场腔以后，其在均匀磁场的作用下偏转飞行，并沿飞行路径将工作气体原子电离成正离子和电子对。程序会逐步精细记录带电粒子在飞行路径上的当前的能量损失、径迹长度、发生电离的位置、电离出电子的电荷量以及粒子进入探测器的时间等物理量，并保存在root文件中，作为下一步模拟的输入数据。带电粒子在路径上损失的能量正比于电离所产生电子的数量。图1是Geant4程序模拟不同动量和种类的带电粒子从入射窗进入探测器，并穿过气体介质留下的径迹图。

图1　模拟得到的粒子在TPC中的径迹Fig.1　Simulated tracks of particles at TPC.

3　电子的漂移和读出

Geant4模拟程序最终只是给出了粒子径迹上每个电离位置损失的能量，而损失的能量正比于产生的初级电离电子的数量，要转换成探测器给出的电信号，则要经过电子的漂移和读出盘读出过程。编写程序，定义 CEE-TPC阳极读出板的大小为100cm×80 cm的一个长方形板，再将读出板均分为10000块大小为10 mm×8 mm的长方形小读出块（pad），每块读出块有一个唯一的编号，在程序运行时独立地读出电荷信号。

不同的气体组成和不同电场强度下，电子的漂移速度不同［34］。根据我们 TPC样机的实际测量结果，在程序中定义电子在漂移场腔中的漂移速度为5 cm·μs-1。程序在运行时，电子首先从原初电离的位置漂移到阳极读出盘附近，这个过程中漂移电子簇会按照给定的扩散系数发散。之后电子簇被阳极板前方的电子倍增器件［35］雪崩放大，被放大后的电子簇到达阳极被各个pad收集，形成电信号。相应的每块pad收集的电量、电子到达pad的时间、pad号等信息被储存在root文件中，作为下一步粒子径迹重建的输入文件。图2是阳极读出盘（X-Z平面）收集的电荷分布信息，右侧灰度条代表电荷量。

4　粒子径迹重建与分析

在该模拟工作中，对带电粒子三维径迹的重建采用了 Kalman Filter程序［36］来实现。该程序由STAR-TPC合作组开发及维护，可对探测器内同时产生的数十条径迹进行重建，既可以处理蒙特卡罗模拟的数据，又可以处理真实的实验数据。在运行时，程序根据粒子在阳极读出块上的电荷量分布和电子漂移时间信息，通过重心法，推算出粒子径迹在三维空间位置上的点阵，再利用复杂的寻迹算法进行径迹重建。算法原理是从TPC边缘开始寻找，然后逐点逆推，将属于同一条径迹上的点全部确定下来，排除来自其它径迹的点的干扰，并且拟合出最佳径迹。寻迹的过程是一个逐步拟合、逐步修正的过程，判断某一点是不是属于当前所寻找的径迹遵循设定的判断标准。程序将多条三维径迹同时重建出来，最后再通过径迹在均匀磁场中的偏转半径、粒子在探测器内的能损等信息，计算得到粒子的动量及误差。

运行程序，将CEE-TPC的Geant4模拟的结果输入到Kalman Filter程序中进行径迹重建，经测试整套程序对不同粒子的重建效率约在40%，对重建的结果进行分析。图3为模拟得到的各粒子的动量分辨率，横坐标为粒子在水平面（X-Z平面）上的初始动量，纵坐标即为粒子的动量分辨率。图4为模拟得到的π介子、质子、氘核的粒子鉴别图，横坐标为重建得到的粒子的初始动量，纵坐标为单位长度上粒子在探测器内的能量损失。

图2　阳极读出盘上的电荷分布Fig.2　Charge distribution on anode pad.

图3　粒子在X-Z平面的动量分辨率Fig.3　Momentum resolution of particles on X-Z plane.

图4　粒子鉴别图Fig.4　Identification of different particles.

由图 3可以看出，对于氘核、质子和π介子在500 MeV·c-1以下，CEE-TPC典型的动量分辨率在5%。而对于较高动量的粒子，由于粒子动量增大会导致偏转半径增大，从而使探测器对半径的测量误差增大，所以动量分辨率会逐渐变差。低能的氘核和质子在200 MeV·c-1附近由于多重散射较大，动量分辨会变差。值得注意的是，由于Geant4模拟过程未加入噪声击中（hit），而真实的实验数据会存在噪声击中点，在Kalman Filter程序进行径迹重建的时候，重建效率和重建速度都会受到一定的影响。

如果减小探测器读出板的尺寸，或者降低均匀磁场的强度，亦或者增加读出块的尺寸，都会使CEE-TPC对动量测量的误差增大。反之，增大均匀磁场的强度及面积，减小读出块的尺寸，可以减小动量测量的误差，但这些受到大接收度超导磁铁尺寸和全波形采样电子学和数据获取系统通道数量的限制。

5　结语

利用Geant4和Kalman Filter等蒙特卡罗程序模拟大型时间投影室 CEE-TPC对带电粒子的探测，最终得到了π介子、质子和氘核的动量分辨率，其典型的动量分辨率在5%左右，满足CEE-TPC进行各种重核碰撞物理实验的技术要求。

1Attie D. TPC review［J］. Nuclear Instruments and Methods in Physics Research A， 2009， 598（1）：89-93. DOI：10.1016/j.nima.2008.08.114

2程文静， 张钊， 易晗， 等. 用于多丝漂移室读出的扇入延迟前端电子学［J］. 核技术， 2016， 39（4）：040403. DOI：10.11889/j.0253-3219.2016.hjs.39.040403

CHENG Wenjing， ZHANG Zhao， YI Han， et al. Development of a fan-in and delay module for MWDC based on Flash-ADC acquisition scheme［J］. Nuclear Techniques， 2016， 39（4）：040403. DOI：10.11889/ j.0253-3219.2016.hjs.39.040403

3Yi H， Zhang Z， Xiao Z G， et al. Prototype studies on the forward MWDC tracking array of external target experiment at HIRFL-CSR［J］. Chinese Physics C， 2014，38（12）：39-43. DOI：10.1088/1674-1137/38/12/126002

4Anderson M， Berkovitz J， Betts W， et al. The STAR time projection chamber：a unique tool for studying high multiplicity events at RHIC［J］. Nuclear Instruments and Methods in Physics Research A， 2003， 499：659-678. DOI：10.1016/S0168-9002（02）01964-2

5Shane R， McIntosh A B， ISobe T， et al. SπRIT：a time-projection chamber for symmetry-energy studies［J］. Nuclear Instruments and Methods in Physics Research A，2015， 784：513-517. DOI：10.1016/j.nima.2015.01.026

6Zhan W L， Xu H S， Sun Z Y， et al. Present status of HIRFL in Lanzhou［J］. International Journal of Modern Physics E， 2006， 15：1941-1956. DOI：10.1142/ S0218301306005526

7Xia J W， Zhan W L， Wei B W， et al. The heavy ion cooler-storage-ring project （HIRFL-CSR） at Lanzhou［J］. Nuclear Instruments and Methods in Physics Research A，2002， 488：11-25. DOI：10.1016/S0168-9002（02）00475-8

8Xiao Z G， Chen L W. Nuclear matter at a HIRFL-CSR energy regime［J］. Journal of Physics G， 2009， 36：064040. DOI：10.1088/0954-3899/36/6/064040

9Lu C G， Duan L M， Xu H S， et al. Test and simulation of a MICROMEGAS detector［J］. Chinese Physics C， 2011， 35：1033-1036. DOI：10.1088/1674-1137/35/11/010

10 Ji C S， Shao M， Zhang H， et al. Prospects for searching the η→e+e- rare decay at the CSR［J］. Chinese Physics C，2013， 37：046201. DOI：10.1088/1674-1137/37/4/046201

11 Xiao Z G， Yong G C， Chen L W， et al. Probing nuclear symmetry energy at high densities using pion， kaon， eta and photon productions in heavy-ion collisions［J］. European Physical Journal A， 2014， 50（2）：1-10. DOI：10.1140/epja/ i2014-14037-6

12 汪鹏飞， 李占奎， 李海霞， 等. HIRFL-CSR外靶终端上硅微条探测器阵列的搭建［J］. 原子核物理评论， 2014， 31（1）：63-68. DOI：10.11804/NuclPhysRev.31.01.063

WANG Pengfei， LI Zhankui， LI Haixia， et al. Build-up of the silicon micro-strip detector array in ETF of HIRFL-CSR［J］. Nuclear Physics Review， 2014， 31（1）：63-68. DOI：10.11804/NuclPhysRev.31.01.063

13 赵兴文， 苏弘， 千奕， 等. 基于 ASIC芯片的硅微条阵列探测器多通道前端电子学模块设计［J］. 原子核物理评 论 ， 2014， 31（4）：499-504. DOI：10.11804/ NuclPhysRev.31.04.499

ZHAO Xingwen， SU Hong， QIAN Yi， et al. Development of a multi-channel front-end electronics module based on ASIC for silicon strip array detectors［J］. Nuclear Physics Review， 2014， 31（4）：499-504. DOI：10.11804/ NuclPhysRev.31.04.499

14 康龙飞， 赵雷， 李敏， 等. HIRFL-CSR 外靶实验读出电子学预研系统［J］. 原子能科学技术， 2015， 49（1）：154-161. DOI：10.7538/yzk.2015.49.01.0154

KANG Longfei， ZHAO Lei， LI Min， et al. Prototype readout electronics system of external experiment in HIRFL-CSR［J］. Atomic Energy Science and Technology， 2015， 49（1）：154-161. DOI：10.7538/yzk.2015.49.01. 0154

15 Zhao L， Kang L F， Zhou J W， et al. A 16-channel high-resolution time and charge measurement module for the external target experiment in the CSR of HIRFL［J］.

Nuclear Science and Techniques， 2014， 25（1）：010401. DOI：10.13538/j.1001-8042/nst.25.010401

16 Kang L F， Zhao L， Zhou J W， et al. A 128-channel high precision time measurement module［J］. Metrology and Measurement Systems， 2013， 20（2）：275-286

17 Liu G H， Ma Y G， Cai X Z， et al. Azimuthal asymmetry of direct photons in intermediate energy heavy-ion collisions［J］. Physics Letters B， 2008， 663：312-316. DOI：10.1016/j.physletb.2008.04.037

18 Li B A， Chen L W， Ko C M. Recent progress and new challenges in isospin physics with heavy-ion reactions［J］. Physics Reports， 2008， 464：113-281. DOI：10.1016/ j.physrep.2008.04.005

19 Danielewicz P， Lacey R， Lynch W G. Determination of the equation of state of dense matter［J］. Science， 2002， 298：1952-1956. DOI：10.1126/science.1078070

20 Lattimer J M， Parkash M. The physics of neutron stars［J］. Science， 2004， 304：536-542. DOI：10.1126/science. 1090720

21 López J A， Ramírez-Homs E. Effect of an electron gas on a neutron-rich nuclear pasta［J］. Nuclear Science and Techniques， 2015， 26：S20502. DOI：10.13538/ j.1001-8042/nst.26.S20502

22 Ma Y G. Application of information theory in nuclear liquid gas phase transition［J］. Physical Review Letters，1999， 83：3617-3620. DOI：10.1103/PhysRevLett.83. 3617

23 Ma Y G， Fang D Q， Sun X Y， et al. Different mechanism of two-proton emission from proton-rich nuclei23Al and22Mg［J］. Physics Letters B， 2015， 743：306-309. DOI：10.1016/j.physletb.2015.02.066

24 Ma Y G， Natowitz J B， Wada R， et al. Critical behavior in light nuclear systems：experimental aspects［J］. Physical Review C， 2005， 71：054606. DOI：10.1103/PhysRevC. 71.054606

25 He W B， Ma Y G， Cao X G， et al. Giant dipole resonance as a fingerprint of α clustering configurations in12C and16O［J］. Physics Review Letter， 2014， 113：032506. DOI：10.1103/PhysRevLett.113.032506

26 Kanada-En'yo Y， Kimura M， Kobayashi F， et al. Cluster structures in stable and unstable nuclei［J］. Nuclear Science and Techniques， 2015， 26：S20501. DOI：10.13538/j.1001-8042/nst.26.S20501

27 李阔昂， 叶沿林. 放射性束物理实验发展现状［J］. 核技术， 2014， 37（10）：100501. DOI：10.11889/j.0253-3219. 2014.hjs.37.100501

LI Kuoang， YE Yanlin. Recent development in experimental RIB physics［J］. Nuclear Techniques， 2014，37（10）：100501. DOI：10.11889/j.0253-3219.2014.hjs. 37.100501

28 张肇文， 申虹. 非均匀核物质状态方程研究［J］. 核技术，2014， 37（10）：100514. DOI：10.11889/j.0253-3219.2014. hjs.37.100514

ZHANG Zhaowen， SHEN Hong. Relativistic equation of state for non-uniform nuclear matter［J］. Nuclear Techniques， 2014， 37（10）：100514. DOI：10.11889/ j.0253-3219.2014.hjs.37.100514

29 Wieman H H， Anderssen E， Greiner L， et al. STAR PIXEL detector mechanical design［J］. JINST4， 2009， 4（5）：5-15. DOI：10.1088/1748-0221/4/05/P05015

30 Wang Y， Cheng J P， Li Y J， et al. An MRPC for fast neutron detection［J］. Chinese Physics C， 2010， 34（1）：88-91. DOI：10.1088/1674-1137/34/1/016

31 Li B A， Ko C M. Formation of superdense hadronic matter in high energy heavy ion collisions［J］. Physical Review C， 1995， 52：20-37. DOI：10.1103/PhysRevC. 52.2037

32 Antcheva I， Ballintijn M， Bellenot B， et al. ROOT - a C++ framework for petabyte data storage statistical analysis andvisualization［J］. Computer Physics Communications，2009， 180（12）：2499-2512. DOI：10.1016/j.cpc.2009.08. 005

33 Gostinelli S， Allison J， Amako K， et al. Geant4 - a simulation toolkit［J］．Nuclear Instruments and Methods in Physics Research A， 2003， 506：250-303. DOI：10.1016/ S0168-9002（03）01368-8

34 曹良俊， 李玉兰， 来永芳， 等. 基于 GEM读出的TPC工作气体的研究［J］. 高能物理与核物理， 2007， 31（5）：475-480. DOI：10.3321/j.issn：0254-3052.2007.05.011

CAO Liangjun， LI Yulan， LAI Yongfang， et al. Study of gaspropertiesforGEM-basedTPC［J］.High Energy Physics and Nuclear Physics， 2007， 31：475-480. DOI：10.3321/j.issn：0254-3052.2007.05.011

35 Wang B L， Liu Q， Liu H B， et al. Ion transportation study for thick gas electron multipliers［J］. Chinese Physics Letters， 2014， 31（12）：35-39. DOI：10.1088/0256-307X/ 31/12/122901

36 Fruhwirth R. Application of Kalman Filtering to track and vertex fitting［J］. Nuclear Instruments and Methods in Physics Research A， 1987， 262：444-450. DOI：10.1016/ 0168-9002（87）90887-4

Simulation of momentum resolution of the CEE-TPC in HIRFL

LI He1，2ZHANG Song1LU Fei1ZHONG Chen1MA Yugang1
1（Shanghai Institute of Applied Physics， Chinese Academy of Sciences， Jiading Campus， Shanghai 201800， China）
2（University of Chinese Academy of Sciences， Beijing 100049， China）

Background: Time Projection Chamber （TPC） is a kind of high resolution gas detectors developed on the basis of Multi Wire Proportional Chamber （MWPC） and Multi Wire Drift Chamber （MWDC）. TPC has an excellent position resolution that can achieve high momentum resolution. Prior to the CEE （Cooler-Storage-Ring （CSR）External-target Experiment）-TPC construction at Heavy Ion Research Facility in Lanzhou （HIRFL）， it is necessary to simulate and estimate the performance of this detector. Purpose: This study aims to estimate the momentum resolution of the TPC via a simulation according to designed parameters. Methods: A series of simulation codes including Geant4， Kalman Filter and other programs were applied to the evaluation of momentum resolution for CEE-TPC using corresponding parameters. Three particles， i.e.， π meson， proton and deuteron， were used for simulation study and other influence factors on momentum resolution were analyzed. Results: It is found that the typical momentum resolution of these three particles is about 5%. Conclusion: According to the simulation results，CEE-TPC can meet the requirements in the heavy ion collisions experiment at CSR to identify three typical particles，i.e.， π meson， proton and deuteron. It can provide referable technical parameters for the construction of the future CEE-TPC.

LI He， male， born in 1990， graduated from Lanzhou University in 2009， master student， major in nuclear energy science and engineering

TPC， Momentum resolution， Monte Carlo simulation

O571

10.11889/j.0253-3219.2016.hjs.39.070401

国家自然科学基金（No.U1332129、No.U11322547、No.11421505）、国家重点基础研究发展计划项目（No.2014CB845400）资助

李贺，男，1990年出生，2009年毕业于兰州大学，现为硕士研究生，核能科学与工程专业

卢飞，E-mail：lufei@sinap.ac.cn

Supported by National Natural Science Foundation of China （No.U1332129， No.U11322547， No.11421505）， the Major State Basic Research Development Program in China （No.2014CB845400）

LU Fei， E-mail：lufei@sinap.ac.cn

2016-04-12，

2016-04-27