fitJetResponse.cc File Reference

#include <cmath>
#include <cassert>
#include <iostream>
#include <map>
#include <filesystem>
#include <TString.h>
#include <TFile.h>
#include <TH1.h>
#include <TH2.h>
#include <TH3.h>
#include <TF1.h>
#include <TKey.h>
#include "Core/CommonTools/interface/DASOptions.h"
#include "Core/JEC/interface/resolution.h"
#include "Core/JEC/interface/DoubleCrystalBall.h"
#include "common.h"
#include "fit.h"
#include <darwin.h>

Include dependency graph for fitJetResponse.cc:

Classes
struct	ResponseFit

Namespaces
	DAS
	DAS::JetEnergy

Functions
vector< double >	GetMergedBinning (TH1 *h, const int threshold=30)
pair< float, float >	findDLogExtrema (const unique_ptr< TH1 > &DLog, const double mu, const pair< float, float > Range)
double	round_to (double value, double precision=1.0)
TH1 *	CumulativeResponse (TH1 *res1D)
void	FitResponseByBin (TDirectory *dir, unique_ptr< TH2 > res2D, const int steering)
void	loopDirsFromGetResponse (TDirectory *dIn, TDirectory *dOut, const int steering, const DT::Slice slice)
void	fitJetResponse (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[]
	)

{
    try {
        DT::StandardInit();
 
        fs::path input, output;
 
        auto options = DAS::Options(
                            "Fit jet response in bins of pt and eta. Any 3D histograms "
                            "with axes following the same convention will be parsed in "
                            "the identical way.",
                            DT::split);
        options.input ("input" , &input , "input ROOT file (from `getJetResponse`)")
               .output("output", &output, "output ROOT file");
 
        const auto& config = options(argc, argv);
        const auto& slice = options.slice();
        const int steering = options.steering();
 
        DAS::JetEnergy::fitJetResponse(input, output, config, steering, slice);
    }
    catch (boost::exception& e) {
        DE::Diagnostic(e);
        return EXIT_FAILURE;
    }
    return EXIT_SUCCESS;
}