Peak center and area estimation in gamma-ray energy spectra using a Mexican-hat wavelet

Abstract

Wavelet analysis is commonly used to detect and localize peaks within a signal，such as in Gamma-ray energy spectra. This paper presents a peak area estimation method based on a new wavelet analysis. Another Mexican Hat Wavelet Signal (MHWS) named after the new MHWS is obtained with the convolution of a Gaussian signal and a MHWS. During the transform, the overlapping background on the Gaussian signal caused by Compton scattering can be subtracted because the impulse response function MHWS is a second-order smooth function, and the amplitude of the maximum within the new MHWS is the net height corresponding to the Gaussian signal height, which can be used to estimate the Gaussian peak area. Moreover, the zero-crossing points within the new MHWS contain the information of the Gaussian variance whose valve should be obtained when the Gaussian peak area is estimated. Further, the new MHWS center is also the Gaussian peak center. With that distinguishing feature, the channel address of a characteristic peak center can be accurately obtained which is very useful in the stabilization of airborne Gamma energy spectra. In particular, a method for determining the correction coefficient k is given, where the peak area is calculated inaccurately because the value of the scale factor in wavelet transform is too small. The simulation and practical applications show the feasibility of the proposed peak center and area estimation method.