code: def best_fit_likelihood(self):
          
          kwargs_lens, kwargs_source, kwargs_lens_light, kwargs_ps, kwargs_cosmo = self.best_fit(bijective=False)
          param_class = self._param_class
          likelihoodModule = self.likelihoodModule
          logL, _ = likelihoodModule.logL(param_class.kwargs2args(kwargs_lens, kwargs_source, kwargs_lens_light, kwargs_ps,
                                                               kwargs_cosmo))
          return logL
docs: returns the log likelihood of the best fit model of the current state of this class
  
          :return: log likelihood, float