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)