def run_test(campaign):
"""Are entries in the target list of a given campaign on silicon?"""
fov = getKeplerFov(campaign)
targetlist_fn = os.path.join(TESTDIR,
"data",
"K2Campaign{0}targets.csv".format(campaign))
targetlist = np.genfromtxt(targetlist_fn,
delimiter=",", dtype=None, names=True)
ra, dec = targetlist["RA_J2000_deg"], targetlist["Dec_J2000_deg"]
mask = ~np.isnan(ra) & ~np.isnan(dec)
for idx in np.where(mask)[0][::500]: # Speed-up: check every 500th
# All the sources in the target list should be on silicon
assert(onSiliconCheck(ra[idx], dec[idx], fov))
# A source +20 deg away should not be on silicon
assert(not onSiliconCheck(ra[idx] + 20, dec[idx], fov))
def check_target(target,
first=0,
last=LAST_CAMPAIGN,
verbose=True,
create_plot=False,
print_speed=True):
"""
Returns
-------
visible_campaigns : list of int
List of campaign numbers in which the object is visible.
"""
K2fov.logger.disabled = True # Disable warnings about prelim fields
ephem = get_ephemeris_dataframe(target, first, last)
visible_campaigns = []
for c in range(first, last + 1):
visible = False
fovobj = K2fov.getKeplerFov(c)
campaign_ephem = ephem.loc[K2fov.getFieldInfo(c)["start"]:K2fov.
getFieldInfo(c)["stop"]]
if create_plot:
import matplotlib.pyplot as pl
from K2fov import projection
ph = projection.PlateCaree()
pl.figure()
fovobj.plotPointing(ph, showOuts=False)
pl.plot(campaign_ephem["ra"],
campaign_ephem["dec"],
color="black",
lw=3)
#plot_padding = 5
#pl.xlim([campaign_ephem["ra"].max() + plot_padding,
# campaign_ephem["ra"].min() - plot_padding])
#pl.ylim([campaign_ephem["dec"].min() - plot_padding,
# campaign_ephem["dec"].max() + plot_padding])
pl.xlabel('R.A. [degrees]')
pl.ylabel('Declination [degrees]')
pl.gca().invert_xaxis()
pl.minorticks_on()
pl.title("Visibility of {0} in K2 C{1}".format(target, c))
plot_fn = "{}-c{}.png".format(target.replace("/", "_"), c)
print("Writing {}".format(plot_fn))
pl.savefig(plot_fn)
pl.close()
for idx, row in campaign_ephem.iterrows():
if onSiliconCheck(row["ra"], row["dec"], fovobj):
logging.debug("{} is visible in C{}.".format(target, c))
logging.debug(campaign_ephem[["dra", "ddec"]])
visible = True
break
if visible:
visible_campaigns.append(c)
if verbose:
mag = campaign_ephem['mag']
if mag.dtype == float: # can also be the string "n/a"
min_mag, max_mag = mag.min(), mag.max()
else:
min_mag, max_mag = -99.9, -99.9
print("Object '{}' is visible in C{} ("
"mag {:.1f}..{:.1f}; "
"{:.1f}..{:.1f}\"/h; "
"ra {:.3f}..{:.3f}; "
"dec {:.3f}..{:.3f}).".format(
target, c, min_mag, max_mag,
campaign_ephem['motion'].min(),
campaign_ephem['motion'].max(),
campaign_ephem['ra'].min(),
campaign_ephem['ra'].max(),
campaign_ephem['dec'].min(),
campaign_ephem['dec'].max()))
continue
K2fov.logger.disabled = False
return visible_campaigns
def check_target(target, first=0, last=LAST_CAMPAIGN, verbose=True,
create_plot=False, print_speed=True):
"""
Returns
-------
visible_campaigns : list of int
List of campaign numbers in which the object is visible.
"""
K2fov.logger.disabled = True # Disable warnings about prelim fields
ephem = get_ephemeris_dataframe(target, first, last)
visible_campaigns = []
for c in range(first, last + 1):
visible = False
fovobj = K2fov.getKeplerFov(c)
campaign_ephem = ephem.loc[K2fov.getFieldInfo(c)["start"]:K2fov.getFieldInfo(c)["stop"]]
if create_plot:
import matplotlib.pyplot as pl
from K2fov import projection
ph = projection.PlateCaree()
pl.figure()
fovobj.plotPointing(ph, showOuts=False)
pl.plot(campaign_ephem["ra"], campaign_ephem["dec"],
color="black", lw=3)
#plot_padding = 5
#pl.xlim([campaign_ephem["ra"].max() + plot_padding,
# campaign_ephem["ra"].min() - plot_padding])
#pl.ylim([campaign_ephem["dec"].min() - plot_padding,
# campaign_ephem["dec"].max() + plot_padding])
pl.xlabel('R.A. [degrees]')
pl.ylabel('Declination [degrees]')
pl.gca().invert_xaxis()
pl.minorticks_on()
pl.title("Visibility of {0} in K2 C{1}".format(target, c))
plot_fn = "{}-c{}.png".format(target.replace("/", "_"), c)
print("Writing {}".format(plot_fn))
pl.savefig(plot_fn)
pl.close()
for idx, row in campaign_ephem.iterrows():
if onSiliconCheck(row["ra"], row["dec"], fovobj):
logging.debug("{} is visible in C{}.".format(target, c))
logging.debug(campaign_ephem[["dra", "ddec"]])
visible = True
break
if visible:
visible_campaigns.append(c)
if verbose:
mag = campaign_ephem['mag']
if mag.dtype == float: # can also be the string "n/a"
min_mag, max_mag = mag.min(), mag.max()
else:
min_mag, max_mag = -99.9, -99.9
print("Object '{}' is visible in C{} ("
"mag {:.1f}..{:.1f}; "
"{:.1f}..{:.1f}\"/h; "
"ra {:.3f}..{:.3f}; "
"dec {:.3f}..{:.3f}).".format(
target,
c,
min_mag,
max_mag,
campaign_ephem['motion'].min(),
campaign_ephem['motion'].max(),
campaign_ephem['ra'].min(),
campaign_ephem['ra'].max(),
campaign_ephem['dec'].min(),
campaign_ephem['dec'].max())
)
continue
K2fov.logger.disabled = False
return visible_campaigns