unfold.cc File Reference

#include <cstdlib>
#include <cassert>
#include <iostream>
#include <vector>
#include <map>
#include <limits>
#include <TString.h>
#include <TFile.h>
#include <TH1.h>
#include <TLegend.h>
#include <TKey.h>
#include <TMatrixD.h>
#include <TVectorD.h>
#include <TDecompSVD.h>
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/Objects/interface/Format.h"
#include "Core/Unfolding/interface/MyTUnfoldDensity.h"
#include "Core/Unfolding/interface/toolbox.h"
#include "Math/VectorUtil.h"
#include <darwin.h>

Include dependency graph for unfold.cc:

This graph shows which files directly or indirectly include this file:

Namespaces
	DAS
	DAS::Unfolding
	DAS::Unfolding::Migrations
	DAS::Unfolding::Checks
	DAS::Unfolding::Tikhonov

Functions
void	CorrectInefficiencies (unique_ptr< TH1 > &unf, const map< TString, unique_ptr< TH1 >> &miss)
double	Condition (const unique_ptr< TH2 > &RM)
void	NegativeBins (const unique_ptr< TH1 > &unf)
unique_ptr< TH2 >	GetPM (const unique_ptr< TH2 > &RM)
void	folding (const unique_ptr< TH1 > &genIn, const unique_ptr< TH1 > &recIn, const unique_ptr< TH2 > &PMin, const map< TString, unique_ptr< TH1 >> &miss, const map< TString, unique_ptr< TH1 >> &fake)
void	BLT (const unique_ptr< TH1 > &dataDist, const unique_ptr< TH2 > &dataCov, const unique_ptr< TH1 > &MC, int rebin=1)
void	CovMargin (unique_ptr< TH1 > &h, unique_ptr< TH2 > &cov, TString name, int ibin)
void	PrepareRegularisation (TString MCvar, const char *fname, unique_ptr< MyTUnfoldDensity > &density)
void	RedoUnfold (unique_ptr< MyTUnfoldDensity > &density, double newtau, TString name, TString subtitle, unique_ptr< TFile > &fOut, const map< TString, unique_ptr< TH1 >> &miss)
double	ScanLcurve (unique_ptr< MyTUnfoldDensity > &density, TDirectory *dOut, const pt::ptree &config)
void	unfold (const fs::path &inputData, const fs::path &inputMC, const fs::path &output, const pt::ptree &config, const int steering, const DT::Slice slice={1, 0})
int	main (int argc, char *argv[])

Variables
static const auto	deps = numeric_limits<double>::epsilon()

Function Documentation

main()

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

{
    try {
        DT::StandardInit();
        DT::MetaInfo::versions["TUnfold"] = TUnfold::GetTUnfoldVersion();
 
        fs::path inputData, inputMC, output;
 
        auto options = DAS::Options(
                            "Run the unfolding with the TUnfold package. The default approach "
                            "is to run least-square minimisation or pure matrix inversion. "
                            "Unlike all other executables, the options (including those for the "
                            "Tikhonov regularisation, if one wishes to try it) may only be "
                            "provided via a configuration file (note: unlike all other commands, "
                            "the `-e` option does not work here). In case no config file "
                            "is provided, the code will run with the default settings. "
                            "Finally, the parallelisation works slightly differently: instead "
                            "of splitting by ranges of events (since we have no events), "
                            "the running is split by the variations, which may be necessary for "
                            "the unfolding of distributions with a large number of bins (> 100).",
                            DT::config | DT::syst | DT::split);
        options.input("inputData", &inputData, "input Data ROOT file (from `getUnfHist*`)")
               .input("inputMC"  , &inputMC  , "input MC ROOT file (from `getUnfHist*`)")
               .output("output", &output, "output ROOT file with unfolded cross section ");
 
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::Unfolding::unfold(inputData, inputMC, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}

Variable Documentation

deps

const auto deps = numeric_limits<double>::epsilon()

static