mgrid = np.log10(M) zgrid = old_div((z_edges[1:]+z_edges[:-1]),2.) zz = np.arange(0.1,2.01,0.05) MMexp = np.arange(13.5,15.71,0.1) MM = 10**MMexp hscgrid = np.loadtxt("data/HSC_DeltalnM_z0_z2_17_04_04.txt") sngrid = old_div(1.,hscgrid) pgrid = np.rot90(sngrid) pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$",ylabel="$z$",ftsize=14) pl.plot2d(pgrid,extent=[MMexp.min(),MMexp.max(),zz.min(),zz.max()],levels=[3.0,5.0],labsize=14,aspect="auto") pl.done(outDir+"origHSCgrid.png") print(hscgrid.shape) #outmerr = interpolateGrid(hscgrid,MM,zz,M,zgrid,regular=False,kind="cubic",bounds_error=False,fill_value=np.inf) outmerr = interpolateGrid(hscgrid,MM,zz,M,zgrid,regular=False,kind="cubic",bounds_error=False) sngrid = old_div(1.,outmerr) pgrid = np.rot90(sngrid) pl = Plotter(xlabel="$\\mathrm{log}_{10}(M)$",ylabel="$z$",ftsize=14) pl.plot2d(pgrid,extent=[mgrid.min(),mgrid.max(),zgrid.min(),zgrid.max()],levels=[3.0,5.0],labsize=14,aspect="auto") pl.done(outDir+"interpHSCgrid.png")
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], labsize=14, aspect="auto", lim=[0., 10.]) pl.done(outDir + "massSense.png") rn[rn < 1] = np.nan mmin = mgrid.min() mmax = mgrid.max() zmin = zgrid.min() zmax = zgrid.max() pgrid = np.rot90((rn)) pl = Plotter(labelX="$\\mathrm{log}_{10}(M)$", labelY="$z$", ftsize=14) pl.plot2d(pgrid, extent=[mmin, mmax, zmin, zmax], labsize=14, aspect="auto") pl.done(outDir + "rn.png") #rn[rn<1]=np.nan
jointgrid = np.sqrt(old_div(1.,jointgridsqinv)) snjoint = old_div(1.,jointgrid) for ind in mindicesList: pl.add(outzgrid,snjoint[ind,:].ravel(),ls="-.") print(mexpgrid[ind]) pl.legendOn(loc='upper right',labsize=10) pl.done(outDir+"slice"+cmbtype+".pdf") #pl.legendOn(loc='upper right',labsize=8) #pl.done(outDir+"slice.pdf") #from orphics.io import Plotter pgrid = np.rot90(old_div(1.,jointgrid)) pl = Plotter(labelX="$\\mathrm{log}_{10}(M)$",labelY="$z$",ftsize=14) pl.plot2d(pgrid,extent=[mmin,mmax,zmin,zmax],levels=[1.0,3.0,5.0],labsize=14) pl.done(outDir+"joint"+cmbtype+".png") #savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgrid_S4-1.0-0.4_grid-default_CMB_all_joint_v0.5.pkl" #savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgridRayUp_S4-1.0-0.4_grid-default_CMB_all_joint_v0.5.pkl" #savefile = "/astro/astronfs01/workarea/msyriac/data/SZruns/v0.6/lensgridRayDn_S4-1.0-0._grid-default_CMB_all_joint_v0.5.pkl" #pickle.dump((outmedges,outzedges,jointgrid),open(savefile,'wb')) save_file = save_func(bigDataDir,expName,cmbGrid,saveCalName,version) pickle.dump((outmedges,outzedges,jointgrid),open(save_file,'wb'))