def test_plot_integrate_ionization_recomb_BH(xHe=0.0):
cosmo = {'omega_b_0' : 0.0463,
'omega_M_0' : 0.26,
'omega_lambda_0' : 1.-0.26,
'omega_k_0' : 0.0,
'h' : 0.72,
'n' : 0.95,
'sigma_8' : 0.85,
#'tau' : 0.09,
#'tau_err' : 0.03,
'omega_n_0' : 0.0,
'N_nu': 0,
'Y_He': 0.24,
'baryonic_effects' : False,
}
nIonDot2 = lambda z1: cr.ionization_from_luminosity(z1,
nIonDotBHmodel2,
xHe=xHe,
**cosmo)
zi = numpy.linspace(3.,15.,500.)
w2i = nIonDot2(zi)
ifunc2 = scipy.interpolate.interp1d(zi, w2i)
ifunc2b = scipy.interpolate.interp1d(zi, w2i * 2.0)
clump_fact_func = lambda z1: 2.0
z = numpy.linspace(15,4,100.)
fIon2, w2, t2 = cr.integrate_ion_recomb(z,
ifunc2,
clump_fact_func,
xHe=xHe,
**cosmo)
fIon2b, w2b, t2b = cr.integrate_ion_recomb(z,
ifunc2b,
clump_fact_func,
xHe=xHe,
**cosmo)
pylab.plot(z, fIon2, "--", label='Model 2')
pylab.plot(z, fIon2b, "-", label='Model 2b')
pylab.plot(z, w2, "0.75", label='M2 no rec.')
pylab.xlim(4,15)
pylab.ylim(0,1.0)
pylab.legend(loc='best')
def test_plot_integrate_ionization_recomb_BH(xHe=0.0):
cosmo = {'omega_b_0' : 0.0463,
'omega_M_0' : 0.26,
'omega_lambda_0' : 1.-0.26,
'omega_k_0' : 0.0,
'h' : 0.72,
'n' : 0.95,
'sigma_8' : 0.85,
#'tau' : 0.09,
#'tau_err' : 0.03,
'omega_n_0' : 0.0,
'N_nu': 0,
'Y_He': 0.24,
'baryonic_effects' : False,
}
nIonDot2 = lambda z1: cr.ionization_from_luminosity(z1,
nIonDotBHmodel2,
xHe=xHe,
**cosmo)
zi = numpy.linspace(3.,15.,500.)
w2i = nIonDot2(zi)
ifunc2 = scipy.interpolate.interp1d(zi, w2i)
ifunc2b = scipy.interpolate.interp1d(zi, w2i * 2.0)
clump_fact_func = lambda z1: 2.0
z = numpy.linspace(15,4,100.)
fIon2, w2, t2 = cr.integrate_ion_recomb(z,
ifunc2,
clump_fact_func,
xHe=xHe,
**cosmo)
fIon2b, w2b, t2b = cr.integrate_ion_recomb(z,
ifunc2b,
clump_fact_func,
xHe=xHe,
**cosmo)
pylab.plot(z, fIon2, "--", label='Model 2')
pylab.plot(z, fIon2b, "-", label='Model 2b')
pylab.plot(z, w2, "0.75", label='M2 no rec.')
pylab.xlim(4,15)
pylab.ylim(0,1.0)
pylab.legend(loc='best')
def integrate_ion_recomb(self):
"""Integrate the recombination rate.
Returns and stores results.
"""
self.make_clump_fact_func()
ion_func = self.make_ion_func()
x, w, t = cr.integrate_ion_recomb(self.z,
ion_func=ion_func,
clump_fact_func =
self.clump_fact_func,
**self.cosmo)
self.xr = x
self.w = w
self.t = t
return x, w, t
