def __init__(self, name='logemission'):
self.fwhm = Parameter(name,
'fwhm',
100.,
tinyval,
hard_min=tinyval,
units="km/s")
self.pos = Parameter(name,
'pos',
5000.,
tinyval,
frozen=True,
units='angstroms')
self.flux = Parameter(name, 'flux', 1.)
self.skew = Parameter(name, 'skew', 1., tinyval, frozen=True)
self.limit = Parameter(name,
'limit',
4.,
alwaysfrozen=True,
hidden=True)
ArithmeticModel.__init__(
self, name,
(self.fwhm, self.pos, self.flux, self.skew, self.limit))
def __init__(self, name='blackbody'):
self.refer = Parameter(name,
'refer',
5000.,
tinyval,
hard_min=tinyval,
frozen=True,
units="angstroms")
self.ampl = Parameter(name,
'ampl',
1.,
tinyval,
hard_min=tinyval,
units="angstroms")
self.temperature = Parameter(name,
'temperature',
3000.,
tinyval,
hard_min=tinyval,
units="Kelvin")
self._argmin = 1.0e-3
self._argmax = 1000.0
ArithmeticModel.__init__(self, name,
(self.refer, self.ampl, self.temperature))
def __init__(self, name='brokenpowerlaw'):
self.refer = Parameter(name, 'refer', 5000., tinyval, hard_min=tinyval, frozen=True, units="angstroms")
self.ampl = Parameter(name, 'ampl', 1., tinyval, hard_min=tinyval, units="angstroms")
self.index1 = Parameter(name, 'index1', 0.1, -10.0, 10.0)
self.index2 = Parameter(name, 'index2', -0.1, -10.0, 10.0)
ArithmeticModel.__init__(self, name, (self.refer, self.ampl, self.index1, self.index2))
def __init__(self, name='shell2d'):
self.xpos = Parameter(name, 'xpos', 0) # p[0]
self.ypos = Parameter(name, 'ypos', 0) # p[1]
self.ampl = Parameter(name, 'ampl', 1) # p[2]
self.r0 = Parameter(name, 'r0', 1, 0) # p[3]
self.width = Parameter(name, 'width', 0.1, 0)
ArithmeticModel.__init__(self, name, (self.xpos, self.ypos, self.ampl, self.r0, self.width))
def __init__(self, name='IC'):
self.TFIR = Parameter(name, 'TFIR', 30, min=0, frozen=True, units='K')
self.uFIR = Parameter(name,
'uFIR',
0.0,
min=0,
frozen=True,
units='eV/cm3'
) # , 0.2eV/cm3 typical in outer disk
self.TNIR = Parameter(name, 'TNIR', 3000, min=0, frozen=True, units='K')
self.uNIR = Parameter(name,
'uNIR',
0.0,
min=0,
frozen=True,
units='eV/cm3'
) # , 0.2eV/cm3 typical in outer disk
# add ECPL params
super(InverseCompton, self).__init__(name=name)
# Initialize model
ArithmeticModel.__init__(self, name,
(self.index, self.ref, self.ampl, self.cutoff,
self.beta, self.TFIR, self.uFIR, self.TNIR,
self.uNIR, self.distance, self.verbose))
self._use_caching = True
self.cache = 10
def __init__(self, name='normbeta1d'):
self.pos = Parameter(name, 'pos', 0)
self.width = Parameter(name, 'width', 1, tinyval, hard_min=tinyval)
self.index = Parameter(name, 'index', 2.5, 0.5, 1000, 0.5)
self.ampl = Parameter(name, 'ampl', 1, 0)
ArithmeticModel.__init__(self, name,
(self.pos, self.width, self.index, self.ampl))
def __init__(self, name='recombination'):
self.refer = Parameter(name, 'refer', 5000., tinyval, hard_min=tinyval, frozen=True, units="angstroms")
self.ampl = Parameter(name, 'ampl', 1., tinyval, hard_min=tinyval, units="angstroms")
self.temperature = Parameter(name, 'temperature', 3000., tinyval, hard_min=tinyval, units="Kelvin")
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval, units="km/s")
ArithmeticModel.__init__(self, name, (self.refer, self.ampl, self.temperature, self.fwhm))
def __init__(self, name='edge'):
self.space = Parameter(name, 'space', 0, 0, 1, 0, 1,
'0 - energy | 1 - wave')
self.thresh = Parameter(name, 'thresh', 1, 0, hard_min=0)
self.abs = Parameter(name, 'abs', 1, 0)
ArithmeticModel.__init__(self, name,
(self.space, self.thresh, self.abs))
def __init__(self, name='emissionvoigt'):
self.center = Parameter(name,
'center',
5000.,
tinyval,
hard_min=tinyval,
frozen=True,
units="angstroms")
self.flux = Parameter(name, 'flux', 1.)
self.fwhm = Parameter(name,
'fwhm',
100.,
tinyval,
hard_min=tinyval,
units="km/s")
self.lg = Parameter(name,
'lg',
1.,
tinyval,
hard_min=tinyval,
frozen=True)
# Create core and wings from Gaussian and Lorentz
self._core = EmissionGaussian()
self._wings = EmissionLorentz()
ArithmeticModel.__init__(self, name,
(self.center, self.flux, self.fwhm, self.lg))
def __init__(self, name='linebroad'):
self.ampl = Parameter(name, 'ampl', 1, 0, hard_min=0)
self.rest = Parameter(name, 'rest', 1000, tinyval, hard_min=tinyval)
self.vsini = Parameter(name, 'vsini', tinyval, tinyval,
hard_min=tinyval)
ArithmeticModel.__init__(self, name,
(self.ampl, self.rest, self.vsini))
def __init__(self, name="IC"):
self.name = name
self.TFIR = Parameter(name, "TFIR", 30, min=0, frozen=True, units="K")
self.uFIR = Parameter(
name, "uFIR", 0.0, min=0, frozen=True, units="eV/cm3"
) # , 0.2eV/cm3 typical in outer disk
self.TNIR = Parameter(
name, "TNIR", 3000, min=0, frozen=True, units="K"
)
self.uNIR = Parameter(
name, "uNIR", 0.0, min=0, frozen=True, units="eV/cm3"
) # , 0.2eV/cm3 typical in outer disk
# add ECPL params
super(InverseCompton, self).__init__(name=name)
# Initialize model
ArithmeticModel.__init__(
self,
name,
(
self.index,
self.ref,
self.ampl,
self.cutoff,
self.beta,
self.TFIR,
self.uFIR,
self.TNIR,
self.uNIR,
self.distance,
self.verbose,
),
)
self._use_caching = True
self.cache = 10
def __init__(self, name='beta1d'):
self.r0 = Parameter(name, 'r0', 1, tinyval, hard_min=tinyval)
self.beta = Parameter(name, 'beta', 1, 1e-05, 10, 1e-05, 10)
self.xpos = Parameter(name, 'xpos', 0, 0, frozen=True)
self.ampl = Parameter(name, 'ampl', 1, 0)
ArithmeticModel.__init__(self, name,
(self.r0, self.beta, self.xpos, self.ampl))
def __init__(self, name='thetestmodel'):
self.notseen = Parameter(name, 'notseen', 10, min=5, max=15, hidden=True)
self.p1 = Parameter(name, 'P1', 1, min=0, max=10)
self.accent = Parameter(name, 'accent', 2, min=0, max=10, frozen=True)
self.norm = Parameter(name, 'norm', 100, min=0, max=1e4)
ArithmeticModel.__init__(self, name, (self.notseen, self.p1, self.accent, self.norm))
def __init__(self, name='disk2d'):
self.xpos = Parameter(name, 'xpos', 0) # p[0]
self.ypos = Parameter(name, 'ypos', 0) # p[1]
self.ampl = Parameter(name, 'ampl', 1) # p[2]
self.r0 = Parameter(name, 'r0', 1, 0) # p[3]
ArithmeticModel.__init__(self, name,
(self.xpos, self.ypos, self.ampl, self.r0))
def __init__(self, name='lorentz2d'):
self.fwhm = Parameter(name, 'fwhm', 10, tinyval, hard_min=tinyval)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name,
'ellip',
0,
0,
0.999,
0,
0.9999,
frozen=True)
self.theta = Parameter(name,
'theta',
0,
-2 * numpy.pi,
2 * numpy.pi,
-2 * numpy.pi,
4 * numpy.pi,
'radians',
frozen=True)
self.ampl = Parameter(name, 'ampl', 1)
ArithmeticModel.__init__(self, name,
(self.fwhm, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl))
self.cache = 0
def __init__(self, name='linebroad'):
self.ampl = Parameter(name, 'ampl', 1, 0, hard_min=0)
self.rest = Parameter(name, 'rest', 1000, tinyval, hard_min=tinyval)
self.vsini = Parameter(name, 'vsini', tinyval, tinyval,
hard_min=tinyval)
ArithmeticModel.__init__(self, name,
(self.ampl, self.rest, self.vsini))
def __init__(self, name="IC"):
self.name = name
self.TFIR = Parameter(name, "TFIR", 30, min=0, frozen=True, units="K")
self.uFIR = Parameter(
name, "uFIR", 0.0, min=0, frozen=True, units="eV/cm3"
) # , 0.2eV/cm3 typical in outer disk
self.TNIR = Parameter(
name, "TNIR", 3000, min=0, frozen=True, units="K"
)
self.uNIR = Parameter(
name, "uNIR", 0.0, min=0, frozen=True, units="eV/cm3"
) # , 0.2eV/cm3 typical in outer disk
# add ECPL params
super().__init__(name=name)
# Initialize model
ArithmeticModel.__init__(
self,
name,
(
self.index,
self.ref,
self.ampl,
self.cutoff,
self.beta,
self.TFIR,
self.uFIR,
self.TNIR,
self.uNIR,
self.distance,
self.verbose,
),
)
self._use_caching = True
self.cache = 10
def __init__(self, name="Bremsstrahlung"):
self.name = name
self.n0 = Parameter(
name, "n0", 1, min=0, max=1e20, frozen=True, units="1/cm3"
)
self.weight_ee = Parameter(
name, "weight_ee", 1.088, min=0, max=10, frozen=True
)
self.weight_ep = Parameter(
name, "weight_ep", 1.263, min=0, max=10, frozen=True
)
# add ECPL params
super().__init__(name=name)
# Initialize model
ArithmeticModel.__init__(
self,
name,
(
self.index,
self.ref,
self.ampl,
self.cutoff,
self.beta,
self.n0,
self.weight_ee,
self.weight_ep,
self.distance,
self.verbose,
),
)
self._use_caching = True
self.cache = 10
def __init__(self, name='DuleTelluric', same_b=False, same_d=False):
self.Cst_Cont = Parameter(name,
'Cst_Cont',
1.,
frozen=False,
min=0.95,
max=1.05)
self.t1_lam_0 = Parameter(name,
't1_lam_0',
5000.,
frozen=False,
min=0.0)
self.t1_b = Parameter(name, 't1_b', 3.5, frozen=False, min=1e-12)
self.t1_d = Parameter(name, 't1_d', 0.0005, frozen=False, min=0)
self.t1_tau_0 = Parameter(name, 't1_tau_0', 0.1, frozen=False, min=0.0)
self.t2_lam_0 = Parameter(name,
't2_lam_0',
5000.,
frozen=False,
min=0.0)
self.t2_b = Parameter(name, 't2_b', 3.5, frozen=False, min=1e-12)
self.t2_d = Parameter(name, 't2_d', 0.0005, frozen=False, min=0)
self.t2_tau_0 = Parameter(name, 't2_tau_0', 0.1, frozen=False, min=0.0)
self.resolution = 80000
if same_b: self.t1_b = self.t2_b
if same_d: self.t1_d = self.t2_d
ArithmeticModel.__init__(
self, name,
(self.Cst_Cont, self.t1_lam_0, self.t1_b, self.t1_d, self.t1_tau_0,
self.t2_lam_0, self.t2_b, self.t2_d, self.t2_tau_0))
def __init__(self, name='beta1d'):
self.r0 = Parameter(name, 'r0', 1, tinyval, hard_min=tinyval)
self.beta = Parameter(name, 'beta', 1, 1e-05, 10, 1e-05, 10)
self.xpos = Parameter(name, 'xpos', 0, 0, frozen=True)
self.ampl = Parameter(name, 'ampl', 1, 0)
ArithmeticModel.__init__(self, name,
(self.r0, self.beta, self.xpos, self.ampl))
def __init__(self, name="Bremsstrahlung"):
self.name = name
self.n0 = Parameter(
name, "n0", 1, min=0, max=1e20, frozen=True, units="1/cm3"
)
self.weight_ee = Parameter(
name, "weight_ee", 1.088, min=0, max=10, frozen=True
)
self.weight_ep = Parameter(
name, "weight_ep", 1.263, min=0, max=10, frozen=True
)
# add ECPL params
super(Bremsstrahlung, self).__init__(name=name)
# Initialize model
ArithmeticModel.__init__(
self,
name,
(
self.index,
self.ref,
self.ampl,
self.cutoff,
self.beta,
self.n0,
self.weight_ee,
self.weight_ep,
self.distance,
self.verbose,
),
)
self._use_caching = True
self.cache = 10
def __init__(self, name='recombination'):
self.refer = Parameter(name,
'refer',
5000.,
tinyval,
hard_min=tinyval,
frozen=True,
units="angstroms")
self.ampl = Parameter(name,
'ampl',
1.,
tinyval,
hard_min=tinyval,
units="angstroms")
self.temperature = Parameter(name,
'temperature',
3000.,
tinyval,
hard_min=tinyval,
units="Kelvin")
self.fwhm = Parameter(name,
'fwhm',
100.,
tinyval,
hard_min=tinyval,
units="km/s")
ArithmeticModel.__init__(
self, name, (self.refer, self.ampl, self.temperature, self.fwhm))
def __init__(self, name='edge'):
self.space = Parameter(name, 'space', 0, 0, 1, 0, 1,
'0 - energy | 1 - wave')
self.thresh = Parameter(name, 'thresh', 1, 0, hard_min=0)
self.abs = Parameter(name, 'abs', 1, 0)
ArithmeticModel.__init__(self, name,
(self.space, self.thresh, self.abs))
def __init__(self, name='normbeta1d'):
self.pos = Parameter(name, 'pos', 0)
self.width = Parameter(name, 'width', 1, tinyval, hard_min=tinyval)
self.index = Parameter(name, 'index', 2.5, 0.5, 1000, 0.5)
self.ampl = Parameter(name, 'ampl', 1, 0)
ArithmeticModel.__init__(self, name,
(self.pos, self.width, self.index, self.ampl))
def __init__(self, name='absorptionedge'):
self.edgew = Parameter(name, 'edgew', 5000., tinyval, frozen=True, units='angstroms')
self.tau = Parameter(name, 'tau', 0.5)
self.index = Parameter(name, 'index', 3.0, alwaysfrozen=True,
hidden=True)
ArithmeticModel.__init__(self, name,
(self.edgew, self.tau, self.index))
def __init__(self, name='absorptionlorentz'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.ewidth = Parameter(name, 'ewidth', 1.)
ArithmeticModel.__init__(self, name, (self.fwhm, self.pos, self.ewidth))
def __init__(self, name='logabsorption'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.tau = Parameter(name, 'tau', 0.5)
ArithmeticModel.__init__(self, name, (self.fwhm, self.pos,
self.tau))
def __init__(self, name='blackbody'):
self.refer = Parameter(name, 'refer', 5000., tinyval, hard_min=tinyval, frozen=True, units="angstroms")
self.ampl = Parameter(name, 'ampl', 1., tinyval, hard_min=tinyval, units="angstroms")
self.temperature = Parameter(name, 'temperature', 3000., tinyval, hard_min=tinyval, units="Kelvin")
self._argmin = 1.0e-3
self._argmax = 1000.0
ArithmeticModel.__init__(self, name, (self.refer, self.ampl, self.temperature))
def __init__(self, name='bpl1d'):
self.gamma1 = Parameter(name, 'gamma1', 0, -10, 10)
self.gamma2 = Parameter(name, 'gamma2', 0, -10, 10)
self.eb = Parameter(name, 'eb', 100, 0, 1000, 0)
self.ref = Parameter(name, 'ref', 1, frozen=True)
self.ampl = Parameter(name, 'ampl', 1, 1e-20)
ArithmeticModel.__init__(self, name,
(self.gamma1, self.gamma2,
self.eb, self.ref, self.ampl))
def __init__(self, name='bpl1d'):
self.gamma1 = Parameter(name, 'gamma1', 0, -10, 10)
self.gamma2 = Parameter(name, 'gamma2', 0, -10, 10)
self.eb = Parameter(name, 'eb', 100, 0, 1000, 0)
self.ref = Parameter(name, 'ref', 1, frozen=True)
self.ampl = Parameter(name, 'ampl', 1, 1e-20)
ArithmeticModel.__init__(self, name,
(self.gamma1, self.gamma2,
self.eb, self.ref, self.ampl))
def __init__(self, name='Sync'):
self.B = Parameter(name, 'B', 1, min=0, max=10, frozen=True, units='G')
# add ECPL params
super(Synchrotron, self).__init__(name=name)
# Initialize model
ArithmeticModel.__init__(self, name, (self.index, self.ref, self.ampl,
self.cutoff, self.beta, self.B,
self.distance, self.verbose))
self._use_caching = True
self.cache = 10
def __init__(self, name='pp'):
self.nh = Parameter(name, 'nH', 1, min=0, frozen=True, units='1/cm3')
# add ECPL params
super(PionDecay, self).__init__(name=name)
# Initialize model
ArithmeticModel.__init__(self, name, (self.index, self.ref, self.ampl,
self.cutoff, self.beta, self.nh,
self.distance, self.verbose))
self._use_caching = True
self.cache = 10
def __init__(self, name='emissionlorentz'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.flux = Parameter(name, 'flux', 1.)
self.kurt = Parameter(name, 'kurt', 2., frozen=True)
ArithmeticModel.__init__(self, name, (self.fwhm, self.pos,
self.flux, self.kurt))
def __init__(self, name='accretiondisk'):
self.ref = Parameter(name, 'ref', 5000., frozen=True, units='angstroms')
self.beta = Parameter(name, 'beta', 0.5, -10, 10)
self.ampl = Parameter(name, 'ampl', 1.)
self.norm = Parameter(name, 'norm', 20000.0, tinyval, alwaysfrozen=True,
hidden=True)
ArithmeticModel.__init__(self, name,
(self.ref, self.beta, self.ampl, self.norm))
def __init__(self, name='renamedpars'):
self.period = Parameter(name, 'period', 1, 1e-10, 10, tinyval)
self.offset = Parameter(name, 'offset', 0, 0, hard_min=0)
self.ampl = Parameter(name,
'ampl',
1,
1e-05,
hard_min=0,
aliases=['norm'])
ArithmeticModel.__init__(self, name,
(self.period, self.offset, self.ampl))
def __init__(self, name='fm'):
self.ebv = Parameter(name, 'ebv', 0.5) # E(B-V)
self.x0 = Parameter(name, 'x0', 4.6) # Position of Drude bump
self.width = Parameter(name, 'width', 0.06) # Width of Drude bump
self.c1 = Parameter(name, 'c1', 0.2)
self.c2 = Parameter(name, 'c2', 0.1)
self.c3 = Parameter(name, 'c3', 0.02)
self.c4 = Parameter(name, 'c4', 0.1)
ArithmeticModel.__init__(self, name,
(self.ebv, self.x0, self.width, self.c1,
self.c2, self.c3, self.c4))
def __init__(self,name='pp'):
self.index = Parameter(name , 'index' , 2.1 , min=-10 , max=10)
self.ref = Parameter(name , 'ref' , 60 , min=0 , frozen=True , units='TeV')
self.ampl = Parameter(name , 'ampl' , 100 , min=0 , max=1e60 , hard_max=1e100 , units='1e30/eV')
self.cutoff = Parameter(name , 'cutoff' , 0 , min=0 , frozen=True , units='TeV')
self.beta = Parameter(name , 'beta' , 1 , min=0 , max=10 , frozen=True)
self.nh = Parameter(name , 'nH' , 1 , min=0 , frozen=True , units='1/cm3')
self.verbose = Parameter(name , 'verbose' , 0 , min=0 , frozen=True)
ArithmeticModel.__init__(self,name,(self.index,self.ref,self.ampl,self.cutoff,self.beta,self.nh,self.verbose))
self._use_caching = True
self.cache = 10
def __init__(self, name='seaton'):
self.ebv = Parameter(name, 'ebv', 0.5)
self._xtab = numpy.array([0.0, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5,
1.6, 1.7, 1.8, 1.9, 2.0, 2.1,
2.2, 2.3, 2.4, 2.5, 2.6, 2.7])
self._extab = numpy.array([0.0, 1.36, 1.64, 1.84, 2.04, 2.24, 2.44,
2.66, 2.88, 3.14, 3.36, 3.56, 3.77,
3.96, 4.15, 4.26, 4.40, 4.52, 4.64])
ArithmeticModel.__init__(self, name, (self.ebv,))
def __init__(self, name='emissionvoigt'):
self.center = Parameter(name, 'center', 5000., tinyval, hard_min=tinyval, frozen=True, units="angstroms")
self.flux = Parameter(name, 'flux', 1.)
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval, units="km/s")
self.lg = Parameter(name, 'lg', 1., tinyval, hard_min=tinyval, frozen=True)
# Create core and wings from Gaussian and Lorentz
self._core = EmissionGaussian()
self._wings = EmissionLorentz()
ArithmeticModel.__init__(self, name, (self.center, self.flux, self.fwhm, self.lg))
def __init__(self, name='absorptiongaussian'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.ewidth = Parameter(name, 'ewidth', 1.)
self.limit = Parameter(name, 'limit', 4., alwaysfrozen=True,
hidden=True )
ArithmeticModel.__init__(self, name, (self.fwhm, self.pos,
self.ewidth, self.limit))
def __init__(self, name='fm'):
self.ebv = Parameter(name, 'ebv', 0.5) # E(B-V)
self.x0 = Parameter(name, 'x0', 4.6) # Position of Drude bump
self.width = Parameter(name, 'width', 0.06) # Width of Drude bump
self.c1 = Parameter(name, 'c1', 0.2)
self.c2 = Parameter(name, 'c2', 0.1)
self.c3 = Parameter(name, 'c3', 0.02)
self.c4 = Parameter(name, 'c4', 0.1)
ArithmeticModel.__init__(self, name, (self.ebv, self.x0,
self.width, self.c1, self.c2,
self.c3, self.c4))
def __init__(self, name='opticalgaussian'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.tau = Parameter(name, 'tau', 0.5)
self.limit = Parameter(name, 'limit', 4., alwaysfrozen=True,
hidden=True )
ArithmeticModel.__init__(self, name, (self.fwhm, self.pos,
self.tau, self.limit))
def __init__(self, name='hubblereynolds'):
self.r0 = Parameter(name, 'r0', 10, 0, hard_min=0)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name, 'ellip', 0, 0, 0.999, 0, 0.9999)
self.theta = Parameter(name, 'theta', 0, 0, 2*numpy.pi, -2*numpy.pi,
4*numpy.pi, 'radians')
self.ampl = Parameter(name, 'ampl', 1)
ArithmeticModel.__init__(self, name,
(self.r0, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl))
self.cache = 0
def __init__(self, name='hubblereynolds'):
self.r0 = Parameter(name, 'r0', 10, 0, hard_min=0)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name, 'ellip', 0, 0, 0.999, 0, 0.9999)
self.theta = Parameter(name, 'theta', 0, -2*numpy.pi, 2*numpy.pi, -2*numpy.pi,
4*numpy.pi, 'radians')
self.ampl = Parameter(name, 'ampl', 1)
ArithmeticModel.__init__(self, name,
(self.r0, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl))
self.cache = 0
def __init__(self, name="voigtabsorptionline"):
self.lam_0 = Parameter(name, "lam_0", 5000.0, frozen=False, min=0.0)
self.b = Parameter(name, "b", 3.5, frozen=False, min=1e-12)
self.d = Parameter(name, "d", 0.0005, frozen=False, min=0)
self.N = Parameter(name, "N", 999, frozen=True, hidden=True, min=0.0)
self.f = Parameter(name, "f", 999, frozen=True, hidden=True, min=0.0)
self.tau_0 = Parameter(name, "tau_0", 0.1, frozen=False, min=0.0)
ArithmeticModel.__init__(
self, name,
(self.lam_0, self.b, self.d, self.N, self.f, self.tau_0))
def __init__(self, name='parametercase'):
self.period = Parameter(name, 'Period', 1, 1e-10, 10, tinyval)
self.offset = Parameter(name, 'Offset', 0, 0, hard_min=0)
self.ampl = Parameter(name, 'Ampl', 1, 1e-05, hard_min=0, aliases=["NORM"])
with warnings.catch_warnings(record=True) as warn:
warnings.simplefilter("always", DeprecationWarning)
pars = (self.perioD, self.oFFSEt, self.NORM)
validate_warning(warn)
self._basemodel = Sin()
ArithmeticModel.__init__(self, name, pars)
def __init__(self, name='VoigtLine'):
self.lam_0 = Parameter(name, 'lam_0', 5000., frozen=False, min=0.0)
self.b = Parameter(name, 'b', 3.5, frozen=False, min=1e-12)
self.d = Parameter(name, 'd', 0.0005, frozen=False, min=0)
self.N = Parameter(name, 'N', 999, frozen=True, hidden=True, min=0.0)
self.f = Parameter(name, 'f', 999, frozen=True, hidden=True, min=0.0)
self.tau_0 = Parameter(name, 'tau_0', 0.1, frozen=False, min=0.0)
ArithmeticModel.__init__(
self, name,
(self.lam_0, self.b, self.d, self.N, self.f, self.tau_0))
def __init__(self, name='parametercase'):
self.period = Parameter(name, 'Period', 1, 1e-10, 10, tinyval)
self.offset = Parameter(name, 'Offset', 0, 0, hard_min=0)
self.ampl = Parameter(name, 'Ampl', 1, 1e-05, hard_min=0, aliases=["NORM"])
with warnings.catch_warnings(record=True) as warn:
warnings.simplefilter("always", DeprecationWarning)
pars = (self.perioD, self.oFFSEt, self.NORM)
validate_warning(warn)
self._basemodel = Sin()
ArithmeticModel.__init__(self, name, pars)
def __init__(self, name='polynomial'):
pars = []
for i in xrange(6):
pars.append(Parameter(name, 'c%d' % i, 0, frozen=True))
pars[0].val = 1
pars[0].frozen = False
for p in pars:
setattr(self, p.name, p)
self.offset = Parameter(name, 'offset', 0, frozen=True)
pars.append(self.offset)
ArithmeticModel.__init__(self, name, pars)
def __init__(self, name='seaton'):
self.ebv = Parameter(name, 'ebv', 0.5)
self._xtab = numpy.array([
0.0, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1,
2.2, 2.3, 2.4, 2.5, 2.6, 2.7
])
self._extab = numpy.array([
0.0, 1.36, 1.64, 1.84, 2.04, 2.24, 2.44, 2.66, 2.88, 3.14, 3.36,
3.56, 3.77, 3.96, 4.15, 4.26, 4.40, 4.52, 4.64
])
ArithmeticModel.__init__(self, name, (self.ebv, ))
def __init__(self, name='sersic2d'):
self.r0 = Parameter(name, 'r0', 10, 0, hard_min=0)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name, 'ellip', 0, 0, 0.999, 0, 0.9999)
self.theta = Parameter(name, 'theta', 0, -2*numpy.pi, 2*numpy.pi, -2*numpy.pi,
4*numpy.pi, 'radians')
self.ampl = Parameter(name, 'ampl', 1)
self.n = Parameter(name,'n', 1, .1, 10, 0.01, 100, frozen=True )
ArithmeticModel.__init__(self, name,
(self.r0, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl, self.n))
self.cache = 0
def __init__(self, name='polynomial'):
pars = []
for i in xrange(6):
pars.append(Parameter(name, 'c%d' % i, 0, frozen=True))
pars[0].val = 1
pars[0].frozen = False
for p in pars:
setattr(self, p.name, p)
self.offset = Parameter(name, 'offset', 0, frozen=True)
pars.append(self.offset)
ArithmeticModel.__init__(self, name, pars)
def __init__(self, name='logemission'):
self.fwhm = Parameter(name, 'fwhm', 100., tinyval, hard_min=tinyval,
units="km/s")
self.pos = Parameter(name, 'pos', 5000., tinyval, frozen=True, units='angstroms')
self.flux = Parameter(name, 'flux', 1.)
self.skew = Parameter(name, 'skew', 1., tinyval, frozen=True)
self.limit = Parameter(name, 'limit', 4., alwaysfrozen=True,
hidden=True )
ArithmeticModel.__init__(self, name,
(self.fwhm, self.pos, self.flux,
self.skew, self.limit))
def __init__(self, name='bbody'):
self.space = Parameter(name,
'space',
0,
0,
1,
0,
1,
'0 - energy | 1 - wave',
alwaysfrozen=True)
self.kT = Parameter(name, 'kT', 1, 0.1, 1000, 0, 1e+10, 'keV')
self.ampl = Parameter(name, 'ampl', 1, 1e-20, 1e+20, 1e-20, 1e+20)
ArithmeticModel.__init__(self, name, (self.space, self.kT, self.ampl))
def __init__(self, name='sersic2d'):
self.r0 = Parameter(name, 'r0', 10, 0, hard_min=0)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name, 'ellip', 0, 0, 0.999, 0, 0.9999)
self.theta = Parameter(name, 'theta', 0, 0, 2*numpy.pi, -2*numpy.pi,
4*numpy.pi, 'radians')
self.ampl = Parameter(name, 'ampl', 1)
self.n = Parameter(name,'n', 1, .1, 10, 0.01, 100, frozen=True )
ArithmeticModel.__init__(self, name,
(self.r0, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl, self.n))
self.cache = 0
def __init__(self, name='beta2d'):
self.r0 = Parameter(name, 'r0', 10, tinyval, hard_min=tinyval)
self.xpos = Parameter(name, 'xpos', 0)
self.ypos = Parameter(name, 'ypos', 0)
self.ellip = Parameter(name, 'ellip', 0, 0, 0.999, 0, 0.9999,
frozen=True)
self.theta = Parameter(name, 'theta', 0, -2*numpy.pi, 2*numpy.pi, -2*numpy.pi,
4*numpy.pi, 'radians', True)
self.ampl = Parameter(name, 'ampl', 1)
self.alpha = Parameter(name, 'alpha', 1, -10, 10)
ArithmeticModel.__init__(self, name,
(self.r0, self.xpos, self.ypos, self.ellip,
self.theta, self.ampl, self.alpha))
self.cache = 0
def __init__(self, name='sm'):
self.ebv = Parameter(name, 'ebv', 0.5)
self._xtab = numpy.array([
0.00, 0.29, 0.45, 0.80, 1.11, 1.43, 1.82, 2.27, 2.50, 2.91, 3.65,
4.00, 4.17, 4.35, 4.57, 4.76, 5.00, 5.26, 5.56, 5.88, 6.25, 6.71,
7.18, 8.00, 8.50, 9.00, 9.50, 10.00
])
self._extab = numpy.array([
0.00, 0.16, 0.38, 0.87, 1.50, 2.32, 3.10, 4.10, 4.40, 4.90, 6.20,
7.29, 8.00, 8.87, 9.67, 9.33, 8.62, 8.00, 7.75, 7.87, 8.12, 8.15,
8.49, 9.65, 10.55, 11.55, 12.90, 14.40
])
ArithmeticModel.__init__(self, name, (self.ebv, ))
def __init__(self, name="Known Velocity Line"):
self.v_cloud = Parameter(name,
"v_cloud",
0.0,
frozen=False,
min=-200,
max=200)
self.b = Parameter(name, "b", 3.5, frozen=False, min=1e-12)
self.d = Parameter(name, "d", 0.0005, frozen=False, min=1e-12)
self.N = Parameter(name, "N", 999, frozen=False, hidden=False, min=0)
self.f = Parameter(name, "f", 999, frozen=True, hidden=False, min=0)
self.lab_lam_0 = Parameter(name, "lab_lam_0", 5000, frozen=True)
ArithmeticModel.__init__(
self, name,
(self.v_cloud, self.b, self.d, self.N, self.f, self.lab_lam_0))
def __init__(self, name='jdpileup'):
self.alpha = Parameter(name, 'alpha', 0.5, 0, 1, 0, 1)
self.g0 = Parameter(name, 'g0', 1, tinyval, 1, tinyval, 1, frozen=True)
self.f = Parameter(name, 'f', 0.95, 0.9, 1, 0, 1)
self.n = Parameter(name, 'n', 1, tinyval, 100, tinyval, 2048,
alwaysfrozen=True)
self.ftime = Parameter(name, 'ftime', 3.241, tinyval, 5, tinyval, 100,
'sec', alwaysfrozen=True)
self.fracexp = Parameter(name, 'fracexp', 0.987, 0, 1, 0, 1,
alwaysfrozen=True)
self.nterms = Parameter(name, 'nterms', 30, 1, 100, 1, 100,
alwaysfrozen=True)
self._results = None
ArithmeticModel.__init__(self, name,
(self.alpha, self.g0, self.f, self.n,
self.ftime, self.fracexp, self.nterms))
