Bayesian full-spectrum analysis of high-resolution gamma-ray spectra with energy scale correction and MCMC uncertainty quantification: Dataset and software

##################################################################################################

This repository contains  MATLAB code and the utilized gamma-spectrometric data for the paper:

Bayesian full-spectrum analysis of high-resolution gamma-ray spectra with energy scale correction and MCMC uncertainty quantification
Gerd Wübbeler, Matthias Stein, Henrike Fleischhack, Stefan Röttger, Anja Honig, and Clemens Elster

DOI:   Wuebbeler et al 2025 Meas. Sci. Technol. https://doi.org/10.1088/1361-6501/adf138

##################################################################################################

Run "bayes_gamma_fsa.m" for an example calculation of the posterior distribution for the isotope activities for a given mixture spectrum.
The calculation may take several hours. The computing time benefits from using the MATLAB Parallel Computing Toolbox.

The script "import_gamma_fsa_data.m" can be used to import the gamma-ray data stored in the CSV files.
The generated data structures are saved in the file "PTB_Gamma_FSA_2025.mat" for convenience.

Required MATLAB Toolboxes
- Optimization Toolbox
- Parallel Computing Toolbox (optional, improves run time)

The code utilizes tools from the MCMC Diagnostics Toolbox available at https://users.aalto.fi/~ave/code/mcmcdiag/. (see 'mcmcdiag' subdirectory)

The subdirectory 'PTB_Gamma_data_2025' contains Excel and CSV files with the gamma-ray spectra and additional information such as livetimes, activities etc..
See the file 's2_description_of_dataset.pdf' in that directory for a detailed description of the files.
The file 's1_calculation_of_single_activities.pdf' in the same subdirectory provides a detailed description of the peak evalulation method applied at PTB.

Contact

In case of problems or questions please contact  `gerd.wuebbeler@ptb.de`.