I need to do a simultaneous fit of the following spectral data. I need to fit skewed gaussians to a spectrum A, take the first and second derivative of this fitted skewed_gaussian spectrum and fit that two spectrums B and C each with their own coefficients. I would like to determine the residuals between the actual and fitted data and use this as a Cost function to minimize. I am trying to find the optimal way to do this fit. I thought of using minimize from scipy but I know it can be a bit slow. What is the best way to approach this? Continue reading...
