MakeEfficiencyHistos.h

Go to the documentation of this file.

#include <algorithm>
#include <cstdlib>
#include <functional>
#include <vector>
 
#include "Core/CommonTools/interface/binnings.h"
#include "Core/CommonTools/interface/toolbox.h"
 
#include "ExtractBinning.h"
 
#include <TEfficiency.h>
#include <TH3.h>
#include <TString.h>
 
using namespace std;
 
namespace DAS::Normalisation {
 
unique_ptr<TH1F> makeTriggerEffHist (const unique_ptr<TEfficiency> &t_eff, int iHLT, int iEta)
{
    auto  h3D = t_eff->GetTotalHistogram();  // used only to extract binning
 
    // Extracting HLT bins (thresholds)
    auto h3D_HLT = unique_ptr<TH1>( dynamic_cast<TH1*>(h3D->Clone("h3D_HLT")) );  // create one clone for HLT bins
    vector<double> HLTBins = extractBinning(h3D_HLT, "X");
 
    // Extracting resolution bins
    auto h3D_res = unique_ptr<TH1>( dynamic_cast<TH1*>(h3D->Clone("h3D_res")) );  // a separate clone for Res bins
    vector<double> ResBins = extractBinning(h3D_res, "Z");
    const int nResBins = ResBins.size()-1;
 
    TString name = "efficiency";
    TString title = Form("HLT%.0f_eta%d", HLTBins[iHLT-1], iEta+1);
    unique_ptr<TH1F> h = make_unique<TH1F>(name, title, nResBins, 0, 2*HLTBins[iHLT-1]);
    for (int ipt = 1; ipt <= nResBins; ++ipt){        // loop over pt bins (z axis of TEfficiency)
        h->SetBinContent( ipt, t_eff->GetEfficiency       (t_eff->GetGlobalBin(iHLT, iEta, ipt)) );
        h->SetBinError  ( ipt, t_eff->GetEfficiencyErrorUp(t_eff->GetGlobalBin(iHLT, iEta, ipt)) );
    } // end of pt loop
 
    h->SetDirectory(nullptr);
    return h;
}
 
} // end of DAS::Normalisation namespace