maskok=~mask ip=np.arange(12*nside**2) ipok=ip[~mask] #ellbins=[0,50,100,150,200,250,300,350,3*nside] ellbins = [0,50,100,150,200] nbins=len(ellbins)-1 ellmin=np.array(ellbins[0:nbins]) ellmax=np.array(ellbins[1:nbins+1])-1 ellval=(ellmin+ellmax)/2 binspec=pyquad.binspectrum(spectra,ellmin,ellmax) reload(qml) ll=np.arange(int(np.max(ellbins))+1) bl=np.exp(-0.5*ll**2*(fwhmrad/2.35)**2) ds_dcb=qml.compute_ds_dcb_par(ellbins,nside,ipok,bl,polar=True,temp=False) signoise=0.05 mapi,mapq,mapu=hp.synfast(spectra[1:],nside,fwhm=fwhmrad,pixwin=True,new=True) themapi=mapi.copy() themapi[mask]=0 themapi[~mask]+=np.random.randn(len(themapi[~mask]))*signoise themapi[~mask]-=np.mean(themapi[~mask]) themapq=mapq.copy() themapq[mask]=0 themapq[~mask]+=np.random.randn(len(themapq[~mask]))*signoise themapq[~mask]-=np.mean(themapq[~mask]) themapu=mapu.copy() themapu[mask]=0 themapu[~mask]+=np.random.randn(len(themapu[~mask]))*signoise themapu[~mask]-=np.mean(themapu[~mask])