def loadClipAndCompFits(data, row):
epic = data[row, 0]
clip = dpc.loadClipboard("clips/c%09i-05.clip" %(epic))
clip = pl.serveTask(clip)
mp.figure(2)
mp.clf()
dpp.plotData(clip)
mp.title('Epic %.0f' %(epic))
mp.figure(3)
mp.clf()
compareFits(clip)
mp.title('Epic %.0f' %(epic))
mp.figure(1)
print "**************"
print row, epic
print clip.bls
flux = clip.detrend.flux_frac
flags = clip.detrend.flags
noi = noise.computeSgCdpp_ppm(flux[~flags]) * 1e-6
print "BLS SNR= %.2f" %(clip.bls.depth/noi)
def loadClipAndCompFits(data, row):
epic = data[row, 0]
clip = dpc.loadClipboard("clips/c%09i-05.clip" % (epic))
clip = pl.serveTask(clip)
mp.figure(2)
mp.clf()
dpp.plotData(clip)
mp.title('Epic %.0f' % (epic))
mp.figure(3)
mp.clf()
compareFits(clip)
mp.title('Epic %.0f' % (epic))
mp.figure(1)
print "**************"
print row, epic
print clip.bls
flux = clip.detrend.flux_frac
flags = clip.detrend.flags
noi = noise.computeSgCdpp_ppm(flux[~flags]) * 1e-6
print "BLS SNR= %.2f" % (clip.bls.depth / noi)
def displayOne(clip):
mp.figure(1)
mp.clf()
dpp.plotData(clip)
mp.figure(2)
mp.clf()
cb.compareFits(clip)
def displayOne(clip):
mp.figure(1)
mp.clf()
dpp.plotData(clip)
mp.figure(2)
mp.clf()
cb.compareFits(clip)
def summaryPlot1(output):
"""
Plot the data and some info from the clipboard
output is a clipboard.
"""
epicid=str(output['value'])
blsper=output['bls.period']
trapsnr=output['trapFit.snr']
trapper=output['trapFit.period_days']
trapdur=output['trapFit.period_days']
logTlpp = np.log10(output['lpp.TLpp'])
plt.clf()
plt.subplot(2,2,(1,2))
pp.plotData(output)
titlewords="%s P=%.2f" %(epicid,blsper)
plt.title(titlewords)
plt.subplot(223)
pp.plotFolded(output)
titlewords="P=%.2f d dur=%.2f h SNR=%.1f " % (trapper, trapdur,trapsnr)
plt.title(titlewords)
plt.xlim((0,trapper))
plt.subplot(224)
plt.cla()
try:
plt.plot(np.arange(1,142),output['lpp.binnedFlux'],'bo')
lab="logLPP=%.2f" % (logTlpp)
plt.xlim((-.9,141.9))
plt.title(lab)
except:
pass
try:
plt.figtext(.15,.96,output['disposition.fluxVet.disp'],color='r',fontsize=14)
plt.figtext(.12,.93,output['disposition.reasonForFail'],color='r',fontsize=10)
except:
pass
def plot_multipages(outfile, clip, intext, figtype="png"):
"""Create a line of text for the exporter
Inputs:
------------
outfile
A name for the pdf file to be written to
clip
A clipboard object
intext
A string containing information to be put on the first page.
Takes a clipboard, clip, and create plots
put these plots all into one multi paged document
specifieed by outfile
"""
plt.ioff()
dotperinch = 120
figuresize = (11, 8)
# The PDF document
pdfout = outfile + '.pdf'
pdf_pages = PdfPages(pdfout)
# Create a figure instance (ie. a new page)
#pdf_pages.attach_note(('KIC %u [%u]' % (clip.value, clip.disposition.isCandidate)),positionRect=[100,200,10,400])
fig = plt.figure(1, figsize=figuresize, dpi=dotperinch)
plt.figtext(0.2, 0.85, intext, color='r', fontsize=15)
plt.figtext(0.15, 0.2, clip.disposition, color='b', fontsize=14)
plt.title('Disposition Information for EPIC %u' % clip['value'])
if clip.disposition.isSignificantEvent:
steltxt = clip.get('stellar', defaultValue='No Stellar Available')
plt.figtext(0.7, 0.3, steltxt)
plantxt = clip.get('planet', defaultValue='No Planet Param. Avail.')
plt.figtext(0.7, 0.2, plantxt)
fig.patch.set_visible(False)
plt.gca().axis('off')
plt.savefig(pdf_pages, format='pdf')
outname = "%s-%s.%s" % (outfile, "info", figtype)
plt.savefig(outname)
plt.close()
fig = plt.figure(1, figsize=figuresize, dpi=dotperinch)
# Plot whatever you wish to plot
pp.summaryPlot1(clip)
# Done with the page
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "summary", figtype)
fig.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
# Plot whatever you wish to plot
pp.plotData(clip, nPanel=3)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "data", figtype)
fig.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.indivTransitPlot(clip, 6)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "indiv", figtype)
plt.savefig(outname)
plt.close()
try:
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.blsPlot(clip)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "bls", figtype)
plt.savefig(outname)
plt.close()
except AttributeError:
pass
try:
#Plot centroid plots
(fig2, fig1) = dip.plotWrapper(clip)
except ValueError:
fig1 = plt.figure()
fig2 = plt.figure()
fig1.set_size_inches(figuresize)
pdf_pages.savefig(fig1, dpi=dotperinch)
outname = "%s-%s.%s" % (outfile, "cent1", figtype)
fig1.savefig(outname)
fig2.set_size_inches(figuresize)
pdf_pages.savefig(fig2, dpi=dotperinch)
outname = "%s-%s.%s" % (outfile, "cent2", figtype)
fig2.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.lppDiagnostic(clip)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "lpp", figtype)
plt.savefig(outname)
plt.close()
# Write the PDF document to the disk
pdf_pages.close()
plt.close()
def plot_multipages(outfile, clip, intext, figtype="png"):
"""Create a line of text for the exporter
Inputs:
------------
outfile
A name for the pdf file to be written to
clip
A clipboard object
intext
A string containing information to be put on the first page.
Takes a clipboard, clip, and create plots
put these plots all into one multi paged document
specifieed by outfile
"""
plt.ioff()
dotperinch = 120
figuresize = (11, 8)
# The PDF document
pdfout = outfile + ".pdf"
pdf_pages = PdfPages(pdfout)
# Create a figure instance (ie. a new page)
# pdf_pages.attach_note(('KIC %u [%u]' % (clip.value, clip.disposition.isCandidate)),positionRect=[100,200,10,400])
fig = plt.figure(1, figsize=figuresize, dpi=dotperinch)
plt.figtext(0.2, 0.85, intext, color="r", fontsize=15)
plt.figtext(0.15, 0.2, clip.disposition, color="b", fontsize=14)
plt.title("Disposition Information for EPIC %u" % clip["value"])
if clip.disposition.isSignificantEvent:
steltxt = clip.get("stellar", defaultValue="No Stellar Available")
plt.figtext(0.7, 0.3, steltxt)
plantxt = clip.get("planet", defaultValue="No Planet Param. Avail.")
plt.figtext(0.7, 0.2, plantxt)
fig.patch.set_visible(False)
plt.gca().axis("off")
plt.savefig(pdf_pages, format="pdf")
outname = "%s-%s.%s" % (outfile, "info", figtype)
plt.savefig(outname)
plt.close()
fig = plt.figure(1, figsize=figuresize, dpi=dotperinch)
# Plot whatever you wish to plot
pp.summaryPlot1(clip)
# Done with the page
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "summary", figtype)
fig.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
# Plot whatever you wish to plot
pp.plotData(clip, nPanel=3)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "data", figtype)
fig.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.indivTransitPlot(clip, 6)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "indiv", figtype)
plt.savefig(outname)
plt.close()
try:
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.blsPlot(clip)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "bls", figtype)
plt.savefig(outname)
plt.close()
except AttributeError:
pass
try:
# Plot centroid plots
(fig2, fig1) = dip.plotWrapper(clip)
except:
fig1 = plt.plot()
fig2 = plt.plot()
fig1.set_size_inches(figuresize)
pdf_pages.savefig(fig1, dpi=dotperinch)
outname = "%s-%s.%s" % (outfile, "cent1", figtype)
fig1.savefig(outname)
fig2.set_size_inches(figuresize)
pdf_pages.savefig(fig2, dpi=dotperinch)
outname = "%s-%s.%s" % (outfile, "cent2", figtype)
fig2.savefig(outname)
plt.close()
fig = plt.figure(figsize=figuresize, dpi=dotperinch)
pp.lppDiagnostic(clip)
pdf_pages.savefig()
outname = "%s-%s.%s" % (outfile, "lpp", figtype)
plt.savefig(outname)
plt.close()
# Write the PDF document to the disk
pdf_pages.close()
plt.close()
