'*', '*', '*', '*', '*', '*', '*', '*', '*', '*', '*', '*'
]
data_J_b = [
'*', '*', '*', 4, '*', '*', '*', '*', 9, 10, '*', '*', '*', '*', 15, '*',
'*', '*', 19, '*', '*', 22, '*', 24, 25, 26, '*'
]
#import 850um data
ID = []
RA = []
DEC = []
flux_850 = []
for line in cat_850.readlines()[2:]:
tmp = line.split()
ID.append(tmp[0])
RA.append(co.convHMS(tmp[1]))
DEC.append(co.convDMS(tmp[2]))
flux_850.append(float(tmp[3]))
#import MIPS IDs
ID_MIPS = []
RA_MIPS = []
DEC_MIPS = []
for line in cat_MIPS.readlines()[2:]:
tmp = line.split()
if tmp[1][0:2] == '17':
ID_MIPS.append(tmp[0])
RA_MIPS.append(co.convHMS(tmp[1]))
DEC_MIPS.append(co.convDMS(tmp[2]))
else:
pass
sub2_1700 = data_1700[2] #but don't use these or you'll double count
Rband_1700 = data_1700[
3] #undetected in 24um have estimates of their 24um flux from the R-band mag
all_1700 = data_1700[4] #everything detected at 24um
width_1700 = np.logspace(0.05 + 0.15, 3.5 + 0.15, len(
data_1700[0])) - bins_1700
#850um data
ID_850 = []
RA_850 = []
DEC_850 = []
cat_850 = open('../h1549/4sig.gaia', 'r')
for line in cat_850.readlines()[2:]:
tmp = line.split()
ID_850.append(tmp[0])
RA_850.append(co.convHMS(tmp[1]))
DEC_850.append(co.convDMS(tmp[2]))
ID_850 = np.array(ID_850)
RA_850 = np.array(RA_850)
DEC_850 = np.array(DEC_850)
cat_850.close()
#proto-z LBG data for HS1549
ID = []
RA = []
DEC = []
RA_1549 = co.convHMS('15:51:52.604')
DEC_1549 = co.convDMS('19:11:22.90')
flux = []
z = []
cat = open('../h1549/MIPS/h1549_24um_proto_cat.gaia', 'r')
S) - 25452.43730 * erf(.6324555320 * np.sqrt(S)) + 25870
# From Geach et al. 2017 S2-CLS paper
def boost(snr):
return 1.0 + 0.2 * (snr / 5.0)**(-2.3)
def S2CLS_core(S, a, b, c):
return a * (S / b)**(c) * np.exp(-S / b)
#############
# HS1549 DATA
#############
RA_map_1549 = co.convHMS('15:51:53.582')
DEC_map_1549 = co.convDMS('19:10:50.30')
data_1549 = open('../h1549/flux_4sig.gaia', 'r')
ID_1549 = []
RA_1549 = []
DEC_1549 = []
flux_1549 = []
flux_1549_deboosted = []
snr_1549 = []
for line in data_1549.readlines()[2:]:
tmp = line.split()
if tmp[0] == '1': #SMA resolved source
ID_1549.append('1.1')
ID_1549.append('1.2')
ID_1549.append('1.3')
#sources without multi-lambda data
data_24_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*','*']#29???
data_ch4_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*',29]
data_ch2_b = ['*','*','*','*',5,6,7,'*','*','*','*','*','*',14,'*','*','*',18,19,20,'*','*','*','*','*','*','*','*',29]
data_J_b = ['*','*','*','*',5,6,7,8,'*','*','*',12,'*',14,'*',16,17,18,19,20,'*','*',23,'*',25,'*','*',28,29]
#import 850um data
ID = []
RA = []
DEC = []
flux_850 = []
for line in cat_850.readlines()[2:]:
tmp = line.split()
ID.append(tmp[0])
RA.append(co.convHMS(tmp[1]))
DEC.append(co.convDMS(tmp[2]))
flux_850.append(float(tmp[3]))
#import IRAC IDs
ID_IRAC = []
RA_IRAC = []
DEC_IRAC = []
for line in cat_IRAC.readlines()[2:]:
tmp = line.split()
if tmp[1][0:2] == '15':
ID_IRAC.append(tmp[0])
RA_IRAC.append(co.convHMS(tmp[1]))
DEC_IRAC.append(co.convDMS(tmp[2]))
else:
pass
import xlrd
import matplotlib.patches as patches
from PIL import Image
# Define a Gaussian function
def Gaussian(x, A, mu, sigma):
return A*np.exp(-0.5*((x-mu)/sigma)**2)
# Physical properties from Ned Wright cosmo calc.
DL = 101.6*1e6 #pc
scale = 0.471 #kpc/"
# Coma Cluster properties from Zwicky+1942
z_coma = 0.023100
v_coma = 6925.0 #km/s
sigma_coma = 1000.0 #km/s
RA_coma=co.convHMS('12:59:48.7') #deg
DEC_coma=co.convDMS('+27:58:50') #deg
################
# Importing Data
################
loc = ('coma.xlsx')
wb=xlrd.open_workbook(loc)
sheet=wb.sheet_by_index(0)
# Sorting data and converting to useful units
ID = []
RA = []
DEC = []
type = []
import coords as co
mom0_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom0.fits'
mom1_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom1.fits'
mom2_map = '../moment_maps/spt0348_C+_dirty_contsub_briggs_robust05.image.mom2.fits'
############
#MOMENT MAPS
############
#set up figure
fig = pl.figure()
pl.rc('font', size=30)
pl.rc('mathtext', default='regular')
RA_c = co.convHMS('03:48:42.271')
DEC_c = co.convDMS('-62:20:50.85')
#moment-0 map
mom0 = ap.FITSFigure(mom0_map, figure=fig, subplot=[0, 0, 0.75, 1])
mom0.tick_labels.set_xformat('hh:mm:ss')
mom0.tick_labels.set_yformat('dd:mm:ss')
mom0.hide_xaxis_label()
mom0.hide_xtick_labels()
mom0.hide_yaxis_label()
mom0.hide_ytick_labels()
mom0.recenter(RA_c, DEC_c, width=5. / 3600, height=5. / 3600)
mom0.show_colorscale(cmap='jet', vmax=12, vmin=0.0, interpolation='nearest')
mom0.add_colorbar(axis_label_text=r'S$_{[CII]}$ (Jy beam$^{-1}$ km s$^{-1}$)',
axis_label_rotation=-90,
axis_label_pad=50)