def myBBth(ell, r):
clBB = qc.get_Dl_fromlib(ell,
r,
lib=binned_camblib,
unlensed=True,
specindex=2)[1]
return clBB
def get_cmb(self, coverage):
okpix = coverage > (np.max(coverage) * float(0))
maskpix = np.zeros(12 * self.nside**2)
maskpix[okpix] = 1
Namaster = nam.Namaster(maskpix,
lmin=40,
lmax=2 * self.nside - 1,
delta_ell=30)
ell = np.arange(2 * self.nside - 1)
global_dir = os.getcwd()
binned_camblib = qc.bin_camblib(Namaster,
global_dir + '/camblib.pkl',
self.nside,
verbose=False)
cls = qc.get_Dl_fromlib(ell,
self.r,
lib=binned_camblib,
unlensed=False)[0]
mycls = qc.Dl2Cl_without_monopole(ell, cls)
np.random.seed(self.seed)
maps = hp.synfast(mycls.T, self.nside, verbose=False, new=True)
return maps
def myclth(ell, r):
clth = qc.get_Dl_fromlib(ell, r, lib=binned_camblib,
unlensed=True)[1]
return clth