def get_weights():
import homogenise_nz
res = load_data()
list_snr_bins_cats = []
for isb,snr_bin in enumerate(config['snr_bins']):
select = (res['SNR'] > snr_bin[0]) & (res['SNR'] < snr_bin[1])
res_bin = res[select]
res_bin_z = (res_bin['photoz'], res_bin['w'])
list_snr_bins_cats.append(res_bin_z)
label = 'snr.%s' % config['method']
list_weights = homogenise_nz.get_weights(list_snr_bins_cats,target_nz_index=0,label=label,photoz_min=0.3,photoz_max=1.3,plots=True)
import cPickle as pickle
filename_pickle = 'weights.%s.pp2' % (config['method'])
pickle.dump(list_weights,open(filename_pickle,'w'))
print 'saved ', filename_pickle
import pdb; pdb.set_trace()
pl.show()
def get_weights():
import homogenise_nz
res = load_data()
list_snr_bins_cats = []
for isb, snr_bin in enumerate(config['snr_bins']):
select = (res['SNR'] > snr_bin[0]) & (res['SNR'] < snr_bin[1])
res_bin = res[select]
res_bin_z = (res_bin['photoz'], res_bin['w'])
list_snr_bins_cats.append(res_bin_z)
label = 'snr.%s' % config['method']
list_weights = homogenise_nz.get_weights(list_snr_bins_cats,
target_nz_index=0,
label=label,
photoz_min=0.3,
photoz_max=1.3,
plots=True)
import cPickle as pickle
filename_pickle = 'weights.%s.pp2' % (config['method'])
pickle.dump(list_weights, open(filename_pickle, 'w'))
print 'saved ', filename_pickle
import pdb
pdb.set_trace()
pl.show()
ra_modest_errorflag_infoflag_spte_zcut_airmass_upper = ra_modest_errorflag_infoflag_spte_zcut[include_upper]
dec_modest_errorflag_infoflag_spte_zcut_airmass_upper = dec_modest_errorflag_infoflag_spte_zcut[include_upper]
e1_modest_errorflag_infoflag_spte_zcut_airmass_upper = e1_modest_errorflag_infoflag_spte_zcut[include_upper]
e2_modest_errorflag_infoflag_spte_zcut_airmass_upper = e2_modest_errorflag_infoflag_spte_zcut[include_upper]
c1_modest_errorflag_infoflag_spte_zcut_airmass_upper = c1_modest_errorflag_infoflag_spte_zcut[include_upper]
c2_modest_errorflag_infoflag_spte_zcut_airmass_upper = c2_modest_errorflag_infoflag_spte_zcut[include_upper]
m_modest_errorflag_infoflag_spte_zcut_airmass_upper = m_modest_errorflag_infoflag_spte_zcut[include_upper]
w_modest_errorflag_infoflag_spte_zcut_airmass_upper = w_modest_errorflag_infoflag_spte_zcut[include_upper]
z_modest_errorflag_infoflag_spte_zcut_airmass_upper = z_modest_errorflag_infoflag_spte_zcut[include_upper]
a = plt.hist(z_modest_errorflag_infoflag_spte_zcut_airmass_lower,200); X = np.zeros([200,2]); binsx = a[1]; bins = (binsx[0:-1]+binsx[1:])/2.0; X[:,0] = bins; X[:,1] = a[0]/np.sum(a[0]); np.savetxt('/Users/drgk/DES/data/systematics/nz_v7_r_airmass_lower_zmin03_zmax13.dat',X,delimiter =' ')
a = plt.hist(z_modest_errorflag_infoflag_spte_zcut_airmass_upper,200); X = np.zeros([200,2]); binsx = a[1]; bins = (binsx[0:-1]+binsx[1:])/2.0; X[:,0] = bins; X[:,1] = a[0]/np.sum(a[0]); np.savetxt('/Users/drgk/DES/data/systematics/nz_v7_r_airmass_upper_zmin03_zmax13.dat',X,delimiter =' ')
list_split = [ (z_modest_errorflag_infoflag_spte_zcut,w_modest_errorflag_infoflag_spte_zcut), (z_modest_errorflag_infoflag_spte_zcut_airmass_upper,w_modest_errorflag_infoflag_spte_zcut_airmass_upper) , (z_modest_errorflag_infoflag_spte_zcut_airmass_lower,w_modest_errorflag_infoflag_spte_zcut_airmass_lower) ]
list_nz_weights = homogenise_nz.get_weights(list_split,target_nz_index=0,photoz_min=0.3,photoz_max=1.3,photoz_nbins=50)
nz_weight_bin0 = list_nz_weights[0] # these should be all one (or very close) as we used this bin as a target n(z)
nz_weight_bin1 = list_nz_weights[1]
nz_weight_bin2 = list_nz_weights[2]
# ## Upper Airmass Catalogue
f = open('/Users/drgk/DES/SV_tests/athena_cats/im3shape_v7_r_shears_airmass_r_upper.dat','w')
g = open('/Users/drgk/DES/SV_tests/athena_cats/im3shape_v7_r_m_airmass_r_upper.dat','w')
for i in range(0,len(ra_modest_errorflag_infoflag_spte_zcut_airmass_upper)):
print str(100*(i+1)/float(len(ra_modest_errorflag_infoflag_spte_zcut_airmass_upper)))+' percent complete: shear Airmass Upper'
ra = ra_modest_errorflag_infoflag_spte_zcut_airmass_upper[i]
dec = dec_modest_errorflag_infoflag_spte_zcut_airmass_upper[i]
e1 = e1_modest_errorflag_infoflag_spte_zcut_airmass_upper[i] - c1_modest_errorflag_infoflag_spte_zcut_airmass_upper[i]
e2 = e2_modest_errorflag_infoflag_spte_zcut_airmass_upper[i] - c2_modest_errorflag_infoflag_spte_zcut_airmass_upper[i]