def METHOD_NAME(self, params, T, Xs):
    alpha_params = params["alpha_"]
    alpha_ = safe_exp(np.dot(Xs["alpha_"], alpha_params))
    beta_params = params["beta_"]
    beta_ = np.exp(np.dot(Xs["beta_"], beta_params))
    return np.logaddexp(beta_ * (np.log(np.clip(T, 1e-25, np.inf)) - np.log(alpha_)), 0)