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()