In [10]:

    

%load_ext autoreload
%autoreload 2
from orphics import maps,cosmology,io,stats
from enlib import enmap
import numpy as np


    

    




The autoreload extension is already loaded. To reload it, use:
  %reload_ext autoreload

In [11]:

    

lc = cosmology.LimberCosmology(lmax=2000,pickling=True,skipPower=False,low_acc=True)


    

    




Generating theory Cls...
initializing power...
Initializing CMB window..

In [12]:

    

lc.addStepNz("g1",0.43,0.51,bias=2)
lc.addStepNz("g2",0.51,0.57,bias=2)
lc.addStepNz("g3",0.57,0.7,bias=2)
lc.addStepNz("gall",0.43,0.7,bias=2)


    

    




('Initializing galaxy window for ', 'g1', ' ...')
('Initializing galaxy window for ', 'g2', ' ...')
('Initializing galaxy window for ', 'g3', ' ...')
('Initializing galaxy window for ', 'gall', ' ...')

In [13]:

    

ellrange = np.arange(0,2000,1)
lc.generateCls(ellrange)
cl11 = lc.getCl("g1","g1")
cl22 = lc.getCl("g2","g2")
cl33 = lc.getCl("g3","g3")
cl00 = lc.getCl("gall","gall")


    

In [15]:

    

pl = io.Plotter()
pl.add(ellrange,cl11)
pl.add(ellrange,cl22)
pl.add(ellrange,cl33)
pl.add(ellrange,cl00,ls="--")
pl.done()