randomraN=rnd.field('RA')
randomdecN=rnd.field('DEC')
randomzN=rnd.field('z')
random_wfkpN=rnd.field('WEIGHT_FKP')
random_wN=random_wfkpN
##### Homogeneity
# define r bins
rmin=1.
rmax=200.
nbins=50
r,dd,rr,dr=galtools.paircount_data_random(dataraS,datadecS,datazS,randomraS,randomdecS,randomzS,cosmo,rmin,rmax,nbins,log=1,file='dr10South_pairs.txt',nproc=12)
r,dd,rr,dr=galtools.paircount_data_random(dataraN,datadecN,datazN,randomraN,randomdecN,randomzN,cosmo,rmin,rmax,nbins,log=1,file='dr10North_pairs.txt',nproc=12)
#weighted
r,dd,rr,dr=galtools.paircount_data_random(dataraS,datadecS,datazS,randomraS,randomdecS,randomzS,cosmo,rmin,rmax,nbins,log=1,file='dr10South_pairs_weighted.txt',nproc=2,wdata=data_wS,wrandom=random_wS)
r,dd,rr,dr=galtools.paircount_data_random(dataraN,datadecN,datazN,randomraN,randomdecN,randomzN,cosmo,rmin,rmax,nbins,log=1,file='dr10North_pairs_weighted.txt',nproc=14,wdata=data_wN,wrandom=random_wN)
#### z bins
# define r bins
rmin=1.
rmax=200.
rmax=200.
nbinsr=50
bins=np.array(zeros(nbins),dtype='|S20')
for i in arange(nbins):
bins[i]='z_'+str(zedge[i])+'_'+str(zedge[i+1])
# South
for i in arange(np.size(zedge)-1):
print(' ')
print(bins[i])
mini=zedge[i]
maxi=zedge[i+1]
wdata=where((datazS >= mini) & (datazS <= maxi))
wrandom=where((randomzS >= mini) & (randomzS <= maxi))
r,dd,rr,dr=galtools.paircount_data_random(dataraS[wdata],datadecS[wdata],datazS[wdata],randomraS[wrandom],randomdecS[wrandom],randomzS[wrandom],cosmo,rmin,rmax,nbinsr,log=1,file='dr10South_'+bins[i]+'_pairs_weighted.txt',nproc=8,wdata=data_wS,wrandom=random_wS)
# North
for i in arange(np.size(zedge)-1):
print(' ')
print(bins[i])
mini=zedge[i]
maxi=zedge[i+1]
wdata=where((datazN >= mini) & (datazN <= maxi))
wrandom=where((randomzN >= mini) & (randomzN <= maxi))
r,dd,rr,dr=galtools.paircount_data_random(dataraN[wdata],datadecN[wdata],datazN[wdata],randomraN[wrandom],randomdecN[wrandom],randomzN[wrandom],cosmo,rmin,rmax,nbinsr,log=1,file='dr10North_'+bins[i]+'_pairs_weighted.txt',nproc=16,wdata=data_wS,wrandom=random_wS)
randomzN=rnd.field('z')
random_wfkpN=rnd.field('WEIGHT_FKP')
random_wN=random_wfkpN
# define r bins
rmin=0.
rmax=200.
nbins=50
nproc=8
#weighted
r,dd,rr,dr=galtools.paircount_data_random(dataraS,datadecS,datazS,randomraS,randomdecS,randomzS,cosmo,rmin,rmax,nbins,log=None,file='dr10_V6_South_pairs_forxi.txt',nproc=nproc,wdata=data_wS,wrandom=random_wS)
r,dd,rr,dr=galtools.paircount_data_random(dataraN,datadecN,datazN,randomraN,randomdecN,randomzN,cosmo,rmin,rmax,nbins,log=None,file='dr10_V6_North_pairs_forxi.txt',nproc=nproc,wdata=data_wN,wrandom=random_wN)
##### read north and south for xi
rs,dds,rrs,drs,ngs,nrs=galtools.read_pairs('dr10_V6_South_pairs_forxi.txt')
rn,ddn,rrn,drn,ngn,nrn=galtools.read_pairs('dr10_V6_North_pairs_forxi.txt')
ra,dda,rra,dra,nga,nra=galtools.combine_regions('dr10_V6_South_pairs_forxi.txt','dr10_V6_North_pairs_forxi.txt')
r=rn
# define r bins
rmin = 0.
rmax = 200.
nbins = 50
nproc = 8
#weighted
r, dd, rr, dr = galtools.paircount_data_random(
dataraS,
datadecS,
datazS,
randomraS,
randomdecS,
randomzS,
cosmo,
rmin,
rmax,
nbins,
log=None,
file='dr10_V6_South_pairs_forxi.txt',
nproc=nproc,
wdata=data_wS,
wrandom=random_wS)
r, dd, rr, dr = galtools.paircount_data_random(
dataraN,
datadecN,
datazN,
randomraN,
randomdecN,
randomzN,
wok=where(rnd.field('z')<maxz)
rnd=rnd[wok]
rnd=rnd[0:size(dataraS)*3]
randomraS=rnd.field('RA')
randomdecS=rnd.field('DEC')
randomzS=rnd.field('z')
random_wfkpS=rnd.field('WEIGHT_FKP')
random_wS=random_wfkpS
# define r bins
rmin=0.
rmax=200.
nbins=50
r,dd,rr,dr=galtools.paircount_data_random(dataraS,datadecS,datazS,randomraS,randomdecS,randomzS,cosmo,rmin,rmax,nbins,log=None,file='dr10South_pairs_forxi_weighted.txt',nproc=8,wdata=data_wS,wrandom=random_wS)
r,dd,rr,dr=galtools.paircount_data_random(dataraS,datadecS,datazS,randomraS,randomdecS,randomzS,cosmo,rmin,rmax,nbins,log=None,file='dr10South_pairs_forxi.txt',nproc=8)
south=np.loadtxt('dr10South_pairs_forxi.txt')
rs_nw=south[:,0]
dds_nw=south[:,1]
rrs_nw=south[:,2]
drs_nw=south[:,3]
south=np.loadtxt('dr10South_pairs_forxi_weighted.txt')
rs=south[:,0]
dds=south[:,1]
rrs=south[:,2]
