myNls = NlGenerator(lmap, theory, gradCut=gradCut)
for polComb in ['TT', 'TE', 'EE', 'EB', 'ET', 'TB']:
for beamY in beamRange:
beamell = beamscale(beamY)
for tellmaxY, pellmaxY in [(3000, 5000), (beamell, beamell)]:
for noiseX,beamX,lab,tellminX,tellmaxX,pellminX,pellmaxX in \
[(noiseY,beamY,"sameGrad",tellminY,tellmaxY,pellminY,pellmaxY), \
(30.,7.0,"planckGrad",2,3000,2,3000)]:
noiseTX = noiseX
noisePX = np.sqrt(2.) * noiseTX
noiseTY = noiseY
noisePY = np.sqrt(2.) * noiseTY
kmax = getMax(polComb, tellmaxY, pellmaxY)
i += 1
print((i, tellmaxY, pellmaxY, kmax, "delens:",
delensTolerance))
bin_edges = np.arange(kmin, kmax, dell) + dell
myNls.updateBins(bin_edges)
myNls.updateNoise(beamX,noiseTX,noisePX,tellminX,tellmaxX, \
pellminX,pellmaxX,beamY=beamY,noiseTY=noiseTY, \
noisePY=noisePY,tellminY=tellminY,tellmaxY=tellmaxY, \
pellminY=pellminY,pellmaxY=pellmaxY)
ls, Nls = myNls.getNl(polComb=polComb, halo=halo)
fileName = saveRoot + getFileNameString([
'gradCut', 'polComb', 'beamY', 'noiseY', 'grad',
gradCut = 2000
halo = True
#beamX = 1.0
#beamY = 1.0
noiseTX = 10.0
noisePX = 14.14
noiseTY = 10.0
noisePY = 14.14
tellmin = 2
tellmax = 8000
gradCut = 2000
pellmin = 2
pellmax = 8000
polComb = 'EB'
kmin = 100
kmax = getMax(polComb, tellmax, pellmax)
deg = 10.
px = 0.5
dell = 10
bin_edges = np.arange(kmin, kmax, dell) + dell
theory = loadTheorySpectraFromCAMB(cambRoot,
unlensedEqualsLensed=False,
useTotal=False,
lpad=9000)
lmap = lm.makeEmptyCEATemplate(raSizeDeg=deg,
decSizeDeg=deg,
pixScaleXarcmin=px,
pixScaleYarcmin=px)
myNls = NlGenerator(lmap, theory, bin_edges, gradCut=gradCut)
