#### Mask racenter = 0.0 deccenter = -57.0 maxang = 20.0 center = equ2gal(racenter, deccenter) nsmaskinit = nside veccenter = hp.ang2vec(pi / 2 - np.radians(center[1]), np.radians(center[0])) vecpix = hp.pix2vec(nsmaskinit, np.arange(12 * nsmaskinit ** 2)) cosang = np.dot(veccenter, vecpix) maskok = np.degrees(np.arccos(cosang)) < maxang ### Make Mask Map mapang = XPol.map_ang_from_edges(maskok) maskmap = XPol.apodize_mask(maskok, 2, mapang=mapang) # hp.gnomview(maskmap,rot=[racenter,deccenter],coord=['G','C'],reso=15) wl = hp.anafast(maskmap, regression=False) wl = wl[0 : lmax + 1] maps = hp.synfast(spectra[1:], nside, fwhm=0, pixwin=True, new=True) cls, newl, Mll, MllBinned, MllBinnedInv, p, q, pseudocls = XPol.get_spectra(maps, maskmap, 2 * nside - 1, 20, 20) nbins = len(newl) nbmc = 100 allclsout = np.zeros((nbmc, 6, nbins)) allcls = np.zeros((nbmc, 6, lmax + 1))