MakeResponseHistos.h

Go to the documentation of this file.

#include <algorithm>
#include <cstdlib>
#include <functional>
#include <vector>
 
#include "Core/Objects/interface/Jet.h"
#include "Core/Objects/interface/Event.h"
 
#include "Core/CommonTools/interface/binnings.h"
#include "Core/CommonTools/interface/toolbox.h"
 
#include "Core/JEC/interface/common.h"
 
#include <TString.h>
#include <TH3D.h>
 
namespace DAS::Normalisation {
 
template<typename T = TH3D> std::unique_ptr<T> makeRespHist
    (TString name, std::vector<double> pt_edges = DAS::JetEnergy::pt_JERC_edges,
     TString title = ";p_{T}^{gen};|#eta^{rec}|;#frac{p_{T}^{rec}}{p_{T}^{gen}}")
{
    using namespace std;
    using namespace DAS::JetEnergy;
 
    static int nResBins = 200;
    static auto resBins = MakeBinning<double>({0.,2.}, nResBins);
    int nPtBins = pt_edges.size()-1;
    unique_ptr<T> h = make_unique<T>(name, title,
                                     nPtBins, pt_edges.data(),
                                     nAbsEtaBins, abseta_edges.data(),
                                     nResBins, resBins.data());
 
    h->SetDirectory(nullptr);
    return h;
}
 
} // end of DAS::Normalisation namespace