I’m currently trying to unfold the neutron spectrum using Bonner spheres spectrometry (which is essentially an under-determined problem). I have been using least_squares from scipy.optimize but I’m getting the feeling that there’s a smoothing step which is “flattening” my spectrum. The peaks are being broadened and I feel there is some sort of continuity conditions or smoothing steps that are being used to a achieve a somewhat less optimal results. Given that nature of the problem, there should be some narrow peaks (neutrons at a given energy value) for the number of counts/measurements to become accurate.