L1TGMTClient.cc

Go to the documentation of this file.

#include "DQM/L1TMonitorClient/interface/L1TGMTClient.h"

#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"

#include <sstream>

L1TGMTClient::L1TGMTClient(const edm::ParameterSet& ps) {
    parameters_ = ps;
    initialize();
}

L1TGMTClient::~L1TGMTClient() {
    LogDebug("TriggerDQM") << "[TriggerDQM]: ending... ";
}

//--------------------------------------------------------
void L1TGMTClient::initialize() {

    // base folder for the contents of this job
    monitorName_ = parameters_.getUntrackedParameter<std::string> (
            "monitorName", "");
    LogDebug("TriggerDQM") << "Monitor name = " << monitorName_ << std::endl;

    output_dir_ = parameters_.getUntrackedParameter<std::string> ("output_dir",
            "");
    LogDebug("TriggerDQM") << "DQM output dir = " << output_dir_ << std::endl;

    input_dir_ = parameters_.getUntrackedParameter<std::string> ("input_dir",
            "");
    LogDebug("TriggerDQM") << "DQM input dir = " << input_dir_ << std::endl;

    m_runInEventLoop = parameters_.getUntrackedParameter<bool> (
            "runInEventLoop", false);
    m_runInEndLumi = parameters_.getUntrackedParameter<bool> ("runInEndLumi",
            false);
    m_runInEndRun = parameters_.getUntrackedParameter<bool> ("runInEndRun",
            false);
    m_runInEndJob = parameters_.getUntrackedParameter<bool> ("runInEndJob",
            false);

}

//--------------------------------------------------------
void L1TGMTClient::dqmEndJob(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {

    // booking histograms in the output_dir_

    ibooker.setCurrentFolder(output_dir_);

    eff_eta_dtcsc = bookClone1DVB(ibooker, igetter, "eff_eta_dtcsc", "efficiency DTCSC vs eta",
            "eta_DTCSC_and_RPC");

    if (eff_eta_dtcsc != nullptr) {
        eff_eta_dtcsc->setAxisTitle("eta", 1);
        if(eff_eta_dtcsc->getTH1F()->GetSumw2N() == 0) eff_eta_dtcsc->getTH1F()->Sumw2();

    }

    eff_eta_rpc = bookClone1DVB(ibooker, igetter, "eff_eta_rpc", "efficiency RPC vs eta",
            "eta_DTCSC_and_RPC");

    if (eff_eta_rpc != nullptr) {
        eff_eta_rpc->setAxisTitle("eta", 1);
        if(eff_eta_rpc->getTH1F()->GetSumw2N() == 0) eff_eta_rpc->getTH1F()->Sumw2();

    }

    eff_phi_dtcsc = bookClone1D(ibooker, igetter, "eff_phi_dtcsc", "efficiency DTCSC vs phi",
            "phi_DTCSC_and_RPC");

    if (eff_phi_dtcsc != nullptr) {
        eff_phi_dtcsc->setAxisTitle("phi (deg)", 1);
        if(eff_phi_dtcsc->getTH1F()->GetSumw2N() == 0) eff_phi_dtcsc->getTH1F()->Sumw2();

    }

    eff_phi_rpc = bookClone1D(ibooker, igetter, "eff_phi_rpc", "efficiency RPC vs phi",
            "phi_DTCSC_and_RPC");

    if (eff_phi_rpc != nullptr) {
        eff_phi_rpc->setAxisTitle("phi (deg)", 1);
        if(eff_phi_rpc->getTH1F()->GetSumw2N() == 0) eff_phi_rpc->getTH1F()->Sumw2();

    }

    eff_etaphi_dtcsc = bookClone2D(ibooker, igetter, "eff_etaphi_dtcsc",
            "efficiency DTCSC vs eta and phi", "etaphi_DTCSC_and_RPC");

    if (eff_etaphi_dtcsc != nullptr) {
        eff_etaphi_dtcsc->setAxisTitle("eta", 1);
        eff_etaphi_dtcsc->setAxisTitle("phi (deg)", 2);
        if(eff_etaphi_dtcsc->getTH2F()->GetSumw2N() == 0) eff_etaphi_dtcsc->getTH2F()->Sumw2();

    }

    eff_etaphi_rpc = bookClone2D(ibooker, igetter, "eff_etaphi_rpc",
            "efficiency RPC vs eta and phi", "etaphi_DTCSC_and_RPC");

    if (eff_etaphi_rpc != nullptr) {
        eff_etaphi_rpc->setAxisTitle("eta", 1);
        eff_etaphi_rpc->setAxisTitle("phi (deg)", 2);
        if(eff_etaphi_rpc->getTH2F()->GetSumw2N() == 0) eff_etaphi_rpc->getTH2F()->Sumw2();

    }

    processHistograms(ibooker, igetter);
}

//--------------------------------------------------------

void L1TGMTClient::dqmEndLuminosityBlock(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const edm::LuminosityBlock& lumiSeg, const edm::EventSetup& evSetup) {
 }

//--------------------------------------------------------
void L1TGMTClient::processHistograms(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter) {

    LogDebug("TriggerDQM") << "L1TGMTClient: processing..." << std::endl;

    makeEfficiency1D(ibooker, igetter, eff_eta_dtcsc, "eta_DTCSC_and_RPC", "eta_RPC_only");
    makeEfficiency1D(ibooker, igetter, eff_eta_rpc, "eta_DTCSC_and_RPC", "eta_DTCSC_only");

    makeEfficiency1D(ibooker, igetter, eff_phi_dtcsc, "phi_DTCSC_and_RPC", "phi_RPC_only");
    makeEfficiency1D(ibooker, igetter, eff_phi_rpc, "phi_DTCSC_and_RPC", "phi_DTCSC_only");

    makeEfficiency2D(ibooker, igetter, eff_etaphi_dtcsc, "etaphi_DTCSC_and_RPC", "etaphi_RPC_only");
    makeEfficiency2D(ibooker, igetter, eff_etaphi_rpc, "etaphi_DTCSC_and_RPC", "etaphi_DTCSC_only");
}

void L1TGMTClient::makeRatio1D(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, MonitorElement* mer, std::string h1Name, std::string h2Name) {

    igetter.setCurrentFolder(output_dir_);

    TH1F* h1 = get1DHisto(input_dir_ + "/" + h1Name, igetter);
    TH1F* h2 = get1DHisto(input_dir_ + "/" + h2Name, igetter);

    if (mer == nullptr) {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::makeRatio1D: monitoring element zero, not able to retrieve histogram"
                << std::endl;
        return;
    }

    TH1F* hr = mer->getTH1F();

    if (hr && h1 && h2) {
        hr->Divide(h1, h2, 1., 1., " ");
    }
}

void L1TGMTClient::makeEfficiency1D(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, MonitorElement* meeff, std::string heName, std::string hiName) {

    igetter.setCurrentFolder(output_dir_);

    TH1F* he = get1DHisto(input_dir_ + "/" + heName, igetter);
    TH1F* hi = get1DHisto(input_dir_ + "/" + hiName, igetter);

    if (meeff == nullptr) {
        LogDebug("TriggerDQM")
                << "L1TGMTClient::makeEfficiency1D: monitoring element zero, not able to retrieve histogram"
                << std::endl;
        return;
    }

    TH1F* heff = meeff->getTH1F();

    if (heff && he && hi) {
        TH1F* hall = (TH1F*) he->Clone("hall");
        hall->Add(hi);
        heff->Divide(he, hall, 1., 1., "B");
        delete hall;
    }
}

void L1TGMTClient::makeEfficiency2D(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, MonitorElement* meeff, std::string heName, std::string hiName) {

    igetter.setCurrentFolder(output_dir_);

    TH2F* he = get2DHisto(input_dir_ + "/" + heName, igetter);
    TH2F* hi = get2DHisto(input_dir_ + "/" + hiName, igetter);

    if (meeff == nullptr) {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::makeEfficiency2D: monitoring element zero, not able to retrieve histogram"
                << std::endl;
        return;
    }

    TH2F* heff = meeff->getTH2F();

    if (heff && he && hi) {
        TH2F* hall = (TH2F*) he->Clone("hall");
        hall->Add(hi);
        heff->Divide(he, hall, 1., 1., "B");
        delete hall;
    }
}

TH1F* L1TGMTClient::get1DHisto(std::string meName, DQMStore::IGetter &igetter) {

    MonitorElement * me_ = igetter.get(meName);

    if (!me_) {
        LogDebug("TriggerDQM") << "\nL1TGMTClient: " << meName << " NOT FOUND.";
        return nullptr;
    }

    return me_->getTH1F();
}

TH2F* L1TGMTClient::get2DHisto(std::string meName, DQMStore::IGetter &igetter) {
    MonitorElement * me_ = igetter.get(meName);

    if (!me_) {
        LogDebug("TriggerDQM") << "\nL1TGMTClient: " << meName << " NOT FOUND.";
        return nullptr;
    }
    return me_->getTH2F();
}

MonitorElement* L1TGMTClient::bookClone1D(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string& name, const std::string& title, const std::string& hrefName) {

    MonitorElement* me;

    TH1F* href = get1DHisto(input_dir_ + "/" + hrefName, igetter);

    if (href) {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone1D: booking histogram "
                << hrefName << std::endl;
        const unsigned nbx = href->GetNbinsX();
        const double xmin = href->GetXaxis()->GetXmin();
        const double xmax = href->GetXaxis()->GetXmax();
    ibooker.setCurrentFolder(output_dir_);
        me = ibooker.book1D(name, title, nbx, xmin, xmax);
    } else {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone1D: not able to clone histogram "
                << hrefName << std::endl;
        me = nullptr;
    }

    return me;
}

MonitorElement* L1TGMTClient::bookClone1DVB(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string& name, const std::string& title, const std::string& hrefName) {

    MonitorElement* me;

    TH1F* href = get1DHisto(input_dir_ + "/" + hrefName, igetter);

    if (href) {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone1DVB: booking histogram "
                << hrefName << std::endl;
        int nbx = href->GetNbinsX();
        if (nbx > 99)
            nbx = 99;
        float xbins[100];
        for (int i = 0; i < nbx; i++) {
            xbins[i] = href->GetBinLowEdge(i + 1);
        }
        xbins[nbx] = href->GetXaxis()->GetXmax();

        ibooker.setCurrentFolder(output_dir_);
        me = ibooker.book1D(name, title, nbx, xbins);

    } else {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone1DVB: not able to clone histogram "
                << hrefName << std::endl;
        me = nullptr;
    }

    return me;
}

MonitorElement* L1TGMTClient::bookClone2D(DQMStore::IBooker &ibooker, DQMStore::IGetter &igetter, const std::string& name,
        const std::string& title, const std::string& hrefName) {

    MonitorElement* me;

    TH2F* href = get2DHisto(input_dir_ + "/" + hrefName, igetter);

    if (href) {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone2D: booking histogram "
                << hrefName << std::endl;
        const unsigned nbx = href->GetNbinsX();
        const double xmin = href->GetXaxis()->GetXmin();
        const double xmax = href->GetXaxis()->GetXmax();
        const unsigned nby = href->GetNbinsY();
        const double ymin = href->GetYaxis()->GetXmin();
        const double ymax = href->GetYaxis()->GetXmax();
    ibooker.setCurrentFolder(output_dir_);
        me = ibooker.book2D(name, title, nbx, xmin, xmax, nby, ymin, ymax);
    } else {
        LogDebug("TriggerDQM")
                << "\nL1TGMTClient::bookClone2D: not able to clone histogram "
                << hrefName << std::endl;
        me = nullptr;
    }

    return me;
}