Title :
Nonparametric bayesian inference in nuclear spectrometry
Author :
Barat, ?‰ric ; Dautremer, Thomas ; Montagu, Thierry
Abstract :
We address the problem of X/gamma-ray spectra estimation in the fields of nuclear physics. Bayesian estimation of experimental backgrounds has been studied in [1] involving splines. Since Dirichlet processes (DP) sit on discrete measures, they provide an appealing prior for photopeaks. On the other hand, in order to tackle the complexity of experimental backgrounds, we consider a Polya Tree Mixture (PTM) - with suitable parameters yielding distribution continuity - for which predictive densities exhibit better smoothness properties than a single Polya Tree. Furthermore, it is easy to introduce some physical Compton line approximation formula (e.g. Klein-Nishina) in the base measure of the Polya Tree, or some physically driven local modifications of the PTM prior parameters. As backgrounds depend on photopeaks locations, we propose a hierarchical model where the PTM is conditioned on the DP. We use a beta prior for the mixing proportion between the DP and the PTM. Energies are not directly observed due to detection devices noises which introduce a convolution of both discrete and continuous measures by an assumed gaussian kernel whose variance is an unknown linear function of energy. Thus, the proposed semiparametric model for experimental data becomes a hierarchical Polya Tree-Dirichlet mixture of normal kernels. Besides, observed energies are binned in an histogram introducing additional quantification noise. The quantities of interest are usually posterior functionals of the DP mixing distribution. This implies an inverse problem which is carried out in the framework of finite stick-breaking representation. Thanks to conjugacy, draws from the posterior DP and PTM are easily obtained. The approach yields to a global peaks/background separation while offering spectrum resolution enhancement. The method is illustrated on experimental HPGe spectra.
Keywords :
Gaussian processes; X-ray detection; X-ray spectroscopy; gamma-ray detection; gamma-ray spectroscopy; germanium radiation detectors; inverse problems; solid scintillation detectors; trees (mathematics); DP mixing distribution; Dirichlet process; HPGe spectra; Klein-Nishina; LaBr3 scintillators; Polya Tree-Dirichlet mixture; X-ray spectra; detection device noises; discrete measures; finite stick-breaking representation; gamma-ray spectra; gaussian kernel; global background separation; global peaks separation; histogram; inverse problem; linear energy function; nonparametric Bayesian inference; nuclear spectrometry; physical Compton line approximation formula; quantification noise; semiparametric model; smoothness properties; splines; Bayesian methods; Convolution; Energy measurement; Gaussian noise; Histograms; Inverse problems; Kernel; Noise measurement; Nuclear physics; Spectroscopy; Bayesian nonparametrics; Dirichlet mixture; Pólya tree; density estimation; nuclear physics; stick breaking priors; x/γ-ray spectra;
Conference_Titel :
Nuclear Science Symposium Conference Record, 2007. NSS '07. IEEE
Print_ISBN :
978-1-4244-0922-8
Electronic_ISBN :
1095-7863
DOI :
10.1109/NSSMIC.2007.4436469