Example notebook

Lorem ipsum dolor sit amet, consectetur adipiscing elit. Vivamus porttitor blandit eros quis tempor. Donec vel leo tincidunt, egestas odio sed, euismod nunc. Nullam odio dui, consequat sit amet enim non, fringilla iaculis est. Fusce varius tellus et dolor cursus tempor. Pellentesque dapibus convallis iaculis. Suspendisse tincidunt lacinia dui, ac rhoncus leo tempor quis. Aliquam erat volutpat. Mauris mauris arcu, malesuada at tellus eu, feugiat tincidunt tellus. Nam vel tempor lorem. Fusce bibendum orci eleifend, dictum elit efficitur, consectetur lacus. Etiam placerat, odio quis venenatis sollicitudin, turpis turpis euismod ante, eu dictum augue quam at felis. Proin dapibus condimentum fermentum. Nulla in maximus magna. Aenean at dolor tristique, pulvinar leo ut, bibendum sapien. Mauris sollicitudin condimentum pretium. Sed porta nisi nec condimentum rutrum.

Suspendisse lacus nibh, efficitur a tortor nec, consequat posuere turpis. Nullam a porta quam. Phasellus eleifend lacus quis arcu tristique, sit amet molestie leo venenatis. In hac habitasse platea dictumst. Lorem ipsum dolor sit amet, consectetur adipiscing elit. Aenean vehicula accumsan nulla, quis laoreet tortor lobortis suscipit. Morbi id turpis et eros viverra pretium id vel velit. Morbi erat mauris, auctor quis orci vitae, vehicula maximus erat. Donec convallis massa sit amet dictum porta. Lorem ipsum dolor sit amet, consectetur adipiscing elit.



In [25]:

    
%matplotlib inline

from astropy.io import fits as pyfits
import matplotlib.pyplot as plt
import seaborn as sns
import numpy as np

sns.set(style='white')



In [12]:

    
def wd2jd(wd):
    return (wd / 86400.) + 2453005.5



In [18]:

    
class FitsFile(object):
    def __init__(self, filename):
        with pyfits.open(filename) as infile:
            self.data = infile[1].data

        self.jd = wd2jd(self.data['tmid'])
        self.flux = self.data['tamflux2']
        
    def plot(self, axis, period, epoch):
        phase = ((self.jd - epoch) / period) % 1
        phase[phase > 0.8] -= 1.0
        
        axis.plot(phase, self.flux, '.', alpha=0.5)

WASP-12 b

We plot the phase folded lightcurve below:



In [26]:

    
f = FitsFile('example.fits')

fig, axis = plt.subplots()
f.plot(axis, 1.09142245, 2454508.97682)

Now a distribution analysis:



In [42]:

    
x = np.random.standard_normal(size=1024)
y = np.random.standard_exponential(size=1024)

sns.jointplot(x, y, kind='kde', stat_func=None)









    Out[42]:





<seaborn.axisgrid.JointGrid at 0x11591a590>