Unfolding analysis of LaBr₃:Ce gamma spectrum with a detector response matrix constructing algorithm based on energy resolution calibration

Unfolding analysis of LaBr 3 :Ce gamma spectrum with a detector response matrix constructing algorithm based on energy resolution calibration 1 2 3 first-author 2 3 corresp Corresponding author, tuoxg@cdut.edu.cn Corresponding author, tuoxg@cdut.edu.cn tuoxg@cdut.edu.cn 1 1 1 3 3 Fundamental Science on Nuclear Wastes and Environmental Safety Laboratory, Southwest University of Science and Technology, Mianyang 621010, China Fundamental Science on Nuclear Wastes and Environmental Safety Laboratory, Southwest University of Science and Technology, Mianyang 621010, China Sichuan University of Science & Engineering, Zigong 643000, China Sichuan University of Science & Engineering, Zigong 643000, China State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China With respect to the gamma spectrum, the energy resolution improves with increase in energy. The counts of full energy peak change with energy, and this approximately complies with the Gaussian distribution. This study mainly examines a method to deconvolve the LaBr 3 :Ce gamma spectrum with a detector response matrix constructing algorithm based on energy resolution calibration. In the algorithm, the full width at half maximum (FWHM) of full energy peak was calculated by the cubic spline interpolation algorithm and calibrated by a square root of a quadratic function that changes with the energy. Additionally, the detector response matrix was constructed to deconvolve the gamma spectrum. Furthermore, an improved SNIP algorithm was proposed to eliminate the background. In the experiment, several independent peaks of 152 Eu, 137 Cs, and 60 Co sources were detected by a LaBr 3 :Ce scintillator that were selected to calibrate the energy resolution. The Boosted Gold algorithm was applied to deconvolve the gamma spectrum. The results showed that the peak position difference between the experiment and the deconvolution was within ±2 channels and the relative error of peak area was approximately within 0.96%–6.74%. Finally, a 133 Ba spectrum was deconvolved to verify the efficiency and accuracy of the algorithm in unfolding the overlapped peaks. Detector response matrix Energy resolution calibration LaBr 3 :Ce scintillator SNIP background elimination Boosted Gold deconvolution algorithm 1. Introduction 1 Gamma spectrum analysis technology is an important method for qualitative and quantitative analyses of the gamma radiation nuclide. It plays a significant role in fields including nuclear technology and radiation protection. Generally, the gamma-ray spectrometers mainly include scintillator detectors and semiconductor detectors. The LaBr 3 :Ce scintillator was developed recently [ 1 1 ] when compared with the conventional NaI(Tl) scintillator that displays a higher energy resolution and higher photo yield [ 2 2 ]. However, the LaBr 3 :Ce scintillator still possess a few disadvantages including intrinsic radioactivity that is caused by the decay of 138 La. Conversely, the 138 La emits gamma-ray with energy corresponding to 788.7 keV through β - decay. In contrast, it yields gamma-rays with energy corresponding to 1435.7 keV and X-rays with energy corresponding to 32 keV by K-electron capture that potentially results in a sum-peak with energy corresponding to 1467.7 keV [ 3 3 ]. These features enhance the difficulty of gamma spectrum analysis. Additionally, with respect to the limited resolution of LaBr 3 :Ce scintillator, the overlapping peak is a common situation. Thus, peak decomposition is essential. Presently, the principle unfolding method of commercial gamma analysis software corresponds to peak function fitting, which involves fitting the spectrum data with a function and obtaining the peak area by integrating the function [ 4 4 , 5 5 ]. In terms of multiplets, the method depends on the form of peak function and initial values for iteration fitting. The program is significantly more complicated in most cases. From the signaling system viewpoint, the spectrum involves convolution by inputting an impulse signal with the detector response function. Therefore, the input radiation signal can be inverted through deconvolution. Thus, the function fitting is a forward method, and the deconvolution is a retrieval method. The earliest deconvolution algorithm corresponded to the direct demodulation method that was used for image restoration in high-energy astronomy by Li Tibei in 1993 [ 6 6 ] in which the Gaussian Seidel iteration method was applied. However, the algorithm caused large oscillations and produced negative values, and this is unacceptable in a gamma spectrum with positive data. A few other deconvolution algorithms were proposed and include Richardson–Lucy algorithm [ 7 7 ] and maximizes a posteriori deconvolution algorithm [ 8 8 ]. Nevertheless, they are unable to completely decompose the cluster of peaks [ 9 9 ]. A study by Miroslav Morháč [ 10 10 ] was the first to introduce the Gold deconvolution algorithm for gamma spectrum decomposition in 1997. The Gold algorithm is an iteration algorithm that was proposed by R. Gold in 1964 [ 11 11 ]. It is a positive and stable algorithm that was studied by Miroslav Morháč for several years [ 12 12 ]. However, with respect to the computer running theory, the Gold deconvolution algorithm leads to slow convergence and involves significant memory consumption, which is more obvious under the circumstances of big data and multi-dimensions [ 9 9 ]. In order to address these problems, Miroslav Morháč et al. in 2011 presented a boosted Gold deconvolution algorithm in which a boosting coefficient was introduced to enhance the operating efficiency [ 8 8 ]. A few researchers also adopted this algorithm for the gamma spectrum [ 13 13 , 14 14 ]. Essentially, the Gold deconvolution algorithm involves solving an equation set with an ill-posed problem through iterations [ 7 7 ], and the establishment of the response matrix is the crucial factor that affects the accuracy of the deconvolution results. Generally, the response matrix is a two-dimensional matrix that is constructed by a set of standard sources with numerous independent peaks, and the normalization of the peak data is treated as a response function in this energy region [ 15 15 ]. This method depends on the experimental peak data, and it is not affected by the form of response function although it requires several standard sources, and this is not satisfied in many laboratories. The Monte Carlo (MC) simulation can overcome the defects of the experiment and can save considerable time. However, it must accurately model the same because the veracity of input parameters and description for interaction mechanism influence the results, and thus, the MC results must be verified by performing an experiment [ 16 16 ]. Specifically, with respect to the gamma spectrum, the energy resolution changes with the energy. The Gaussian function is a good response function for the gamma spectrum detector, and the standard deviation of Gaussian function is determined by the energy resolution [ 17 17 ]. In the study, the detector response matrix was constructed by a series of Gaussian response functions that were established based on energy resolution calibration with a low amount of sources, and this was followed by performing the unfolding analysis of the LaBr 3 :Ce scintillator gamma spectrum. 2. Materials and methods 1 2.1 Gold deconvolution algorithm 2 The gamma spectrum y ( x ) is a convolution in which the gamma ray function f ( x ) is input with the detector response function h ( x ), and the f ( x ) is a δ -function (impulse function) in ideal conditions. The convolution is an integration that is expressed as follows: <math display="block" id="M1"><mrow><mi>y</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>h</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>*</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mrow><mo>+</mo><mi>∞</mi></mrow></msubsup><mi>h</mi></mrow></mstyle><mrow><mo>(</mo><mrow><mi>x</mi><mo>−</mo><mi>t</mi></mrow><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mtext>d</mtext><mi>t</mi><mo>.</mo></mrow></math> With respect to a discrete system, it is converted into a summation as follows: <math display="block" id="M2"><mrow><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>−</mo><mn>1</mn></mrow></msubsup><mi>h</mi></mstyle><mrow><mo>(</mo><mrow><mi>i</mi><mo>−</mo><mi>k</mi></mrow><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo> </mo><mo>,</mo><mtext> </mtext><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>…</mo><mo>,</mo><mtext> </mtext><mi>N</mi><mo>−</mo><mn>1</mn><mo>,</mo></mrow></math> where N denotes the length of f vector. Based on the convolution theory of discrete signals [ 18 18 ], if the length of f and h correspond to N , then the length of y corresponds to 2 N -1. Thus, the matrix form is as follows: <math display="block" id="M3"><mrow><mi>y</mi><mo>=</mo><mi>H</mi><mi>f</mi><mo>,</mo></mrow></math> where H denotes a (2 N -1)× N response matrix in which each column consists of h , y is a 2 N -1 column vector with zero padding, and f is a N column vector. Therefore, the matrix form is expressed as follows: <math display="block" id="M4"><mrow><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mi>N</mi><mo>−</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mo>=</mo><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mn>0</mn><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>3</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mo>.</mo></mrow></math> It is assumed that the response function and the output vector (spectrum) of the system given above are known. The deconvolution represents the solution of the overdetermined system with linear equations, i.e., the calculation of the vector f . The goal of the deconvolution method is to improve the resolution in the spectrum and decompose the multiplets. Based on all the methods that were studied, the Gold deconvolution is considered as the most stable method [ 19 19 ]. It is suitable to process positive definite data. Specifically, the method is described as follows: <math display="block" id="M5"><mtable columnalign="left"><mtr><mtd><mi>y</mi><mo>=</mo><mi>H</mi><mi>f</mi><mo>,</mo></mtd></mtr><mtr><mtd><msup><mi>H</mi><mi>T</mi></msup><mi>y</mi><mo>=</mo><msup><mi>H</mi><mi>T</mi></msup><mi>H</mi><mi>f</mi><mo>,</mo></mtd></mtr><mtr><mtd><mo> </mo><mtext>A</mtext><mo>=</mo><msup><mi>H</mi><mi>T</mi></msup><mi>H</mi><mo>,</mo><mtext> </mtext><msup><mi>y</mi><mo>'</mo></msup><mo>=</mo><msup><mi>H</mi><mi>T</mi></msup><mi>y</mi><mo>,</mo></mtd></mtr><mtr><mtd><msubsup><mi>f</mi><mi>i</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo>=</mo><mfrac><mrow><msubsup><mi>y</mi><mi>i</mi><mo>'</mo></msubsup></mrow><mrow><msubsup><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>m</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>−</mo><mn>1</mn></mrow></msubsup><msub><mi>A</mi><mrow><mi>i</mi><mi>m</mi></mrow></msub><msubsup><mi>f</mi><mi>m</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup></mrow></mfrac><msubsup><mi>f</mi><mi>i</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup><mo>,</mo><mtext> </mtext><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>⋯</mo><mo>,</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo>,</mo></mtd></mtr></mtable></math> where n denotes the iterations times, and the initial values <math id="M6"><mrow><msup><mi>f</mi><mrow><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mn>1</mn><mo>,</mo><mn>1</mn><mo>,</mo><mo>⋯</mo><mo>,</mo><mn>1</mn></mrow><mo>]</mo></mrow></mrow><mi>T</mi></msup></mrow></math> . Eq. (5) is the Gold deconvolution algorithm. The advantage of the Gold method is that the solution f is positive. However, the method does not concentrate a peak into the δ-function. Miroslav Morháč et al. [ 8 8 ] presented a high-resolution Boosted Gold deconvolution algorithm in which a boosting coefficient p was introduced to completely concentrate the area of the peak into a single channel. In simple terms, a loop was added, and p power of the results of Gold iteration was considered as the initial value into the next Gold iteration as follows: <math display="block" id="M7"><mrow><msubsup><mi>f</mi><mi>i</mi><mrow><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></mrow></msubsup><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><msubsup><mi>f</mi><mi>i</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup></mrow><mo>]</mo></mrow></mrow><mi>p</mi></msup><mo>.</mo></mrow></math> 2.2 Response matrix construction based on resolution calibration 2 The detector response function is simplified as a full-energy peak function after eliminating the background in which the Compton scattering background is included. The change in the counts of full energy peak with energy approximately comply with the Gaussian distribution, and thus the Gaussian function is used as a detector response function [ 17 17 ]. The standard deviation of Gaussian function is determined by the energy resolution of the gamma spectrum. The energy resolution is defined as the full width at half maximum (FWHM) of the peak in the pulse height spectrum divided by its energy as follows: <math display="block" id="M8"><mrow><mi>η</mi><mo>=</mo><mfrac><mrow><mi>F</mi><mi>W</mi><mi>H</mi><mi>M</mi><mrow><mo>(</mo><mrow><mtext>keV</mtext></mrow><mo>)</mo></mrow></mrow><mi>E</mi></mfrac><mo>*</mo><mn>100</mn><mi>%</mi><mo>,</mo></mrow></math> where η denotes the energy resolution corresponding to the ability of the spectrometer to separate the different energy particles, and E denotes the energy of the gamma-rays in units keV. The energy resolution increases with decreases in the η value. The energy resolution improves with increases in the energy although the FWHM of the peak also increases [ 20 20 , 21 21 ], and thus, the FWHM reflects the energy resolution. The energy resolution is the most important indicator for scintillator gamma spectrometer and is mainly determined by four factors [ 22 22 , 23 23 ]: the statistics of the scintillator photon creation process (η ph ), fluctuation of the multiplication coefficient of the photomultiplier (η PMT ), uncertainty of photoelectric conversion efficiency (η phelect ), and the noise of the electronics (η noise ). Additionally, it is not possible to detect other factors (η other ), such as pulse accumulation, incident window thickness and amplitude-analyzer width per channel, in experiments, and they are treated as constants [ 23 23 ]. The first factor is inherent to the detector material and is reciprocally proportional to the square root of energy <math id="M9"><mrow><msub><mi>η</mi><mrow><mtext>ph</mtext></mrow></msub><mo>∝</mo><mn>1</mn><mo>/</mo><msqrt><mi>E</mi></msqrt></mrow></math> [ 24 24 ]. The second factor is experimentally determined based on the measured number of photoelectrons, and it depends on the light output of the crystal and efficiency of photoelectron collection at the first dynode and is approximately reciprocally proportional to the energy <math id="M10"><mrow><msub><mi>η</mi><mrow><mtext>PMT</mtext></mrow></msub><mo>∝</mo><mn>1</mn><mo>/</mo><mi>E</mi></mrow></math> [ 24 24 ]. Finally, both the uncertainty of photoelectric conversion efficiency and the noise of the electronics are constants [ 23 23 ]. These factors are independent of each other, and thus, the energy resolution of the scintillator is expressed based on the cascading event as follows: <math display="block" id="M11"><mrow><mtext>η</mtext><mo>=</mo><msqrt><mrow><msubsup><mi>η</mi><mrow><mtext>ph</mtext></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>η</mi><mrow><mtext>PMT</mtext></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>η</mi><mrow><mtext>phelect</mtext></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>η</mi><mrow><mtext>noise</mtext></mrow><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>η</mi><mrow><mtext>other</mtext></mrow><mn>2</mn></msubsup></mrow></msqrt><mo>=</mo><msqrt><mrow><mi>a</mi><mo>+</mo><mfrac><mi>b</mi><mi>E</mi></mfrac><mo>+</mo><mfrac><mi>c</mi><mrow><msup><mi>E</mi><mn>2</mn></msup></mrow></mfrac><mtext> </mtext></mrow></msqrt><mo>.</mo></mrow></math> This corresponds to the energy resolution calibration where a , b , and c denote constants. Eq. (7) is combined with Eq. (8), and the relationship of FWHM with energy E is expressed as follows: <math display="block" id="M12"><mrow><mi>F</mi><mi>W</mi><mi>H</mi><mi>M</mi><mo>=</mo><msqrt><mrow><mi>a</mi><mo>⋅</mo><msup><mi>E</mi><mn>2</mn></msup><mo>+</mo><mi>b</mi><mo>⋅</mo><mi>E</mi><mo>+</mo><mi>c</mi></mrow></msqrt><mo>.</mo></mrow></math> In practice, the standard deviation σ of the Gaussian function exhibits the following relationship with FWHM : <math display="block" id="M13"><mrow><mi>σ</mi><mrow><mo>(</mo><mi>E</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><msqrt><mrow><mn>2</mn><mtext>ln</mtext><mn>2</mn></mrow></msqrt></mrow></mfrac><mi>F</mi><mi>W</mi><mi>H</mi><mi>M</mi><mo>.</mo></mrow></math> Hence, the standard deviation is not a constant because it is a function of the photon energy. With respect to the discrete gamma spectrum data, each channel corresponds to different energies. Thus, the detector response function with Gaussian function is expressed as follows: <math display="block" id="M14"><mrow><msub><mi>h</mi><mi>m</mi></msub><mrow><mo>(</mo><mrow><mi>i</mi><mo>,</mo><msub><mi>σ</mi><mi>m</mi></msub></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mi>π</mi></mrow></msqrt><msub><mi>σ</mi><mi>m</mi></msub></mrow></mfrac><mi>exp</mi><mrow><mo>(</mo><mrow><mo>−</mo><mfrac><mrow><msup><mrow><mrow><mo>(</mo><mrow><mi>i</mi><mo>−</mo><mi>u</mi></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup></mrow><mrow><mn>2</mn><msubsup><mi>σ</mi><mi>m</mi><mn>2</mn></msubsup></mrow></mfrac></mrow><mo>)</mo></mrow><mo>,</mo><mtext> </mtext><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>⋯</mo><mo>,</mo><mi>N</mi><mo>−</mo><mn>1.</mn><mtext> </mtext><mi>m</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>⋯</mo><mo>,</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo>,</mo></mrow></math> where σ m denotes the standard deviation of m channel, and u denotes the position of gamma peak. Equation (11) indicates that different energies correspond to different response functions. Thus, Eq. (4) is re-expressed as follows: <math display="block" id="M15"><mrow><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mi>N</mi><mo>−</mo><mn>3</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>y</mi><mo stretchy="false">(</mo><mn>2</mn><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mo>=</mo><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mn>0</mn><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mn>0</mn><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr><mtr><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>3</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mi>h</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mo>⋯</mo></mrow></mtd></mtr><mtr><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd><mtd><mo>⋮</mo></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mrow><mo>[</mo><mrow><mtable><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>0</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mo>⋅</mo></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>2</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr><mtd><mrow><mi>f</mi><mo stretchy="false">(</mo><mi>N</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow><mo>]</mo></mrow><mo>.</mo></mrow></math> in which the response matrix consists of different response functions based on the resolution calibration. 2.3 Background elimination 2 The background affects the accuracy of the deconvolution results unless the background can be removed from the measured spectrum [ 13 13 ]. The statistics sensitive nonlinear iterative peak-clipping (SNIP) algorithm is a baseline background elimination algorithm that was proposed by C.G. Ryan in 1988 [ 25 25 ]. It was proven as an efficient method to eliminate the complete background of a spectrum [ 26 26 ]. The vectors from <math id="M16"><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></msubsup></mrow></math> , <math id="M17"><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></msubsup></mrow></math> to <math id="M18"><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup></mrow></math> are calculated in a step by step manner in which n denotes a user-adjustable iteration parameter. Specifically, it denotes a width factor of the filter window. The new value in the channel i in the n -th iteration is obtained as follows: <math display="block" id="M19"><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup><mo>=</mo><mi>min</mi><mrow><mo>{</mo><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mi>n</mi></mrow><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup></mrow><mo>]</mo></mrow></mrow><mo>}</mo></mrow><mo>.</mo></mrow></math> After the iteration process, we obtain the resulting baseline spectrum. The parameter n strongly influences the shape of the estimated background [ 26 26 ]. If a lower value is set for n , then the background is excessively eliminated, and the background is eliminated incompletely if a high value is set for n . Miroslav Morháč improved the SNIP algorithm with a clipping window that is adaptive to peak regions widths [ 27 27 ]. It constructed a vector <math display="block" id="M20"><mrow><mi>r</mi><mo>=</mo><mrow><mo>[</mo><mrow><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mo>…</mo><mn>0</mn><mo>,</mo><msub><mi>w</mi><mn>1</mn></msub><mo>,</mo><mo>…</mo><msub><mi>w</mi><mn>1</mn></msub><mo>,</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>0</mn><mo>,</mo><msub><mi>w</mi><mn>2</mn></msub><mo>,</mo><msub><mi>w</mi><mn>2</mn></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>w</mi><mn>2</mn></msub><mo>,</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>0</mn><mo>,</mo><msub><mi>w</mi><mi>k</mi></msub><mo>,</mo><msub><mi>w</mi><mi>k</mi></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>w</mi><mi>k</mi></msub><mo>,</mo><mn>0</mn><mo>,</mo><mn>0</mn><mo>,</mo><mo>…</mo><mn>0</mn></mrow><mo>]</mo></mrow></mrow></math> of the same length with the spectrum in which <math id="M21"><mrow><msub><mi>w</mi><mi>k</mi></msub></mrow></math> represents the width of the k -th peak region. Then it is assumed that: <math display="block" id="M22"><mrow><mi>n</mi><mo>=</mo><mtext>max</mtext><mrow><mo>{</mo><mrow><msub><mi>w</mi><mi>j</mi></msub></mrow><mo>}</mo></mrow><mo>,</mo><mtext> </mtext><mi>j</mi><mo>∈</mo><mrow><mo>(</mo><mrow><mn>1</mn><mo>,</mo><mi>k</mi></mrow><mo>)</mo></mrow><mo>.</mo></mrow></math> The improved SNIP algorithm is defined as follows: <math display="block" id="M23"><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow></mrow></msubsup><mo>=</mo><mrow><mo>{</mo><mrow><mtable><mtr><mtd><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mrow><mo>{</mo><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo>,</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>[</mo><mrow><msubsup><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mi>n</mi></mrow><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo>+</mo><msubsup><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mi>n</mi></mrow><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup></mrow><mo>]</mo></mrow></mrow><mo>}</mo></mrow><mi>i</mi><mi>f</mi><mtext> </mtext><mi>n</mi><mo>≤</mo><mi>r</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mi>i</mi><mrow><mrow><mo>(</mo><mrow><mi>n</mi><mo>−</mo><mn>1</mn></mrow><mo>)</mo></mrow></mrow></msubsup><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mi>o</mi><mi>t</mi><mi>h</mi><mi>e</mi><mi>r</mi><mi>w</mi><mi>i</mi><mi>s</mi><mi>e</mi></mrow></mtd></mtr></mtable></mrow></mrow><mo>.</mo></mrow></math> This method is used to select the n adaptive to the peak width, although it evidently depends on the peak searching and peak-region width determination. Additionally, with respect to the overlapping peak-region, the width value w is not suitable for each peak. However, with respect to gamma spectrum, the width of peak changes with energy, and the scanning width is fixed as twice of the gamma peak FWHM of channel i [ 25 25 ]. Subsequently, the vector r is changed as follows: <math display="block" id="M24"><mrow><mi>r</mi><mo>=</mo><mi>r</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><mrow><mn>2</mn><mo>*</mo><mi>F</mi><mi>W</mi><mi>H</mi><msub><mi>M</mi><mi>i</mi></msub></mrow><mo>)</mo></mrow><mo>,</mo></mrow></math> where round denotes a function in MATLAB for round numbers, and factor 2 indicates the relation between the peak-width and FWHM based on a previous study [ 25 25 ]. Therefore, we estimate the precise value of n based on the energy resolution calibration and perform background subtraction that does not involve the peak-region. 3. Results and discussion 1 3.1 Experimental 2 The LaBr 3 :Ce scintillator exhibits high-energy resolution (≤3.2% @ 662 keV), high light yield (~ 60000 photons/MeV), and fast principal decay constant (≤ 25 ns), which is increasingly focused on in the field of radiation detection and displays potential as opposed to a NaI(Tl) detector [ 2 2 ]. However, the energy resolution of the LaBr 3 :Ce scintillator is still very low when compared with the HPGe semiconductor detector, and thus, the overlapping peaks commonly exist. In this study, the experiment was performed by a LaBr 3 :Ce scintillator detector that comprised of a Φ7.62 cm×7.62 cm LaBr 3 :Ce crystal with the BrilLanCeTM 380 series and a Φ7.62 cm PMT with a type of XP5300 manufactured by Saint-Gobain that corresponds to a multi-channel analyzer (MCA) with a type of DSP-jr2.0 manufactured by ORTEC and a lead room. The 152 Eu, 137 Cs, 60 Co, and 133 Ba gamma sources were detected. The 152 Eu, 137 Cs, 60 Co sources were used for resolution calibration because of independent peaks and the energy uniform distribution on the full spectrum. 133 Ba was used to test the unfolding capability. In this study, the procedures of unfolding analysis of LaBr 3 :Ce gamma spectrum are described as shown in Fig. 1. Fig. 1. 3.2 Deconvolution analysis and discussion 2 Several independent peaks of 152 Eu, 137 Cs, and 60 Co sources were selected to calibrate the energy resolution as shown in Table 1 Table 1 . As described in Section 2.2, the experimental data of FWHM was calibrated by <math id="M25"><mrow><mi>F</mi><mi>W</mi><mi>H</mi><mi>M</mi><mo>=</mo><msqrt><mrow><mi>c</mi><mo>+</mo><mi>b</mi><mi>E</mi><mo>+</mo><mi>a</mi><msup><mi>E</mi><mn>2</mn></msup></mrow></msqrt></mrow></math> . This function exhibits a friendly mathematical form although it is difficult to directly fit E and FWHM . We fitted the square of energy resolution η 2 and reciprocal of the energy 1/ E with a quadratic polynomial by using Eq. (8) to determine the constant coefficients a , b , and c . In this study, the fit results are shown in Fig. 2. Fig. 2. We obtained the following coefficients: a=0.0004, b=0.0603, and c=172.86. 1 1 Energy (keV) 1 1 Nuclide 1 1 FWHM(channel) 1 1 Experimental peak position 1 1 Deconvolution peak position 1 1 Experimental peak area 1 1 Deconvolution peak area 1 1 Peak area error (%) 1 1 121.78 1 1 152 Eu 1 1 13.66 1 1 106 1 1 107 1 1 47291 1 1 46480 1 1 -1.72 1 1 244.70 1 1 152 Eu 1 1 14.15 1 1 215 1 1 216 1 1 8820 1 1 9069 1 1 2.82 1 1 344.28 1 1 152 Eu 1 1 15.46 1 1 303 1 1 304 1 1 28837 1 1 28561 1 1 -0.96 1 1 443.97 1 1 152 Eu 1 1 18.61 1 1 391 1 1 392 1 1 1989 1 1 2123 1 1 6.74 1 1 661.66 1 1 137 Cs 1 1 20.03 1 1 583 1 1 583 1 1 55671 1 1 54456 1 1 -2.18 1 1 778.90 1 1 152 Eu 1 1 20.09 1 1 685 1 1 687 1 1 8116 1 1 7887 1 1 -2.82 1 1 964.08 1 1 152 Eu 1 1 24.08 1 1 849 1 1 848 1 1 7311 1 1 7682 1 1 5.07 1 1 1173.23 1 1 60 Co 1 1 23.68 1 1 1029 1 1 1030 1 1 6081 1 1 6228 1 1 2.42 1 1 1332.49 1 1 60 Co 1 1 33.27 1 1 1169 1 1 1169 1 1 5178 1 1 5408 1 1 4.44 After the energy resolution calibration, the response function shape changes with energy, as shown in Fig. 3 (a) Fig. 3 (a) . It is observed that the peak-width increases with the energy. Based on the convolution theory of discrete signal, the response vector was extended to the response matrix based on the energy resolution calibration as shown in Fig. 3 (b) Fig. 3 (b) . With respect to the independent peaks of 152 Eu, 137 Cs, and 60 Co, the results of deconvolution are shown in Fig. 4. Fig. 4. It is observed that the characteristic peaks efficiently converged to a single channel, namely the δ-function (impulse spectrum). Subsequently, two indicators of the peak position and peak area were examined to illustrate the effectiveness of qualitative and quantitative analyses, respectively. The results are shown in Table 1 Table 1 , and the experimental peak area is calculated by total peak area method that is expressed as follows: <math display="block" id="M26"><mrow><msub><mi>A</mi><mtext>n</mtext></msub><mo>=</mo><munderover><mstyle mathsize="140%" displaystyle="true"><mo>∑</mo></mstyle><mrow><mi>i</mi><mo>=</mo><mi>L</mi></mrow><mi>R</mi></munderover><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>−</mo><mrow><mo>(</mo><mrow><msub><mi>y</mi><mtext>R</mtext></msub><mo>+</mo><msub><mi>y</mi><mtext>L</mtext></msub></mrow><mo>)</mo></mrow><mtext>*</mtext><mfrac><mrow><mi>R</mi><mo>−</mo><mi>L</mi><mo>+</mo><mn>1</mn></mrow><mn>2</mn></mfrac><mo>,</mo></mrow></math> where <math id="M27"><mrow><msub><mi>A</mi><mtext>n</mtext></msub></mrow></math> denotes the net peak area, R and L denote the right boundary channel and left boundary channel, respectively, as determined by the first-order derivative peak searching algorithm, and y i denotes the counts of channel i , <math id="M28"><mrow><msub><mi>y</mi><mtext>R</mtext></msub></mrow></math> denotes the counts of channel R , and <math id="M29"><mrow><msub><mi>y</mi><mtext>L</mtext></msub></mrow></math> denotes the counts for channel L . From the results shown in Table 1 Table 1 , it is observed that the peak position difference between the deconvolved and the experimental is within ±2 channels. It demonstrated that the deconvolution does not significantly impact the peak position. The relative error of the deconvolution area with the experimental was approximately within 0.96%–6.74%. It is observed that the deconvolution peak area does not significantly differ from the experimental peak area. Subsequently, we illustrated the unfolding capability of overlapping peaks by using this algorithm. The 133 Ba source was detected in the same conditions by the LaBr 3 :Ce scintillator, and the response matrix was constructed based on the resolution calibration with the sources of 152 Eu, 137 Cs, and 60 Co. The 133 Ba source mainly has five obvious gamma peaks with energies corresponding 81.00 keV, 276.40 keV, 302.84 keV, 356.01 keV, and 383.85 keV. However, the energies corresponding to 276.40 keV and 302.84 keV and 356.01 keV and 383.85 keV overlap. The results of deconvolution are shown in Fig. 5 Fig. 5 and Table 2 Table 2 . It is observed that all the overlapping peaks were unfolded into an impulse spectrum and the peak position was accurate. The total peak area method (Eq. (18)) cannot be used for the overlapping peaks area calculation. The experimental peak area was calculated by the Gaussian function least-squares fitting method. The results show that the relative error of the deconvolution area with the fitting area was approximately within 0.81%–3.97%. Thus, we conclude that the deconvolution is efficient and accurate in unfolding the overlapping peaks. 1 1 Energy (keV) 1 1 FWHM(channel) 1 1 Experimental peak position 1 1 Deconvolved peak position 1 1 Fittingpeak area 1 1 Deconvolved peak area 1 1 Peak area error (%) 1 1 276.40 1 1 13.43 1 1 240 1 1 240 1 1 4291 1 1 4400 1 1 2.54 1 1 302.84 1 1 14.66 1 1 264 1 1 263 1 1 12871 1 1 12986 1 1 0.89 1 1 356.01 1 1 15.85 1 1 310 1 1 310 1 1 43229 1 1 42881 1 1 -0.81 1 1 383.85 1 1 14.60 1 1 335 1 1 335 1 1 6604 1 1 6866 1 1 3.97 In contrast, it should be noted that the peak of 1435.7 keV and the sum peak of 1467.7 keV emitted from 138 La in LaBr 3 :Ce materials are completely unfolded and distinguished from Fig. 4 Fig. 4 and Fig. 5. Fig. 5. This is important for gamma nuclide analysis of LaBr 3 :Ce detector when the nuclide’s gamma energy is close to 1467.7 keV such as 40 K (1460.82 keV). Additionally, as shown in Fig. 4 Fig. 4 and Fig. 5 Fig. 5 , several miscellaneous peaks are present in the deconvolution spectrum as opposed to characteristic peaks, and this was caused by the noise of the spectrum, although it did not evidently impact the characteristic peaks and it does not constitute the focus of the present study. 4. Conclusion 1 In this study, we presented a detector response matrix constructing algorithm to deconvolve the LaBr3:Ce gamma spectrum based on energy resolution calibration and an improved SNIP background elimination algorithm. The 152 Eu, 137 Cs, and 60 Co sources were detected to calculate the FWHM, and it was satisfactorily calibrated by the square root of a quadratic function. The detector response matrix was constructed based on the energy resolution calibration. The Boosted Gold algorithm was used to deconvolve the gamma spectrum. The results when combined with the detector response matrix constructing algorithm indicate that the Boosted Gold deconvolution algorithm is efficient in terms of dealing with the independent peaks and also with the overlapping peaks. It displays advantages in converging the peak area to a single channel and does not significantly influence the peak position. Furthermore, after the construction of the detector response matrix, it is not necessary for the deconvolution process to follow the fitting algorithm with respect to discriminating the overlapping peaks’ positions. The results demonstrate that the response matrix construction method is simplified and accurate based on the energy resolution calibration. References E. V. D. van Loef , P. Dorenbos , C. W. E. van Eijk et al ., High-energy-resolution scintillator: Ce3+ activated LaBr3 . Appl. Phys. Lett . 79 , 1573 - 1575 ( 2001 ). doi: 10.1063/1.1385342 http://doi.org/10.1063/1.1385342 High-energy-resolution scintillator: Ce3+ activated LaBr3 K. S. Shah , J. Glodo , M. Klugerman et al ., LaBr3: Ce scintillators for gamma-ray spectroscopy . IEEE Trans. Nucl. Sci . 50 , 2410 - 2413 ( 2003 ). doi: 10.1109/TNS.2003.820614 http://doi.org/10.1109/TNS.2003.820614 LaBr3: Ce scintillators for gamma-ray spectroscopy R. Nicolini , F. Camera , N. Blasi et al ., Investigation of the properties of a 1″×1″ LaBr3: Ce scintillator . Nucl. Instrum. Meth. A 582 , 554 - 561 ( 2007 ). doi: 10.1016/j.nima.2007.08.221 http://doi.org/10.1016/j.nima.2007.08.221 Investigation of the properties of a 1″×1″ LaBr3: Ce scintillator Z. Li , Y. W. Zhang , S. F. Sun et al ., Study on Response Function Model for γ-ray Scintillation Detector . Atom. Energy Sci. Tech . 49 , 1354 - 1358 ( 2015 ). doi: 10.7538/yzk.2015.49.08.1354 http://doi.org/10.7538/yzk.2015.49.08.1354 (in Chinese) Study on Response Function Model for γ-ray Scintillation Detector R. Shi , X. G. Tuo , H. L. Li et al ., 239Pu alpha spectrum analysis based on PIPS detector response function and variations with vacuum and distance . Nucl. Sci. Tech . 28 , 4 - 10 ( 2017 ). doi: 10.1007/s41365-016-0163-x http://doi.org/10.1007/s41365-016-0163-x 239Pu alpha spectrum analysis based on PIPS detector response function and variations with vacuum and distance T. P. Li , M. Wu , Reconstruction of objects by direct demodulation . Astrophys. Space Sci . 215 , 213 - 227 ( 1994 ). doi: 10.1007/BF00660079 http://doi.org/10.1007/BF00660079 Reconstruction of objects by direct demodulation L. B. Lucy , An iterative technique for the rectification of observed distributions . Astron. J . 79 , 745 - 754 ( 1974 ). doi: 10.1086/111605 http://doi.org/10.1086/111605 An iterative technique for the rectification of observed distributions M. Morháč , V. Matoušek , High-resolution boosted deconvolution of spectroscopic data . J. Comput. Appl. Math . 235 , 1629 - 1640 ( 2011 ). doi: 10.1016/j.cam.2010.09.005 http://doi.org/10.1016/j.cam.2010.09.005 High-resolution boosted deconvolution of spectroscopic data M. Morháč , V. Matoušek , Complete positive deconvolution of spectrometric data . Digit. Signal Process . 19 , 372 - 392 ( 2009 ). doi: 10.1016/j.dsp.2008.06.002 http://doi.org/10.1016/j.dsp.2008.06.002 Complete positive deconvolution of spectrometric data M. Morháč , J. Kliman , V. Matoušek et al ., Efficient one-and two-dimensional Gold deconvolution and its application to γ-ray spectra decomposition . Nucl. Instrum. Meth. A 401 , 385 - 408 ( 1997 ). doi: 10.1016/S0168-9002(97)01058-9 http://doi.org/10.1016/S0168-9002(97)01058-9 Efficient one-and two-dimensional Gold deconvolution and its application to γ-ray spectra decomposition R. Gold , An iterative unfolding method for response matrices . Math. Comput. Res. Develop. Rep . 6984 , 1 - 39 ( 1964 ). doi: 10.2172/4634295 http://doi.org/10.2172/4634295 An iterative unfolding method for response matrices M. Morháč , Deconvolution methods and their applications in the analysis of γ-ray spectra . Nucl. Instrum. Meth. A 559 , 119 - 123 ( 2006 ). doi: 10.1016/j.nima.2005.11.129 http://doi.org/10.1016/j.nima.2005.11.129 Deconvolution methods and their applications in the analysis of γ-ray spectra L. Li , X. G. Tuo , M. Z. Liu et al ., High-resolution boosted reconstruction of γ-ray spectra . Nucl. Sci. Tech . 25 , 18 - 24 ( 2014 ). doi: 10.13538/j.1001-8042/nst.25.050202 http://doi.org/10.13538/j.1001-8042/nst.25.050202 High-resolution boosted reconstruction of γ-ray spectra X. Y. Ai , Y. X. Wei , W. Y. Xiao , Direct Demodulation method for γ spectra analysis detected by CZT detectors . J. Tsinghua Univ. (Sci. & Tech.) . 46 , 821 - 824 ( 2006 ). doi: 10.3321/j.issn:1000-0054.2006.06.018 http://doi.org/10.3321/j.issn:1000-0054.2006.06.018 (in Chinese) Direct Demodulation method for γ spectra analysis detected by CZT detectors M. Jandel , M. Morháč , J. Kliman et al ., Decomposition of continuum γ-ray spectra using synthesized response matrix . Nucl. Instrum. Meth. A 516 , 172 - 183 ( 2004 ). doi: 10.1016/j.nima.2003.07.047 http://doi.org/10.1016/j.nima.2003.07.047 Decomposition of continuum γ-ray spectra using synthesized response matrix C. M. Salgado , C. C. Conti , P. H. B. Becker , Determination of HPGe detector response using MCNP5 for 20-150keV X-rays . Appl. Radiat. Isot . 64 , 700 - 705 ( 2006 ). doi: 10.1016/j.apradiso.2005.12.011 http://doi.org/10.1016/j.apradiso.2005.12.011 Determination of HPGe detector response using MCNP5 for 20-150keV X-rays Z. Li , X. G. Tuo , R. Shi et al ., A statistical approach to fit Gaussian part of full-energy peaks from Si (PIN) and SDD X-ray spectrometers . Sci. China Technol. Sc . 57 , 19 - 24 ( 2014 ). doi: 10.1007/s11431-013-5427-7 http://doi.org/10.1007/s11431-013-5427-7 A statistical approach to fit Gaussian part of full-energy peaks from Si (PIN) and SDD X-ray spectrometers M. Morháč , V. Matoušek , J. Kliman , Efficient algorithm of multidimensional deconvolution and its application to nuclear data processing . Digit. Signal Process . 13 , 144 - 171 ( 2003 ). doi: 10.1016/S1051-2004(02)00011-8 http://doi.org/10.1016/S1051-2004(02)00011-8 Efficient algorithm of multidimensional deconvolution and its application to nuclear data processing P. Bandžuch , M. Morháč , J. Krištiak , Study of the Van Cittert and Gold iterative methods of deconvolution and their application in the deconvolution of experimental spectra of positron annihilation . Nucl. Instrum. Meth. A 384 , 506 - 515 ( 1997 ). doi: 10.1016/S0168-9002(96)00874-1 http://doi.org/10.1016/S0168-9002(96)00874-1 Study of the Van Cittert and Gold iterative methods of deconvolution and their application in the deconvolution of experimental spectra of positron annihilation Z. Li , X. G. Tuo , R. Shi et al ., Analytic fitting and simulation methods for characteristic X-ray peaks from Si-PIN detector . Nucl. Sci. Tech . 24 , 060206 ( 2013 ). doi: 10.13538/j.1001-8042/nst.2013.06.007 http://doi.org/10.13538/j.1001-8042/nst.2013.06.007 Analytic fitting and simulation methods for characteristic X-ray peaks from Si-PIN detector Z. Li , X. G. Tuo , J. B. Yang et al ., Statistical distribution based detector response function of a Si (PIN) detector for Kα and Kβ X-ray . Chinese Phys. C 37 , 018202 ( 2013 ). doi: 10.1088/1674-1137/37/1/018202 http://doi.org/10.1088/1674-1137/37/1/018202 Statistical distribution based detector response function of a Si (PIN) detector for Kα and Kβ X-ray P. Dorenbos , J. T. M. de Haas , C. W. E. van Eijk , Non-proportionality in the scintillation response and the energy resolution obtainable with scintillation crystals . IEEE Trans. Nucl. Sci . 42 , 2190 - 2202 ( 1995 ). doi: 10.1109/23.489415 http://doi.org/10.1109/23.489415 Non-proportionality in the scintillation response and the energy resolution obtainable with scintillation crystals M. Moszyński , W. Czarnacki , W. Klamra et al ., Intrinsic energy resolution of pure NaI studied with large area avalanche photodiodes at liquid nitrogen temperatures . Nucl. Instrum. Meth. A 505 , 63 - 67 ( 2003 ). doi: 10.1016/S0168-9002(03)01021-0 http://doi.org/10.1016/S0168-9002(03)01021-0 Intrinsic energy resolution of pure NaI studied with large area avalanche photodiodes at liquid nitrogen temperatures M. Moszynski , J. Zalipska , M. Balcerzyk et al ., Intrinsic energy resolution of NaI(Tl) . Nucl. Instrum. Meth. A 484 , 259 - 269 ( 2002 ). doi: 10.1016/S0168-9002(01)01964-7 http://doi.org/10.1016/S0168-9002(01)01964-7 Intrinsic energy resolution of NaI(Tl) C. G. Ryan , E. Clayton , W. L. Griffin et al ., SNIP: a statistics-sensitive background treatment for the quantitative analysis of PIXE spectra in geoscience applications . Nucl. Instrum. Meth. B 34 , 396 - 402 ( 1988 ). doi: 10.1016/0168-583X(88)90063-8 http://doi.org/10.1016/0168-583X(88)90063-8 SNIP: a statistics-sensitive background treatment for the quantitative analysis of PIXE spectra in geoscience applications M. Morháč , V. Matoušek , Peak clipping algorithms for background estimation in spectroscopic data . Appl. Spectrosc . 62 , 91 - 106 ( 2008 ). doi: 10.1366/000370208783412762 http://doi.org/10.1366/000370208783412762 Peak clipping algorithms for background estimation in spectroscopic data M. Morháč , An algorithm for determination of peak regions and baseline elimination in spectroscopic data . Nucl. Instrum. Meth. A 600 , 478 - 487 ( 2009 ). doi: 10.1016/j.nima.2008.11.132 http://doi.org/10.1016/j.nima.2008.11.132 An algorithm for determination of peak regions and baseline elimination in spectroscopic data 10.1007/s41365-017-0340-6 This study was supported by the National Natural Science Foundation of China (Nos. 41374130 and 41604154). 2017-03-30 2017-05-20 2017-06-11 2017-12-19 2018-01-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-017-0340-6 http://dx.doi.org/10.1007/s41365-017-0340-6 . 2018 Nuclear Science and Techniques 01 29