def ionization(self, z, maglim=None): xH = cr.ionization_from_luminosity( z, self._iPhotFunc, rate_is_tfunc=self.extrap_var == 't', ratedensityfunc_args=[maglim], **self.cosmo) return xH
def ionization(self, z, maglim=None): xH = cr.ionization_from_luminosity(z, self._iPhotFunc, rate_is_tfunc = self.extrap_var == 't', ratedensityfunc_args=[maglim], **self.cosmo) return xH
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')