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))