How does CosmoMC find the best-fit values for the cosmological parameters? I understand that this can be done with 'action = 2', but what actually happens inside? I was looking for corresponding code to "action_maxlike=2" in source/GeneralSetup.f90 but I couldn't find anything. I'm also trying to look for error estimates to those best-fit values.