def guess(self, dep, *args, **kwargs):
pos = get_position(dep, *args)
param_apply_limits(pos, self.eb, **kwargs)
ref = guess_reference(self.ref.min, self.ref.max, *args)
param_apply_limits(ref, self.ref, **kwargs)
norm = guess_amplitude_at_ref(self.ref.val, dep, *args)
param_apply_limits(norm, self.ampl, **kwargs)
def get_xspec_position(y, x, xhi=None):
if xhi is not None:
if x[0] > x[-1] and xhi[0] > xhi[-1]:
lo = _hc / xhi
hi = _hc / x
x, xhi = lo, hi
else:
if x[0] > x[-1]:
x = _hc / x
return get_position(y, x, xhi)
def guess(self, dep, *args, **kwargs):
ampl = guess_amplitude(dep, *args)
pos = get_position(dep, *args)
fwhm = guess_fwhm(dep, *args)
param_apply_limits(pos, self.pos, **kwargs)
norm = fwhm['val'] * numpy.sqrt(numpy.pi) * \
numpy.exp(lgam(self.index.val - 0.5) - lgam(self.index.val))
for key in ampl.keys():
ampl[key] *= norm
param_apply_limits(ampl, self.ampl, **kwargs)
def guess(self, dep, *args, **kwargs):
pos = get_position(dep, *args)
fwhm = guess_fwhm(dep, *args)
param_apply_limits(pos, self.pos, **kwargs)
param_apply_limits(fwhm, self.fwhm, **kwargs)
norm = guess_amplitude(dep, *args)
if fwhm != 10:
aprime = norm['val'] * self.fwhm.val * numpy.pi / 2.
ampl = {'val': aprime,
'min': aprime / _guess_ampl_scale,
'max': aprime * _guess_ampl_scale}
param_apply_limits(ampl, self.ampl, **kwargs)
else:
param_apply_limits(norm, self.ampl, **kwargs)
