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)