notebook.community

Edit and run

COSMOS is a deep, wide area, multi-wavelength survey aimed at measuring the evolution of galaxies on scales from a few Kpc to 10's of Mpc.

The COSMOS survey is centered at (J2000): RA +150.11916667 (10:00:28.600) DEC +2.20583333 (+02:12:21.00)



In [2]:

    
%pylab inline









    



Populating the interactive namespace from numpy and matplotlib



In [3]:

    
import imaginglss
decals = imaginglss.DECALS('/project/projectdirs/m779/imaginglss/dr2.conf.py')



In [4]:

    
dr = decals.datarelease



In [5]:

    
RA = 150.11916667
DEC = 2.20583333
cosmos = dr.create_footprint((RA -0.7, RA + 0.7, DEC-0.7, DEC+0.7 ))



In [6]:

    
cosmos_DR2 = cosmos.intersect(dr.footprint)



In [7]:

    
print cosmos_DR2









    



Footprint: len(bricks)=42 , area=2.62287 degrees, range=FootPrintRange(ramin=149.5, ramax=151.0, decmin=1.375, decmax=3.125, area=2.6228742650495449)



In [8]:

    
print cosmos









    



Footprint: len(bricks)=42 , area=2.62287 degrees, range=FootPrintRange(ramin=149.5, ramax=151.0, decmin=1.375, decmax=3.125, area=2.6228742650495449)



In [9]:

    
cat = dr.create_catalogue(cosmos_DR2)



In [10]:

    
from imaginglss.analysis import targetselection
from imaginglss.analysis import completeness



In [11]:

    
sigma = {'r':5, 'g':5, 'z': 5}



In [12]:

    
LRG = completeness.LRG(sigma)(targetselection.LRG(cat))









    



/global/project/projectdirs/m779/yfeng1/source/imaginglss/imaginglss/utils/npyquery.py:436: RuntimeWarning: invalid value encountered in power
  r = node.function(*ops)
/global/project/projectdirs/m779/yfeng1/source/imaginglss/imaginglss/utils/npyquery.py:436: RuntimeWarning: invalid value encountered in greater
  r = node.function(*ops)



In [13]:

    
plot(LRG['RA'], LRG['DEC'], ', ')









    Out[13]:





[<matplotlib.lines.Line2D at 0x7f2f19fedf10>]



In [ ]: