A fast digital filter algorithm for gamma-ray spectroscopy with double-exponential decaying scintillators

Abstract

Scintillators like CsI(Na), having double-exponential decay times, typically cannot be used in high-count rate applications due to the complicated pulse shapes created by the convolution of their light output decay curves with the decay constant of charge integrating preamplifiers. We present here a novel digital filtering algorithm that is capable of using CsI(Na) at input count rates exceeding 250 kcps, while still achieving good energy resolution. We used a 2.54 cm diameter and 2.54 cm long CsI(Na) crystal, whose scintillation light can be best described by a short component with a 550 ns decay time and a long component with a 4 s decay time. The crystal was coupled to a 2.86 cm diameter photomultiplier tube. The digital filtering algorithm was implemented in XIA's all-digital Polaris spectrometer, in which five running sums were captured from each digitized scintillation pulse and the Polaris's on-board DSP read these sums and used a set of precomputed coefficients to reconstruct the pulse's total light output as a measure of the deposited energy. The algorithm was tested at different input count rates, ranging from 19 kcps to 270 kcps using a 1 mCi 137 Cs source. The energy resolution (full-width at half-maximum) at 662 keV was 10.7% at 19 kcps and 11.7% at 270 kcps with a filter rise time of 1.0 s, and improved to 7.0% and 8.4%, respectively, with a filter rise time of 3.2 s. The energy peak shifted by less than 0.3% for input count rates below the maximum throughput point. Output count rates of 65.3 and 17.8 kcps were obtained with filter rise time of 1.0 and 3.2 s, respectively, at an input count rate of 270 kcps. This algorithm can be easily adapted to other double-exponential decaying scintillators by changing the decay times used in the energy reconstruction formula.Index Terms-Author, please supply your own keywords or send a blank e-mail to keywords@ieee.org to receive a list of suggested keywords.

A general-purpose data acquisition system and a waveform analysis algorithm based on digitization

“…Pixie-16 模块采用数字梯形滤波器 [8,10] 将电荷灵敏型前置放大器产生的典型指数衰减信号转 换为梯形脉冲，用于快时间和脉冲幅度测量。 每个滤波器有两个主要参数，分别是梯形滤波的上升 时间和平台时间。对于快时间滤波 (fast filter)，滤波器采用较小的时间参数，此时梯形脉冲幅度与 输入信号上升沿的斜率成正比。数字恒比定时滤波器 (CFD filter)利用经典的恒比定时算法 [8,11] ，提 取快时间信号的与幅度无关的定时信息。对于脉冲幅度滤波(slow filter)，采用较大的时间参数，此 时梯形脉冲的幅度与输入信号的幅度成正比。在幅度测量中还需要第三个参数来校正前置放大器输 出信号的指数衰减 [10] 1 和 k = 10；(c) 拟合得到的α能谱。 [13] Figure 7. (Color online) Decomposition method of the pile-up pulse.…”

A segmented, enriched N-type germanium detector for neutrinoless double beta-decay experiments

“…Typically a pulser is used to decouple the additional contributing factors to the measured peak width, but since SEGA currently is not equipped with a pulser we placed a 60 Co source above the crystal and recorded data with very long traces (40 µs). Having 40 µs traces allowed us to apply an offline digital trapezoidal filter outlined in [31], which is the same filter implemented in the XIA module, while varying the shaping parameters, eliminating the need for multiple data sets for which the shaping parameters were varied from run to run. FWHM values were measured by fitting the 1332 keV photopeak in the trapezoidal filter output spectra using a Gaussian function with a linear background.…”

“…Energy measurements were either made using the trapezoidal filter in the XIA, DGF-4C module [31], or in some special cases, directly from trace data. The XIA trapezoidal filter is designed to be relatively insensitive to the signal rate and uses a large running sample of the baseline to accurately determine its contributions to the filter sums.…”

“…It was also necessary to avoid spurious triggers and incorrect integrated energy values reported by the XIA DGF-4C module when large image charge was present. To address these issues the trace from each channel was used to determine the energy using an offline trapezoidal filter [31]. This filter is similar to the trapezoidal filter the XIA DGF-4C module uses to calculate the energy, but results in lower resolution because the trace samples used for this particular study were 3 µs long and included only about 300 ns of baseline data.…”

An X-ray image intensifier for microsecond time-resolved experiments

“…The X-ray image intensifier has phosphor in the input window (CsI:Na + ) as well as in the output window. Since the major decay component of CsI:Na + is known to be 550 ns [12], it is expected that the persistence of CsI:Na + dominates the decay of the image intensifier when YAG is used in the output window. However, since CsI:Na + in the image intensifier is a needle-shaped crystal sandwiched between beryllium and photocathode, its persistence may be affected by this structure.…”