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