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])
