def systematics_552():
# Bad1 http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/04/17/f2/c6/p13/v1/PTF_201004172696_i_p_scie_t062817_u011575385_f02_p110002_c06.fits?center=127.403,19.6476deg&size=100px
# Bad2 http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/04/17/f2/c6/p13/v1/PTF_201004173180_i_p_scie_t073759_u011575280_f02_p110002_c06.fits?center=127.403,19.6476deg&size=100px
# Good http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/04/07/f2/c6/p13/v1/PTF_201004071439_i_p_scie_t032714_u011539562_f02_p110002_c06.fits?center=127.403,19.6476deg&size=100px
mjd_offset = 54832
lc1 = session.query(LightCurve).filter(LightCurve.objid == 552).one()
bad_light_curve = PTFLightCurve.fromDBLightCurve(lc1)
light_curves = session.query(LightCurve).filter(func.q3c_radial_query(LightCurve.ra, LightCurve.dec, lc1.ra, lc1.dec,30/3600.)).all()
bad_obs1 = 55303.26964
bad_obs2 = 55303.31804
good_obs = 55293.14391
# seeing, airmass, filename, mjd
imlist = np.genfromtxt("data/aas_552_imagelist.txt", skiprows=4, usecols=[11,12,20,25], dtype=[("seeing", float), ("airmass", float), ("filename", "|S100"), ("mjd", float)])
idx_sort = np.argsort(imlist["mjd"])
imlist = imlist[idx_sort]
print "Bad1:", imlist["filename"][imlist["mjd"] == 55303.26964]
print "Bad2:", imlist["filename"][imlist["mjd"] == 55303.31804]
print "Good:", imlist["filename"][imlist["mjd"] == 55293.14391]
return
plt.plot(imlist["mjd"], imlist["seeing"], "r.")
plt.show()
def systematics_9347():
# http://kanaloa.ipac.caltech.edu/ibe/search/ptf/dev/process?POS=129.568,19.6232
# one http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/05/15/f2/c6/p13/v1/PTF_201005152355_i_p_scie_t053906_u011486277_f02_p110004_c06.fits?center=129.568,19.6232deg&size=150px
# two http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/04/25/f2/c6/p13/v1/PTF_201004251929_i_p_scie_t043750_u011578017_f02_p110004_c06.fits?center=129.568,19.6232deg&size=150px
mjd_offset = 54832
lc1 = session.query(LightCurve).filter(LightCurve.objid == 9347).one()
bad_light_curve = PTFLightCurve.fromDBLightCurve(lc1)
light_curves = session.query(LightCurve).filter(func.q3c_radial_query(LightCurve.ra, LightCurve.dec, lc1.ra, lc1.dec,30/3600.)).all()
bad_light_curve.plot()
print [x.ra for x in light_curves]
def systematics_9347():
# http://kanaloa.ipac.caltech.edu/ibe/search/ptf/dev/process?POS=129.568,19.6232
# one http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/05/15/f2/c6/p13/v1/PTF_201005152355_i_p_scie_t053906_u011486277_f02_p110004_c06.fits?center=129.568,19.6232deg&size=150px
# two http://kanaloa.ipac.caltech.edu/ibe/data/ptf/dev/process/proc/2010/04/25/f2/c6/p13/v1/PTF_201004251929_i_p_scie_t043750_u011578017_f02_p110004_c06.fits?center=129.568,19.6232deg&size=150px
mjd_offset = 54832
lc1 = session.query(LightCurve).filter(LightCurve.objid == 9347).one()
bad_light_curve = PTFLightCurve.fromDBLightCurve(lc1)
light_curves = session.query(LightCurve).filter(func.q3c_radial_query(LightCurve.ra, LightCurve.dec, lc1.ra, lc1.dec,30/3600.)).all()
bad_light_curve.plot()
print [x.ra for x in light_curves]
print [x.dec for x in light_curves]
print lc1.ra, lc1.dec
return
