Skip to content

Commit

Permalink
Merge 76474b9 into 0b21274
Browse files Browse the repository at this point in the history
  • Loading branch information
profxj committed Oct 26, 2017
2 parents 0b21274 + 76474b9 commit ff2d906
Showing 1 changed file with 6 additions and 4 deletions.
10 changes: 6 additions & 4 deletions linetools/guis/spec_widgets.py
Original file line number Diff line number Diff line change
Expand Up @@ -385,13 +385,15 @@ def on_key(self, event):
sign=1
# Amplitude
Aguess = np.max(self.spec.flux[pix]-lconti[pix])
Cguess = np.mean(self.spec.wavelength[pix])
Cguess = np.mean(self.spec.wavelength[pix].value)
sguess = 0.1*np.abs(self.adict['wv_1']-self.adict['wv_2'])
# Fit
g_init = models.Gaussian1D(amplitude=Aguess, mean=Cguess, stddev=sguess)
fitter = fitting.LevMarLSQFitter()
parm = fitter(g_init, self.spec.wavelength[pix].value,
sign*(self.spec.flux[pix]-lconti[pix]))
#QtCore.pyqtRemoveInputHook()
#pdb.set_trace()
#QtCore.pyqtRestoreInputHook()
parm = fitter(g_init, self.spec.wavelength[pix].value, sign*(self.spec.flux[pix]-lconti[pix]))
# Error
var = [fitter.fit_info['param_cov'][ii,ii] for ii in range(3)]
sig = np.sqrt(var) # amplitude, mean, stddev
Expand Down Expand Up @@ -420,7 +422,7 @@ def on_key(self, event):
# Message
mssg = 'Gaussian Fit: '
mssg = mssg+' :: Mean={:g}, Amplitude={:g}, sigma={:g}, flux={:g}'.format(
parm.mean.value, parm.amplitude.value, parm.stddev.value, flux)
parm.mean.value*self.spec.wavelength.unit, parm.amplitude.value, parm.stddev.value*self.spec.wavelength.unit, flux)
mssg = mssg+' :: sig(Mean)={:g}, sig(Amplitude)={:g}, sig(sigma)={:g}, sig(flux)={:g}'.format(
sig_dict['mean'], sig_dict['amplitude'], sig_dict['stddev'], min(sig_flux1,sig_flux2))
mssg = mssg+' :: EW ={:g}'.format(EW)
Expand Down

0 comments on commit ff2d906

Please sign in to comment.