基于波形面积的雪崩光探测器增益测量方法
【技术领域】
[0001] 本发明属于微弱光探测方法技术领域，具体涉及一种基于波形面积的雪崩光探测 器增益测量方法。
【背景技术】
[0002] 所谓光探测器的增益是指光探测器的输出信号幅度和初始信号幅度之比；信号幅 度包括有电荷量、电压或电流等。
[0003] 现有光探测器增益的测量方法主要有以下三种：
[0004] 第一种方法是利用电压表或电流表直接测量输出信号幅度和初始信号幅度， 两者的比值即为增益；该种方法适用于增益较小的探测器（参考文献：R.H.Haitz， A. Goetzberger, R. M. Scarlett and ff. Shockley, Avalanche effects in Silicon p-n junctions. I. Localized Photomultiplication studies on microplasmas, Journal of applied physics, Vol.34, No. (6), 1963, pp. 1581-1590)；
[0005] 第二种方法是使用电荷数字转换器（即QDC，Charge to Digital Converter)或模 数转换器（即 ADC，Analog to Digital Converter)或多道分析器（即 MCA，Multi-Channel Analyzer)将光探测器输出的信号电荷量转化为数字量进行增益测量；通过分析相邻两光 电子峰的道数之差，乘以仪器的电荷转换率，即可算出探测器的增益；（参考文献：Paolo Finocchiaro,Alfio Pappalardo, Luigi Cosentino,etal. , Characterization of a Novel 100-Channel Silicon Photomultiplier-Part II ：Charge and Time,IEEE TRANSACTIONS ON ELECTRON DEVICES，VOL. 55, NO. 10, OCTOBER，（2008) pp. 2765-2773);
[0006] 第三种方法是利用放射性元素辐射的具有特征能量的射线来标定探测器的增益， (参考文南犬：G. F. Knoll，Radiation Detection and measurement，2nd edition, Wiley, New York, 2001, pp395. ；A. N. Otte, J. Hose, R. Mirzoyan, A. Romaszkiewicz, M. Teshima, A.Theaa, A measurement of the photon detection efficiency of silicon photomultipliers, Nuclear Instruments and Methods in Physics Research A, 567(2006)360 - 363)。
[0007] 第一类方法只适用于增益较小的探测器，对于高增益的探测器（>105)，探测器已 具有单光子响应能力，采用这种方法测量的结构不准确。
[0008] 第二类增益测量方法，其对应的系统较为复杂，所需的高速QDC、门产生器 等仪器造价昂贵，且属于专用测量装置，通用性不强，详见参考文献：(M.Gr〇dZiCka， M. Moszyr nski，T. Szcz，etal.，Energy resolution of small scintillation detectors with SiPM light readout，JINST，VOL. 8, P02017,（2013)1-17)。另外在利用该方法进行 测量之前，需要事先对QDC的电荷数字转化率进行定标，所以测量步骤较为繁琐。
[0009] 第三种方法需要在完备的核辐射保护装置下进行测量，否则对人体有伤害，给测 量造成了不便。
【发明内容】
[0010] 本发明的目的在于提供一种基于波形面积的雪崩光探测器增益测量方法，解决了 现有雪崩光探测器增益测量方法中存在的系统复杂、通用性不强及测量便捷性差的问题。
[0011] 本发明所采用的技术方案是，基于波形面积的雪崩光探测器增益测量方法，具体 按照以下步骤实施：
[0012] 步骤1、由稳压电源进行供电，将多像素光子探测器输出的雪崩信号输入到通用 数字示波器中，由通用数字示波器显示出波形电压值；
[0013] 若待测多像素光子探测器的信号幅度太小，就要先将多像素光子探测器输出的雪 崩信号通入高速放大器内进行放大，再输入通用数字示波器中；
[0014] 步骤2、经步骤1后，将通用数字示波器显示的波形电压值对时间进行积分，求取 出波形面积值；
[0015] 步骤3、预先在计算机内嵌入计算模块，利用计算机内的计算模型对经步骤2得到 的波形面积值A进行处理，获取增益值；
[0016] 步骤4、为了提高测量的准确性，要重复步骤1~步骤3至少100次，获得多个增益 值；
[0017] 步骤5、对经步骤4获取的多个增益值进行直方图统计，完成基于波形面积的雪崩 光探测器增益的测量。
[0018] 本发明的特点还在于：
[0019] 步骤2中波形面积值A通过如下算法获得：
[0021] 式中，VU1)是在时间点h测量的电脉冲波形信号的电压值，At是采样时间间隔， Ta、Tb分别是采样时间的起始时刻和终止时刻。
[0022] 步骤3具体按照以下步骤实施：
[0023] 步骤3. 1、利用计算机内嵌的计算模型将经步骤2得到的波形面积值A换算为电荷 量Q，具体按照以下算法实施：
[0025] 在式中，&是负载电阻值，A是波形面积值，也是电压波形的积分值；
[0026] 步骤3. 2、将经步骤3. 1换算得到的电荷量Q除以测量时间内初始载流子的电量 e，即得到多像素光子探测器的增益值Gain ;
[0028] 式中，Gain为探测器的增益值；
[0029] 若步骤1中采用了高速放大器，在进行增益换算时，只要将电荷量Q再除以高速放 大器的增益系数η即可，具体按照以下算法实施： CN 105181153 A 说明书 3/6 页
[0031] 多像素光子探测器采用MPPC Cl 1209-10型多像素光子探测器。
[0032] 高速放大器采用HSA-Y-2-40放大器。
[0033] 步骤5具体按照以下方式实施：
[0034] 经步骤4获得多个增益值后，利用计算机内置的直方图统计软件根据多个增益值 绘制出统计直方图，该统计直方图的峰值对应的增益值即为多像素光子探测器的平均增 益，从而完成基于波形面积的雪崩光探测器增益的测量。
[0035] 通用数字示波器为DP04102B-L型数字示波器。
[0036] 稳压电源为IT6235型精密稳压电源，用于提供电力供应。
[0037] 本发明的有益效果在于：
[0038] (1)本发明基于波形面积的雪崩光探测器增益测量方法，无需利用其它定标手段 定标即可实现待测探测器的增益测量，测量便捷，符合当今人们对于便捷性的需求。
[0039] (2)本发明基于波形面积的雪崩光探测器增益测量方法，只需使用实验室通用的 测试仪器及示波器进行测量，有效的降低了测量成本，解决了现有方法存在的测试系统价 格高昂，测量成本高的问题。
[0040] (3)本发明基于波形面积的雪崩光探测器增益测量方法通用性强，适用范围广泛， 具有非常大的应用前景。
【附图说明】
[0041] 图1是本发明基于波形面积的雪崩光探测器增益测量方法中涉及的测量系统的 工作状态示意图；
[0042] 图2是本发明基于波形面积的雪崩光探测器增益测量方法的流程图；
[0043] 图3是利用本发明基于波形面积的雪崩光探测器增益测量方法测得的多像素光 子探测器增益随偏压的变化趋势图。
[0044] 图中，1.多像素光子探测器，2.高速放大器，3.通用数字示波器，4.计算机，5.稳 压电源。
【具体实施方式】
[0045] 下面结合附图和【具体实施方式】对本发明进行详细说明。
[0046] 本发明基于波形面积的雪崩光探测器增益测量方法中应用的测量装置，其结构如 图1所示，包括有多像素光子探测器1，多像素光子探测器1外接稳压电源5,多像素光子探 测器1依次与高速放大器2、通用数字示波器3及计算机4连接；计算机4内嵌有计算模型， 该计算模型能求取波形面积，并将波形面积换算为增益值，还能对增益值进行统计；最终通 过对获取的所有增益值进行统计，实现本发明基于波形面积的雪崩光探测器增益的测量方 法。
[0047] 本发明基于波形面积的雪崩光探测器增益测量方法，如图2所示，具体按照以下 步骤实施：
[0048] 步骤1、由稳压电源5进行供电，将多像素光子探测器1输出的雪崩信号输入到通 用数字示波器3中，由通用数字示波器3显示出波形电压值；
[0049] 若待测多像素光子探测器1的信号幅度太小，就要先将多像素光子探测器1输出 的雪崩信号通入高速放大器2内进行放大，再输入通用数字示波器3中；
[0050] 其中的多像素光子探测器1采用MPPC Cl 1209-10型多像素光子探测器；
[0051] 当然，也可以用单点探测器或探测器阵列来替代多像素光子探测器1 ;
[0052] 高速放大器 2 为 HSA-Y-2-40 放大器，3 dB band width IOkHz-L 9GHz， noisefigure 4. 9dB, voltage Gain 40dB(100X)；
[0053] 通用数字示波器3为DP04102B-L型数字示波器；利用通用数字示波器3可以记录 和传输多像素光子探测器1输出的雪崩信号波形数据；
[0054] 稳压电源5为IT6235型精密稳压电源，用于提供电力供应。
[0055] 步骤2、经步骤1后，将通用数字示波器3显示的波形电压值对时间进行积分，求取 出波形面积值A ;
[0056] 波形面积值A按照如下算法获得：
[0058] 式中，VU1)是在时间点