for q in quadrants:
filename = '/home/ntejos/J2218/reduction/mos/%s/%s/catalog.txt' % (p,
q)
colnames = 'object,ra,dec,z,zf,temp,mag,mag_err'
try:
catalog = readtxt(filename, names=colnames)
except:
continue
print p, q
for i, obj in enumerate(catalog.object):
ra.append(catalog.ra[i])
dec.append(catalog.dec[i])
z.append(catalog.z[i])
zf.append(catalog.zf[i])
mag.append(catalog.mag[i])
ra = np.array(ra)
dec = np.array(dec)
z = np.array(z)
zf = np.array(zf)
mag = np.array(mag)
#clean duplicates
for i in range(len(ra)):
cond = (np.fabs(ra[i] - ra[i + 1:]) <
ratol) & (np.fabs(dec[i] - dec[i + 1:]) < dectol)
if np.sum(cond) > 0:
z[i] = -99.
writetxt('/home/ntejos/catalogs/Q2218/VIMOS.txt', [ra, dec, z, mag])
from barak.pyvpfit import readf26
from astro.io import writetxt
q0107a = readf26('/home/ntejos/catalogs/Q0107/A.f26')
q0107b = readf26('/home/ntejos/catalogs/Q0107/B.f26')
q0107c = readf26('/home/ntejos/catalogs/Q0107/C.f26')
cond = q0107a.lines.name == 'H I '
a = q0107a.lines[cond]
cond = q0107b.lines.name == 'H I '
b = q0107b.lines[cond]
cond = q0107c.lines.name == 'H I '
c = q0107c.lines[cond]
a.name[:] = 'HI'
b.name[:] = 'HI'
c.name[:] = 'HI'
writetxt('/home/ntejos/catalogs/Q0107/A_HI.txt',[a.name,a.z,a.zsig,a.b,a.bsig,a.logN,a.logNsig])
writetxt('/home/ntejos/catalogs/Q0107/B_HI.txt',[b.name,b.z,b.zsig,b.b,b.bsig,b.logN,b.logNsig])
writetxt('/home/ntejos/catalogs/Q0107/C_HI.txt',[c.name,c.z,c.zsig,c.b,c.bsig,c.logN,c.logNsig])
from barak.pyvpfit import readf26
from astro.io import writetxt
q1022 = readf26('/home/ntejos/COS/q1022/FUV/q1022_fuv_all.f26')
cond = q1022.lines.name == 'H I '
a = q1022.lines[cond]
a.name[:] = 'HI'
writetxt('/home/ntejos/catalogs/Q1022/HI.txt',
[a.name, a.z, a.zsig, a.b, a.bsig, a.logN, a.logNsig])
rows = fh.readlines()
fh.close()
if rows[0].startswith('#'):
print name, 'already converted, skipping'
return None
header = []
i = 0
while not rows[i].startswith('*****'):
header.append('# ' + rows[i].strip())
i += 1
i +=1
wa, tr = [], []
for r in rows[i:]:
if not r.strip():
continue
w, t = map(float, r.split())
wa.append(w)
tr.append(t / 100.)
return adict(wa=wa, tr=tr, hd='\n'.join(header) + '\n')
if 1:
names = glob('*.txt')
for name in names:
f = readfilt(name)
if f is None:
continue
writetxt(name, [f.wa, f.tr], header=f.hd)
from barak.pyvpfit import readf26
from astro.io import writetxt
q1022 = readf26('/home/ntejos/COS/q1022/FUV/q1022_fuv_all.f26')
cond = q1022.lines.name == 'H I '
a = q1022.lines[cond]
a.name[:] = 'HI'
writetxt('/home/ntejos/catalogs/Q1022/HI.txt',[a.name,a.z,a.zsig,a.b,a.bsig,a.logN,a.logNsig])
for p in pointings:
for q in quadrants:
filename = '/home/ntejos/J1022/reduction/mos/%s/%s/catalog.txt' %(p,q)
colnames = 'object,ra,dec,z,zf,temp,mag,mag_err'
try:
catalog = readtxt(filename,names=colnames)
except:
continue
print p,q
for i,obj in enumerate(catalog.object):
ra.append(catalog.ra[i])
dec.append(catalog.dec[i])
z.append(catalog.z[i])
zf.append(catalog.zf[i])
mag.append(catalog.mag[i])
ra = np.array(ra)
dec = np.array(dec)
z = np.array(z)
zf = np.array(zf)
mag = np.array(mag)
#clean duplicates
for i in range(len(ra)):
cond = (np.fabs(ra[i]-ra[i+1:]) < ratol) & (np.fabs(dec[i] - dec[i+1:])<dectol)
if np.sum(cond)>0:
z[i]=-99.
writetxt('/home/ntejos/catalogs/Q1022/VIMOS.txt',[ra,dec,z,mag])
for p in pointings:
for q in quadrants:
filename = "/home/ntejos/J2218/reduction/mos/%s/%s/catalog.txt" % (p, q)
colnames = "object,ra,dec,z,zf,temp,mag,mag_err"
try:
catalog = readtxt(filename, names=colnames)
except:
continue
print p, q
for i, obj in enumerate(catalog.object):
ra.append(catalog.ra[i])
dec.append(catalog.dec[i])
z.append(catalog.z[i])
zf.append(catalog.zf[i])
mag.append(catalog.mag[i])
ra = np.array(ra)
dec = np.array(dec)
z = np.array(z)
zf = np.array(zf)
mag = np.array(mag)
# clean duplicates
for i in range(len(ra)):
cond = (np.fabs(ra[i] - ra[i + 1 :]) < ratol) & (np.fabs(dec[i] - dec[i + 1 :]) < dectol)
if np.sum(cond) > 0:
z[i] = -99.0
writetxt("/home/ntejos/catalogs/Q2218/VIMOS.txt", [ra, dec, z, mag])
isc[i] = 1
zg.append(z)
zg = np.array(zg)[isgal]
isd = isd[isgal]
isv = isv[isgal]
isg = isg[isgal]
isc = isc[isgal]
ra = M[isgal].ALPHA_J2000
dec = M[isgal].DELTA_J2000
# index into the master catalogue
ind = M[isgal].NUMBER - 1
mag = M[isgal].MAG_AUTO
template = M[isgal].VIMOS_TEMPLATE
rec = np.rec.fromarrays([ra,dec,zg,isd,isv,isg,isc,mag,template],
names='ra,dec,z,dei,vim,gmos,cfht,mag,template')
return rec
gals = read_gal()
vimos = gals[gals.vim==1]
writetxt('/home/ntejos/catalogs/Q0107/VIMOS.txt',[vimos.ra,vimos.dec,vimos.z,vimos.mag,vimos.template])
deimos = gals[gals.dei==1]
writetxt('/home/ntejos/catalogs/Q0107/DEIMOS.txt',[deimos.ra,deimos.dec,deimos.z,deimos.mag])
cfht = gals[gals.cfht==1]
writetxt('/home/ntejos/catalogs/Q0107/CFHT.txt',[cfht.ra,cfht.dec,cfht.z,cfht.mag])
gmos = gals[gals.gmos==1]
writetxt('/home/ntejos/catalogs/Q0107/GMOS.txt',[gmos.ra,gmos.dec,gmos.z,gmos.mag])
from barak.pyvpfit import readf26
from astro.io import writetxt
q0107a = readf26('/home/ntejos/catalogs/Q0107/A.f26')
q0107b = readf26('/home/ntejos/catalogs/Q0107/B.f26')
q0107c = readf26('/home/ntejos/catalogs/Q0107/C.f26')
cond = q0107a.lines.name == 'H I '
a = q0107a.lines[cond]
cond = q0107b.lines.name == 'H I '
b = q0107b.lines[cond]
cond = q0107c.lines.name == 'H I '
c = q0107c.lines[cond]
a.name[:] = 'HI'
b.name[:] = 'HI'
c.name[:] = 'HI'
writetxt('/home/ntejos/catalogs/Q0107/A_HI.txt',
[a.name, a.z, a.zsig, a.b, a.bsig, a.logN, a.logNsig])
writetxt('/home/ntejos/catalogs/Q0107/B_HI.txt',
[b.name, b.z, b.zsig, b.b, b.bsig, b.logN, b.logNsig])
writetxt('/home/ntejos/catalogs/Q0107/C_HI.txt',
[c.name, c.z, c.zsig, c.b, c.bsig, c.logN, c.logNsig])
isd = isd[isgal]
isv = isv[isgal]
isg = isg[isgal]
isc = isc[isgal]
ra = M[isgal].ALPHA_J2000
dec = M[isgal].DELTA_J2000
# index into the master catalogue
ind = M[isgal].NUMBER - 1
mag = M[isgal].MAG_AUTO
template = M[isgal].VIMOS_TEMPLATE
rec = np.rec.fromarrays([ra, dec, zg, isd, isv, isg, isc, mag, template],
names='ra,dec,z,dei,vim,gmos,cfht,mag,template')
return rec
gals = read_gal()
vimos = gals[gals.vim == 1]
writetxt('/home/ntejos/catalogs/Q0107/VIMOS.txt',
[vimos.ra, vimos.dec, vimos.z, vimos.mag, vimos.template])
deimos = gals[gals.dei == 1]
writetxt('/home/ntejos/catalogs/Q0107/DEIMOS.txt',
[deimos.ra, deimos.dec, deimos.z, deimos.mag])
cfht = gals[gals.cfht == 1]
writetxt('/home/ntejos/catalogs/Q0107/CFHT.txt',
[cfht.ra, cfht.dec, cfht.z, cfht.mag])
gmos = gals[gals.gmos == 1]
writetxt('/home/ntejos/catalogs/Q0107/GMOS.txt',
[gmos.ra, gmos.dec, gmos.z, gmos.mag])
fh.close()
if rows[0].startswith('#'):
print name, 'already converted, skipping'
return None
header = []
i = 0
while not rows[i].startswith('*****'):
header.append('# ' + rows[i].strip())
i += 1
i += 1
wa, tr = [], []
for r in rows[i:]:
if not r.strip():
continue
w, t = map(float, r.split())
wa.append(w)
tr.append(t / 100.)
return adict(wa=wa, tr=tr, hd='\n'.join(header) + '\n')
if 1:
names = glob('*.txt')
for name in names:
f = readfilt(name)
if f is None:
continue
writetxt(name, [f.wa, f.tr], header=f.hd)
