def growth(self):
"The instantiated growth model"
if np.issubclass_(self.growth_model, gf.GrowthFactor):
return self.growth_model(self.cosmo, **self.growth_params)
else:
return get_model(self.growth_model, "hmf.growth_factor", cosmo=self.cosmo,
**self.growth_params)
def filter(self):
"""
Instantiated model for filter/window functions.
"""
if issubclass_(self.filter_model, Filter):
return self.filter_model(self.k,self._unnormalised_power, **self.filter_params)
elif isinstance(self.filter_model, basestring):
return get_model(self.filter_model, "hmf.filters", k=self.k,
power=self._unnormalised_power, **self.filter_params)
def growth(self):
"The instantiated growth model"
if np.issubclass_(self.growth_model, gf.GrowthFactor):
return self.growth_model(self.cosmo, **self.growth_params)
else:
return get_model(self.growth_model,
"hmf.growth_factor",
cosmo=self.cosmo,
**self.growth_params)
def transfer(self):
"""
The instantiated transfer model
"""
if np.issubclass_(self.transfer_model, tm.TransferComponent):
return self.transfer_model(self.cosmo, **self.transfer_params)
elif isinstance(self.transfer_model, basestring):
return get_model(self.transfer_model, "hmf.transfer_models", cosmo=self.cosmo,
**self.transfer_params)
def transfer(self):
"""
The instantiated transfer model
"""
if np.issubclass_(self.transfer_model, tm.TransferComponent):
return self.transfer_model(self.cosmo, **self.transfer_params)
elif isinstance(self.transfer_model, basestring):
return get_model(self.transfer_model,
"hmf.transfer_models",
cosmo=self.cosmo,
**self.transfer_params)
def wdm(self):
"""
The instantiated WDM model.
Contains quantities relevant to WDM.
"""
if np.issubclass_(self.wdm_transfer, WDM):
return self.wdm_transfer(self.wdm_mass, self.cosmo,self.z,
**self.wdm_params)
elif isinstance(self.wdm_transfer, basestring):
return get_model(self.wdm_model, __name__, mx=self.wdm_mass, cosmo=self.cosmo,
z=self.z,**self.wdm_params)
def filter(self):
"""
Instantiated model for filter/window functions.
"""
if issubclass_(self.filter_model, Filter):
return self.filter_model(self.k, self._unnormalised_power,
**self.filter_params)
elif isinstance(self.filter_model, basestring):
return get_model(self.filter_model,
"hmf.filters",
k=self.k,
power=self._unnormalised_power,
**self.filter_params)
def hmf(self):
"""
Instantiated model for the hmf fitting function.
"""
if issubclass_(self.hmf_model, ff.FittingFunction):
return self.hmf_model(m=self.m, nu2=self.nu, z=self.z,
delta_halo=self.delta_halo, omegam_z=self.cosmo.Om(self.z),
delta_c=self.delta_c, n_eff=self.n_eff,
** self.hmf_params)
elif isinstance(self.hmf_model, basestring):
return get_model(self.hmf_model, "hmf.fitting_functions",
m=self.m, nu2=self.nu, z=self.z,
delta_halo=self.delta_halo, omegam_z=self.cosmo.Om(self.z),
delta_c=self.delta_c, n_eff=self.n_eff,
** self.hmf_params)
def dndm(self):
"""
The number density of haloes in WDM, ``len=len(m)`` [units :math:`h^4 M_\odot^{-1} Mpc^{-3}`]
"""
dndm = super(MassFunctionWDM, self).dndm
if self.alter_dndm is not None:
if np.issubclass_(self.alter_dndm, WDMRecalibrateMF):
alter = self.alter_dndm(M=self.M, dndm0=dndm, wdm=self.wdm,
**self.alter_params)
else:
alter = get_model(self.alter_dndm, __name__, M=self.M,
dndm0=dndm, wdm=self.wdm, **self.alter_params)
dndm = alter.dndm_alter()
return dndm
def wdm(self):
"""
The instantiated WDM model.
Contains quantities relevant to WDM.
"""
if np.issubclass_(self.wdm_model, WDM):
return self.wdm_model(self.wdm_mass, self.cosmo, self.z,
**self.wdm_params)
elif isinstance(self.wdm_transfer, basestring):
return get_model(self.wdm_model,
__name__,
mx=self.wdm_mass,
cosmo=self.cosmo,
z=self.z,
**self.wdm_params)
def dndm(self):
"""
The number density of haloes in WDM, ``len=len(m)`` [units :math:`h^4 M_\odot^{-1} Mpc^{-3}`]
"""
dndm = super(MassFunctionWDM, self).dndm
if self.alter_dndm is not None:
if np.issubclass_(self.alter_dndm, WDMRecalibrateMF):
alter = self.alter_dndm(M=self.M,
dndm0=dndm,
wdm=self.wdm,
**self.alter_params)
else:
alter = get_model(self.alter_dndm,
__name__,
M=self.M,
dndm0=dndm,
wdm=self.wdm,
**self.alter_params)
dndm = alter.dndm_alter()
return dndm
def hmf(self):
"""
Instantiated model for the hmf fitting function.
"""
if issubclass_(self.hmf_model, ff.FittingFunction):
return self.hmf_model(m=self.m,
nu2=self.nu,
z=self.z,
delta_halo=self.delta_halo,
omegam_z=self.cosmo.Om(self.z),
delta_c=self.delta_c,
n_eff=self.n_eff,
**self.hmf_params)
elif isinstance(self.hmf_model, basestring):
return get_model(self.hmf_model,
"hmf.fitting_functions",
m=self.m,
nu2=self.nu,
z=self.z,
delta_halo=self.delta_halo,
omegam_z=self.cosmo.Om(self.z),
delta_c=self.delta_c,
n_eff=self.n_eff,
**self.hmf_params)
