SmearChain.cpp File Reference

Allows you to smear contour. For example after performing sets of study one finds out that used sets of uncertainty doesn't fully cover analysis need. Then one can smear additionally contour. More...

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

Include dependency graph for SmearChain.cpp:

Go to the source code of this file.

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

Detailed Description

Allows you to smear contour. For example after performing sets of study one finds out that used sets of uncertainty doesn't fully cover analysis need. Then one can smear additionally contour.

Definition in file SmearChain.cpp.

Function Documentation

main()

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

Definition at line 32 of file SmearChain.cpp.

                                 {
  SetMaCh3LoggerFormat();
  if (argc != 3)
  {
    MACH3LOG_ERROR("How to use: {} MCMC_Output.root config", argv[0]);
    throw MaCh3Exception(__FILE__ , __LINE__ );
  }
  MACH3LOG_INFO("Producing single fit output");
  std::string filename = argv[1];
  std::string config = argv[2];
  SmearChain(filename, config);
 
  return 0;
}

SmearChain()

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

Main function creating MCMCProcessor and calling Smear Chain.

Parameters

inputFile	MCMC Chain
config	Config file with settings

Definition at line 12 of file SmearChain.cpp.

{
  MACH3LOG_INFO("File for study: {}", inputFile);
 
  YAML::Node Settings = M3OpenConfig(config);
 
  // Make the processor
  auto Processor = std::make_unique<MCMCProcessor>(inputFile);
  Processor->SetOutputSuffix("_Smear_MCMC");
  Processor->Initialise();
 
  const auto& Smear = Settings["SmearChain"];
 
  std::vector<std::string> Names = Get<std::vector<std::string>>(Smear["Smear"][0], __FILE__, __LINE__);
  std::vector<double> ErrorValue = Get<std::vector<double>>(Smear["Smear"][1], __FILE__, __LINE__);
 
  bool SaveUnsmearedBranch = GetFromManager<bool>(Smear["SaveUnsmearedBranch"], false);
  Processor->SmearChain(Names, ErrorValue, SaveUnsmearedBranch);
}