def make_sdss_color_grids(filename,l, b, delta_l, delta_b): wanted_columns = ["U","G","R","I","Z"] photometry = segue.select_stars_in_area(filename,l,b,delta_l,delta_b,wanted_columns) good_data = np.isfinite(photometry["U"]) & np.isfinite(photometry["G"]) & np.isfinite(photometry["R"]) & np.isfinite(photometry["I"]) & np.isfinite(photometry["Z"])
print >> sys.stderr, sys.argv[0].split("/")[-1] + ": " + str(err.msg)
print >> sys.stderr, "\t for help use --help"
return 2
## If no -i input set use the first argument as the fits filename
if not filename and len(args)>0:
filename = args[0]
else:
print "Bad things happened"
sys.exit(1) ## this should be improved
## Main Program Starts Here
## Read list of clusters
clusterfile = "../ClassMaterial/data/sdss_clusters.txt"
clusterdata = ascii.read(clusterfile, delimiter=",", data_start=0, names=["Name", "Width"])
for cluster in clusterdata:
coordinates = coord.GalacticCoordinates.from_name(cluster["Name"])
l = (coordinates.l.degrees+180.)*u.deg ## l returned by GalacticCoordinates.from_name are -180 to 180, so add 180 to get a 0-360 value
b = coordinates.b.degrees*u.deg
dl = 5.0*u.deg
db = 5.0*u.deg
data = segue.select_stars_in_area(filename, l, b, dl, db, "FEH_ADOP", "RV_ADOP", "DIST_ADOP", "L", "B")
print cluster["Name"], l, b, len(data)
if __name__ == "__main__":
sys.exit(main())
