def METHOD_NAME(self, minimize_kwargs=dict(tol=0.01, method=None)):
    """
    Minimize the death rate and select the best performing model.
    Parameters
    ----------
    tol: tolerance to pass to the optimizer.
    Returns
    -------
    minimization_results: dict
        Results dict from scipy.optimize.minimize.
    """
    self.minimization_results = minimize(
        self._loss_function, x0=self.x0, bounds=self.optimization_bounds, **minimize_kwargs
    )
    self.best_model = self.seir_model_class(
        **self.seir_model_args,
        suppression_policy=self.parametric_policy(
            self.minimization_results["x"],
            t_list=self.seir_model_args["t_list"],
            **self.parametric_policy_kwargs
        )
    )
    self.best_model.METHOD_NAME()
    return self.minimization_results