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'