Exemple #1
0
	def __init__(self, planet=None, star=None, instrument=None, directory=None,depthassumedforplotting=None, **kwargs):
		'''Initialize the parameters of the model.'''

		# setup the talker
		Talker.__init__(self)

		# create an empty array of parameters for eb
		self.ebparams = np.zeros(eb.NPAR, dtype=np.double)

		# keep track of a depth for plotting, if necessary
		self.depthassumedforplotting=depthassumedforplotting


		# load the model, if possible, and if a directory was given
		self.directory = directory
		if directory is not None and planet is None:
			self.load(directory)
		else:
			# define the subsets of the parameters, maybe by falling to defaults
			if planet is None:
				planet = Planet()
			if star is None:
				star = Star()
			if instrument is None:
				instrument = Instrument()
			self.planet = planet
			self.star = star
			self.instrument = instrument
Exemple #2
0
    def __init__(self, filename='test_marples.txt'):
        """initialize a Marples object, by reading it from a file"""

        # initialize the talker
        Talker.__init__(self)

        # read in the raw marple text file
        self.read(filename=filename)

        # populate a uniform grid (which will be mostly zeros) with definied marples
        self.gridify()
Exemple #3
0
	def __init__(self, filename=None, **kwargs):
		"""general initialization for PDF instances"""

		# initialize the talker
		Talker.__init__(self)

		# by default, PDF starts out uninitialized
		self._initialized = False

		# if possible, load the PDF from a file
		if filename is not None:
			self.load(filename)
Exemple #4
0
    def __init__(self, filename, **kwargs):
        '''Initialize a Relation object.'''

        # decide if it should be chatty
        Talker.__init__(self, pithy=True, **kwargs)

        # store the data filename

        self.filename = filename

        # load the data table
        self.load()
Exemple #5
0
	def __init__(self, directory=None, **kwargs):
		Talker.__init__(self)
		if directory is not None:
			self.load(directory)
		else:
			self.registered = []
			for key, value in kwargs.iteritems():
				# set up as a parameter
				self.__dict__[key] = Parameter(key, value)

				# fix everything initially (probably unnecessary?)
				self.__dict__[key].fix()

				# register this key in the parameters list
				self.registered.append(key)
Exemple #6
0
	def __init__(self,  name=None, ra=None, dec=None, **kw):
		'''Initialize the star, gathering its data from Simbad;
			if you supply...
				ra and dec (with astropy units attached)
				pmra and pmdec (in sky-projected mas/year) and epoch
			they will be used instead of the Simbad values.'''

		Talker.__init__(self)

		# decide how to initialize: use coordinates? or use simbad? or else?
		if ra is not None and dec is not None:
			# create the star from the input coordinates
			self.fromCoords(ra=ra, dec=dec, name=name, **kw)
		elif name is not None:
			# create the star from Simbad, searching on name
			self.fromSimbad(name, **kw)
		else:
			raise ValueError('please call {0}() either with "name" or "ra" + "dec" defined'.format(self.__class__.__name__))
Exemple #7
0
    def __init__(self, name=None, ra=None, dec=None, **kw):
        '''Initialize the star, gathering its data from Simbad;
			if you supply...
				ra and dec (with astropy units attached)
				pmra and pmdec (in sky-projected mas/year) and epoch
			they will be used instead of the Simbad values.'''

        Talker.__init__(self)

        # decide how to initialize: use coordinates? or use simbad? or else?
        if ra is not None and dec is not None:
            # create the star from the input coordinates
            self.fromCoords(ra=ra, dec=dec, name=name, **kw)
        elif name is not None:
            # create the star from Simbad, searching on name
            self.fromSimbad(name, **kw)
        else:
            raise ValueError(
                'please call {0}() either with "name" or "ra" + "dec" defined'.
                format(self.__class__.__name__))
Exemple #8
0
    def __init__(self, planet='WASP94Ab'):
        Talker.__init__(self)
        data = astropy.io.ascii.read(tablefile)
        match = (data['NAME'] == planet).nonzero()[0]
        assert (len(match == 1))
        self.row = data[match]
        self.mass = self.row['MASS'][0]  # in Jupiter
        self.radius = self.row['R'][0]  # in Jupiter
        self.period = self.row['PER'][0]
        self.a_over_r = self.row['AR'][0]
        self.stellar_teff = self.row['TEFF']
        self.stellar_mass = self.row['MSTAR']
        self.stellar_radius = self.row['RSTAR']

        self.teq = self.stellar_teff / np.sqrt(2 * self.a_over_r)
        self.mu = 2.2
        self.surfacegravity = u.G * self.mass * u.Mearth / (self.radius *
                                                            u.Rearth)**2
        self.scaleheight = u.k_B * self.teq / u.mp / self.mu / self.surfacegravity

        print self.__dict__
Exemple #9
0
 def __init__(self, planet, i, phasefrommidtransit=0.0):
     Talker.__init__(self)
     self.planet = planet
     self.midtransit = planet.period*(i + phasefrommidtransit) + planet.transit_epoch
     self.i = i
Exemple #10
0
    def __init__(self, marple, output='untrimmed'):
        Talker.__init__(self)
        self.marple = marple

        self.directory = output + '/'
        zachopy.utils.mkdir(self.directory)
Exemple #11
0
 def __init__(self, model, **kwargs):
     Talker.__init__(self)
     self.model = model
     self.model.lastfit = self
Exemple #12
0
    def __init__(self):

        Talker.__init__(self)
        self.period_max = 20.0
        self.read()
        self.gridify()
Exemple #13
0
    def __init__(self, filename='/Users/zkbt/Desktop/mo00113182+5908400_marples.txt'):

        Talker.__init__(self)
        self.period_max = 20.0
        self.read(filename=filename)
        self.gridify()
Exemple #14
0
 def __init__(self, pops=None):
     Talker.__init__(self)
     self.named = []
     self.pops = pops
Exemple #15
0
 def __init__(self, remake=False):
     Talker.__init__(self, line=80)
     self.db = star_to_files.db
     self.load(remake=remake)
Exemple #16
0
 def __init__(self):
     Talker.__init__(self)
     self.named = []
Exemple #17
0
 def __init__(self, *args, **kwargs):
     emcee.EnsembleSampler.__init__(self, *args, **kwargs)
     Talker.__init__(self)
     self.nametag = 'emcee'
     self.labels = None