tellmax,
pellmin,
pellmax,
beamY=beamY,
noiseTY=noiseTY,
noisePY=noisePY)
ls, Nls = myNls.getNl(polComb=polComb, halo=halo)
#ls,Nls = np.loadtxt("../SZ_filter/data/LA_pol_Nl.txt",unpack=True,delimiter=',')
sn, kappa = NFWMatchedFilterSN(cc,
Mexp,
c,
z,
ells=ls,
Nls=Nls,
kellmax=kellmax,
overdensity=overdensity,
critical=critical,
atClusterZ=atClusterZ,
returnKappa=True)
# io.quickPlot2d(kappa,outDir+"kappa.png")
snM, kappa = NFWMatchedFilterSN(cc,
Mexp,
c,
z,
ells=ls,
Nls=Nls,
kellmax=kellmax,
overdensity=overdensity,
ellkk = np.arange(2, 9000, 1)
Clkk = theory.gCl("kk", ellkk)
pl = Plotter(scaleY='log', scaleX='log')
pl.add(ellkk, 4. * Clkk / 2. / np.pi)
pl.add(ls, 4. * Nls / 2. / np.pi)
pl.legendOn(loc='lower left', labsize=10)
pl.done("output/nl.png")
overdensity = 200.
critical = False
atClusterZ = False
# overdensity=180.
# critical=False
# atClusterZ=False
kellmax = 8000
sn, k, std = NFWMatchedFilterSN(cc,
Mexp,
c,
z,
ells=ls,
Nls=Nls,
kellmax=kellmax,
overdensity=overdensity,
critical=critical,
atClusterZ=atClusterZ)
print sn * np.sqrt(1000.)