-
Notifications
You must be signed in to change notification settings - Fork 0
/
long_periods.py
85 lines (64 loc) · 4.01 KB
/
long_periods.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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
"""
This is a script whose purpose is to identify long periodic candidates
and let me study them.
"""
from __future__ import division
import os
import numpy as np
import matplotlib.pyplot as plt
from periodic_selector import long_periodic_selector, best_long_period
from orion_plot import OrionStarData
from tablemate_comparisons import ukvar_spread
from official_star_counter import maxvars_spread_per
from variables_data_filterer import variables_photometry
from lightcurve_maker import lc_and_phase_and_colors
from table_maker import make_megeath_class_column
# These can be broken into two groups:
# 1. Anyone with a period longer than 20 days from the official periodics (shorter than 20 days)
# 2. Long periods identified by hand greater than 50 days.
megeath_class_column = make_megeath_class_column()
path = os.path.expanduser("~/Dropbox/Bo_Tom/lightcurve_book/long_periods/")
path2 = os.path.expanduser("~/Dropbox/Bo_Tom/lightcurve_book/short_periods/")
def find_20_50_day_periodics_and_save_to_file():
blp = best_long_period(long_periodic_selector(maxvars_spread_per, max_period=100), max_period=100)
ukvar_blp = blp.where(np.in1d(blp.SOURCEID, ukvar_spread.SOURCEID))
print "There are {0} stars with nominal periods between 50 and 100.".format(len(ukvar_blp))
for sid, period in zip(ukvar_blp.SOURCEID, ukvar_blp.best_period):
sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0])
fig = lc_and_phase_and_colors(sd, period)
fig.ax_j_phase.set_title("ONCvar {}".format(sd.name))
fig.ax_j_lc.set_title(sid)
fig.savefig("{0}/{1}/{2}.png".format(path, "50_100", sd.name))
plt.close()
blp2 = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=20, max_period=50), min_period=20, max_period=50)
ukvar_blp2 = blp2.where(np.in1d(blp2.SOURCEID, ukvar_spread.SOURCEID))
print "There are {0} stars with nominal periods between 20 and 50.".format(len(ukvar_blp2))
for sid, period in zip(ukvar_blp2.SOURCEID, ukvar_blp2.best_period):
sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0])
fig = lc_and_phase_and_colors(sd, period)
fig.ax_j_phase.set_title("ONCvar {}".format(sd.name))
fig.ax_j_lc.set_title(sid)
fig.savefig("{0}/{1}/{2}.png".format(path, "20_50", sd.name))
plt.close()
blp3 = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=10, max_period=20), min_period=10, max_period=20)
ukvar_blp3 = blp3.where(np.in1d(blp3.SOURCEID, ukvar_spread.SOURCEID))
print "There are {0} stars with nominal periods between 10 and 20.".format(len(ukvar_blp3))
for sid, period in zip(ukvar_blp3.SOURCEID, ukvar_blp3.best_period):
sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0])
fig = lc_and_phase_and_colors(sd, period)
fig.ax_j_phase.set_title("ONCvar {}".format(sd.name))
fig.ax_j_lc.set_title(sid)
fig.savefig("{0}/{1}/{2}.png".format(path, "10_20", sd.name))
plt.close()
def find_min_max_day_periodics_and_write(min_period=2, max_period=10, path=path2):
blp = best_long_period(long_periodic_selector(maxvars_spread_per, min_period=min_period, max_period=max_period),
min_period=min_period, max_period=max_period)
ukvar_blp = blp.where(np.in1d(blp.SOURCEID, ukvar_spread.SOURCEID))
print "There are {0} stars with nominal periods between {1} and {2}.".format(len(ukvar_blp), min_period, max_period)
for sid, period in zip(ukvar_blp.SOURCEID, ukvar_blp.best_period):
sd = OrionStarData(variables_photometry, sid, ukvar_spread.UKvar_ID[ukvar_spread.SOURCEID==sid][0])
fig = lc_and_phase_and_colors(sd, period)
fig.ax_j_phase.set_title("ONCvar {} ({})".format(sd.name, megeath_class_column[ukvar_spread.SOURCEID==sid][0]))
fig.ax_j_lc.set_title("{0}, Q={1}".format(sid, int((ukvar_spread.autovar+ukvar_spread.strict)[[ukvar_spread.SOURCEID==sid]][0])))
fig.savefig("{0}/{1}_{2}/{3}.png".format(path, min_period, max_period, sd.name))
plt.close()