getUnfPhysDist.cc File Reference

#include <cstdlib>
#include <cassert>
#include <iostream>
#include <vector>
#include <array>
#include <filesystem>
#include <functional>
#include <algorithm>
#include <type_traits>
#include <TString.h>
#include <TFile.h>
#include <TH1.h>
#include <TKey.h>
#include <TVectorT.h>
#include "Math/VectorUtil.h"
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/CommonTools/interface/toolbox.h"
#include "Core/Unfolding/interface/Observable.h"
#include <darwin.h>

Include dependency graph for getUnfPhysDist.cc:

Classes
struct	Dist
struct	Matrix

Namespaces
	DAS
	DAS::Unfolding
	DAS::Unfolding::ExtractHistogram

Functions
template<typename T >
void	nodeLoop (TDirectory *dOut, TUnfoldBinning *binning, T t, ostream &cout, TString name, int depth)
template<typename T >
void	axisLoop (TDirectory *dOut, TUnfoldBinning *node, T t, ostream &cout, vector< int > ibins, int shift, TString name="", TString title="")
void	getUnfPhysDist (const fs::path &input, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
int	main (int argc, char *argv[])

Function Documentation

main()

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

Todo:

automate this by looking up the metainfo

{
    try {
        DT::StandardInit();
        DT::MetaInfo::versions["TUnfold"] = TUnfold::GetTUnfoldVersion();
 
        fs::path input, output;
 
        auto options = DAS::Options(
                            "Extract distributions with physical binning. The structure "
                            "of the input file is preserved, with one directory per "
                            "systematic variation and further subdirectories. For TH1s: "
                            "`variation/distribution/observable/histograms`. For TH2s:"
                            "`variation/distribution/observable/observable/histograms`.",
                            DT::split | DT::config | DT::syst);
        options.input ("input" , &input , "input ROOT file")
               .output("output", &output, "output ROOT file")
               .args("observables", "unfolding.observables", "list of observables " 
                                                             "(consistent with input of `getUnfHist`)");
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::Unfolding::getUnfPhysDist(input, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}