def PlotcmbWeights(self,outfile):
#plot weights
pl = Plotter()
for ii in range(len(self.freq)):
pl.add(self.evalells,self.W_ll_cmb[:,ii],label=str(self.freq[ii])+' GHz')
pl.legend(loc='lower left',labsize=10)
pl.done(outfile)
N2 = Nz[np.logical_and(zrange > zleft, zrange <= zright)].sum()
currentAxis.add_patch(
Rectangle((zcent - xerr + pad, 0),
2 * xerr - pad / 2.,
N,
facecolor=col)) #,alpha=0.5))
#currentAxis.add_patch(Rectangle((zcent - xerr+pad+pad/3., 0), 2*xerr-pad/2., N2, facecolor=col))
pl.add([0, 0], [0, 0], ls='-', linewidth=4, label=expName, color=col)
massSense = lndM #*100./np.sqrt(Nmz)
massSense = interpolate_grid(
massSense, masses, zrange, 10**mexp_new, z_new,
regular=True) #,kind="cubic",bounds_error=False,fill_value=np.inf)
print((massSense.shape), testcount)
fsense = massSense / np.sqrt(rn)
pl.legend(labsize=9, loc='upper right')
pl._ax.set_ylim(1, 5.e4) # fsky
pl._ax.set_xlim(0., 3.)
pl.done(outDir + "clNofz.pdf")
fsense[fsense > 10.] = np.nan
from orphics.io import Plotter
import os
mmin = mgrid.min()
mmax = mgrid.max()
zmin = zgrid.min()
zmax = zgrid.max()
pgrid = np.rot90((fsense))
pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$", ylabel="$z$", ftsize=14)
pl.plot2d(pgrid,
extent=[mmin, mmax, zmin, zmax],
fclbbTot = lambda x: fdlbb(x)+fclbb(x)*fgPer/100.
r = rSigma(fsky,ellBBRange,fnBBSmall,fdCls,fclbbTot,fflbb)
cprint("sigma(r) with delensing: "+ str(r*1e4)+"e-4",color="green",bold=True)
rdelens.append(r)
outDir = os.environ['WWW']+"plots/fixtime/"+saveName+"_"
pl = Plotter(xlabel="$f_{\\mathrm{sky}}$",ylabel="delensing %")
pl.add(fskyList,efficiencies)
pl.done(outDir + "efficiencies.png")
# pl = Plotter(xlabel="$f_{\\mathrm{sky}}$",ylabel="sig(mnu)")
# pl.add(fskyList,mnus)
# pl.done(outDir + "mnus.png")
pl = Plotter(xlabel="$f_{\\mathrm{sky}}$",ylabel="Clkk S/N")
pl.add(fskyList,sns)
pl.done(outDir + "sns.png")
pl = Plotter(xlabel="$f_{\\mathrm{sky}}$",ylabel="sig(r)",xscale='log',yscale='log')
pl.add(fskyList,rs,ls="--",label="no delensing")
pl.add(fskyList,rdelens,label="with delensing")
#pl._ax.set_ylim(1.e-4,1.e-2)
pl.legend(loc='upper right',labsize=12)
pl.done(outDir + "rs.png")