In [ ]:

    

import matplotlib.pyplot as plt
from lightcurve_pipeline.utils.utils import SETTINGS
from lightcurve_pipeline.database.database_interface import engine
from lightcurve_pipeline.database.database_interface import session
from lightcurve_pipeline.database.database_interface import Stats


    

In [ ]:

    

# Make plot of 'poisson factor' for all lightcurves
stats = session.query(Stats).all()
poisson_factor_all = [item.poisson_factor for item in stats if item.poisson_factor is not None]
poisson_factor = [pf for pf in poisson_factor_all if pf < 10]
fig = plt.figure()
ax = fig.add_subplot(111)
ax.set_xlim((0, 10))
ax.hist(poisson_factor, bins=100)
ax.set_xlabel(r'stdev / $\mu$')
ax.set_title('Poisson Factor (all)')
ax.text(4, 4000, '{} lightcuves > 10 (max={})'.format(len(poisson_factor_all) - len(poisson_factor), max(poisson_factor_all)), size=10)
plt.savefig('poisson_histogram.png')
plt.close()


    

In [ ]:

    

# Make plot of 'poisson factor' for only composite lightcurves
stats = session.query(Stats).all()
stats = session.query(Stats).all()
poisson_factor_all = [item.poisson_factor for item in stats if item.poisson_factor is not None and 'composite' in item.lightcurve_path]
poisson_factor = [pf for pf in poisson_factor_all if pf < 10]
fig = plt.figure()
ax = fig.add_subplot(111)
ax.set_xlim((0, 10))
ax.hist(poisson_factor, bins=100)
ax.set_xlabel(r'stdev / $\mu$')
ax.set_title('Poisson Factor (composite only)')
ax.text(4, 1000, '{} lightcuves > 10 (max={})'.format(len(poisson_factor_all) - len(poisson_factor), max(poisson_factor_all)), size=10)
plt.savefig('poisson_histogram_composite.png')
plt.close()


    

In [ ]: