beamPX = beamTY = beamPY = beamNow
beamY = beamTY
nTX, nPX, nTY, nPY = myNls.updateNoiseAdvanced(
beamTX, noiseTX, beamPX, noisePX, tellminX, tellmaxX,
pellminX, pellmaxX, beamTY, noiseTY, beamPY, noisePY,
tellminY, tellmaxY, pellminY, pellmaxY, (lkneeTX, lkneePX),
(alphaTX, alphaPX), (lkneeTY, lkneePY), (alphaTY, alphaPY),
None, None, None, None, None, None, None, None, fgFuncX,
fgFuncY, None, None, None, None, None, None, None, None)
ls, Nls, ells, dclbb, efficiency = myNls.getNlIterative(
pols,
kellmin,
kellmax,
tellmax,
pellmin,
pellmax,
dell=dell,
halo=True)
ls = ls[1:-1]
Nls = Nls[1:-1]
clkk_binned = clfunc(ls)
Nls += clkk_binned
if miscenter:
ray = old_div(beamNow, 2.)
else:
decSizeDeg=arc / 60.,
pixScaleXarcmin=px,
pixScaleYarcmin=px)
myNls = NlGenerator(lmap, theory, bin_edges, gradCut=gradCut)
myNls.updateNoise(beam, noiseT, noiseP, tellmin, tellmax, pellmin, pellmax)
polCombList = ['TT', 'EE', 'ET', 'EB', 'TB']
colorList = ['red', 'blue', 'green', 'orange', 'purple']
ellkk = np.arange(2, 9000, 1)
Clkk = theory.gCl("kk", ellkk)
lsmv, Nlmv, ells, dclbb, efficiency = myNls.getNlIterative(polCombList,
kellmin,
kellmax,
tellmax,
pellmin,
pellmax,
dell=10,
halo=True,
plot=True)
print(efficiency)
pl = Plotter(scaleY='log', scaleX='log')
pl.add(ellkk, 4. * Clkk / 2. / np.pi)
# CHECK THAT NORM MATCHES HU/OK
for polComb, col in zip(polCombList, colorList):
ls, Nls = myNls.getNl(polComb=polComb, halo=halo)
# try:
# huFile = 'data/hu_'+polComb.lower()+'.csv'
