diff --git a/abacusnbody/hod/abacus_hod.py b/abacusnbody/hod/abacus_hod.py index 09910bf..7f21716 100644 --- a/abacusnbody/hod/abacus_hod.py +++ b/abacusnbody/hod/abacus_hod.py @@ -153,8 +153,8 @@ def __init__(self, sim_params, HOD_params, clustering_params = None, chunk=-1, n self.shearbins = np.linspace(-0.5, 0.5, 101) self.halo_mass_func, edges = np.histogramdd( - np.vstack((np.log10(self.halo_data['hmass']), - self.halo_data.get('hdeltac', np.zeros(len(self.halo_data['hmass']))), + np.vstack((np.log10(self.halo_data['hmass']), + self.halo_data.get('hdeltac', np.zeros(len(self.halo_data['hmass']))), self.halo_data.get('hfenv', np.zeros(len(self.halo_data['hmass']))))).T, bins = [self.logMbins, self.deltacbins, self.fenvbins], weights = self.halo_data['hmultis']) @@ -168,9 +168,9 @@ def __init__(self, sim_params, HOD_params, clustering_params = None, chunk=-1, n assert 'hdeltac' in self.halo_data.keys() if self.want_shear: assert 'hshear' in self.halo_data.keys() - + self.halo_mass_func_wshear, edges = np.histogramdd( - np.vstack((np.log10(self.halo_data['hmass']), + np.vstack((np.log10(self.halo_data['hmass']), self.halo_data.get('hdeltac', np.zeros(len(self.halo_data['hmass']))), self.halo_data.get('hfenv', np.zeros(len(self.halo_data['hmass']))), self.halo_data.get('hshear', np.zeros(len(self.halo_data['hmass']))))).T,