def nbar_g(lnm, hmf, HODparams):
integrand = hmf * \
(N_cen(np.exp(lnm), HODparams[:2]) + \
N_sat(np.exp(lnm), HODparams))
return inty(integrand, lnm)
def I_2(self, z):
"""
Integral 2-halo term in the galaxy power spectrum
"""
ilist = []
for i in range(len(self.universe.k)):
ilist.append((1. / self.get_nbar(z)) * \
inty(self._I2inty(z, i), self.lnm))
return np.array(ilist)
def P_1h(self, z1, z2):
"""
1-halo term in the galaxy power spectrum
"""
if z1 != z2:
return 0.
else:
plist = []
for i in range(len(self.universe.k)):
plist.append((1. / self.get_nbar(z1) ** 2) * \
inty(self._1h_int(z1, i), self.lnm))
return np.array(plist)
