common.h

Go to the documentation of this file.

#pragma once
 
#include <iostream>
#include <optional>
#include <regex>
#include <string>
#include <sstream>
 
#include "Core/CommonTools/interface/binnings.h"
#include "Core/CommonTools/interface/toolbox.h"
#include "Core/Objects/interface/Jet.h"
#include "Core/Objects/interface/Event.h"
 
#include <darwin.h>
 
namespace DAS::JetEnergy {
 
inline static const std::vector<double> pt_JERC_edges = {15, 17, 20, 23, 27, 30, 35, 40, 45, 57, 72, 90, 120, 150, 200, 300, 400, 550, 750, 1000, 1500, 2000, 2500, 3000, 3500, 4000, 4500, 5000, 5500, 6000, 6500, 7000, 7500, 8000, 8500, 9000, 9500};
inline static const int nPtJERCbins = pt_JERC_edges.size()-1;
 
inline static const std::vector<double> abseta_edges = { 0.000, 0.261, 0.522, 0.783, 1.044, 1.305, 1.566, 1.740, 1.930, 2.043, 2.172, 2.322, 2.500, 2.650, 2.853, 2.964, 3.139, 3.489, 3.839, 5.191 };
inline static const int nAbsEtaBins = abseta_edges.size()-1;
inline static const std::vector<TString> absetaBins = MakeTitle(abseta_edges, "|#eta|", false, true, [](double v) { return Form("%.3f", v);});
 
inline static const std::map<int,std::vector<double>> rho_edges = {
    { 2016, {0, 6.69, 12.39, 18.09, 23.79, 29.49, 35.19, 40.9, 999} }, // EOY16
    { 2017, {0, 7.47, 13.49, 19.52, 25.54, 31.57, 37.59, 999} }, // UL17
    { 2018, {0, 7.35, 13.26, 19.17, 25.08, 30.99, 36.9, 999} }, // UL18
 
    { 2024, {0.0, 8.23, 14.49, 20.75, 27.01, 33.27, 39.53, 45.79, 52.05, 999} }, // Winter24 (copy of Summer23)
};
 
inline int GetNRhoBins (int year) { return rho_edges.at(year).size()-1; }
inline std::vector<TString> rhoBins (int year)
{
    return MakeTitle(rho_edges.at(year), "#rho", false, false,
                        [](double v) { return Form("%.2f", v); });
}
 
inline static const std::vector<double> eta_unc_edges = MakeBinning({-5.4, 5.4}, 0.2, BinningOptions::step);
inline static const int nEtaUncBins = eta_unc_edges.size()-1;
 
inline static const std::vector<std::string> JES_variations {"AbsoluteStat", "AbsoluteScale",  "AbsoluteMPFBias", "Fragmentation", "SinglePionECAL", "SinglePionHCAL", "FlavorQCD", "RelativeJEREC1", "RelativeJEREC2", "RelativeJERHF", "RelativePtBB", "RelativePtEC1", "RelativePtEC2", "RelativePtHF", "RelativeBal", "RelativeSample", "RelativeFSR", "RelativeStatFSR", "RelativeStatEC", "RelativeStatHF", "PileUpDataMC", "PileUpPtRef", "PileUpPtBB", "PileUpPtEC1", "PileUpPtEC2", "PileUpPtHF"}; 
 
} // end of DAS::JetEnergy namespace