applyJERsmearing.cc File Reference

#include <cassert>
#include <cstdlib>
#include <algorithm>
#include <functional>
#include <filesystem>
#include <iostream>
#include <limits>
#include <random>
#include <vector>
#include <TChain.h>
#include <TDirectory.h>
#include <TH2.h>
#include <TFile.h>
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/CommonTools/interface/toolbox.h"
#include "Core/CommonTools/interface/ControlPlots.h"
#include "Core/Objects/interface/Event.h"
#include "Core/Objects/interface/Format.h"
#include "Core/Objects/interface/Jet.h"
#include "Core/JEC/interface/JMEmatching.h"
#include <correction.h>
#include <darwin.h>
#include "common.h"

Include dependency graph for applyJERsmearing.cc:

Namespaces
	DAS
	DAS::JetEnergy

Functions
void	applyJERsmearing (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[]
	)

{
    try {
        DT::StandardInit();
 
        vector<fs::path> inputs;
        fs::path output;
 
        auto options = DAS::Options(
            "Apply jet energy smearing correction. In practice, the correction factors "
            "are stored in a dedicated vector in the reconstructed jets. To apply the "
            "pure stochastic (scaling) smearing, use a null (infinite) core size; hybrid "
            "smearing (scaling in core and stochastic in tails) is recommended. Note: JME "
            "recommends 3 sigmas, but Run 2 studies suggest not more than 2 sigmas.",
            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>("tables", "corrections.JER.tables",
                              "JSON file containing the corrections")
               .arg<string>("campaign", "corrections.JER.campaign",
                            "campaign (e.g. `Summer19UL18_RunA`)")
               .arg<float>("core_width", "corrections.JER.core_width",
                           "width of the Gaussian core in units of sigma "
                           "(for the matching algorithm)");
 
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::JetEnergy::applyJERsmearing(inputs, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
 
    return EXIT_SUCCESS;
}