def process_acs_field(brick, fieldnum): """Compute ACS SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_acs', 'f475w', 'f814w') try: A, header = compute_area("f475w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f814w", brick, fieldnum) except NoWeightmap: raise NoDataError sb475, e475 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb814, e814 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return { "name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f475w": sb475, "f814w": sb814, "f475w_err": e475, "f814w_err": e814 }
def process_uv_field(brick, fieldnum): """Compute UV SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_uv', 'f275w', 'f336w') try: A, header = compute_area("f275w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f336w", brick, fieldnum) except NoWeightmap: raise NoDataError sb275, e275 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb336, e336 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return { "name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f275w": sb275, "f336w": sb336, "f275w_err": e275, "f336w_err": e336 }
def process_ir_field(brick, fieldnum): """Compute IR SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_ir', 'f110w', 'f160w') try: A, header = compute_area("f110w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f160w", brick, fieldnum) except NoWeightmap: raise NoDataError sb110, e110 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb160, e160 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return { "name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f110w": sb110, "f160w": sb160, "f110w_err": e110, "f160w_err": e160 }
def process_uv_field(brick, fieldnum): """Compute UV SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_uv', 'f275w', 'f336w') try: A, header = compute_area("f275w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f336w", brick, fieldnum) except NoWeightmap: raise NoDataError sb275, e275 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb336, e336 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return {"name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f275w": sb275, "f336w": sb336, "f275w_err": e275, "f336w_err": e336}
def process_field(fieldname): """Compute SB for a single Brown HST field.""" data = load_photometry(fieldname, 'brown', 'f606w', 'f814w') image_path = brown_image_path(fieldname, "f606w") msk_path = brown_phot_path(fieldname, kind="msk") header = fits.getheader(image_path) A = compute_area(fits.getdata(msk_path), header) sb606, e606 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb814, e814 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) log.info(fieldname) log.info("R {:.4f} kpc".format(rad.kpc)) log.info("mu_606: {:.6f}".format(sb606)) log.info("mu_814: {:.6f}".format(sb814)) return {"name": fieldname, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f606w": sb606, "f814w": sb814, "f606w_err": e606, "f814w_err": e814}
def process_acs_field(brick, fieldnum): """Compute ACS SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_acs', 'f475w', 'f814w') try: A, header = compute_area("f475w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f814w", brick, fieldnum) except NoWeightmap: raise NoDataError sb475, e475 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb814, e814 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return {"name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f475w": sb475, "f814w": sb814, "f475w_err": e475, "f814w_err": e814}
def process_ir_field(brick, fieldnum): """Compute IR SB in this brick+field.""" fieldname = "phat_b{0:02d}_f{1:02d}".format(brick, fieldnum) data = load_photometry(fieldname, 'phat_ir', 'f110w', 'f160w') try: A, header = compute_area("f110w", brick, fieldnum) except NoWeightmap: try: A, header = compute_area("f160w", brick, fieldnum) except NoWeightmap: raise NoDataError sb110, e110 = compute_sb(data['m1'], data['e1'], data['cfrac'], A) sb160, e160 = compute_sb(data['m2'], data['e2'], data['cfrac'], A) ra0, dec0, rad = compute_field_coord(header) return {"name": fieldname, "brick": brick, "field": fieldnum, "ra": ra0, "dec": dec0, "radius": rad.kpc, "f110w": sb110, "f160w": sb160, "f110w_err": e110, "f160w_err": e160}