Unfolding code
Bonner spheres are widely used in neutron spectrometry.
Getting back the neutron energy distribution requires unfolding codes.
Many codes based on different approaches exist, e.g.
- UMG: Unfolding with MAXED (MAXimum Entropy Deconvolution) and GRAVEL (a program based on SAND, Spectrum Analysis by Neutron Detectors)
- MITOM (MIguel TOMà¡s unfolding code)
- FRUIT (FRascati Unfolding Interactive Tool)
Unfortunately, most of them are black boxes.
So I decided to write a plain open source one in order to deepen my understanding of unfolding.
It is written in Root/C++ and brings in the following features:
- Educated guess spectrum
- combination of the analytical neutron components used in FRUIT (thermal, epithermal and fast)
- seed parameters derived from the IAEA neutron spectra compendium
- Minimization performed by Minuit
The normalization in FRUIT is a bit tricky.
An explanation is given here (look for the string "logbook deconvolution").
Results obtained so far are very satisfying, being in good agreement with other codes.
----------------------------------------------------------------------------
IMG ----------------------------------------------------------------------------
The file should soon be available...