Correlation function of DR12Q SDSS Catalog bins

In [80]:

    

import healpix_util as hu
import astropy as ap
import numpy as np
from astropy.io import fits
from astropy.table import Table
import astropy.io.ascii as ascii
from astropy.io import fits
from astropy.constants import c
import matplotlib.pyplot as plt
import math as m
from math import pi
import scipy.special as sp
from astroML.decorators import pickle_results
from scipy import integrate
import warnings
from sklearn.neighbors import BallTree
import pickle
import pymangle
import cython_metric
import lcdmmetric
import pyfits
#from astroML.datasets import fetch_sdss_specgals
#from astroML.correlation import bootstrap_two_point_angular
%matplotlib inline


    

In [81]:

    

fdat=ascii.read("./output/DR12Qbin1.dat")


    

In [82]:

    

fdat


    

    
Out[82]:




<Table length=48309>

z
rar
decr
float64
float64
float64
1.360358
9.3e-05
-0.035487
1.134544
0.000682
0.243272
0.835325
0.000758
-0.142169
0.752729
0.000782
-0.099872
0.896221
0.000918
-0.082249
0.882033
0.001199
0.042834
0.702284
0.001299
0.240805
1.301611
0.001425
0.021841
0.704229
0.001539
0.236405
0.840221
0.001691
0.339623
...
...
...
0.761227
6.281947
0.510654
0.864865
6.282053
-0.055125
0.854921
6.282082
0.454166
0.752217
6.282573
0.223546
0.813033
6.282579
-0.087908
1.172096
6.282584
-0.027623
0.979325
6.282665
0.288813
0.786509
6.282744
0.260439
0.835873
6.28277
0.15974
1.197253
6.283072
0.369546

In [83]:

    

z=fdat['z']
ra=fdat['rar']
dec=fdat['decr']


    

In [84]:

    

Ez = lambda x: 1/m.sqrt(0.3*(1+x)**3+0.7)

np.vectorize(Ez)
#Calculate comoving distance of a data point using the Redshift - This definition is based on the cosmology model we take. Here the distance for E-dS universe is considered. Also note that c/H0 ratio is cancelled in the equations and hence not taken.

def DC_LCDM(z):
  return integrate.quad(Ez, 0, z)[0]
DC_LCDM=np.vectorize(DC_LCDM)


    

In [85]:

    

dr12qbin1 = open("./output/dr12qbin1srarf.dat",'w')
dr12qbin1.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin1.write("%f\t " %z[i])
    dr12qbin1.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin1.write("%f\n " %DC_LCDM(z[i]))
dr12qbin1.close()


    

In [86]:

    

dr12qbin1dat=ascii.read("./output/dr12qbin1srarf.dat")


    

In [87]:

    

NSIDE=512
dr12qbin1pix=hu.HealPix("ring",NSIDE)


    

In [88]:

    

pixdata = open("./output/pixdatadr12qbin1.dat",'w')
pixdata.write("z\t pix \n")

for i in range(0,len(dr12qbin1dat)):
    pixdata.write("%f\t" %dr12qbin1dat['z'][i])
    pixdata.write("%d\n" %dr12qbin1pix.eq2pix(dr12qbin1dat['rar'][i]*180.0/pi,dr12qbin1dat['decr'][i]*180.0/pi))
pixdata.close()


    

In [89]:

    

pixdata = ascii.read("./output/pixdatadr12qbin1.dat")
hpixdata=np.array(np.zeros(hu.nside2npix(NSIDE)))
for j in range(len(pixdata)):
    hpixdata[pixdata[j]['pix']]+=1


    

In [90]:

    

hpixdata


    

    
Out[90]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.])

In [91]:

    

hu.mollview(hpixdata,rot=180)
plt.savefig("./plots/dr12qbin1moll.pdf")


    

In [92]:

    

hu.orthview(hpixdata)
plt.savefig("./plots/dr12qbin1orth.pdf")


    

In [93]:

    

fdat=ascii.read("./output/DR12Qbin2.dat")


    

In [94]:

    

fdat


    

    
Out[94]:




<Table length=42650>

z
rar
decr
float64
float64
float64
1.618846
7.1e-05
0.084296
1.576134
0.000149
0.605257
1.710103
0.000157
0.266244
1.833563
0.000239
0.19133
2.001412
0.000373
0.150584
2.0892
0.000401
-0.096075
1.640517
0.00042
-0.075666
2.004653
0.000488
0.529287
1.545614
0.000513
0.105756
1.615499
0.000626
-0.178954
...
...
...
1.869809
6.282173
0.574029
1.605174
6.282174
0.534793
2.078732
6.282216
0.587197
1.852452
6.282393
0.301937
1.700308
6.282434
0.289916
1.762858
6.282488
0.189745
1.639977
6.282548
0.424658
1.621909
6.282646
-0.073565
1.569409
6.282703
-0.077083
1.849221
6.28282
0.033848

In [95]:

    

z=fdat['z']
ra=fdat['rar']
dec=fdat['decr']


    

In [96]:

    

dr12qbin2 = open("./output/dr12qbin2srarf.dat",'w')
dr12qbin2.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin2.write("%f\t " %z[i])
    dr12qbin2.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin2.write("%f\n " %DC_LCDM(z[i]))
dr12qbin2.close()


    

In [97]:

    

dr12qbin2dat=ascii.read("./output/dr12qbin2srarf.dat")


    

In [98]:

    

NSIDE=512
dr12qbin2pix=hu.HealPix("ring",NSIDE)


    

In [99]:

    

pixdata = open("./output/pixdatadr12qbin2.dat",'w')
pixdata.write("z\t pix \n")

for i in range(0,len(dr12qbin2dat)):
    pixdata.write("%f\t" %dr12qbin2dat['z'][i])
    pixdata.write("%d\n" %dr12qbin2pix.eq2pix(dr12qbin2dat['rar'][i]*180.0/pi,dr12qbin2dat['decr'][i]*180.0/pi))
pixdata.close()


    

In [100]:

    

pixdata = ascii.read("./output/pixdatadr12qbin2.dat")
hpixdata=np.array(np.zeros(hu.nside2npix(NSIDE)))
for j in range(len(pixdata)):
    hpixdata[pixdata[j]['pix']]+=1


    

In [101]:

    

hpixdata


    

    
Out[101]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.])

In [102]:

    

hu.mollview(hpixdata,rot=180)
plt.savefig("./plots/dr12qbin2moll.pdf")


    

In [103]:

    

hu.orthview(hpixdata)
plt.savefig("./plots/dr12qbin2orth.pdf")


    

In [ ]:

    




    

In [104]:

    

fdat=ascii.read("./output/DR12Qbin3.dat")


    

In [105]:

    

fdat


    

    
Out[105]:




<Table length=139641>

z
rar
decr
float64
float64
float64
2.309097
3.3e-05
0.31021
2.497941
4.8e-05
0.261357
2.332655
0.0001
-0.023126
2.373299
0.000119
0.532663
2.542303
0.000128
-0.130654
2.164559
0.000141
-0.004206
2.543471
0.00016
0.204148
2.651057
0.000171
0.385902
2.354148
0.000231
0.021868
2.164761
0.000257
0.533554
...
...
...
2.759595
6.283047
0.565462
2.733773
6.283068
0.087709
2.552077
6.283114
-0.147293
2.398304
6.28312
0.049193
2.181301
6.28314
0.519839
2.439179
6.283143
-0.171497
2.187777
6.283164
0.152742
2.415508
6.283166
0.176571
2.451015
6.28317
0.505355
2.397667
6.283184
0.606452

In [106]:

    

z=fdat['z']
ra=fdat['rar']
dec=fdat['decr']


    

In [107]:

    

dr12qbin3 = open("./output/dr12qbin3srarf.dat",'w')
dr12qbin3.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin3.write("%f\t " %z[i])
    dr12qbin3.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin3.write("%f\n " %DC_LCDM(z[i]))
dr12qbin3.close()


    

In [108]:

    

dr12qbin3dat=ascii.read("./output/dr12qbin3srarf.dat")


    

In [109]:

    

NSIDE=512
dr12qbin3pix=hu.HealPix("ring",NSIDE)


    

In [110]:

    

pixdata = open("./output/pixdatadr12qbin3.dat",'w')
pixdata.write("z\t pix \n")

for i in range(0,len(dr12qbin3dat)):
    pixdata.write("%f\t" %dr12qbin3dat['z'][i])
    pixdata.write("%d\n" %dr12qbin3pix.eq2pix(dr12qbin3dat['rar'][i]*180.0/pi,dr12qbin3dat['decr'][i]*180.0/pi))
pixdata.close()


    

In [111]:

    

pixdata = ascii.read("./output/pixdatadr12qbin3.dat")
hpixdata=np.array(np.zeros(hu.nside2npix(NSIDE)))
for j in range(len(pixdata)):
    hpixdata[pixdata[j]['pix']]+=1


    

In [112]:

    

hpixdata


    

    
Out[112]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.])

In [113]:

    

hu.mollview(hpixdata,rot=180)
plt.savefig("./plots/dr12qbin3moll.pdf")


    

In [114]:

    

hu.orthview(hpixdata)
plt.savefig("./plots/dr12qbin3orth.pdf")


    

In [115]:

    

fdat=ascii.read("./output/DR12Qbin4.dat")


    

In [116]:

    

fdat


    

    
Out[116]:




<Table length=42945>

z
rar
decr
float64
float64
float64
3.088839
0.000103
0.34928
2.906448
0.000252
0.280095
2.879372
0.000361
-0.004858
3.009077
0.000378
0.239244
3.123795
0.000428
0.594248
3.250361
0.000465
-0.120452
2.828129
0.000544
0.204236
3.089968
0.000897
0.398611
3.278693
0.000911
0.556924
3.22843
0.001024
0.026139
...
...
...
2.974564
6.282025
0.069631
3.128922
6.282221
0.212078
3.18261
6.282488
0.390924
2.81591
6.282607
-0.061368
2.976977
6.282751
0.199639
3.421283
6.282824
0.528812
3.128951
6.282996
-0.02931
3.166911
6.283033
0.624115
3.064749
6.283085
0.007521
3.113391
6.283173
0.605993

In [117]:

    

z=fdat['z']
ra=fdat['rar']
dec=fdat['decr']


    

In [118]:

    

dr12qbin4 = open("./output/dr12qbin4srarf.dat",'w')
dr12qbin4.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin4.write("%f\t " %z[i])
    dr12qbin4.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin4.write("%f\n " %DC_LCDM(z[i]))
dr12qbin4.close()


    

In [119]:

    

dr12qbin4dat=ascii.read("./output/dr12qbin4srarf.dat")


    

In [120]:

    

NSIDE=512
dr12qbin4pix=hu.HealPix("ring",NSIDE)


    

In [121]:

    

pixdata = open("./output/pixdatadr12qbin4.dat",'w')
pixdata.write("z\t pix \n")

for i in range(0,len(dr12qbin4dat)):
    pixdata.write("%f\t" %dr12qbin4dat['z'][i])
    pixdata.write("%d\n" %dr12qbin4pix.eq2pix(dr12qbin4dat['rar'][i]*180.0/pi,dr12qbin4dat['decr'][i]*180.0/pi))
pixdata.close()


    

In [122]:

    

pixdata = ascii.read("./output/pixdatadr12qbin4.dat")
hpixdata=np.array(np.zeros(hu.nside2npix(NSIDE)))
for j in range(len(pixdata)):
    hpixdata[pixdata[j]['pix']]+=1


    

In [123]:

    

hpixdata


    

    
Out[123]:





array([ 0.,  0.,  0., ...,  0.,  0.,  0.])

In [124]:

    

hu.mollview(hpixdata,rot=180)
plt.savefig("./plots/dr12qbin4moll.pdf")


    

In [125]:

    

hu.orthview(hpixdata)
plt.savefig("./plots/dr12qbin4orth.pdf")


    

In [126]:

    

def DC_LC(z):
    return np.log(1+z)


    

In [127]:

    

dr12qbin1=ascii.read("./output/dr12qbin1srarf.dat")


    

In [128]:

    

z=dr12qbin1['z']
ra=dr12qbin1['rar']
dec=dr12qbin1['decr']


    

In [129]:

    

dr12qbin1 = open("./output/dr12qbin1srarfLC.dat",'w')
dr12qbin1.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin1.write("%f\t " %z[i])
    dr12qbin1.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin1.write("%f\n " %DC_LC(z[i]))
dr12qbin1.close()


    

In [130]:

    

dr12qbin2=ascii.read("./output/dr12qbin2srarf.dat")


    

In [131]:

    

z=dr12qbin2['z']
ra=dr12qbin2['rar']
dec=dr12qbin2['decr']


    

In [132]:

    

dr12qbin2 = open("./output/dr12qbin2srarfLC.dat",'w')
dr12qbin2.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin2.write("%f\t " %z[i])
    dr12qbin2.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin2.write("%f\n " %DC_LC(z[i]))
dr12qbin2.close()


    

In [133]:

    

dr12qbin3=ascii.read("./output/dr12qbin3srarf.dat")


    

In [134]:

    

z=dr12qbin3['z']
ra=dr12qbin3['rar']
dec=dr12qbin3['decr']


    

In [135]:

    

dr12qbin3 = open("./output/dr12qbin3srarfLC.dat",'w')
dr12qbin3.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin3.write("%f\t " %z[i])
    dr12qbin3.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin3.write("%f\n " %DC_LC(z[i]))
dr12qbin3.close()


    

In [136]:

    

dr12qbin4=ascii.read("./output/dr12qbin4srarf.dat")


    

In [137]:

    

z=dr12qbin4['z']
ra=dr12qbin4['rar']
dec=dr12qbin4['decr']


    

In [138]:

    

dr12qbin4 = open("./output/dr12qbin4srarfLC.dat",'w')
dr12qbin4.write("z\t rar\t decr\t s\n")

for i in range(0,len(z)):
    dr12qbin4.write("%f\t " %z[i])
    dr12qbin4.write("%f\t %f\t " %(ra[i],dec[i]))
    dr12qbin4.write("%f\n " %DC_LC(z[i]))
dr12qbin4.close()


    

In [ ]: