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