-
Notifications
You must be signed in to change notification settings - Fork 0
/
c4_plots.py
executable file
·63 lines (49 loc) · 1.86 KB
/
c4_plots.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
import logging, itertools, os
from datetime import date
import astropy.io.ascii as at
import matplotlib.pyplot as plt
from k2spin.config import *
from k2spin import plot
today = date.today().isoformat()
def plot_list(results_list):
"""
"""
res = at.read(base_path+"tables/"+results_list)
f = open("/home/stephanie/my_papers/hyadesk2/figure_sets/f8.tbl","w")
count = 1
for i, epic in enumerate(res["EPIC"]):
logging.info(epic)
outfilename = "ktwo{0}-c04_lc_analysis.png".format(epic)
plot.paper_lcs(epic,res[i])
plt.savefig(base_path+"plot_outputs/"+outfilename,bbox_inches="tight")
if ((epic==2107361051) or (epic==2107361050) or
(epic==210963067) or (epic==2109630670) or
(epic==210675409)):
# Use lc from smaller centroiding box for 210735105
# but daofind lc for 210963067
# 210675409 is too bright but still in my list somehow
# note I never ran 211037886
continue
elif epic==2109630671:
save_epic = 210963067
else:
save_epic = epic
figsetname = "f8_{0}.eps".format(count)
f.write("{0} & EPIC {1}\n".format(figsetname,save_epic))
plt.savefig("/home/stephanie/my_papers/hyadesk2/figure_sets/"+figsetname,bbox_inches="tight")
plt.close("all")
count += 1
f.close()
if __name__=="__main__":
plot_list("c4_lcs_aps_results_2015-12-18_comments.csv")
"""
lc_file = "ktwo210408563-c04.csv"
epic = "210408563"
ap = 5
res = at.read(base_path+"tables/c4_lcs_aps_results_2015-12-18.csv")
plot.paper_lcs(epic,res[4])
plt.savefig("/home/stephanie/my_papers/hyadesk2/sample_lc.eps",
bbox_inches="tight")
plt.savefig("/home/stephanie/Dropbox/plots_for_sharing/sample_lc.png",
bbox_inches="tight")
"""