/portal/howtos/man_howto_tutorial_sourcelist/man_howto_tutorial_sourcelist.shtmlprevious | /portal/howtos/man_howto_tutorial_introduction/man_howto_tutorial_introduction.shtmlhome | /portal/howtos/man_howto_tutorial_surfphot/man_howto_tutorial_surfphot.shtmlnext |
How many raw twilight flat frames satisfy the above criteria?
awe> qinstr = RawTwilightFlatFrame.instrument.name == 'OMEGACAM' awe> qfilter = RawTwilightFlatFrame.filter.name == 'OCAM_g_SDSS' awe> qchip = RawTwilightFlatFrame.chip.name == 'ESO_CCD_#92' awe> qim = RawTwilightFlatFrame.imstat.median > 1000 awe> qqf = RawTwilightFlatFrame.quality_flags == 0 awe> query = qinstr & qfilter & qchip & qim & qqf awe> print(len(query)) 2314
awe> t = [rtf.MJD_OBS for rtf in query] awe> y = [rtf.imstat.median/rtf.EXPTIME for rtf in query] or, quicker, awe> result = [(rtf.MJD_OBS, rtf.imstat.median/rtf.EXPTIME) for rtf in query] awe> t, y = zip(*result) awe> pylab.scatter(t, y)
awe> tfrac = [k - math.floor(k) for k in t] awe> pylab.clf() awe> pylab.scatter(tfrac, y)
awe> query = (RawBiasFrame.instrument.name == 'OMEGACAM') &\ (RawBiasFrame.DATE_OBS > datetime.datetime(2011,8,1)) &\ (RawBiasFrame.DATE_OBS < datetime.datetime(2012,8,1)) &\ (RawBiasFrame.chip.name == 'ESO_CCD_#96') &\ (RawBiasFrame.quality_flags == 0) &\ (RawBiasFrame.is_valid > 0) awe> res = [(b.DATE_OBS, b.imstat.median) for b in query] awe> dat, val = zip(*res) awe> pylab.clf() awe> pylab.scatter(dat, val)
/portal/howtos/man_howto_tutorial_sourcelist/man_howto_tutorial_sourcelist.shtmlprevious | /portal/howtos/man_howto_tutorial_introduction/man_howto_tutorial_introduction.shtmlhome | /portal/howtos/man_howto_tutorial_surfphot/man_howto_tutorial_surfphot.shtmlnext |