def bin_edges_from_config(Config, section):
from orphics.tools.stats import npspace
spacing = Config.get(section, "spacing")
minim = Config.getfloat(section, "left_edge")
maxim = Config.getfloat(section, "right_edge")
num = Config.getint(section, "num_bins")
return npspace(minim, maxim, num, scale=spacing)
cc = ClusterCosmology(cosmoDict,
constDict,
clTTFixFile="../SZ_filter/data/cltt_lensed_Feb18.txt")
M500 = 5.e14
c500 = 1.18
z = 0.5
delta = 500
R500 = cc.rdel_c(M500, z, delta)
# R500_alt = cc.rdel_c_alt(M500,c500,delta)
# print R500, R500_alt
# sys.exit()
# R500 = R500_alt
Rrange = npspace(0.1, 10. * R500, 100, "log")
rhofunc = halos.rho_nfw(M500, c500, R500)
rhos = rhofunc(Rrange)
pl = Plotter(scaleY='log',
scaleX='log',
labelX="$R$ (Mpc/h)",
labelY="$\\rho (h^2 M_{\\odot}/{\\mathrm{Mpc}^3})$")
pl.add(Rrange, rhos)
pl.done("output/rhos.png")
comL = cc.results.comoving_radial_distance(z) * cc.h
thetaS = R500 / c500 / comL
arcmin = 0.1
arcmax = 50.