DiagMCMC.cpp

Go to the documentation of this file.

#include "Fitters/MCMCProcessor.h"
#include "Manager/Manager.h"
 
 
void DiagMCMC(const std::string& inputFile, const std::string& config)
{
  MACH3LOG_INFO("File for study: {}", inputFile);
 
  YAML::Node Settings = M3OpenConfig(config);
 
  // Make the processor
  auto Processor = std::make_unique<MCMCProcessor>(inputFile);
  Processor->SetOutputSuffix("_MCMC_Diag");
  //KS:Turn off plotting detector and some other setting
  Processor->SetExcludedTypes(GetFromManager<std::vector<std::string>>(Settings["DiagMCMC"]["ExcludedTypes"], {}));
  Processor->SetExcludedNames(GetFromManager<std::vector<std::string>>(Settings["DiagMCMC"]["ExcludedNames"], {}));
  Processor->SetExcludedGroups(GetFromManager<std::vector<std::string>>(Settings["DiagMCMC"]["ExcludedGroups"], {}));
  Processor->SetPlotRelativeToPrior(GetFromManager<bool>(Settings["DiagMCMC"]["PlotRelativeToPrior"], false));
  //KS: Use 20 batches for batched means
  Processor->SetnBatches(GetFromManager<int>(Settings["DiagMCMC"]["nBatches"], 20));
  Processor->SetnLags(GetFromManager<int>(Settings["DiagMCMC"]["nLags"], 25000));
  Processor->SetPrintToPDF(GetFromManager<bool>(Settings["PrintToPDF"], true));
  Processor->Initialise();
  if(Settings["MaxEntries"]) {
    Processor->SetEntries(Get<int>(Settings["MaxEntries"], __FILE__, __LINE__));
  }
  //KS: finally call main method
  Processor->DiagMCMC();
}
 
int main(int argc, char *argv[]) {
  SetMaCh3LoggerFormat();
  if (argc != 3)
  {
    MACH3LOG_ERROR("How to use: DiagMCMC MCMC_Output.root config");
    throw MaCh3Exception(__FILE__ , __LINE__ );
  }
  MACH3LOG_INFO("Producing single fit output");
  std::string filename = argv[1];
  std::string config = argv[2];
  DiagMCMC(filename, config);
 
  return 0;
}