Optimal parameter choice of CR-RCm digital filter in nuclear pulse processing

NUCLEAR ELECTRONICS AND INSTRUMENTATION

Optimal parameter choice of CR-RC^m digital filter in nuclear pulse processing

Huai-Qiang Zhang，

Zhuo-Dai Li，

Bin Tang，

He-Xi Wu

Nuclear Science and Techniques

Vol.30, No.7

Article number 108

Published in print 01 Jul 2019

Available online 05 Jun 2019

DOI：10.1007/s41365-019-0638-7

63700

CR-RC^m filters are widely used in nuclear energy spectrum measurement systems. The choice of parameters of a CR-RC^m digital filter directly affects its performance in terms of energy resolution and pulse count rate in digital nuclear spectrometer systems. A numerical recursive model of a CR differential circuit and RC integration circuit is derived, which shows that the shaping result of CR-RC^m is determined by the adjustment parameter (k, it determines the shaping time of the shaper) and the integral number (m). Furthermore, the amplitude-frequency response of CR-RC^m is analyzed, which shows that it is a bandpass filter; the larger the shaping parameters (k and m), the narrower is the frequency band. CR-RC^m digital Gaussian shaping is performed on the actual sampled nuclear pulse signal under different shaping parameters. The energy spectrum of ¹³⁷Cs is measured based on the LaBr₃(Ce) detector under different parameters. The results show that the larger the shaping parameters (m and k), the closer the shaping result is to Gaussian shape, the wider is the shaped pulse, the higher is the energy resolution, and the lower is the pulse count rate. For the same batch of pulse signals, the energy resolution is increased from 3.8% to 3.5%, and the full energy peak area is reduced from 7815 to 6503. Thus, the optimal shaping parameters are m=3 and k=0.95. These research results can provide a design reference for the development of digital nuclear spectrometer measurement systems.

CR-RCm digital filterNuclear pulse signalGaussian shapingEnergy resolution

1 Introduction

In analog nuclear spectrometer systems, to improve the signal-to-noise ratio (SNR) of nuclear pulse processing, a CR-RC^m circuit is often used to filter the nuclear pulse signal and shape it into a quasi-Gaussian shape [1-3]. In digital nuclear spectrometer systems, based on a high-speed analog-to-digital converter (ADC) and a high-performance field-programmable gate array (FPGA), the filtering of a nuclear pulse signal is often realized using a digital algorithm. Triangular filters, trapezoidal filters, and cusp-like filters are widely used in digital nuclear pulse processing [4-11]. Owing to the simple implementation and high SNR of a Gaussian filter, it is an important digital filter shaping method. The method mainly includes a time domain based on the digital Sallen–Key [12-16], and a transfer function based on CR-RC^m [17]. The digital Gaussian filter based on CR-RC^m is often used in digital nuclear spectrometer systems, owing to its superior performance and simple shaping principle and method. Based on the basic principles of CR and RC circuits, a recursive model of the digital filter shaping algorithm is derived. Furthermore, the actual sampled nuclear signal is processed using the CR-RC^m filter, and the effect of the digital filter is analyzed based on the amplitude-frequency response characteristics. The performance is compared and analyzed in terms of the energy resolution and pulse count rate under different forming parameters.

2 Filter shaping principle of CR-RC^m

The CR-RC^m filter shaping circuit used in nuclear energy spectrum measurement systems is shown in Fig. 1. In this figure, V_in is the input of the circuit, CR is the differential circuit (high-pass filter circuit), RC is the integral circuit (low-pass filter circuit), and V_out is the output of the circuit. The input nuclear pulse signal is shaped into a quasi-Gaussian shape by the circuit.

Fig. 1.

Circuit of a CR-RC^m filter

As shown in Fig. 1, R₁C₁=R₂C₂=R_mC_m=τ. V_in (as a function of the time t) is the input signal. When it is a step signal, it is represented by Eq. (1) in the time and complex frequency domains, where Q is a charge and Ci is an input capacitor.

V_{in} (t) = \frac{Q}{C_{i}} u (t) \Rightarrow V_{in} (s) = \frac{Q}{C_{i}} \cdot \frac{1}{S} .

(1)

The transfer function of CR differential shaping is

H_{CR} (s) = \frac{S τ}{1 + S τ} .

(2)

The transfer function of RC integral shaping is

H_{RC} (s) = \frac{1}{1 + S τ} .

(3)

From Eq. (2) and Eq. (3), the output is

\begin{array}{l} V_{out} (s) = \frac{Q}{C_{i}} \frac{τ}{{(1 + S τ)}^{m + 1}} \\ \Rightarrow & V_{out} (t) = \frac{Q}{C_{i} m!} {(\frac{t}{τ})}^{m} e^{- \frac{t}{τ}} \end{array}

(4)

If the input signal V_in is a negative exponential signal, it is necessary to add a pole-zero cancellation circuit before the differential shaping to achieve Gaussian shaping of the nuclear pulse signal [18].

3 Derivation of algorithm of CR-RC^m digital filter

A CR-RC^m filter consists of a CR high-pass filter and m RC low-pass filters. To digitize the nuclear signal more effectively, the CR and RC circuits are separately digitized [19, 20].

(1) The differential equation for CR high-pass circuits according to Kirchhoff’s current law (KCL) is shown in Eq. (5).

\frac{d V_{i} (t)}{d t} = \frac{V_{o} (t)}{C R} + \frac{d V_{o} (t)}{d t}

(5)

where V_i(t) and Vo(_t) are the input and output signals of the CR circuit, respectively. For the discrete signal sequence digitized by the high-speed ADC, V_i(t) is converted into a digital sequence X[n], and V_o(t) is converted into a digital sequence Y[n]. The derivation can yield the digital recursive formula shown in Eq. (6).

Y [n] = k (X [n] - X [n - 1]) + k \cdot Y [n - 1]

(6)

As shown in Eq. (6), k = CR/(CR+T), where T is the sampling period of the ADC. Thus, the transfer function in the z-domain of the CR high-pass filter and its frequency response function are [20, 21] expressed as follows:

\begin{array}{l} H_{CR} (z) = k \frac{1 - z^{- 1}}{1 - k \cdot z^{- 1}} \\ \Rightarrow & H_{CR} (e^{j w}) = k \frac{1 - e^{- j w}}{1 - k e^{- j w}} \end{array}

(7)

In Eq. (7), w is the angular frequency of the component signal, and Euler’s formula transformation for Eq. (7) is shown in Eq. (8).

H_{CR} (e^{j w}) = k \frac{1 - \cos (w) + j \cdot \sin (w)}{1 - k (\cos (w) + j \cdot \sin (w))}

(8)

The modulo of Eq. (8) is shown in Eq. (9).

| H_{CR} (e^{j w}) | = k \frac{\sqrt{2 - 2 \cos (w)}}{\sqrt{1 + k^{2} - 2 k \cdot \cos (w)}}

(9)

(2)The digital recursive equation, frequency response function, and modes of the RC low-pass circuit according to the same principle are shown in Eq. (10), Eq. (11), and Eq. (12), respectively.

Y [n] = (1 - k) \cdot X [n] + k \cdot Y [n - 1],

(10)

H_{RC} (e^{j w}) = \frac{1 - k}{1 - k \cdot \cos (w) + j \cdot k \cdot \sin (w)},

(11)

| H_{RC} (e^{j w}) | = \frac{1 - k}{\sqrt{1 + k^{2} - 2 k \cdot \cos (w)}} .

(12)

As shown in Eq. (10), k =RC/(RC+T), where T is the sampling period of the ADC.

According to Eq. (8) and Eq. (11), the frequency response function of the CR-RC^m filter is shown in Eq. (13).

H (e^{j w}) = \frac{{(1 - k)}^{m} \cdot k \cdot (1 - k \cos (w) + j \cdot k \cdot \sin (w))}{1 - k \cdot \cos (w) + j \cdot k \cdot \sin (w)}

(13)

The modulo of Eq. (13) is shown in Eq. (14).

| H (e^{j w}) | = \frac{{(1 - k)}^{m} \cdot k \cdot \sqrt{2 - 2 \cos (w)}}{{(\sqrt{1 + k^{2} - 2 k \cdot \cos (w)})}^{m + 1}}

(14)

Based on the SNR coefficient of the CR-RC^m filter $F = \frac{m! e^{m}}{m^{m}} {(2 m - 1)}^{\frac{1}{4}} {[\frac{(2 m - 3)!!}{2 (2 m)!!}]}^{\frac{1}{2}}$ , for m=1, F=1.36; for m=2, F=1.22; for m=3, F=1.18; for m=4, F=1.17; and for m=5, F=1.16. Thus, with the increase in the parameter m, the change in the parameter F is small, and the filtering system becomes more complicated. Therefore, in practical applications, the parameter m is generally 3 or 4. For the value of the parameter k, the sampling period of the high-speed ADC is generally 100 ns, 50 ns, and 20 ns in the digital nuclear energy spectrum. When RC=τ =1 μs, the corresponding k values are k=0.91, k=0.95, and k=0.98. According to Eq. (9) and Eq. (12), the frequency response curve of the CR and RC filters is shown in Fig. 2.

Fig. 2.

Frequency response curve of CR and RC filters (different k values)

As shown in Fig. 2, the amplitude-frequency response curve of the CR filter shows an incremental increase in frequency. As the value of k increases, the response curve rises faster, and it approaches 1 sooner. The amplitude-frequency response curve of the RC filter shows a decline. As the value of k increases, the response curve falls faster, and it approaches zero sooner. When m=3, k=0.91, k=0.95, and k=0.98, the frequency response curve of CR-RC^m is shown in Fig. 3. When k=0.95, m=2, m=3, and m=4,the frequency response curve of CR-RC^m is shown in Fig. 4.

Fig. 3.

Frequency response curve of the CR-RC³ filter (different k values)

Fig. 4.

Frequency response curve of the CR-RC^m filter (different m values)

As shown in Fig. 3 and Fig. 4, the frequency response curve of the digitized CR-RC^m filter shows that it is a bandpass filter. The parameters m and k can be used to adjust the bandpass region of the filter, and the larger the values of m and k, the narrower is the bandpass region.

4 Algorithm implementation

¹³⁷Cs is measured using the CH249-02 NaI detector (Hamamatsu). The high voltage is +560 V, and the output pulse of the signal conditioning circuit is approximately 10μs. USB-5133 (National Instruments Corporation) is used to obtain the original nuclear pulse signal (the sampling rate is 50 MHz). According to the numerical recursive formula of the CR-RC^m filter, for m=3, the CR-RC^m shaping results under different k values are as shown in Fig. 5.

Fig. 5.

Actual nuclear pulse signal of CR-RC^m shaping

As shown in Fig. 5, with the increase of parameter k, the rising edge of the CR-RC^m is gentler, and the closer to a Gaussian shape, the lower the pulse amplitude; however, for longer pulse shaping time, the generalization of the pile-up occurs, and the original signal is a negative exponentially decreasing signal. After CR differential shaping, a certain degree of zero-crossing signal is present. For a pulse with a small k value or for a long interval of the original pulse, there is no pulse pile-up. The undershoot signal of zero has no effect on the subsequent pulse amplitude extraction. The k value is too large, and the longer the shaping time, the higher is the pile-up probability. Furthermore, the undershoot of the zero crossing directly affects the amplitude extraction of subsequent adjacent pulses [22].

5 Test and analysis of CR-RC^m filter performance

The LaBr₃(Ce) detector (Beijing Huakailong’s 2" ×2") is used to measure the ¹³⁷Cs energy spectrum. The ADC sampling rate is 20 MHz. For the same batch of sampling pulse data, when m=3, k=0.91, k=0.95,and k=0.98, the energy spectrum is shown in Fig. 6.

Fig. 6.

Energy spectrum of ¹³⁷Cs (partial channels)

As shown in Fig. 6, to demonstrate the energy spectrum results under different k values more clearly, the data of partial channels are displayed. Similarly, when k=0.95, m=2, m=3, and m=4,the statistical analysis of the energy resolution and the full energy peak area of the energy spectrum is shown in Table 1.

Performance comparison of different shaping parameters

Shaping parameter	Energy resolution	Full energy Peak area
CR-RC³ k=0.91	3.8%	7815
CR-RC³ k=0.95	3.6%	7362
CR-RC³ k=0.98	3.5%	6537
CR-RC² k=0.95	3.7%	7471
CR-RC⁴ k=0.95	3.6%	6503

As shown in Table 1, for the same batch of pulses, after the CR-RC^m digital filter shaping process under different parameters, with the increase in k and m values, the energy resolution is increased. This is because the width of the shaping pulse is too large, and the full energy peak area of the energy spectrum is lowered owing to the pile-up identification [23].

From Table 1, when m=3 and k=0.91, the energy resolution and the full energy peak area are, respectively, 3.8% and 7815; when m=3 and k=0.95, the energy resolution and the full energy peak area are, respectively, 3.6% and 7362; when m=3 and k=0.98, the energy resolution and the full energy peak area are, respectively, 3.5% and 6537; when k=0.95 and m=2, the energy resolution and pulse counting rates are, respectively, 3.7% and 7471; when k=0.95 and m=4, the energy resolution and the full energy peak area are, respectively, 3.6% and 6503.

6 Conclusion

CR-RC^m digital filters are widely used in digital nuclear energy spectrum measurement systems. The filter parameters directly affect the energy resolution and pulse count rate of the system. The numerical recursive algorithm of digital filtering of CR and RC is based on KCL. The measured signal is processed by the CR-RC^m digital filter, and the energy resolution and maximum pulse of the system are obtained under different k and m values in the energy spectrum of ¹³⁷Cs. From the comparative analysis, combined with the above indicators, it is determined that if the system requires high energy resolution, the values of k and m can be larger. Moreover, when it requires a high count rate, the values of k and m can be smaller. When the system needs to balance the performance of both, the optimal parameters of the CR-RC^m digital filter are k = 0.95 and m=3. The relevant results can provide a design reference for further use of the CR-RC^m digital filter.

References

[1]

H.Q. Zhang, H.X. Wu, H.T. Ning.

Research on the Simulation and Gaussian Shaping of Nuclear Signal Based on Multisim Software

. Nuclear Electronics & Detection Technology. 35(8), 845-848 (2015). doi: 10.3969/j.issn.0258-0934.2015.08.023 (in Chinese)