RHat.cpp

Go to the documentation of this file.

// MaCh3 includes
#include "Fitters/RHatCalculator.h"
 
 
// *******************
 
// *******************
int main(int argc, char *argv[]) {
// *******************
  SetMaCh3LoggerFormat();
  M3::Utils::MaCh3Welcome();
 
  std::vector<std::string> MCMCFile;
  int Nchains = 0;
  if (argc < 2)
  {
    MACH3LOG_ERROR("Wrong arguments");
    MACH3LOG_ERROR("./RHat NThin MCMCchain_1.root MCMCchain_2.root MCMCchain_3.root ... [how many you like]");
    throw MaCh3Exception(__FILE__ , __LINE__ );
  }
 
  int NThin = atoi(argv[1]);
 
  //KS Gelman suggests to diagnose on more than one chain
  for (int i = 2; i < argc; i++)
  {
    MCMCFile.push_back(std::string(argv[i]));
    MACH3LOG_INFO("Adding file: {}", MCMCFile.back());
    Nchains++;
  }
 
  if(Nchains == 1)
  {
    MACH3LOG_WARN("Gelman is going to be sad :(. He suggested you should use more than one chain (at least 4). Code works fine for one chain, however, estimator might be biased.");
    MACH3LOG_WARN("Multiple chains are more likely to reveal multimodality and poor adaptation or mixing:");
  }
  MACH3LOG_INFO("Diagnosing {} chains", Nchains);
  auto RHatCalc = std::make_unique<RHatCalculator>(false, MCMCFile, NThin);
  RHatCalc->RunDiagnostic();
 
  return 0;
}