CN105572721A - 用于测量传感器增益的装置、设备及方法 - Google Patents

用于测量传感器增益的装置、设备及方法 Download PDF

Info

Publication number
CN105572721A
CN105572721A CN201510953109.1A CN201510953109A CN105572721A CN 105572721 A CN105572721 A CN 105572721A CN 201510953109 A CN201510953109 A CN 201510953109A CN 105572721 A CN105572721 A CN 105572721A
Authority
CN
China
Prior art keywords
energy
clandestine
action part
gain
illicit action
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510953109.1A
Other languages
English (en)
Other versions
CN105572721B (zh
Inventor
许剑锋
赵指向
黄秋
龚政
彭旗宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sino science and Technology (Shenzhen) Co., Ltd.
Original Assignee
Wuhan Zhongpai Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan Zhongpai Technology Co Ltd filed Critical Wuhan Zhongpai Technology Co Ltd
Priority to CN201510953109.1A priority Critical patent/CN105572721B/zh
Publication of CN105572721A publication Critical patent/CN105572721A/zh
Application granted granted Critical
Publication of CN105572721B publication Critical patent/CN105572721B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • G01T7/005—Details of radiation-measuring instruments calibration techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/20—Measuring radiation intensity with scintillation detectors
    • G01T1/2018—Scintillation-photodiode combinations
    • G01T1/20184—Detector read-out circuitry, e.g. for clearing of traps, compensating for traps or compensating for direct hits
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/20—Measuring radiation intensity with scintillation detectors
    • G01T1/208—Circuits specially adapted for scintillation detectors, e.g. for the photo-multiplier section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/24—Measuring radiation intensity with semiconductor detectors
    • G01T1/248—Silicon photomultipliers [SiPM], e.g. an avalanche photodiode [APD] array on a common Si substrate

Landscapes

  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Measurement Of Radiation (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

本发明提供一种用于测量传感器增益的装置、设备及方法。该装置包括电流检测电路和处理电路。电流检测电路的输入端连接传感器单元的输出端。电流检测电路用于检测传感器单元输出的电流信号并生成相应的检测信号。处理电路的输入端连接电流检测电路的输出端。处理电路用于根据检测信号计算在传感器单元中发生的暗事件的能量，生成暗事件的能量谱图，并基于能量谱图计算传感器单元的增益。根据本发明提供的用于测量传感器增益的装置、设备及方法，通过传感器单元中的暗事件的能量来确定传感器单元的增益，可以快速、准确、高效地确定传感器单元的增益，有利于后续对传感器单元的增益变化进行补偿。

Description

用于测量传感器增益的装置、设备及方法
技术领域
本发明涉及电路领域，具体地，涉及一种用于测量传感器增益的装置、设备及方法。
背景技术
在高能光子(X射线、伽玛光子等)测量系统中，经常采用诸如硅光电倍增管(SiliconPhotomultipliers，SiPM)的光电传感器。SiPM是一种基于硅的光电传感器。SiPM由边长10～100微米左右的小的传感器微元(cell)组成。每个传感器微元都是工作在盖革(Geiger)模式下的雪崩式光电二极管。每个传感器微元每次都只能检测一个可见光子。成百上千的传感器微元组成传感器单元(pixel)。传感器单元的尺寸通常为1平方毫米至几十平方毫米。很多传感器单元组合在一起，又可以组成更大的传感器阵列(例如16x16个3毫米x3毫米的传感器单元组成的阵列)。SiPM阵列和闪烁晶体阵列通过光导层耦合在一起，就构成了基于SiPM的前端检测器。基于SiPM的前端检测器广泛应用于高能光子的检测。相比于传统的光电倍增管(Photomultipliers，PMT)，SiPM具有尺寸小、偏置电压低、时间分辨率高、与核磁共振(MRI)磁场兼容等优点。SiPM的一个缺点是其增益受温度影响较大。SiPM增益的变化，对基于SiPM的前端检测器的许多性能有很大影响，并且会影响高能光子的测量结果。为了补偿SiPM的增益变化带来的影响，可以首先测量出SiPM的增益，以获知当前SiPM的增益变化情况。
因此，需要提供一种用于测量传感器增益的装置，以至少部分地解决现有技术中存在的上述问题。
发明内容
为了至少部分地解决现有技术中存在的问题，根据本发明的一个方面，提供一种用于测量传感器增益的装置。该装置包括电流检测电路和处理电路。电流检测电路的输入端连接传感器单元的输出端。电流检测电路用于检测传感器单元输出的电流信号并生成相应的检测信号。处理电路的输入端连接电流检测电路的输出端。处理电路用于根据检测信号计算在传感器单元中发生的暗事件的能量，生成暗事件的能量谱图，并基于能量谱图计算传感器单元的增益。
根据本发明的另一方面，提供一种用于测量传感器增益的设备。该设备包括与传感器阵列中的多个传感器单元一一对应的多个如上所述的用于测量传感器增益的装置。
根据本发明的另一方面，提供一种用于测量传感器增益的方法。该方法包括：检测传感器单元输出的电流信号并生成相应的检测信号；根据检测信号计算在传感器单元中发生的暗事件的能量；生成暗事件的能量谱图；以及基于能量谱图计算传感器单元的增益。
根据本发明提供的用于测量传感器增益的装置、设备及方法，通过传感器单元中的暗事件的能量来确定传感器单元的增益，由于传感器单元中的暗事件率很高，因此可以在较短的时间内获得较大的参考数据量，从而可以快速、准确、高效地确定传感器单元的增益，有利于后续对传感器单元的增益变化进行补偿。
在发明内容中引入了一系列简化的概念，这些概念将在具体实施方式部分中进一步详细说明。本发明内容部分并不意味着要试图限定所要求保护的技术方案的关键特征和必要技术特征，更不意味着试图确定所要求保护的技术方案的保护范围。
以下结合附图，详细说明本发明的优点和特征。
附图说明
本发明的下列附图在此作为本发明的一部分用于理解本发明。附图中示出了本发明的实施方式及其描述，用来解释本发明的原理。在附图中，
图1示出根据本发明一个实施例的传感器单元以及用于测量传感器增益的装置的示意性框图；
图2示出根据本发明一个实施例的电流检测电路的示意性框图；
图3示出根据本发明一个实施例的电流检测电路的电路示意图；
图4示出根据本发明一个实施例的暗事件的能量谱图；
图5示出根据本发明另一个实施例的暗事件的能量谱图；
图6示出根据本发明又一个实施例的暗事件的能量谱图；以及
图7a示出根据图6所示实施例测量得到的某特定型号的SiPM的增益随过压的变化曲线；
图7b示出根据图6所示实施例测量得到的该特定型号的SiPM的暗事件率随过压的变化曲线；
图8示出根据本发明一个实施例的传感器阵列以及用于测量传感器增益的设备的示意性框图；以及
图9示出根据本发明一个实施例的用于测量传感器增益的方法的流程图。
具体实施方式
在下文的描述中，提供了大量的细节以便能够彻底地理解本发明。然而，本领域技术人员可以了解，如下描述仅涉及本发明的较佳实施例，本发明可以无需一个或多个这样的细节而得以实施。此外，为了避免与本发明发生混淆，对于本领域公知的一些技术特征未进行描述。
根据本发明的一个方面，提供一种用于测量传感器增益的装置。图1示出根据本发明一个实施例的传感器单元110以及用于测量传感器增益的装置120的示意性框图。如图1所示，装置120包括电流检测电路122和处理电路124。
电流检测电路122的输入端连接传感器单元110的输出端。电流检测电路122用于检测传感器单元110输出的电流信号并生成相应的检测信号。
电流检测电路122可以是任何合适的能够检测传感器单元110输出的电流信号的电路。例如，电流检测电路122可以包括示波器，用于检测电流信号并对电流信号进行一些期望处理。可以理解的是，电流检测电路122检测的是传感器单元110在增益测量时段内输出的电流信号。在该增益测量时段内，可能发生有效事件或暗事件，也可能未发生任何事件。在没有事件发生的时段内，传感器单元输出的电流信号为0，电流检测电路122生成的检测信号也是0。
在本文中，有效事件是指高能光子(例如伽玛光子等)在与传感器单元相连的闪烁晶体中作用而引起的在传感器单元中产生电流信号的事件，暗事件是指噪声(通常是热电子)引起的在传感器单元中产生电流信号的事件。在发生有效事件或暗事件时，传感器单元110可以输出一个脉冲电流信号。为了描述方便，在本文中，将在发生有效事件时传感器单元110输出的脉冲电流信号称为有效电流信号，将在发生暗事件时传感器单元110输出的脉冲电流信号称为暗电流信号。有效电流信号的能量远大于暗电流信号的能量，前者通常是后者的几十至几千倍。因此，通过分析传感器单元110输出的电流信号的能量可以确定在传感器单元110中发生的是有效事件还是暗事件。
处理电路124的输入端连接电流检测电路122的输出端。处理电路124用于根据检测信号计算在传感器单元中发生的暗事件的能量，生成暗事件的能量谱图，并基于能量谱图计算传感器单元的增益。
光电传感器(尤其是SiPM)的增益(或增益及串扰率(Crosstalkrate))与温度以及施加到光电传感器上的偏置电压相关。在光电传感器正常工作时，偏置电压通常是不变的，因此其增益(或增益及串扰率)主要受温度影响。以SiPM为例，温度升高时，SiPM的增益(或增益及串扰率)降低，单位时间内暗事件的个数(暗电流的大小)增大。可以通过检测单个暗事件的能量来直接确定SiPM的增益变化情况。在一个示例中，可以根据经验或理论预先设定一个标准暗事件能量。例如，可以将在25度温度以及30V偏置电压下在SiPM内发生的一个暗事件的平均能量设定为标准暗事件能量，将该情况下的增益设定为标准增益。处理电路124通过分析暗事件的能量谱图可以获知在增益测量时段内发生的一个暗事件的能量相对于标准暗事件能量的变化，从而获知在增益测量时段内SiPM的增益相对于标准增益的变化。之后，可以基于增益的变化情况制定合适的补偿策略，以补偿增益变化带来的影响。
根据本发明提供的用于测量传感器增益的装置，通过传感器单元中的暗事件的能量来确定传感器单元的增益，由于传感器单元中的暗事件率很高，因此可以在较短的时间内获得较大的参考数据量，从而可以快速、准确、高效地确定传感器单元的增益，有利于后续对传感器单元的增益变化进行补偿。
可选地，检测信号是数字信号。数字信号由持续时间相等的高电平和低电平组成，数字信号中的所有高电平之和与电流信号对时间的积分成正比。在一个示例中，电流检测电路可以实现为图2所示的形式，以生成上述数字信号。图2示出根据本发明一个实施例的电流检测电路200的示意性框图。
电流检测电路200可以包括积分电路210、比较电路220、传输控制电路230和负反馈电路240。
积分电路210连接传感器单元(如图1所示的传感器单元110)的输出端和负反馈电路240的输出端。积分电路210用于接收来自传感器单元的电流信号和来自负反馈电路240的反馈信号，对电流信号和反馈信号的差进行积分并且输出积分信号。
电流检测电路200是包括负反馈环节的电路，反馈信号被输入到积分电路210。同时，积分电路210还接收来自传感器单元的电流信号。电流信号和反馈信号的流动方向是相反的。例如，如果电流信号是从积分电路210流向传感器单元的，则可以将反馈信号设定为从负反馈电路240流向积分电路210。因此，对于积分电路210来说，实际上最终输入的是电流信号与反馈信号的差，积分电路210可以对二者的差进行积分。应该注意，图2中示出的箭头方向是信号的传输方向，而不一定是信号的流动方向。积分电路210可以采用模拟积分电路实现，例如通过电阻和电容组成的滤波器电路实现。
如图2所示，比较电路220的一个输入端连接积分电路210的输出端并且比较电路220的另一输入端接入参考电平。比较电路220用于将积分信号与参考电平进行比较并生成比较结果。例如，当积分信号的幅度大于参考电平时，比较电路220可以输出高电平，当积分信号的幅度等于或小于参考电平时，比较电路220可以输出低电平。因此，比较电路220输出的比较结果中可以只存在高电平和低电平两种状态。通常，传感器单元输出的电流信号是随时间变化的脉冲电流信号，在这种情况下，积分信号也是随时间变化的信号。因此，比较电路220输出的比较结果是随时间变化而在高电平和低电平两种状态之间切换的信号。
传输控制电路230的输入端连接比较电路220的输出端。传输控制电路230用于利用时钟信号控制比较结果的传输以输出数字信号。如上所述，比较结果是随时间变化而在高电平和低电平两种状态之间切换的信号。在该比较结果中，高电平和低电平的持续时间可能是实时变化的，是无法确定的。所以可以通过传输控制电路230对比较结果进行时间上的量化，使得每段连续的高电平或低电平的持续时间都是时钟信号的周期的整数倍。这种时间上的量化相当于模数转换过程中的时间离散化，因此，从功能性上来看，可以将比较电路220和传输控制电路230这二者视作一个1位的ADC。在传输控制电路230输出的数字信号中，每个高电平和每个低电平的持续时间相等，均等于时钟信号的周期。在一个示例中，高电平可以代表逻辑电平“1”，低电平可以代表逻辑电平“0”，则数字信号是由逻辑电平“1”和“0”组成的序列。假设时钟信号的频率为100Hz，即周期为0.01s，则在数字信号中，单个“1”或“0”的持续时间是0.01s。另外，可以理解的是，当多个“1”或多个“0”连续出现时，该多个“1”或多个“0”的持续时间是0.01s的整数倍。传输控制电路230可以是寄存器或受时钟信号控制的开关电路等。
负反馈电路240的输入端连接传输控制电路230的输出端。负反馈电路240用于将数字信号转换为反馈信号并且将反馈信号反馈给积分电路210。负反馈电路240可以包括数模转换器(DAC)，用于对数字信号进行数模转换以将其转换为模拟信号。具体地，该DAC可以是1位的DAC，以将传输控制电路230输出的由“1”和“0”组成的序列转换为模拟信号，例如转换为幅度随时间变化的电压信号。负反馈电路240可以进一步包括电阻。DAC经由该电阻连接到积分电路210的输入端。因此，上述电压信号将产生流经该电阻的电流信号，即反馈信号。反馈信号可以抵消来自传感器单元的电流信号的一部分，避免积分电路210所输出的积分信号过大，以保持电路稳定。此外，由于反馈信号在不断抵消来自传感器单元的电流信号，所以电流信号引发的反馈信号的累加值可以视作来自传感器单元的电流信号的累加值。同时，反馈信号的幅度与数字信号的幅度成正比。因此，当电流信号的持续时间已经结束并且反馈信号的幅度稳定在零(即针对电流信号的负反馈作用已经停止)时，可以利用数字信号来计算光子的能量。应该注意，反馈信号不宜过大或过小。反馈信号过大会导致电流信号的抵消速度过快，使得数字信号包含的误差增大，影响测量精度。相反，反馈信号过小会导致电流信号的抵消速度过慢，使得无法及时减小积分信号的值，从而也会影响测量精度。反馈信号的幅度可以根据实际需要来确定，本发明不对此进行限制。
图3示出根据本发明一个实施例的电流检测电路300的电路示意图。电流检测电路300可以包括积分电路310、比较电路320、传输控制电路330和负反馈电路340。
积分电路310可以采用运算放大器、电阻Rin和电容C来实现。比较电路320可以采用比较器实现，其参考电平是地电平。传输控制电路330可以采用D触发器构成的寄存器来实现，其时钟信号的频率是fs’。负反馈电路340可以采用一个DAC和一个电阻Rf来实现。该DAC可以是现场可编程门阵列(FPGA)的一个引脚。整个电流检测电路300可以由模拟电路部分和FPGA部分来实现，如图3所示。根据需要，电路中也可以加入二极管，或者分压电路来减少或者隔绝电路中可能存在的偏置直流电流。本领域技术人员根据以上对电流检测电路200的描述，可以理解电流检测电路300的原理和工作方式，不再赘述。
生成数字化的检测信号可以方便后续对有效事件或暗事件的能量进行计算。
可选地，处理电路124可以通过以下方式计算在传感器单元110中发生的暗事件的能量：根据数字信号中的高电平和低电平的出现规律确定脉冲事件是否发生；当确定脉冲事件发生时开始利用数字信号计算脉冲事件的能量，直至根据数字信号中的高电平和低电平的出现规律确定脉冲事件结束为止；以及根据脉冲事件的能量判断脉冲事件是有效事件还是暗事件，如果脉冲事件是暗事件，则将脉冲事件的能量确定为暗事件的能量。
脉冲事件是指在传感器单元中产生脉冲电流信号的事件，其可能是有效事件，也可能是暗事件。如上文所述，数字信号是“0”和“1”组成的序列。当传感器单元中未发生有效事件或暗事件时，电流检测电路122、200或300输出的是一系列连续的“0”，而在发生有效事件或暗事件时，电流检测电路122、200或300输出的信号中应当包含若干个“1”。可以总结脉冲事件发生和脉冲事件未发生时的“1”和“0”的出现规律，并根据总结出的规律来区分脉冲事件发生和脉冲事件未发生这两种情况。例如，如果在数字信号中，在很多个(例如10个)连续的“0”后面出现了一个“1”，则说明此时发生了脉冲事件，其有可能是有效事件，也有可能是暗事件。随后，可以开始利用数字信号中的“1”计算脉冲事件的能量。例如，可以对数字信号中的“1”进行累加或对其进行计数。当数字信号中连续出现5个“0”时，说明脉冲事件结束，此时可以停止对“1”进行累加或计数。这样，最后可以获得脉冲事件的相对能量，其可以反映脉冲事件的绝对能量的大小。如上文所述，有效事件的能量与暗事件的能量相差较大，因此可以通过能量来区分两者。例如，如果最后获得的脉冲事件的相对能量是几千，则可以认为脉冲事件是有效事件，此时可以将计算出的能量丢弃。如果最后获得的脉冲事件的相对能量是几百，则可以认为脉冲事件是暗事件，此时可以将计算出的能量作为暗事件的能量并将其用于生成能量谱图。这种测量暗事件能量的方式可以视为是触发方式，即当有暗事件触发时才测量暗事件的能量。
可选地，能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数。处理电路124可以通过以下方式计算传感器单元110的增益：根据能量谱图中的、与单个暗事件的能量相对应的第一谱峰的横坐标计算增益，或根据第一谱峰的横坐标和能量谱图中的、与两个暗事件的能量相对应的第二谱峰的横坐标之间的差计算增益。
根据上文所述的方法，可以在每次有暗事件发生时，都计算获得该暗事件的一个相对能量值。可以在测量时段(例如5秒钟)内，不停地进行暗事件的检测和能量计算，并针对在该测量时段内检测到的所有暗事件，进行谱分析，得到能量谱图。图4示出根据本发明一个实施例的暗事件的能量谱图。图4所示的能量谱图中具有两个较明显的谱峰，左边的谱峰是普通暗事件的能量谱峰(即第一谱峰)，右边的谱峰是存在串扰(crosstalk)的暗事件的能量谱峰(即第二谱峰)。第二谱峰的横轴位置差不多是第一谱峰的两倍。具体地，第一谱峰在能量为大约40的位置处，40代表单个暗事件的能量，因此，第一谱峰与单个暗事件的能量相对应。第二谱峰在能量为大约80的位置处，80代表两个暗事件的能量，因此，第二谱峰与两个暗事件的能量相对应。根据上文描述可知，处理电路124在每次计算能量时，基于从暗事件开始至暗事件结束的时间内的数字信号来计算暗事件的能量，因此其每次获得的基本是单个暗事件的能量。
暗事件能量的增益即代表传感器单元的增益。通过第一谱峰的横坐标或第一谱峰的横坐标与第二谱峰的横坐标之间的差可以获知当前单个暗事件的能量。通过当前单个暗事件的能量与标准暗事件能量的差异可以计算增益。例如，假设标准暗事件的能量是30，标准增益是100。根据如图4所示的谱图，已知当前单个暗事件的能量为接近40。这样，传感器单元的当前增益可以认为是标准增益的接近1.33倍，即当前增益为133。可以理解的是，由于干扰等因素存在，第一谱峰和第二谱峰的位置可能不准确，例如可能向左或向右移动。然而，由于第一谱峰与第二谱峰是同步移动的，因此二者的横坐标之间的差始终是不变的，因此与直接利用第一谱峰的横坐标来计算增益的方式相比，利用第一谱峰的横坐标与第二谱峰的横坐标之间的差来计算增益的方式的准确度更高。
可选地，处理电路124可以通过以下方式计算在传感器单元110中发生的暗事件的能量：根据数字信号中的高电平和低电平的出现规律确定有效事件未发生的时间段；从在有效事件未发生的时间段内生成的数字信号中选择在预定时段内生成的数字信号；以及根据在预定时段内生成的数字信号计算暗事件的能量。
可以首先根据数字信号中的“1”和“0”的出现规律确定有效事件发生的时间段，进而可以确定有效事件未发生的时间段。例如，如果一段数字信号中的“1”出现了5000次，可以认为有效事件发生。在与其相邻的另一段数字信号中的“1”出现了4000次，也可以认为有效事件发生。在这两个有效事件之间，可能发生了多次暗事件。因此，可以在两个有效事件之间，选择一个或多个预定时段，并选择在该预定时段内生成的数字信号。该预定时段可以是任何合适的时间段，例如10微秒，本发明不对此进行限制。例如，如果两个有效事件之间的整个时间段长度是1毫秒，则可以将这1毫秒分为100个10微秒。分别计算每个10微秒内包含的暗事件的能量，可以得到100个暗事件能量值。可以在增益测量时段(例如5秒)内持续进行这样的有效事件检测和暗事件能量计算。在这种方式中，不采用触发方式来启动暗事件的能量测量，而是计算未发生有效事件(仅有暗事件)的预定时段内所有暗事件的能量，以用于生成能量谱图。
可选地，预定时段设定为使得在预定时段内的平均暗事件个数满足第一预设条件。能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数，能量谱图包括多个谱峰并且能量谱图中的第N个谱峰与N个暗事件的能量相对应。处理电路124可以通过以下方式计算传感器单元的增益：根据能量谱图中的、任意两个相邻谱峰的横坐标之间的差计算增益。可选地，第一预设条件是平均暗事件个数小于20。当平均暗事件个数小于20时，可以使得能量谱图中出现多个离散的谱峰，即能量谱图呈现泊松分布。当然，可以理解的是，使能量谱图呈现泊松分布的平均暗事件个数也可以具有其他合适的选值范围，本发明不对此进行限制。下面对此进行详细描述。
根据以上描述可知，处理电路124每次计算预定时段内的暗事件的能量。如果该预定时段设定得合适，则可以使得在每个预定时段内发生多个暗事件，也就使得每次计算得到的暗事件的能量等于多个暗事件的能量之和。图5示出根据本发明另一个实施例的暗事件的能量谱图。如图5所示，该能量谱图中具有比图4更多的谱峰，并且高度最高的谱峰也不再是第一谱峰。这是因为，时间段设定得较长，使得每个时间段内发生多个暗事件的能量的可能性远大于只发生单个暗事件的可能性，从而使得每次计算得到的暗事件的能量总是大于单个暗事件的能量。平均暗事件个数是指在(增益测量时段内的)所选择的所有预定时段内发生的暗事件个数相对预定时段个数的平均。
图5所述的能量谱图包括如下几个谱峰。第一个谱峰的横坐标代表的是(预定时段内发生的)一个暗事件的能量；第二个谱峰的横坐标代表的是(预定时段内发生的)两个暗事件的能量；依此类推，第N个谱峰的横坐标代表的是(预定时段内发生的)N个暗事件的能量。第N个谱峰的横坐标基本是第一个谱峰的N倍。暗事件能量的增益即为传感器单元的增益。当温度变化引起增益变化时，可能导致图5所示的所有谱峰向右侧移动。通过跟踪谱峰中的一个(或者几个)的位置变动，可以获知传感器单元的增益变化。这种增益测量方式的优点如下所述。由于第N个谱峰的横坐标是第一个谱峰的横坐标的N倍，因此温度的变化导致的第N个谱峰的横坐标的移动，比第一个谱峰的横坐标的移动要大N倍。因此，通过测算第N个谱峰的横坐标的移动，来测算温度导致的增益变化，比通过测算第一个谱峰的横坐标的移动，来测算温度导致的增益变化，灵敏度和准确度要高一些。另外，由于第一个谱峰不太明显，因此也可以根据任意两个相邻谱峰的横坐标之间的差计算增益。这种方式与上文所述的通过第一谱峰的横坐标与第二谱峰的横坐标之间的差来计算增益的方式类似，不再赘述。
可选地，预定时段设定为使得在预定时段内的平均暗事件个数满足第二预设条件。能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数。处理电路124可以通过以下方式计算传感器单元的增益：对能量谱图进行高斯拟合以获得能量谱图的方差和均值，并且将方差除以均值以获得增益。可选地，第二预设条件是平均暗事件个数大于20。当平均暗事件个数大于20时，可以使得能量谱图近似呈现高斯分布。当然，可以理解的是，使能量谱图呈现高斯分布的平均暗事件个数也可以具有其他合适的选值范围，本发明不对此进行限制。下面对此进行详细描述。
图6示出根据本发明又一个实施例的暗事件的能量谱图。与图5所示实施例不同的是，图6所示实施例的能量累积时间(即所谓的预定时段)比图5所示实施例长很多，从而导致得到的能量谱图中，不再出现离散的谱峰，整个能量谱图呈现如图6所示的近似正态分布(即高斯分布)。
产生这一现象的基本原理如下：
(1)预定时段内的暗事件个数本质上是一个离散的泊松分布；
(2)当预定时段内的暗事件个数较少时，可能出现离散的谱峰，但由于测量误差，实际的能量谱图中的谱峰会有一定的宽度；
(3)当预定时段内的暗事件个数较多时，泊松分布逐步接近正态分布。当预定时段内的平均暗事件个数λ大于一定的数值(例如20)时，可以认为泊松分布呈正态分布N(λ,λ)，其可以按正态分布处理；
(4)根据本实施例测量得到的能量谱图的横轴不对应于预定时段内的平均暗事件个数λ，而对应于λG(G是传感器单元的增益)。因此，泊松分布呈正态分布N(λ,λ)后，反映在能量谱图上的正态分布为N(λG,λG2)。因此，把这一正态分布的方差除以其均值即可得到传感器单元的增益G。
需要指出的是，这样计算出来的增益实际上是增益G，“1+串扰率”和“1+单光子检测分辨率平方”三项的乘积。在实际情况中，“1+串扰率”和“1+单光子检测分辨率平方”两项的乘积非常接近于1。因此，通常可以认为，增益G、“1+串扰率”和“1+单光子检测分辨率平方”三项的乘积近似等于增益G。另外，需要指出的是，通常情况下，单光子检测分辨率约为0.05左右，“1+单光子检测分辨率平方”近似为1。因此，计算出来的增益也可以近似为增益G和“1+串扰率”的乘积。由于SiPM测量得到的有效事件的能量正比于增益G和“1+串扰率”的乘积，增益补偿需要同时补偿温度对增益G和“1+串扰率”的影响。因此，计算增益G和“1+串扰率”的乘积比单纯地计算真正的增益更有实际价值。
因此，根据本实施例，计算传感器单元的增益的实施方式实际上非常简单和易于实现，基本步骤如下：
(1)对能量谱进行高斯拟合，计算其均值μ及方差σ2。
(2)计算传感器单元的增益G＝σ2/μ。
图7a示出根据图6所示实施例测量得到的某特定型号的SiPM的增益随过压的变化曲线。图7b示出根据图6所示实施例测量得到的该特定型号的SiPM的暗事件率随过压的变化曲线(暗事件率约等于暗电流除以增益)。从图7a和图7b可以看出图6所示实施例的增益测量结果比较准确。
图5所示的实施例一般要求检测暗事件产生的电流(或者电荷)的电流检测电路具备较高的精度。这往往是以牺牲整个电流检测电路的动态范围(或者死区时间)为代价的。图6所示的实施例对电流检测电路的精度要求较低。
可选地，处理电路124可以进一步用于根据检测信号计算在传感器单元中发生的有效事件的能量E0，并且根据以下公式对有效事件的能量E0进行校正：E＝E0*G0/G，其中，E为经校正的能量，G为所述增益，G0为预定常数。
虽然本实施例描述的是利用处理电路124以及电流检测电路122、200或300来测量有效事件的能量的方法，但可以理解的是，也可以利用其他任何合适的装置测量有效事件的能量。
在本文中，传感器单元的增益G可以实时计算获得，其可以是随时间变化的量。G0也可以具有任何合适的值，本文不对此进行限制。
例如，G0可以等于例如上文所述的标准增益。通过将在当前时刻测量得到的传感器单元的增益与标准增益进行比较，可以获知当前时刻测量得到的有效事件的能量与在标准增益下的能量的差异，从而可以将当前时刻测量得到的能量校正为在标准增益下的能量。如果对多个传感器单元的能量均进行这样的校正，则相当于将这多个传感器单元的能量放在同一标准下，因此它们的能量可以进行相互比较。例如，如果标准增益是100，而当前时刻测量得到的增益G为200，说明当前的增益增大了1倍。因此测量得到的能量也会比在标准增益下的能量增大1倍。如果测量得到的能量E0是2000，则经校正的能量E可以是1000。
另外，也可以将传感器单元的能量校正为实际的能量。在这种情况下，G0可以设定为1。假设，如果已根据本文所描述的实施例获知传感器单元的当前增益是1.5，所检测到的某一有效事件的当前能量是3000，则可以将能量3000除以增益1.5，以得到该有效事件的经校正的能量2000。
可以理解的是，处理电路124计算出的有效事件或暗事件的能量可能是相对能量而非绝对能量，因此传感器单元的增益也可以用相对增益而非绝对增益来表示。
例如，返回参考图4。如上文所述，假设标准暗事件的能量是30，标准增益是100。假设在某一增益测量时段(为描述方便，用“第一增益测量时段”来表示)内获得的能量谱图如图4所示，则在第一增益测量时段内的单个暗事件的能量为接近40。这样，在第一增益测量时段内传感器单元的绝对增益可以认为是标准增益的接近1.33倍，即在第一增益测量时段内的绝对增益为133。在第一增益测量时段前后的预设时间范围内，可以认为增益是不变的。在该时间范围内，如果检测到的有效事件的能量为10000，则可以用能量10000除以绝对增益133，获得校正后的能量75。随后，在另一增益测量时段(为描述方便，用“第二增益测量时段”来表示)内，如果能量谱图向右移动，导致计算得出的在第二增益测量时段内的单个暗事件的能量为接近50。在第二增益测量时段内传感器单元的绝对增益可以认为是标准增益的接近1.67倍，即在该第二增益测量时段内的绝对增益为167。在第二增益测量时段前后的预设时间范围内，可以认为增益是不变的。在该时间范围内，如果检测到的有效事件的能量为10000，则可以用能量10000除以绝对增益167，获得校正后的能量60。这样，可以将校正后的能量75与校正后的能量60相比较，可知两个有效事件之间的能量差异。
可选地，也可以直接将40视作在第一增益测量时段内传感器单元的增益。也就是说，无需基于标准暗事件和标准增益来计算绝对增益。在这种情况下，可以用能量10000除以增益40，获得校正后的能量为250。随后，也可以直接将50视作在第二增益测量时段内传感器单元的增益。在这种情况下，可以用能量10000除以增益50，获得校正后的能量为200。可以将校正后的能量250与校正后的能量200相比较，可知两个有效事件之间的能量差异。75/60与250/200的结果是一样的，因此可见无论是计算绝对增益还是相对增益，都可以正确反映各有效事件的能量之间的相对大小。以上两种增益计算方法可以认为是采用了两种不同的归一化方式，它们可以根据需要进行选择。基于增益校正有效事件的能量可以使得最后获得的能量值能够更准确。
可选地，处理电路124可以采用任何合适的数字电路实现，例如采用FPGA、复杂可编程逻辑器件(CPLD)、数字信号处理器(DSP)、微控制单元(MCU)或中央处理单元(CPU)等实现。
本发明具有如下优点：
(1)较低的硬件成本。本发明实施例提供的电路结构简单，电流检测电路可以采用模拟电路实现，处理电路可以采用数字电路实现，整体不需要过多硬件，因此成本较低。
(2)不干扰传感器单元的正常工作。本发明可以在不干扰传感器单元的正常工作状态的情况下，对通常被认为是噪声的暗事件产生的电流(电荷)进行检测和分析，计算传感器单元的增益。
(3)实时和高精度的增益测量和补偿。传感器单元中的暗事件率一般大于300kHz。这使得能够在足够短的时间内采集足够多的数据，保证谱分析的精度，做到实时高精度地测量和补偿温度导致的增益变化。反应速度非常快。
根据本发明的另一方面，提供一种用于测量传感器增益的设备。该设备包括与传感器阵列中的多个传感器单元一一对应的多个上文所述的用于测量传感器增益的装置。图8示出根据本发明一个实施例的传感器阵列以及用于测量传感器增益的设备的示意性框图。
如图8所示，传感器阵列包括N个传感器单元810，即传感器单元1、传感器单元2……传感器单元N。用于测量传感器增益的设备包括与N个传感器单元一一对应的N个用于测量传感器增益的装置820。每个装置820与上文所述的装置120相同，均包括电流检测电路822和处理电路824。本领域技术人员根据以上关于用于测量传感器增益的装置的描述以及结合图1至7b能够理解用于测量传感器增益的设备的结构、运行方式和优点等，不再赘述。
需要注意的是，在用于测量传感器增益的设备中，各处理电路824均可以采用如上所述的诸如FPGA的数字电路来实现。可以理解，用于测量传感器增益的设备中的任意数量的处理电路824可以集成在同一诸如FPGA的数字电路中。
可选地，处理电路824进一步用于根据检测信号计算在对应传感器单元中发生的有效事件的能量E0，并且根据以下公式对有效事件的能量E0进行校正：E＝E0*G0/G，其中，E为经校正的能量，G为对应传感器单元的增益，G0为预定常数。可选地，G0为每个处理电路824的对应传感器单元的标准增益。对某一传感器阵列来说，各传感器单元的标准增益通常相同。
下面结合SiPM进行描述。SiPM单元的增益通常正比于其偏置过压。偏置过压等于偏置电压与击穿电压之差。由于微电子工艺上的原因，不同SiPM单元的击穿电压都有一定的差异。同一系列SiPM单元组成SiPM阵列时，通常采用同一偏置(工作)电压。由于各SiPM单元的击穿电压略有不同，它们的增益也不一样。SiPM阵列中的各SiPM单元之间的增益差异会影响由SiPM阵列组成的前端检测器的性能。通常情况下，需要对组成SiPM阵列的SiPM单元进行测试和筛选，选择击穿电压差异较小的SiPM单元来组成SiPM阵列。这一方面加大了SiPM阵列的制作成本，另一方面SiPM阵列中的各SiPM单元之间的增益也很难做到完全一致。因此，可以实时地测量SiPM阵列中的各SiPM单元的增益，并通过增益补偿的方法，保证各SiPM单元的增益一致。
以上公式E＝E0*G0/G中的“G0/G”实现了对同一传感器阵列中的各传感器单元的增益的归一化，从而可以保证各传感器单元的增益一致。例如，同一传感器阵列中的三个传感器单元的增益可能分别是1，1.5和0.8。因此，对于能量相等的有效事件来说，每个传感器单元检测到的能量将是不同的。然而，所期望的是，这三个传感器单元检测到的能量相同。因此，可以在后续对能量进行校正(也就是对增益变化进行补偿)时，首先利用一个常数G0分别除以1，1.5和0.8，这相当于进行了增益归一化。G0可以根据需要而定。经过能量校正之后，对于整个传感器阵列来说，已经消除了各传感器单元之间的增益差异对能量测量带来的影响。
根据本发明的另一方面，提供一种用于测量传感器增益的方法。图9示出根据本发明一个实施例的用于测量传感器增益的方法900的流程图。如图9所示，方法900包括以下步骤。
步骤S910，检测传感器单元输出的电流信号并生成相应的检测信号。
步骤S920，根据检测信号计算在传感器单元中发生的暗事件的能量。
步骤S930，生成暗事件的能量谱图。
步骤S940，基于能量谱图计算传感器单元的增益。
可选地，检测信号是数字信号，数字信号由持续时间相等的高电平和低电平组成，数字信号中的所有高电平之和与电流信号对时间的积分成正比。
可选地，步骤S920可以包括：根据数字信号中的高电平和低电平的出现规律确定脉冲事件是否发生；当确定脉冲事件发生时开始利用数字信号计算脉冲事件的能量，直至根据数字信号中的高电平和低电平的出现规律确定脉冲事件结束为止；以及根据脉冲事件的能量判断脉冲事件是有效事件还是暗事件，如果脉冲事件是暗事件，则将脉冲事件的能量确定为暗事件的能量。
可选地，能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数。步骤S940可以包括：根据能量谱图中的、与单个暗事件的能量相对应的第一谱峰的横坐标计算增益，或根据第一谱峰的横坐标和能量谱图中的、与两个暗事件的能量相对应的第二谱峰的横坐标之间的差计算增益。
可选地，步骤S920可以包括：根据数字信号中的高电平和低电平的出现规律确定有效事件未发生的时间段；从在有效事件未发生的时间段内生成的数字信号中选择在预定时段内生成的数字信号；以及根据在预定时段内生成的数字信号计算所述暗事件的能量。
可选地，预定时段设定为使得在预定时段内的平均暗事件个数满足第一预设条件，能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数，能量谱图包括多个谱峰并且能量谱图中的第N个谱峰与N个暗事件的能量相对应，步骤S940可以包括：根据能量谱图中的、任意两个相邻谱峰的横坐标之间的差计算增益。
可选地，第一预设条件是平均暗事件个数小于50。
可选地，预定时段设定为使得在预定时段内的平均暗事件个数满足第二预设条件，能量谱图的横轴表示暗事件的能量，能量谱图的纵轴表示计算得到暗事件的能量的次数，步骤S940可以包括：对能量谱图进行高斯拟合以获得能量谱图的方差和均值；以及将方差除以均值以获得增益。
可选地，第二预设条件是平均暗事件个数大于50。
可选地，方法900可以进一步包括：根据检测信号计算在传感器单元中发生的有效事件的能量E0；以及根据以下公式对有效事件的能量E0进行校正：E＝E0*G0/G，其中，E为经校正的能量，G为所述增益，G0为预定常数。
本领域技术人员根据以上关于用于测量传感器增益的装置的描述以及附图1至8，能够理解本文所公开的用于测量传感器增益的方法900的实施方式及其优点等，为了简洁，本文不对此进行赘述。
现有的针对SiPM的增益补偿方法通常包括如下几种：
(1)温度控制的方法，即通过温度测量和控制系统来直接稳定SiPM的温度，进而稳定SiPM的增益。这种方法需要构建恒温系统，成本高昂。很多情况下，由于物理空间的限制，这种方法可行性不高。
(2)温度补偿的方法，即根据温度变化来调控偏置电压，以稳定SiPM的增益。例如，温度升高时，可以提高偏置电压，以稳定增益。理想情况下，这种方法需要对SiPM阵列中的每一个SiPM单元，设计温度测量和电压控制电路。实际应用中，一组物理空间接近的SiPM单元，可以共用一套温度测量和电压控制电路。这种方法应用最为广泛，但成本较高。
(3)谱峰跟踪的方法。这种方法通过寻找和跟踪在有效事件发生时SiPM检测到的能量谱的光峰(photopeak)谱峰位置，来自动校准增益。这种方法的精度，取决于温度漂移的速度，及有效事件率。要在某一温度实现谱峰的精确定位，需要较多的事件来克服统计误差。通常情况下(例如，PET人体成像)，有效事件率非常低，很难实现实时地跟踪photopeak谱峰。
与上述现有的增益补偿方法相比，本发明提供的通过暗电流测量传感器单元的增益进而对增益进行补偿的方式所采用的电路结构简单，成本低，易于实现。并且，由于传感器单元中的暗事件率很高，因此可以在较短的时间内获得较大的参考数据量，从而可以快速、准确、高效地确定传感器单元的增益，有利于后续对传感器单元的增益变化进行补偿。
虽然本文以SiPM为例描述了本发明的原理和应用，但是应该理解的是，本发明并不局限于此。本发明所提供的用于测量传感器增益的装置及方法还可以应用于PMT或任何其他合适的光电传感器。
本发明已经通过上述实施例进行了说明，但应当理解的是，上述实施例只是用于举例和说明的目的，而非意在将本发明限制于所描述的实施例范围内。此外本领域技术人员可以理解的是，本发明并不局限于上述实施例，根据本发明的教导还可以做出更多种的变型和修改，这些变型和修改均落在本发明所要求保护的范围以内。本发明的保护范围由附属的权利要求书及其等效范围所界定。

Claims (22)

1.一种用于测量传感器增益的装置，包括电流检测电路和处理电路，其中：
所述电流检测电路的输入端连接传感器单元的输出端，所述电流检测电路用于检测所述传感器单元输出的电流信号并生成相应的检测信号；
所述处理电路的输入端连接所述电流检测电路的输出端，所述处理电路用于根据所述检测信号计算在所述传感器单元中发生的暗事件的能量，生成所述暗事件的能量谱图，并基于所述能量谱图计算所述传感器单元的增益。
2.根据权利要求1所述的装置，其特征在于，所述检测信号是数字信号，所述数字信号由持续时间相等的高电平和低电平组成，所述数字信号中的所有高电平之和与所述电流信号对时间的积分成正比。
3.根据权利要求2所述的装置，其特征在于，所述处理电路通过以下方式计算在所述传感器单元中发生的暗事件的能量：
根据所述数字信号中的高电平和低电平的出现规律确定脉冲事件是否发生；
当确定所述脉冲事件发生时开始利用所述数字信号计算所述脉冲事件的能量，直至根据所述数字信号中的高电平和低电平的出现规律确定所述脉冲事件结束为止；以及
根据所述脉冲事件的能量判断所述脉冲事件是有效事件还是暗事件，如果所述脉冲事件是暗事件，则将所述脉冲事件的能量确定为所述暗事件的能量。
4.根据权利要求3所述的装置，其特征在于，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，
所述处理电路通过以下方式计算所述传感器单元的增益：
根据所述能量谱图中的、与单个暗事件的能量相对应的第一谱峰的横坐标计算所述增益，或根据所述第一谱峰的横坐标和所述能量谱图中的、与两个暗事件的能量相对应的第二谱峰的横坐标之间的差计算所述增益。
5.根据权利要求2所述的装置，其特征在于，所述处理电路通过以下方式计算在所述传感器单元中发生的暗事件的能量：
根据所述数字信号中的高电平和低电平的出现规律确定有效事件未发生的时间段；
从在有效事件未发生的时间段内生成的数字信号中选择在预定时段内生成的数字信号；以及
根据所述在预定时段内生成的数字信号计算所述暗事件的能量。
6.根据权利要求5所述的装置，其特征在于，所述预定时段设定为使得在所述预定时段内的平均暗事件个数满足第一预设条件，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，所述能量谱图包括多个谱峰并且所述能量谱图中的第N个谱峰与N个暗事件的能量相对应，
所述处理电路通过以下方式计算所述传感器单元的增益：
根据所述能量谱图中的、任意两个相邻谱峰的横坐标之间的差计算所述增益。
7.根据权利要求6所述的装置，其特征在于，所述第一预设条件是所述平均暗事件个数小于20。
8.根据权利要求5所述的装置，其特征在于，所述预定时段设定为使得在所述预定时段内的平均暗事件个数满足第二预设条件，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，
所述处理电路通过以下方式计算所述传感器单元的增益：
对所述能量谱图进行高斯拟合以获得所述能量谱图的方差和均值，并且将所述方差除以所述均值以获得所述增益。
9.根据权利要求8所述的装置，其特征在于，所述第二预设条件是所述平均暗事件个数大于20。
10.根据权利要求1至9任一项所述的装置，其特征在于，所述处理电路进一步用于根据所述检测信号计算在所述传感器单元中发生的有效事件的能量E0，并且根据以下公式对所述有效事件的能量E0进行校正：
E＝E0*G0/G，
其中，E为经校正的能量，G为所述增益，G0为预定常数。
11.一种用于测量传感器增益的设备，包括与传感器阵列中的多个传感器单元一一对应的多个如权利要求1至9任一项所述的用于测量传感器增益的装置。
12.根据权利要求11所述的设备，其特征在于，所述处理电路进一步用于根据所述检测信号计算在对应传感器单元中发生的有效事件的能量E0，并且根据以下公式对所述有效事件的能量E0进行校正：
E＝E0*G0/G，
其中，E为经校正的能量，G为所述对应传感器单元的增益，G0为预定常数。
13.一种用于测量传感器增益的方法，包括：
检测所述传感器单元输出的电流信号并生成相应的检测信号；
根据所述检测信号计算在所述传感器单元中发生的暗事件的能量；
生成所述暗事件的能量谱图；以及
基于所述能量谱图计算所述传感器单元的增益。
14.根据权利要求13所述的方法，其特征在于，所述检测信号是数字信号，所述数字信号由持续时间相等的高电平和低电平组成，所述数字信号中的所有高电平之和与所述电流信号对时间的积分成正比。
15.根据权利要求14所述的方法，其特征在于，所述计算在所述传感器单元中发生的暗事件的能量包括：
根据所述数字信号中的高电平和低电平的出现规律确定脉冲事件是否发生；
当确定所述脉冲事件发生时开始利用所述数字信号计算所述脉冲事件的能量，直至根据所述数字信号中的高电平和低电平的出现规律确定所述脉冲事件结束为止；以及
根据所述脉冲事件的能量判断所述脉冲事件是有效事件还是暗事件，如果所述脉冲事件是暗事件，则将所述脉冲事件的能量确定为所述暗事件的能量。
16.根据权利要求15所述的方法，其特征在于，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，
所述计算所述传感器单元的增益包括：
根据所述能量谱图中的、与单个暗事件的能量相对应的第一谱峰的横坐标计算所述增益，或根据所述第一谱峰的横坐标和所述能量谱图中的、与两个暗事件的能量相对应的第二谱峰的横坐标之间的差计算所述增益。
17.根据权利要求14所述的方法，其特征在于，所述计算在所述传感器单元中发生的暗事件的能量包括：
根据所述数字信号中的高电平和低电平的出现规律确定有效事件未发生的时间段；
从在有效事件未发生的时间段内生成的数字信号中选择在预定时段内生成的数字信号；以及
根据所述在预定时段内生成的数字信号计算所述暗事件的能量。
18.根据权利要求17所述的方法，其特征在于，所述预定时段设定为使得在所述预定时段内的平均暗事件个数满足第一预设条件，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，所述能量谱图包括多个谱峰并且所述能量谱图中的第N个谱峰与N个暗事件的能量相对应，
所述计算所述传感器单元的增益包括：
根据所述能量谱图中的、任意两个相邻谱峰的横坐标之间的差计算所述增益。
19.根据权利要求18所述的方法，其特征在于，所述第一预设条件是所述平均暗事件个数小于20。
20.根据权利要求17所述的方法，其特征在于，所述预定时段设定为使得在所述预定时段内的平均暗事件个数满足第二预设条件，所述能量谱图的横轴表示所述暗事件的能量，所述能量谱图的纵轴表示计算得到所述暗事件的能量的次数，
所述计算所述传感器单元的增益包括：
对所述能量谱图进行高斯拟合以获得所述能量谱图的方差和均值；以及
将所述方差除以所述均值以获得所述增益。
21.根据权利要求20所述的方法，其特征在于，所述第二预设条件是所述平均暗事件个数大于20。
22.根据权利要求13至21任一项所述的方法，其特征在于，所述方法进一步包括：
根据所述检测信号计算在所述传感器单元中发生的有效事件的能量E0；以及
根据以下公式对所述有效事件的能量E0进行校正：
E＝E0*G0/G，
其中，E为经校正的能量，G为所述增益，G0为预定常数。
CN201510953109.1A 2015-12-11 2015-12-17 用于测量传感器增益的装置、设备及方法 Active CN105572721B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510953109.1A CN105572721B (zh) 2015-12-11 2015-12-17 用于测量传感器增益的装置、设备及方法

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN2015109220088 2015-12-11
CN201510922008 2015-12-11
CN201510953109.1A CN105572721B (zh) 2015-12-11 2015-12-17 用于测量传感器增益的装置、设备及方法

Publications (2)

Publication Number Publication Date
CN105572721A true CN105572721A (zh) 2016-05-11
CN105572721B CN105572721B (zh) 2019-03-22

Family

ID=55883062

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510953109.1A Active CN105572721B (zh) 2015-12-11 2015-12-17 用于测量传感器增益的装置、设备及方法

Country Status (3)

Country Link
US (1) US10838088B2 (zh)
CN (1) CN105572721B (zh)
WO (1) WO2017097147A1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017016469A1 (zh) * 2015-07-29 2017-02-02 武汉中派科技有限责任公司 光子测量前端电路
CN106656390A (zh) * 2016-11-15 2017-05-10 武汉中派科技有限责任公司 用于测量光子时间信息的装置及方法
WO2017097147A1 (zh) * 2015-12-11 2017-06-15 武汉中派科技有限责任公司 用于测量传感器增益的装置、设备及方法
CN110967725A (zh) * 2019-12-18 2020-04-07 广州兰泰胜辐射防护科技有限公司 中子探测探头及中子探测芯片

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020208203A1 (en) 2019-04-12 2020-10-15 Arktis Radiation Detectors Ltd. Panel radiation detector

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5453610A (en) * 1994-05-20 1995-09-26 Summit World Trade Corporation Electronic gain control for photomultiplier used in gamma camera
CN101989526A (zh) * 2009-08-04 2011-03-23 思姆菲舍尔科技公司 使用参考辐射的pmt增益控制系统
CN102341727A (zh) * 2009-03-06 2012-02-01 皇家飞利浦电子股份有限公司 用于单光子计数器的温度补偿和控制电路
CN103424768A (zh) * 2012-05-25 2013-12-04 同方威视技术股份有限公司 一种用于探测器系统的增益稳定装置及其控制方法
CN103675886A (zh) * 2012-09-14 2014-03-26 赛默飞世尔科技有限公司 闪烁探测器增益控制
CN103777226A (zh) * 2012-10-19 2014-05-07 株式会社东芝 光传感器增益检测装置及光传感器增益检测方法
CN104570043A (zh) * 2014-12-18 2015-04-29 中国科学院高能物理研究所 硅光电倍增管的增益控制装置、系统及增益控制方法
CN105074500A (zh) * 2013-02-22 2015-11-18 阿雷瓦核废料回收公司 用于控制多像素光子计数设备的增益和零点的方法以及实施所述方法的光测量系统

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6754298B2 (en) * 2002-02-20 2004-06-22 The Regents Of The University Of Michigan Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method
US8694197B2 (en) * 2011-05-26 2014-04-08 GM Global Technology Operations LLC Gain/amplitude diagnostics of NOx sensors
US8907290B2 (en) * 2012-06-08 2014-12-09 General Electric Company Methods and systems for gain calibration of gamma ray detectors
KR101450806B1 (ko) * 2013-05-02 2014-10-15 서강대학교산학협력단 양전자방출 단층촬영기기의 센서 이득 제어 방법 및 시스템
US9541656B2 (en) * 2013-12-20 2017-01-10 General Electric Company System and method for compensating temperature gain variation in radiation detectors
US10168208B2 (en) * 2015-04-03 2019-01-01 Hitachi High-Technologies Corporation Light amount detection device, immune analyzing apparatus and charged particle beam apparatus that each use the light amount detection device
US9734603B2 (en) * 2015-06-30 2017-08-15 General Electric Company Systems and methods for peak tracking and gain adjustment
US10180351B2 (en) * 2015-12-08 2019-01-15 Texas Instruments Incorporated Dark current compensation for photon counting circuit
CN105572721B (zh) * 2015-12-11 2019-03-22 中派科技（深圳）有限责任公司 用于测量传感器增益的装置、设备及方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5453610A (en) * 1994-05-20 1995-09-26 Summit World Trade Corporation Electronic gain control for photomultiplier used in gamma camera
CN102341727A (zh) * 2009-03-06 2012-02-01 皇家飞利浦电子股份有限公司 用于单光子计数器的温度补偿和控制电路
CN101989526A (zh) * 2009-08-04 2011-03-23 思姆菲舍尔科技公司 使用参考辐射的pmt增益控制系统
CN103424768A (zh) * 2012-05-25 2013-12-04 同方威视技术股份有限公司 一种用于探测器系统的增益稳定装置及其控制方法
CN103675886A (zh) * 2012-09-14 2014-03-26 赛默飞世尔科技有限公司 闪烁探测器增益控制
GB2507863A (en) * 2012-09-14 2014-05-14 Thermo Fisher Scientific Inc Scintillation Detector Gain Control
CN103777226A (zh) * 2012-10-19 2014-05-07 株式会社东芝 光传感器增益检测装置及光传感器增益检测方法
CN105074500A (zh) * 2013-02-22 2015-11-18 阿雷瓦核废料回收公司 用于控制多像素光子计数设备的增益和零点的方法以及实施所述方法的光测量系统
CN104570043A (zh) * 2014-12-18 2015-04-29 中国科学院高能物理研究所 硅光电倍增管的增益控制装置、系统及增益控制方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017016469A1 (zh) * 2015-07-29 2017-02-02 武汉中派科技有限责任公司 光子测量前端电路
US10809395B2 (en) 2015-07-29 2020-10-20 Wuhan Zhong Pai Technology Co., Ltd. Photon measurement front-end circuit with integral module and a negative feedback module
WO2017097147A1 (zh) * 2015-12-11 2017-06-15 武汉中派科技有限责任公司 用于测量传感器增益的装置、设备及方法
US10838088B2 (en) 2015-12-11 2020-11-17 Zhongpai S&T (Shenzhen) Co., Ltd Apparatus, device and method for measuring gain of sensor
CN106656390A (zh) * 2016-11-15 2017-05-10 武汉中派科技有限责任公司 用于测量光子时间信息的装置及方法
CN110967725A (zh) * 2019-12-18 2020-04-07 广州兰泰胜辐射防护科技有限公司 中子探测探头及中子探测芯片
CN110967725B (zh) * 2019-12-18 2024-06-07 中国原子能科学研究院 中子探测探头及中子探测芯片

Also Published As

Publication number Publication date
US10838088B2 (en) 2020-11-17
CN105572721B (zh) 2019-03-22
WO2017097147A1 (zh) 2017-06-15
US20180292549A1 (en) 2018-10-11

Similar Documents

Publication Publication Date Title
CN105548848B (zh) 用于测量击穿电压的装置、设备及方法
Frach et al. The digital silicon photomultiplier—System architecture and performance evaluation
EP2989487B1 (en) Detection of radiation quanta using an optical detector pixel array and pixel cell trigger state sensing circuits
CN105572721A (zh) 用于测量传感器增益的装置、设备及方法
CN109581461B (zh) 核脉冲能量测量方法及系统
CN106656390B (zh) 用于测量光子时间信息的装置及方法
JP5800983B2 (ja) シンチレーションパルス情報取得のための方法及び装置
US9910167B2 (en) Method for digitalizing scintillation pulse
US8981307B2 (en) Pulse height analyzer and nuclear medicine diagnosis apparatus provided with the same
WO2017016469A1 (zh) 光子测量前端电路
CN107450092B (zh) 用于测量光子信息的装置
US7511264B2 (en) Method of resolving ambiguity in photon counting-based detectors
CN106706127B (zh) 基于SiPM的多光子探测方法
US11635531B2 (en) Apparatus for measuring photon information and photon measurement device
Wang et al. Comparison of three pre-amplifier circuits for time readout of SiPM in TOF-PET detectors
JP3015697B2 (ja) 放射線エネルギースペクトル測定装置
US4751390A (en) Radiation dose-rate meter using an energy-sensitive counter
Wang et al. A time-driven FPGA-based fast nuclear charge digitization method
Calò et al. A CMOS front-end for timing and charge readout of silicon photomultipliers
JP7057630B2 (ja) 放射線検出装置
CN111257919A (zh) 一种基于移动触发阈方法评价SiPM探测器的方法
Peng et al. Characterization of silicon photomultipliers and validation of the electrical model
Lerche et al. 5 PET Front-End Electronics
CN110398767A (zh) 射线探测电路
JPS6124657B2 (zh)

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20170519

Address after: 518063 Changhong technology building, No. 18 South twelve Road, Nanshan District science and technology, Guangdong, Shenzhen 903

Applicant after: Sino science and Technology (Shenzhen) Co., Ltd.

Address before: The Northeast Industrial Park Dongxin road East Lake Development Zone 430074 Hubei city of Wuhan Province, No. 11, building D C 1-4 layer

Applicant before: WUHAN ZHONGPAI TECHNOLOGY CO., LTD.

GR01 Patent grant
GR01 Patent grant