Commit f30292e4 authored by Alexandre Beelen's avatar Alexandre Beelen

Add Cowie 2017 tab5 (850micron)

parent 201a4e72
This diff is collapsed.
This diff is collapsed.
# Region file format: DS9 astropy/regions
fk5
circle(188.9640,62.3631,0.0041)
circle(189.2997,62.3700,0.0041)
circle(189.3783,62.2163,0.0041)
circle(188.9828,62.3775,0.0041)
circle(188.9443,62.3370,0.0041)
circle(189.0990,62.0596,0.0041)
circle(189.1134,62.1015,0.0041)
circle(189.1392,62.2357,0.0041)
circle(189.5759,62.1494,0.0041)
circle(189.2868,62.3672,0.0041)
circle(188.9535,62.3180,0.0041)
circle(188.9598,62.1783,0.0041)
circle(189.2972,62.2252,0.0041)
circle(189.1331,62.2874,0.0041)
circle(189.2799,62.2356,0.0041)
circle(189.1217,62.0483,0.0041)
circle(188.8121,62.2528,0.0041)
circle(189.0765,62.2640,0.0041)
circle(188.9719,62.2271,0.0041)
circle(189.1455,62.3232,0.0041)
circle(189.3392,62.1190,0.0041)
circle(189.3458,62.0943,0.0041)
circle(189.1606,62.0408,0.0041)
circle(189.1865,62.3275,0.0041)
circle(189.2566,62.1962,0.0041)
circle(189.2168,62.2072,0.0041)
circle(188.9149,62.2122,0.0041)
circle(189.5519,62.2487,0.0041)
circle(189.1362,62.1059,0.0041)
circle(189.5144,62.1088,0.0041)
circle(189.1920,62.2468,0.0041)
circle(189.0689,62.1174,0.0041)
circle(189.0731,62.3877,0.0041)
circle(189.3079,62.3073,0.0041)
circle(189.4215,62.2058,0.0041)
circle(189.0417,62.1129,0.0041)
circle(189.1410,62.0494,0.0041)
circle(189.4087,62.2923,0.0041)
circle(188.9145,62.2197,0.0041)
circle(189.2012,62.3519,0.0041)
circle(189.0422,62.3454,0.0041)
circle(189.0927,62.2724,0.0041)
circle(189.0671,62.2538,0.0041)
circle(188.9710,62.1582,0.0041)
circle(189.1495,62.1187,0.0041)
circle(188.8692,62.2412,0.0041)
circle(189.0853,62.1110,0.0041)
circle(189.5166,62.1446,0.0041)
circle(189.1483,62.2400,0.0041)
circle(189.1887,62.3388,0.0041)
circle(189.0887,62.2857,0.0041)
circle(189.2155,62.0842,0.0041)
circle(189.3020,62.1766,0.0041)
circle(189.3611,62.1398,0.0041)
circle(189.1148,62.2050,0.0041)
circle(188.9280,62.1408,0.0041)
circle(189.3290,62.2049,0.0041)
circle(189.4338,62.1317,0.0041)
circle(189.0710,62.3255,0.0041)
circle(189.5117,62.3150,0.0041)
circle(189.1445,62.0916,0.0041)
circle(189.1543,62.1479,0.0041)
circle(189.1210,62.1793,0.0041)
circle(189.1538,62.1991,0.0041)
circle(189.2440,62.1587,0.0041)
circle(189.5005,62.2714,0.0041)
circle(189.0840,62.0977,0.0041)
circle(189.5599,62.2396,0.0041)
circle(189.2309,62.0717,0.0041)
circle(188.9165,62.2450,0.0041)
circle(189.5300,62.1796,0.0041)
circle(189.3932,62.0910,0.0041)
circle(189.2608,62.2408,0.0041)
circle(189.5373,62.2405,0.0041)
circle(189.3672,62.3222,0.0041)
circle(189.1822,62.3494,0.0041)
circle(189.3002,62.2033,0.0041)
circle(189.2444,62.2475,0.0041)
circle(188.9297,62.1654,0.0041)
circle(189.2358,62.2021,0.0041)
circle(189.5077,62.2169,0.0041)
circle(189.0234,62.1438,0.0041)
circle(189.0965,62.0891,0.0041)
circle(189.2457,62.3709,0.0041)
circle(189.0358,62.2432,0.0041)
circle(189.1450,62.1614,0.0041)
circle(188.9196,62.2052,0.0041)
circle(189.2913,62.1444,0.0041)
circle(188.9585,62.2606,0.0041)
circle(189.1302,62.1660,0.0041)
circle(189.0369,62.1344,0.0041)
circle(189.4765,62.1831,0.0041)
circle(189.2565,62.3402,0.0041)
circle(189.2175,62.1122,0.0041)
circle(189.3315,62.1726,0.0041)
circle(188.9882,62.2903,0.0041)
circle(189.3684,62.2396,0.0041)
circle(189.0348,62.1479,0.0041)
circle(189.4228,62.3176,0.0041)
circle(189.3287,62.1158,0.0041)
circle(189.0040,62.3019,0.0041)
circle(189.3085,62.1990,0.0041)
circle(189.3874,62.1246,0.0041)
circle(189.2991,62.1886,0.0041)
circle(189.3227,62.1346,0.0041)
circle(189.3542,62.1490,0.0041)
circle(189.0526,62.1431,0.0041)
circle(189.0836,62.3191,0.0041)
circle(189.0019,62.1815,0.0041)
circle(189.1438,62.2114,0.0041)
circle(189.2092,62.1394,0.0041)
circle(189.2511,62.1527,0.0041)
circle(189.3482,62.2871,0.0041)
circle(189.2671,62.1319,0.0041)
circle(189.1938,62.1425,0.0041)
circle(189.3985,62.1824,0.0041)
circle(189.1012,62.1714,0.0041)
circle(189.3053,62.2627,0.0041)
circle(189.3095,62.2024,0.0041)
circle(189.2032,62.1878,0.0041)
circle(189.4147,62.2662,0.0041)
circle(189.2219,62.1880,0.0041)
circle(189.1560,62.2966,0.0041)
circle(189.1907,62.3125,0.0041)
circle(189.0558,62.1502,0.0041)
circle(189.0361,62.2142,0.0041)
circle(189.3439,62.3106,0.0041)
circle(189.0945,62.2749,0.0041)
circle(189.3792,62.1719,0.0041)
circle(189.2899,62.1316,0.0041)
circle(189.0796,62.1679,0.0041)
circle(189.3828,62.1465,0.0041)
circle(189.1757,62.2896,0.0041)
circle(189.2607,62.2173,0.0041)
circle(189.1227,62.2536,0.0041)
circle(189.2962,62.2122,0.0041)
circle(189.2312,62.1372,0.0041)
circle(189.0743,62.2355,0.0041)
circle(189.3185,62.2791,0.0041)
circle(189.3829,62.2707,0.0041)
circle(189.2205,62.3156,0.0041)
circle(189.0710,62.2087,0.0041)
circle(189.0480,62.1760,0.0041)
circle(189.0783,62.1502,0.0041)
circle(189.2163,62.2494,0.0041)
circle(189.0325,62.2471,0.0041)
circle(189.1787,62.2625,0.0041)
circle(189.2772,62.2144,0.0041)
circle(189.4225,62.2142,0.0041)
circle(189.4080,62.2328,0.0041)
circle(189.1644,62.2619,0.0041)
circle(189.0257,62.2104,0.0041)
circle(189.1385,62.1430,0.0041)
circle(189.1853,62.2725,0.0041)
circle(189.2270,62.1475,0.0041)
circle(189.0502,62.1957,0.0041)
circle(189.2624,62.1425,0.0041)
circle(189.0525,62.2060,0.0041)
circle(189.3832,62.2157,0.0041)
circle(189.3654,62.2850,0.0041)
circle(189.2647,62.2764,0.0041)
circle(189.0368,62.1955,0.0041)
circle(189.4033,62.2119,0.0041)
circle(189.2047,62.3039,0.0041)
circle(189.0559,62.2435,0.0041)
circle(189.3347,62.1842,0.0041)
circle(189.2611,62.2339,0.0041)
circle(189.1039,62.1582,0.0041)
circle(189.1144,62.2208,0.0041)
circle(189.1646,62.1681,0.0041)
circle(189.3501,62.2749,0.0041)
circle(189.2957,62.2416,0.0041)
circle(189.2779,62.2689,0.0041)
circle(189.2300,62.1733,0.0041)
circle(189.3049,62.1960,0.0041)
circle(189.2027,62.2044,0.0041)
circle(189.0874,62.1871,0.0041)
circle(189.3070,62.2137,0.0041)
circle(189.2397,62.2819,0.0041)
circle(189.2830,62.1911,0.0041)
circle(189.2551,62.2539,0.0041)
circle(189.1890,62.1967,0.0041)
circle(189.2390,62.2356,0.0041)
circle(189.2744,62.1982,0.0041)
circle(189.1455,62.2044,0.0041)
circle(189.2336,62.2150,0.0041)
......@@ -5,7 +5,7 @@ from astropy.coordinates import SkyCoord
import astropy.units as u
from regions import CircleSkyRegion, write_ds9
indir = Path('GOODSN/rlsaw')
indir = Path('GOODSN/raw')
fitsdir = Path('GOODSN/cat')
regdir = Path('GOODSN/reg')
......@@ -211,3 +211,55 @@ reg = [CircleSkyRegion(coord, fwhm) for coord in coords]
write_ds9(reg, regdir / f'{filename}.reg')
data.write(fitsdir / f'{filename}.fits', overwrite=True)
################################################################################
#
# Greve et al. 2008
#
path = indir / 'Cowie2017_tab5.dat'
filename = 'Cowie2017_850'
bibcode = '2017ApJ...837..139C'
name = 'Cowie et at. 2017'
instrument = 'SCUBA2@JCMT'
fwhm = 14.8 * u.arcsec
data = ascii.read(path,
format='fixed_width', comment='#', delimiter=' ', header_start=None, data_start=0,
names=['ID', 'Name',
'R.A. (J2000)', 'Dec. (J2000)',
'f850', 'e850', 'S/N',
'R.A. (J2000 accurate)', 'Dec. (J2000 accurate)',
'f850 (SMA)', 'f1.4GHz', 'K_s', 'z', 'comment_z', 'z_250'])
data.sort('ID')
for col in [key for key in data.keys() if 'f850' in key or 'e850' in key]:
data[col].unit = u.mJy
data['f1.4GHz']].unit = u.uJy
data['K_s'].unit = u.ABmag
# As SkyCoords object do not support item assignment....
ra = data['R.A. (J2000)']
dec = data['Dec. (J2000)']
mask = data['R.A. (J2000 accurate)'].mask
ra[~mask] = data['R.A. (J2000 accurate)'][~mask]
dec[~mask] = data['Dec. (J2000 accurate)'][~mask]
coords = SkyCoord(ra, dec, unit=(u.hourangle, u.deg))
data['_ra'] = coords.ra
data['_dec'] = coords.dec
data.meta['bibcode'] = bibcode
data.meta['name'] = name
data.meta['inst'] = instrument
reg = [CircleSkyRegion(coord, fwhm) for coord in coords]
# Meta is PER region, so no way to include global meta yet
# reg[0].meta = data.meta
write_ds9(reg, regdir / f'{filename}.reg')
data.write(fitsdir / f'{filename}.fits', overwrite=True)
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment