def METHOD_NAME(self, z):
    """
    :param z: redshift
    :return: critical density of the universe at redshift z in physical units [h^2 M_sun Mpc^-3]
    """
    return self.rhoc * (self.cosmo.efunc(z)) ** 2
    # return self.rhoc*(1+z)**3