def xray_flux(mass, redshift, theta):
eps_f = theta[0]
eps_DM = theta[1]
f_star = theta[2]
S_star = theta[3]
A_C = theta[4]
alpha0 = theta[5]
n_nt = theta[6]
beta = theta[7]
gamma_mod0 = theta[8]
gamma_mod_zslope = theta[9]
x_break = theta[10]
x_smooth = theta[11]
n_nt_mod = theta[12]
clump0 = theta[13]
alpha_clump = theta[14]
beta_clump = theta[15]
gamma_clump = theta[16]
xx_power.set_Flender_params(alpha0, n_nt, beta, eps_f * 1e-6, eps_DM,
f_star, S_star, A_C, gamma_mod0,
gamma_mod_zslope, x_break, x_smooth, n_nt_mod,
clump0, alpha_clump, beta_clump, gamma_clump)
flux = xx_power.return_flux(redshift, mass)
m500 = xx_power.Mvir_to_Mdeltac(redshift, mass, 500.0)
return flux, m500
def ysz_m(mass, redshift, theta):
eps_f = theta[0]
eps_DM = theta[1]
f_star = theta[2]
S_star = theta[3]
A_C = theta[4]
alpha0 = theta[5]
n_nt = theta[6]
beta = theta[7]
gamma_mod0 = theta[8]
gamma_mod_zslope = theta[9]
x_break = theta[10]
x_smooth = theta[11]
n_nt_mod = theta[12]
clump0 = theta[13]
alpha_clump = theta[14]
beta_clump = theta[15]
gamma_clump = theta[16]
xx_power.set_Flender_params(alpha0, n_nt, beta, eps_f * 1e-6, eps_DM,
f_star, S_star, A_C, gamma_mod0,
gamma_mod_zslope, x_break, x_smooth, n_nt_mod,
clump0, alpha_clump, beta_clump, gamma_clump)
ysz = []
m500 = []
for m in mass:
ysz.append(yy_power.return_Ysz(redshift, m))
m500.append(xx_power.Mvir_to_Mdeltac(redshift, m, 500.0))
ysz = np.array(ysz)
m500 = np.array(m500)
E = np.sqrt(Omega_m * (1 + redshift)**3 + Omega_l)
ysz = ysz * E**(-2. / 3.)
return ysz, m500
def tx_m(mass, redshift, theta):
eps_f = theta[0]
eps_DM = theta[1]
f_star = theta[2]
S_star = theta[3]
A_C = theta[4]
alpha0 = theta[5]
n_nt = theta[6]
beta = theta[7]
gamma_mod0 = theta[8]
gamma_mod_zslope = theta[9]
x_break = theta[10]
x_smooth = theta[11]
n_nt_mod = theta[12]
clump0 = theta[13]
alpha_clump = theta[14]
beta_clump = theta[15]
gamma_clump = theta[16]
xx_power.set_Flender_params(alpha0, n_nt, beta, eps_f * 1e-6, eps_DM,
f_star, S_star, A_C, gamma_mod0,
gamma_mod_zslope, x_break, x_smooth, n_nt_mod,
clump0, alpha_clump, beta_clump, gamma_clump)
tx = []
m500 = []
for m in mass:
tx.append(xx_power.return_Tx(redshift, m))
m500.append(xx_power.Mvir_to_Mdeltac(redshift, m, 500.0))
tx = np.array(tx)
m500 = np.array(m500)
return tx, m500
def mgas_m(mass, redshift, theta):
eps_f = theta[0]
eps_DM = theta[1]
f_star = theta[2]
S_star = theta[3]
A_C = theta[4]
A_nt = theta[5]
B_nt = theta[6]
gamma_nt = theta[7]
gamma_mod0 = theta[8]
gamma_mod_zslope = theta[9]
x_break = theta[10]
x_smooth = theta[11]
n_nt_mod = theta[12]
clump0 = theta[13]
alpha_clump = theta[14]
beta_clump = theta[15]
gamma_clump = theta[16]
xx_power.set_Flender_params(eps_f * 1e-6, eps_DM, f_star, S_star, A_C,
A_nt, B_nt, gamma_nt, gamma_mod0,
gamma_mod_zslope, x_break, x_smooth, n_nt_mod,
clump0, alpha_clump, beta_clump, gamma_clump)
mgas = []
m500 = []
for m in mass:
mgas.append(xx_power.return_Mgas(redshift, m))
m500.append(xx_power.Mvir_to_Mdeltac(redshift, m, 500.0))
mgas = np.array(mgas)
m500 = np.array(m500)
return mgas, m500