Updates for spectroscopic data reductions

Updates to add iterative reweighting least squares to make the fitting m..ore robust to outliers. Added a gaussian process fit for the sensitivity function. Added fitting the full cross-correlation option to interidentify.
I'm not sure I fully understand your comment. Do you mean that I should put this logic in interfit and just not have interfit call self.fit()?