def get_data(self):
data=files.read_prior_original()
if self.objtype=='gdev':
w,=where(data['n'] > 2)
else:
w,=where(data['n'] < 2)
self.data=data[w]
def load_data(self):
data=files.read_prior_original()
logic = (data['scale'] < 25)
if self.objtype=='gdev':
w,=where((data['n'] > 3) & logic)
else:
w,=where((data['n'] < 2) & logic)
self.data=data[w]
self.set_sigma()
def get_data(self):
data0 = files.read_prior_original()
eta1, eta2 = lensing.util.g1g2_to_eta1eta2(data0["g"][:, 0], data0["g"][:, 1])
data = eu.numpy_util.add_fields(data0, [("eta", "f8", 2)])
data["eta"][:, 0] = eta1
data["eta"][:, 1] = eta2
self.data = data
def get_data(self):
data0=files.read_prior_original()
eta1,eta2 = lensing.util.g1g2_to_eta1eta2(data0['g'][:,0], data0['g'][:,1])
data=eu.numpy_util.add_fields(data0, [('eta','f8',2)])
data['eta'][:,0] = eta1
data['eta'][:,1] = eta2
self.data=data
def set_bindata(self):
from cluster_step.pipe import CLUSTERSTEP_PSF_STAR,CLUSTERSTEP_GAL
p=files.read_prior_original(type='nosplit')
print p.dtype
wm,=where((p['mag'] > 20) & (p['mag'] < 21))
se1=numpy.zeros(wm.size)
se2=numpy.zeros(wm.size)
for i in xrange(wm.size):
se1[i],se2[i] = lensing.util.g1g2_to_e1e2(p['g'][wm[i],0],p['g'][wm[i],1])
etot2=se1**2 + se2**2
mesq = etot2.mean()
R = 1.0 - 0.5*mesq
print R
data=self.data
wstar,=where( (data['uw_flags']==0)
& ((data['uw_type'] & CLUSTERSTEP_PSF_STAR) != 0) )
#s2n_min=700
#s2n_max=1000
s2n_min=70
s2n_max=80
Tadmom=data['am_irr']+data['am_icc']
wgal,=where( (data['uw_flags']==0)
& ((data['uw_type'] & CLUSTERSTEP_GAL) != 0)
& (Tadmom > 12)
& (data['am_s2n'] > s2n_min)
& (data['am_s2n'] < s2n_max))
print 'ngal:',wgal.size
star_isum = data['uw_isum'][wstar].sum()
star_irr = data['uw_irrsum'][wstar].sum()/star_isum
star_irc = data['uw_ircsum'][wstar].sum()/star_isum
star_icc = data['uw_iccsum'][wstar].sum()/star_isum
gal_isum=data['uw_isum'][wgal].sum()
gal_irr = data['uw_irrsum'][wgal].sum()/gal_isum
gal_irc = data['uw_ircsum'][wgal].sum()/gal_isum
gal_icc = data['uw_iccsum'][wgal].sum()/gal_isum
irr=gal_irr-star_irr
irc=gal_irc-star_irc
icc=gal_icc-star_icc
T = irr+icc
e1=(icc-irr)/T
e2=2*irc/T
print 0.5*e1/R, 0.5*e2/R
stop
