DiagMCMC.cpp File Reference

#include "Fitters/MCMCProcessor.h"
#include "Manager/Manager.h"

Include dependency graph for DiagMCMC.cpp:

Go to the source code of this file.

Functions
void	DiagMCMC (const std::string &inputFile, const std::string &config)
	Main function creating MCMCProcessor and calling MCMC Diagnostic. More...
int	main (int argc, char *argv[])

Function Documentation

DiagMCMC()

void DiagMCMC	(	const std::string &	inputFile,
		const std::string &	config
	)

Main function creating MCMCProcessor and calling MCMC Diagnostic.

Parameters

inputFile	MCMC Chain
config	Config file with settings

Definition at line 8 of file DiagMCMC.cpp.

{
  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();
}

main()

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 34 of file DiagMCMC.cpp.

                                 {
  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;
}