quadrants = ['Q1', 'Q2', 'Q3', 'Q4'] pointings = ['J2218_p3'] ratol = 0.001 dectol = 0.001 ra, dec, z, zf, 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)
def gaussian(x,mean,sigma,amp): return amp*np.exp(-(x-mean)**2/2/(sigma**2)) fields = ['J1005'] z_diff_a = np.array([]) z_diff_b = np.array([]) z_diff_ab = np.array([]) z_diff_q = np.array([]) quadrant = np.array([]) for field in fields: print field filename = '/home/ntejos/catalogs/%s/catalog_%s_total.txt' %(field,field) catalog = readtxt(filename, readnames=True) z_diff_a = np.append(z_diff_a,get_z_diff(catalog,label='a')) z_diff_b = np.append(z_diff_b,get_z_diff(catalog,label='b')) z_diff_ab = np.append(z_diff_ab,get_z_diff(catalog,label='ab')) aux1,aux2 = get_z_diff_quad(catalog,label='ab') z_diff_q = np.append(z_diff_q,aux1) quadrant = np.append(quadrant,aux2) #plots if 1: a_std = float(format(np.std(z_diff_a[np.fabs(z_diff_a)<0.005]),'.5f')) b_std = float(format(np.std(z_diff_b[np.fabs(z_diff_b)<0.005]),'.5f')) ab_std = float(format(np.std(z_diff_ab[np.fabs(z_diff_ab)<0.005]),'.5f'))
import numpy as np quadrants = ['Q1','Q2','Q3','Q4'] pointings = ['J1022_p1','J1022_p2a','J1022_p2b','J1022_p3'] ratol = 0.001 dectol = 0.001 ra,dec,z,zf,mag = [],[],[],[],[] 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)
from subprocess import call from VIMOS_util import \ get_1d_2d_spectra, get_1st_order_region, WMIN, WMAX, MASKATMOS np.seterr(all='ignore') #CMAP = pl.cm.gray CMAP = pl.cm.hot WMIN_PLOT = 5500 WMAX_PLOT = 9500 Ckms = 299792.458 # speed of light km/s, exact LINES = readtxt(astro.datapath + 'linelists/galaxy_lines', names='wa,name,select') def measure_nchi2(twa, tfl, wa, fl, er): masked = np.zeros(len(twa), bool) i0 = 0 for w0, w1 in MASKATMOS: i, j = i0 + twa[i0:].searchsorted([w0, w1]) masked[i:j] = True i0 = j fl1 = np.interp(twa, wa, fl) er1 = np.interp(twa, wa, er) resid = (fl1[~masked] - tfl[~masked]) / er1[~masked] return np.dot(resid, resid) / len(twa)
def gaussian(x, mean, sigma, amp): return amp * np.exp(-(x - mean)**2 / 2 / (sigma**2)) fields = ['J1005'] z_diff_a = np.array([]) z_diff_b = np.array([]) z_diff_ab = np.array([]) z_diff_q = np.array([]) quadrant = np.array([]) for field in fields: print field filename = '/home/ntejos/catalogs/%s/catalog_%s_total.txt' % (field, field) catalog = readtxt(filename, readnames=True) z_diff_a = np.append(z_diff_a, get_z_diff(catalog, label='a')) z_diff_b = np.append(z_diff_b, get_z_diff(catalog, label='b')) z_diff_ab = np.append(z_diff_ab, get_z_diff(catalog, label='ab')) aux1, aux2 = get_z_diff_quad(catalog, label='ab') z_diff_q = np.append(z_diff_q, aux1) quadrant = np.append(quadrant, aux2) #plots if 1: a_std = float(format(np.std(z_diff_a[np.fabs(z_diff_a) < 0.005]), '.5f')) b_std = float(format(np.std(z_diff_b[np.fabs(z_diff_b) < 0.005]), '.5f')) ab_std = float(format(np.std(z_diff_ab[np.fabs(z_diff_ab) < 0.005]), '.5f'))
from subprocess import call from VIMOS_util import \ get_1d_2d_spectra, get_1st_order_region, WMIN, WMAX, MASKATMOS np.seterr(all='ignore') #CMAP = pl.cm.gray CMAP = pl.cm.hot WMIN_PLOT = 5500 WMAX_PLOT = 9500 Ckms = 299792.458 # speed of light km/s, exact LINES = readtxt(astro.datapath + 'linelists/galaxy_lines', names='wa,name,select') def measure_nchi2(twa, tfl, wa, fl, er): masked = np.zeros(len(twa), bool) i0 = 0 for w0,w1 in MASKATMOS: i,j = i0 + twa[i0:].searchsorted([w0, w1]) masked[i:j] = True i0 = j fl1 = np.interp(twa, wa, fl) er1 = np.interp(twa, wa, er) resid = (fl1[~masked] - tfl[~masked]) / er1[~masked] return np.dot(resid, resid) / len(twa)
from astro.io import readtxt,writetxt import numpy as np import pylab as pl from barak.coord import unique_radec,match import pylab as pl import scipy as sc VIMOS = readtxt('/home/ntejos/catalogs/J1005/catalog_J1005_total.txt',readnames=True) VVDS = readtxt('/home/ntejos/catalogs/VVDS/VVDS_F10.cat',readnames=True,comment='#') tol = 3. matches = match(VIMOS.RA_MAPPING,VIMOS.DEC_MAPPING,VVDS.ALPHA_J2000,VVDS.DELTA_J2000,tol) VIMOS_ids = np.where(matches.ind>-1)[0] VVDS_ids = matches.ind[matches.ind>-1] z_diff = VIMOS.ZGAL[VIMOS_ids] - VVDS.Z[VVDS_ids] bins = np.arange(-1,1,0.0005) bins = np.linspace(-0.005,0.005,30) cond = (VIMOS.ZGAL[VIMOS_ids]>0) pl.hist(z_diff[cond],bins,histtype='step',label='std='+str(float(format(np.std(z_diff[(np.fabs(z_diff)<0.005)&(np.fabs(z_diff)>0)]),'.5f')))) pl.legend() pl.xlabel(r'$\Delta z$',fontsize=18) pl.ylabel('#',fontsize=18) pl.title('J1005 VVDS vs VIMOS within '+str(tol)+' arcsec') pl.show() #pl.plot(VIMOS.ZGAL[VIMOS_ids],VVDS.Z[VVDS_ids],'bo')
from astro.io import readtxt, writetxt import numpy as np import pylab as pl from barak.coord import unique_radec, match import pylab as pl import scipy as sc VIMOS = readtxt('/home/ntejos/catalogs/J1005/catalog_J1005_total.txt', readnames=True) VVDS = readtxt('/home/ntejos/catalogs/VVDS/VVDS_F10.cat', readnames=True, comment='#') tol = 3. matches = match(VIMOS.RA_MAPPING, VIMOS.DEC_MAPPING, VVDS.ALPHA_J2000, VVDS.DELTA_J2000, tol) VIMOS_ids = np.where(matches.ind > -1)[0] VVDS_ids = matches.ind[matches.ind > -1] z_diff = VIMOS.ZGAL[VIMOS_ids] - VVDS.Z[VVDS_ids] bins = np.arange(-1, 1, 0.0005) bins = np.linspace(-0.005, 0.005, 30) cond = (VIMOS.ZGAL[VIMOS_ids] > 0) pl.hist(z_diff[cond], bins, histtype='step', label='std=' + str( float( format(
"""Combines all diferent catalogs for different pointings and quadrants in the given field. Some targets can have been observed more than once.""" from astro.io import readtxt from astro.io import writetxt import numpy as np prefix = '/home/ntejos/catalogs/J1005/' catalog = readtxt(prefix+'catalog_J1005_total.txt', readnames=True) #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/Q1005/VIMOS.txt',[ra,dec,z,mag])
"""Combines all diferent catalogs for different pointings and quadrants in the given field. Some targets can have been observed more than once.""" from astro.io import readtxt from astro.io import writetxt import numpy as np prefix = '/home/ntejos/catalogs/J1005/' catalog = readtxt(prefix + 'catalog_J1005_total.txt', readnames=True) #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/Q1005/VIMOS.txt',[ra,dec,z,mag])