Digitalization of inverting filter shaping circuit for nuclear pulse signals

NUCLEAR ELECTRONICS AND INSTRUMENTATION

Digitalization of inverting filter shaping circuit for nuclear pulse signals

Huai-Qiang Zhang ，

Hong-Tao Shi，

Zhuo-Dai Li，

Yu-Wen Li

Nuclear Science and Techniques

Vol.31, No.9

Article number 86

Published in print 01 Sep 2020

Available online 28 Aug 2020

DOI：10.1007/s41365-020-00799-2

81801

In the design of filter shaping circuits for nuclear pulse signals, inverting filter shaping circuits perform better than non-inverting filter shaping circuits. Because these circuits facilitate changing the phase of a pulse signal, they are widely used in processing nuclear pulse signals. In this study, the transfer functions of four types of inverting filter shaping circuits, namely, the common inverting filter shaping, improved inverting filter shaping, multiple feedback low-pass filter shaping, and third-order multiple feedback low-pass filter shaping, in the Laplacian domain, are derived. We establish the numerical recursive function models and digitalize the four circuits, obtain the transfer functions in the Z domain, and analyze the filter performance and amplitude-frequency response characteristics in the frequency domain. Based on the actual nuclear pulse signal of the Si-PIN detector, we realize four types of inverting digital shaping. The results show that under the same shaping parameters, the common inverting digital shaping has better amplitude extraction characteristics, the third-order multiple feedback low-pass digital shaping has better noise suppression performance, and the multiple feedback digital shaping takes into account both pulse amplitude extraction and noise suppression performance.

Nuclear pulse signalInverting filter shaping circuitDigital shapingAmplitude-frequency response

1 Introduction

When nuclear pulse signals are processed, the output nuclear pulse signal of the pre-amplified needs to be filtered and shaped to improve the signal-to-noise ratio (SNR) of the system and change the pulse shape for subsequent processing. Several studies have focused on analog non-inverting filter shaping circuits for nuclear pulse signals. Evariste et al. designed and developed a set of amplifier filter shaping circuits for Si (Li), CdZnTe and CsI detectors, which included a fourth-order integration stage and a differentiation stage [1]. Gan et al. designed a low-noise nuclear pulse signal processing system, in which CR-RC was used to realize filter shaping [2]. Gao et al. designed CR-RC and Sallen-Key circuits based on CdZnTe and Si-PIN detectors, which were used to shaped into quasi-Gaussian [3]. Zhang et al. developed a low-noise integrated circuit based on CR-RC², which exhibited excellent filtering effects [4]. The low-pass feedback terminal of the charge-sensitive preamplifier, based on RC low-pass feedback, was replaced by a junction field-effect transistor and an RC low-pass network for the high-resistance feedback resistor. The noise of the feedback terminal was considerably smaller than the thermal noise introduced by the high resistance feedback resistor [5]. In the work of Zeng et al. [6], a CR low-pass filter circuit was used to process the output signal of the high-purity germanium detector to filter out low-frequency interference. This circuit was equipped with a non-inverting follower and an adjustable gain amplifier circuit to amplify and filter the nuclear pulse signal. Wang et al. developed a compact 16-channel integrated charge-sensitive preamplifier to support the large-scale detector array used in modern nuclear physics experiments. A first-order RC filter circuit was employed to reduce the noise generated by the high voltage cable. Meanwhile, the resistor–capacitor feedback network formed a charge integrating and bleeding circuit that determined the energy sensitivity and pulse decay time [7]. The above circuits were designed and analyzed for their performance as analog non-inverting amplification filter shaping circuits for nuclear pulse signals.

We summarize the major works reporting the digital non-inverting filter shaping method of nuclear pulse signals as follows. Zhang et al. used the digital trapezoidal shaping method to process the simulated nuclear signal and the actual sampled nuclear signal to obtain the best digital shaping parameters [8]. Zhang et al. [9] derived the transfer function of the filter shaping circuit based on Sallen-Key in the Laplacian and Z domains and improved its numerical recursive function model at various values of resistance and capacitance. A comparative analysis of performance indicators, such as amplitude extraction and noise suppression, led to the optimal design model and filter shaping parameters. The CR-RC^m digital shaping method of nuclear pulse signal processing was implemented based on FPGA [10]. Zhou et al. derived the Sallen-Key numerical recursive function model using Kirchhoff’s current law. This was also the earliest digital shaping processing of nuclear pulse signals in the temporal domain [11]. Zhou et al. [12] also used the results of earlier studies to improve the digital Sallen-Key recursive function model and introduce the amplitude adjustable parameter. The transfer function of the digital Sallen-Key circuit was derived in the Z domain and the nuclear pulse signal was convolved with the impulse response function. However, the result, when applied to the processing of the double-exponential nuclear pulse signal, did not improve the real-time performance. This limited the further use of this method [13]. Zhang et al. constructed a set of digital processing platforms for nuclear pulse signals based on MATLAB, which realized the digital shaping processing of simulated nuclear signals and actual sampled nuclear signals under different shaping methods and shaping parameters [14]. In a previous study, the quality factor and cut-off frequency of the digital Sallen-Key shaping method were introduced and simultaneously compared with the trapezoidal shaping results [15]. It was observed that digital Gaussian shaping based on the Sallen-Key method exhibited better noise suppression performance, which was beneficial for pulse amplitude extraction. In contrast, digital Gaussian shaping based on CR-RC^m had better counting rate characteristics, which aided pulse pile-up identification [16]. Zhang et al. studied the digitization of the pole-zero cancellation and Sallen-Key filter shaping circuits and obtained the numerical recursive function model and the amplitude-frequency response curve, respectively [17]. Liu et al. applied the CR-RC^m digital shaping method to the γ-spectrum measurement system and obtained the relationship between the shaping parameters and energy resolution [18]. Regadio et al. designed and developed a set of adjustable digital shapers for particle physics, which could adjust the shaping algorithm and shaping parameters according to actual requirements, and implemented the algorithm in FPGA [19]. Valentin used the unfolding-synthesis technique in digital pulse processing systems by digitizing analog signals according to the unfolding system and converting the nuclear signal into the desired pulse shapes [20]. The adaptive digital pulse shaping technique was investigated for real-time signal processing, and a digital Sallen–Key low-pass filter was implemented to enhance the SNR and reduce baseline drifting in this system [21].

In comparison with the non-inverting digital filter shaping method, the inverting digital shaping method exhibits a better SNR, which aids in pulse amplitude extraction under the same shaping parameters. The digital model functions of the inverting filter circuits, which were obtained using the numerical solution of differential equations, confirmed that the inverting digital shaping method exhibited better noise suppression performance [22]. However, the amplitude-frequency response was not analyzed, and hence, the numerical recursive function model of the serialized inverting filter shaping circuits was not analyzed and studied.

In this study, the digitalization of common inverting filter shaping circuits, improved inverting filter shaping circuits, multiple feedback low-pass filter shaping circuits and third-order multiple feedback low-pass filter shaping circuits is realized based on the preceding studies. The transfer function in the Laplacian domain, the amplitude-frequency response curve in the Z domain, and the numerical recursive function model are established. These are highly significant for the digital shaping processing of nuclear pulses and the optimal choice of parameters for nuclear pulse signal inverting filter shaping circuits.

2 Digitalization of inverting filter shaping circuits

Inverting filter shaping circuits for nuclear pulse signal processing include common inverting filter shaping circuits, improved inverting filter shaping circuits, multiple feedback low-pass filter shaping circuits, and third-order multiple feedback low-pass filter shaping circuits.

2.1 Digitalization of the common inverting filter shaping circuit

The common inverting filter shaping circuit is composed of an operational amplifier (OPAMP) A, resistors R₁ and R₂ and capacitor C. The circuit is shown in Fig. 1.

Fig. 1.

Schematic diagram of the common inverting filter shaping circuit.

As shown in Fig. 1, V_i and V_o are the input and output voltages, respectively, of the common inverting filter shaping circuit. The transfer function of the common inverting filter shaping circuit in the Laplacian domain is shown in Eq. (1).

H (s) = \frac{V_{o} (s)}{V_{i} (s)} = \frac{- \frac{1}{R_{1} C}}{s + \frac{1}{R_{2} C}}

(1)

When V_i is converted into xi, V_o is converted into yi, where R₁ = R₂ = R. According to differential numerical analysis, Eq. (1) leads to Eq. (2):

R C \cdot \frac{y_{i} - y_{i - 1}}{Δ t} + y_{i} = - x_{i} .

(2)

When xi is converted into a sequence of numbers, x[n], yi is converted into a sequence of numbers, y[n]. Thus, Eq. (3) can be obtained:

y [n] = \frac{k \cdot y [n - 1] - x [n]}{k + 1} .

(3)

When k = RC / Δt (Δt is sampling interval of ADC), the transfer function of the common inverting filter shaping circuit in the Z domain is as shown in Eq. (4):

H (z) = \frac{- 1}{k + 1 - k \cdot z^{- 1}} .

(4)

The frequency response function is:

H (e^{j w}) = \frac{- 1}{k (1 - e^{- j w}) + 1} .

(5)

The modulus of Eq. (5) is shown in Eq. (6).

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

(6)

Eq. (6) may be used to obtain the amplitude- frequency response curves.

2.2 Digitalization of the improved inverting filter shaping circuit

The improved inverting filter shaping circuit is composed of an OPAMP A, resistors R₁, R₂ and R₃ and capacitors C₁ and C₂. The corresponding circuit is shown in Fig. 2.

Fig. 2.

Schematic diagram of the improved inverting filter shaping circuit.

The transfer function in the Laplacian domain of the improved inverting filter shaping circuit shown in Fig. 2 is given by Eq. (7).

\begin{array}{l} H (s) & ​ = ​ \frac{V_{o} (s)}{V_{i} (s)} \\ ​ = ​ \frac{- \frac{1}{R_{1} C_{1}} s - \frac{1}{R_{1} R_{2} C_{1} C_{2}}}{s^{2} + (\frac{1}{R_{2} C_{2}} + \frac{1}{R_{2} C_{1}} + \frac{1}{R_{3} C_{2}}) s + \frac{1}{R_{2} R_{3} C_{1} C_{2}}} \end{array}

(7)

When R₁ = R₂ = R₃ = R, C₁ = C₂ = C, according to the differential numerical analysis, Eq. (7) changes into Eq. (8).

\begin{array}{l} {(R C)}^{2} \cdot \frac{\frac{y_{i} ​ - ​ y_{i ​ - ​ 1}}{Δ t} ​ - ​ \frac{y_{i ​ - ​ 1} ​ - ​ y_{i ​ - ​ 2}}{Δ t}}{Δ t} ​ + ​ 3 R C \cdot \frac{y_{i} ​ - ​ y_{i ​ - ​ 1}}{Δ t} ​ + ​ y_{i} \\ ​ = ​ & ​ - ​ R C \cdot \frac{x_{i} ​ - ​ x_{i ​ - ​ 1}}{Δ t} ​ - ​ x_{i} \end{array}

(8)

Then Eq. (9) can be obtained as

\begin{array}{l} y [n] = & \frac{(2 k^{2} + 3 k) \cdot y [n - 1] - k^{2} \cdot y [n - 2]}{k^{2} + 3 k + 1} \\ - \frac{- (k + 1) \cdot x [n] - k \cdot x [n - 1]}{k^{2} + 3 k + 1} \end{array}

(9)

When k = RC / Δt, the transfer function of the improved inverting filter shaping circuit in the Z domain is shown in Eq. (10).

H (z) = \frac{- k (1 - z^{- 1}) - 1}{k^{2} {(1 - z^{- 1})}^{2} + 3 k (1 - z^{- 1}) + 1}

(10)

The frequency response function is shown in Eq. (11):

H (e^{j w}) = \frac{- k (1 - e^{- j w}) - 1}{k^{2} {(1 - e^{- j w})}^{2} + 3 k (1 - e^{- j w}) + 1}

(11)

The modulus of Eq. (11) is shown in Eq. (12).

\begin{array}{l} | H (e^{j w}) | & = \sqrt{\frac{1 + 2 k + 2 k^{2} - 2 k (1 + k) \cdot c o s (w)}{A - B \cdot c o s (w) + C \cdot c o s (2 w)}} \\ A = 1 + 6 k + 20 k^{2} + 18 k^{3} + 6 k^{4} \\ B = 2 k (3 + 11 k + 12 k^{2} + 4 k^{3}) \\ C = 2 k^{2} (1 + 3 k + k^{2}) \end{array}

(12)

The amplitude-frequency response curves can be obtained using Eq. (12).

2.3 Digitalization of the multiple feedback low-pass filter shaping circuit

The multiple feedback low-pass filter shaping circuit consists of an OPAMP A, resistors R₁, R₂ and R₃ and capacitors C₁ and C₂. The corresponding circuit is shown in Fig. 3.

Fig. 3.

Schematic diagram of the multiple feedback low-pass filter shaping circuit.

In the Laplacian domain, the transfer function of the multiple feedback low-pass filter shaping circuit shown in Fig. 3 is given by Eq. (13).

\begin{array}{l} H (s) & = \frac{V_{o} (s)}{V_{i} (s)} \\ = \frac{- \frac{1}{R_{1} C_{1} R_{2} C_{2}}}{s^{2} + s \frac{1}{C_{1}} (\frac{1}{R_{1}} + \frac{1}{R_{2}} + \frac{1}{R_{3}}) + \frac{1}{R_{2} C_{1} R_{3} C_{2}}} \end{array}

(13)

When R₁ = R₂ = R₃ = R, C₁ = C₂ = C, according to differential numerical analysis, Eq. (13) can be changed into Eq. (14).

\begin{array}{l} {(R C)}^{2} \cdot \frac{\frac{y_{i} - y_{i - 1}}{Δ t} - \frac{y_{i - 1} - y_{i - 2}}{Δ t}}{Δ t} \\ + 3 R C \cdot \frac{y_{i} - y_{i - 1}}{Δ t} + y_{i} = - x_{i} \end{array}

(14)

Further derivation yields Eq. (15) as

y [n] = \frac{(2 k^{2} + 3 k) \cdot y [n - 1] - k^{2} \cdot y [n - 2] - x [n]}{k^{2} + 3 k + 1} .

(15)

When k = RC / Δt, the transfer function of the multiple feedback low-pass filter shaping circuit in the Z domain is shown in Eq. (16)

H (e^{j w}) = \frac{- 1}{k^{2} {(1 - z^{- 1})}^{2} + 3 k (1 - z^{- 1}) + 1} .

(16)

The frequency response function is

H (e^{j w}) = \frac{- 1}{k^{2} {(1 - e^{- j w})}^{2} + 3 k (1 - e^{- j w}) + 1} .

(17)

The modulus of Eq. (17) is given by Eq. (18).

\begin{array}{l} | H (e^{j w}) | & = \sqrt{\frac{1}{A - B \cdot \cos (w) + C \cdot \cos (2 w)}} \\ A = 1 + 6 k + 20 k^{2} + 18 k^{3} + 6 k^{4} \\ B = 2 k (3 + 11 k + 12 k^{2} + 4 k^{3}) \\ C = 2 k^{2} (1 + 3 k + k^{2}) \end{array}

(18)

The amplitude- frequency response curves can be obtained using Eq. (18).

2.4 Digitalization of the third-order multiple feedback low-pass filter shaping circuit

The third-order multiple feedback low-pass filter shaping circuit consists of an OPAMP A, resistors R₁, R₂, R₃ and R₄ and capacitors C₁, C₂ and C₃. The corresponding circuit is shown in Fig. 4.

Fig. 4.

Schematic diagram of the third-order multiple feedback low-pass filter shaping circuit.

The transfer function in the Laplacian domain of the third-order multiple feedback low-pass filter shaping circuit shown in Fig. 4 is given by Eq.(19).

\begin{array}{l} H (s) = \frac{V_{o} (s)}{V_{i} (s)} \\ = \frac{- \frac{1}{R_{1} R_{2} R_{3} C_{1} C_{2} C_{3}}}{s^{3} + s^{2} (\frac{1}{R_{1} C_{1}} + \frac{1}{R_{2} C_{1}} + \frac{1}{R_{2} C_{2}} + \frac{1}{R_{3} C_{2}} + \frac{1}{R_{4} C_{2}}) + \cdot \cdot \cdot} \\ s (\frac{1}{R_{1} R_{2} C_{1} C_{2}} + \frac{1}{R_{1} R_{3} C_{1} C_{2}} + \frac{1}{R_{2} R_{3} C_{1} C_{2}}) + \cdot \cdot \cdot \\ s (\frac{1}{R_{2} R_{4} C_{1} C_{2}} + \frac{1}{R_{1} R_{4} C_{1} C_{2}} + \frac{1}{R_{3} R_{4} C_{2} C_{3}}) + \cdot \cdot \cdot \\ \frac{R_{1} + R_{2}}{R_{1} R_{2} R_{3} R_{4} C_{1} C_{2} C_{3}} \end{array}

(19)

When R₁ = R₂ = R₃ = R₄ = R,C₁ = C₂ = C₃ = C, Eq. (19) can be changed into Eq. (20) using differential numerical analysis.

\begin{array}{l} {(R C)}^{3} ​ \cdot ​ \frac{\frac{\frac{y_{i} - y_{i - 1}}{Δ t} ​ - ​ \frac{y_{i - 1} - y_{i - 2}}{Δ t}}{Δ t} ​ - ​ \frac{\frac{y_{i - 2} - y_{i - 3}}{Δ t} ​ - ​ \frac{y_{i - 3} - y_{i - 4}}{Δ t}}{Δ t}}{Δ t} + \dots \\ 5 (R C)^{2} ​ \cdot ​ \frac{\frac{y_{i} - y_{i - 1}}{Δ t} - \frac{y_{i - 1} - y_{i - 2}}{Δ t}}{Δ t} + \dots \\ 6 R C ​ \cdot ​ \frac{y_{i} - y_{i - 1}}{Δ t} + y_{i} = - x_{i} \end{array}

(20)

Upon further derivation, Eq. (21) can be obtained.

\begin{array}{l} y [n] ​ = ​ & \frac{(2 k^{3} ​ + ​ 10 k^{2} ​ + ​ 6 k) ​ \cdot ​ y [n - 1] ​ - ​ 5 k^{2} ​ \cdot ​ y [n - 2]}{k^{3} ​ + ​ 5 k^{2} ​ + ​ 6 k ​ + ​ 2} ​ - ​ \dots \\ \frac{- 2 k^{3} \cdot y [n - 3] - k^{3} \cdot y [n - 4] - x [n]}{k^{3} + 5 k^{2} + 6 k + 2} \end{array}

(21)

When k = RC / Δt, the transfer function of the third-order multiple feedback low-pass filter shaping circuit in the Z domain is shown in Eq. (22).

\begin{array}{l} H (z) ​ = ​ & \frac{​ - ​ 1}{k^{3} {(1 ​ - ​ z^{​ - ​ 1})}^{3} (1 ​ + ​ z^{​ - ​ 1}) ​ + ​ 5 k^{2} {(1 ​ - ​ z^{​ - ​ 1})}^{2} ​ - ​ 6 k (1 ​ - ​ z^{​ - ​ 1}) ​ + ​ 2} \end{array}

(22)

The frequency response function is shown in Eq. (23),

\begin{array}{l} H ​ (z) ​ ​ = ​ \frac{​ - ​ 1}{k^{3} {(1 ​ ​ - ​ ​ e^{​ - ​ j w} ​)}^{3} (1 ​ ​ + ​ ​ e^{​ - ​ j w} ​) ​ ​ + ​ ​ 5 k^{2} ​ {(​ 1 ​ ​ - ​ ​ e^{​ - ​ j w} ​)}^{2} ​ ​ - ​ ​ 6 k (1 ​ ​ - ​ ​ e^{​ - ​ j w} ​) ​ ​ + ​ 2} \end{array}

(23)

and its modulus is given by Eq. (24).

\begin{array}{l} ​ | ​ H ​ (​ e^{j w} ​) ​ | ​ & ​ = ​ \frac{1}{\sqrt{2 ​ \cdot ​ (A ​ - ​ ​ B ​ \cdot ​ \cos (w) ​ ​ - ​ C ​ ​ \cdot ​ \cos (2 w) ​ ​ + ​ ​ D ​ ​ \cdot ​ \cos (3 w) ​ ​ - ​ ​ E ​ ​ \cdot ​ \cos (4 w) ​)}} \\ A ​ = ​ 5 k^{6} ​ + ​ 25 k^{5} ​ + ​ 57 k^{4} ​ - ​ 88 k^{3} ​ + ​ 46 k^{2} ​ - ​ 12 k ​ + ​ 2 \\ B ​ = ​ 2 k (2 k^{5} ​ + ​ 10 k^{4} ​ + ​ 41 k^{3} ​ - ​ 58 k^{2} ​ + ​ 28 k ​ - ​ 6) \\ C ​ = ​ k^{2} (4 k^{4} ​ + ​ 20 k^{3} ​ - ​ 37 k^{2} ​ + ​ 30 k ​ - ​ 10) \\ D ​ = ​ 4 k^{6} ​ + ​ 20 k^{5} ​ - ​ 18 k^{4} ​ + ​ 4 k^{3} \\ E ​ = ​ k^{6} ​ + ​ 5 k^{5} ​ - ​ 6 k^{4} ​ + ​ 2 k^{3} \end{array}

(24)

Eq. (24) may be used to obtain the amplitude-frequency response curves.

2.5 Amplitude-frequency response curves of inverting filter shaping circuits

The amplitude-frequency response curves of common inverting filters at different shaping parameter values, according to Eq. (6), are shown in Fig. 5(a). Fig. 5(b) shows the amplitude-frequency response curves of the improved inverting filters at different shaping parameter values according to Eq. (12). The amplitude-frequency response curves of the multiple feedback low-pass filters at different shaping parameter values, according to Eq. (18), are shown in Fig. 5(c). Similarly, Fig. 5(d) shows the amplitude-frequency response curves of third-order multiple feedback low-pass filters at different shaping parameter values according to Eq. (24).

Fig. 5.

(Color online) Frequency response curves at different shaping parameter values for the four types of filter shaping circuits. (a) Amplitude-frequency response curves of the common inverting filter shaping circuit; (b)Amplitude-frequency response curves of the improved inverting filter shaping circuit; (c) Amplitude-frequency response curves of the multiple feedback low-pass filter shaping circuit; (d) Amplitude-frequency response curves of the third-order multiple feedback low-pass filter shaping circuit.

Figure 5(a) shows that the common inverting filter shaping function model functions as a low-pass filter. When k is small, the amplitude-frequency response curve declines relatively slowly. As the value of k increases, the amplitude-frequency response curve declines faster. The closer the value of k is to zero, the smaller the cut-off frequency and the lower the noise suppression performance.

Figure 5(b) shows that the cut-off frequency of the improved inverting filter is lower than that of the common inverting filter. Thus, it performs better as a low-pass filter. Higher values of the shaping parameter k indicate better filtering performance.

From Fig. 5(c), we observe that the multiple feedback low-pass filter shaping function model behaves as a low-pass filter. A comparison with the common inverting filter reveals similar behavior.

As shown in Fig. 5(d), the cut-off frequency of the third-order multiple feedback low-pass filter is lower than that of the multiple feedback low-pass filter, which indicates better low-pass filtering performance. The amplitude of this filter is the lowest among the four filters.

According to Fig. 5, a large RC value, i.e., a large value of k leads to a small cut-off frequency of the amplitude-frequency response curve and better noise suppression performance by digital filter shaping. This is beneficial for amplitude extraction and energy resolution of the system, but it simultaneously widens the digital Gaussian output, produces smaller peaks, and causes a slower time response, which are not conducive for pulse pile-up identification. This affects pulse counting rate the system. On the contrary, a smaller value of k implies a larger cut-off frequency, which defeats the purpose of filter shaping and affects the effective extraction and analysis of pulse amplitudes.

3 Comparative analysis of digital shaping results

The time constant of the front-end circuit of the signal output of the Si-PIN detector system is 10 μs. Digital shaping is performed using different shaping methods and shaping parameters at a sampling rate of 40 MHz. Eqs. (3), (9), (15), and (21) are used to implement the four types of inverting digital filter shaping processes for nuclear pulse signals.

3.1 Shaping results with inverting digital filter

The original signal is shaped using the same parameter value (k = 10). The numerical recursive modes of the common inverting filter shaping, improved inverting filter shaping, multiple feedback low-pass filter shaping, and third-order multiple feedback low-pass filter shaping are used to process nuclear pulse signals, and the results are shown in Fig. 6(a). Considering the noise suppression performance and pulse amplitude extraction, the digital shaping outputs of the multi-feedback low-pass filter shaping circuit and the third-order multi-feedback low-pass filter shaping circuit at different shaping parameter values are selected for comparison and analysis. For the same signal, the digital shaping results of the above two methods when k = 5, 10, 20, are shown in Fig. 6(b) and Fig. 6(c), respectively.

Fig. 6.

(Color online) Shaping results of inverting digital filters. (a) Comparison of the digital shaping results; (b) Digital shaping outputs of the multi-feedback low-pass filter at different parameters; (c) Digital shaping outputs of the third-order multi-feedback low-pass filter at different parameters.

Figure 6(a) shows that at the same shaping parameter values, common inverting digital shaping output exhibits the largest amplitude, but the worst shaping symmetry and noise suppression performance. In comparison, improved inverting digital shaping performs better in terms of noise suppression and waveform symmetry. The third-order multi-feedback low-pass digital shaping output has the smallest pulse amplitude and the best noise suppression performance, whereas the multi-feedback low-pass digital shaping output exhibits the best noise suppression and pulse amplitude comprehensive performance.

Figure 6(b) indicates that the pulse amplitude, noise suppression, and peak position information of the multiple feedback digital shaping output differ at different shaping parameter values. As k increases, the noise suppression performance improves, the pulse amplitude becomes smaller, and the peak position of the shaping output pulse shifts to the right. The above conclusion is consistent with the theoretical analysis and the results of the amplitude-frequency response curve.

As shown in Fig. 6(c), the digital shaping results of the third-order multi-feedback low-pass filter resemble those in Fig. 6(b). However, in comparison with the multiple feedback digital shaping output, the third-order multiple feedback digital shaping output has an additional RC low-pass filter circuit at the front end, and therefore, under the same shaping parameter values, the third-order multiple feedback digital shaping output has better noise suppression performance. Additionally, the presence of the RC filter circuit lowers the amplitude of the output pulse amplitude slightly. Therefore, in actual inverse digital shaping, the multiple feedback low-pass filter shaping circuit is capable of appropriately increasing the shaping parameter values to obtain a better SNR, and the third-order multiple feedback filter shaping circuit can appropriately reduce the shaping parameter value. Thus, the premise of ensuring better noise suppression leads to a higher SNR.

3.2 Performance analysis of digital shaping methods

To compare and analyze the outputs of inverting digital shaping under the same shaping parameter values and different methods, or the same shaping methods and different shaping parameter values, the amplitude, symmetry, and SNR after inverting digital shaping are quantitatively evaluated.

From Fig. 6, we observe that the output pulse amplitude of the inverting digital shaping has different degrees of attenuation for different shaping parameters. We quantitatively evaluate the amplitude attenuation index by using Eq. (25) to evaluate the amplitude attenuation.

d = (1 - \frac{A}{A_{0}}) \times 100 %

(25)

Here, A₀ and A are the pulse amplitudes before and after shaping, respectively, and d is the amplitude attenuation coefficient. A smaller d indicates lower pulse amplitude attenuation.

Eq. (26) is used to quantitatively evaluate the symmetry of the digital shaping output.

δ = \frac{\sum_{i = 1}^{N} {(x_{0 - i} - x_{0 + i})}^{2}}{N}

(26)

Here, i = 1, … N, is an integer, x_0-i is the ith point on the left side of the signal peak position, x_0+i is the ith point on the right side of the signal peak position, and N is the range of symmetry, i. e., the number of points. δ is used to indicate the degree of symmetry and a smaller δ implies better symmetry. For example, the δ of a Gaussian pulse is 0.

Eq. (27) is used to quantitatively evaluate the SNR of the digital shaping output. A higher SNR indicates better noise suppression performance.

S N R = \frac{V_{max}}{\sqrt{(\sum_{i = 1}^{N} {(V_{i} - V_{i^{'}})}^{2}) / N}}

(27)

Here, V_max is the maximum value of the output pulse amplitude, Vi is the pulse amplitude value corresponding to the i^th point, $V_{i^{'}}$ is the pulse amplitude value of the i^th point obtained using polynomial fitting, and N is the number of calculation points.

Table. 1 shows the amplitude, symmetry, and SNR corresponding to Fig. 6. In particular, Tables. 1(a), 1(b), and 1(c) correspond to Figs. 6(a), 6(b), and 6(c), respectively.

The amplitude attenuation is accurately estimated by using the third pulse in Fig. 5 for calculating the amplitude attenuation parameters.

Table. 1(a) shows that the common inverting digital shaping method has the smallest amplitude attenuation index (14.15%), but its SNR (38.69) and symmetry (120.35) are the worst. In comparison, the third-order multiple feedback low-pass digital shaping method has the best SNR (52.30) and symmetry (0.61), but its amplitude attenuation index is poor (37.29%), while the multiple feedback low-pass digital shaping method exhibits the best synthesis performance.

Performance comparison of inverting digital filters.

Shaping methods	d	δ	SNR
[Performance comparison of different shaping methods.]
Common inverting	14.15%	120.35	38.69
Improved inverting	24.68%	9.57	40.30
Multiple feedback	26.01%	3.45	46.23
Third-order multiple	37.29%	0.61	52.30
Shaping parameter	d	δ	SNR
[Performance comparison of the multiple feedback digital shaping method.]
k=5	17.44%	13.09	40.61
k=10	26.01%	3.45	46.23
k=20	38.18%	0.018	53.57
[Performance comparison of the third-order multiple feedback digital shaping method.]
k=5	25.45%	2.67	45.42
k=10	37.29%	0.61	52.30
k=20	48.95%	0.0062	63.02

As evident from Table. 1(b), for the multiple feedback low-pass digital shaping output, the amplitude attenuation index increases from 17.44% to 38.18%, the symmetry index parameter decreases from 13.09 to 0.018, and the SNR index is increased from 40.61 to 53.57, as k increases. The value of k = 10 is recommended to obtain the best comprehensive performance.

Table. 1(c) shows that for the third-order multiple feedback low-pass digital shaping output, the amplitude attenuation index increases from 25.45% to 48.95%, the symmetry index parameter decreases from 2.67 to 0.0062, and the SNR index increases from 45.42 to 63.02 as k increases. After comprehensive consideration, k = 5 is found to be the best choice.

4 Conclusion

The digitalization of the common inverting filter shaping circuit, improved inverting filter shaping circuit, multiple feedback low-pass filter shaping circuit, and third-order multiple feedback low-pass filtering shaping circuits were achieved. The transfer functions of the four types of circuits in the Laplacian and Z domains were obtained. An analysis of the filter performance showed that all the four types of filter shaping circuits exhibited low-pass filter performance. The digital shaping model of the third-order multiple feedback low-pass filter is preferred in terms of noise suppression and waveform symmetry performance. In contrast, in terms of comprehensive performance of amplitude extraction, noise suppression, and waveform symmetry, the digital shaping model of the multiple feedback low-pass filter is preferred.

References

[1]

W.T. Evariste, S.U. Hong, Q. Yi et al.,

Design and simulation of Gaussian shaping amplifier made only with CMOS FET for FEE of particle detector

. Nucl. Sci. Tech. 21(5), 312-315 (2010). https://doi.org/10.13538/j.1001-8042/nst.21.