applyDataNormalisation.cc File Reference

#include <cstdlib>
#include <cassert>
#include <iostream>
#include <vector>
#include <map>
#include <filesystem>
#include "Core/Objects/interface/Jet.h"
#include "Core/Objects/interface/Event.h"
#include "Core/CommonTools/interface/ControlPlots.h"
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/CommonTools/interface/toolbox.h"
#include <TFile.h>
#include <TChain.h>
#include <TH2.h>
#include "Core/Normalisation/interface/LumiUnc.h"
#include "applyDataNormalisation.h"
#include <darwin.h>

Include dependency graph for applyDataNormalisation.cc:

Namespaces
	DAS
	DAS::Normalisation

Functions
void	applyDataNormalisation (const vector< fs::path > &inputs, 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:

make optional with /dev/null

{
    try {
        DT::StandardInit();
 
        vector<fs::path> inputs;
        fs::path output;
 
        auto options = DAS::Options(
                            "Normalise the data samples to the luminosity. Prescaled events may either "
                            "be normalised to the respective trigger luminosities or using trigger precales.\n"
                            "Existing strategies:\n"
                            " - 'pt': assume that the jets are sorted in pt and try to match the leading jet"
                            " - 'eta': sort the jets in eta and try to match the most fwd and second most fwd",
                            DT::config | DT::split | DT::Friend | DT::syst);
        options.inputs("inputs", &inputs, "input ROOT file(s) or directory")
               .output("output", &output, "output ROOT file")
               .arg<fs::path>("luminosities" , "corrections.normalisation.luminosities" , "2-column file with luminosity per trigger"       )
               .arg<fs::path>("uncertainties", "corrections.normalisation.uncertainties", "luminosity uncertaintiy (including correlations)") 
               .arg<fs::path>("turnons"      , "corrections.normalisation.turnons"      , "2-column file with turn-on per trigger"          )
               .arg<fs::path>("efficiencies" , "corrections.normalisation.efficiencies" , "output ROOT file of `getTriggerTurnons`"         )
               .arg<string>("strategy", "corrections.normalisation.strategy", "'eta' or 'pt'"        )
               .arg<string>("method"  , "corrections.normalisation.method"  , "'prescales' or 'lumi'");
 
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::Normalisation::applyDataNormalisation(inputs, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}