HcalOnlineHarvesting.cc

Go to the documentation of this file.

#include "DQM/HcalTasks/interface/HcalOnlineHarvesting.h"

using namespace hcaldqm;
using namespace hcaldqm::constants;

HcalOnlineHarvesting::HcalOnlineHarvesting(edm::ParameterSet const& ps) :
    DQHarvester(ps), _nBad(0), _nTotal(0), _reportSummaryMap(nullptr)
{

    //  NOTE: I will leave Run Summary Generators in place 
    //  just not triggering on endJob!
    _vsumgen.resize(nSummary);
    _vnames.resize(nSummary);
    _vmarks.resize(nSummary);
    for (uint32_t i=0; i<_vmarks.size(); i++)
        _vmarks[i]=false;
    _vnames[fRaw]="RawTask";
    _vnames[fDigi]="DigiTask";
    _vnames[fReco]="RecHitTask";
    _vnames[fTP]="TPTask";
    _vnames[fPedestal]="PedestalTask";

    _vsumgen[fRaw] = new hcaldqm::RawRunSummary("RawRunHarvesting",
        _vnames[fRaw], ps);
    _vsumgen[fDigi] = new hcaldqm::DigiRunSummary("DigiRunHarvesting", 
        _vnames[fDigi],ps);
    _vsumgen[fReco] = new hcaldqm::RecoRunSummary("RecoRunHarvesting",
        _vnames[fReco], ps);
    _vsumgen[fTP] = new hcaldqm::TPRunSummary("TPRunHarvesting",
        _vnames[fTP], ps);
    _vsumgen[fPedestal] = new hcaldqm::PedestalRunSummary("PedestalRunHarvesting",
        _vnames[fPedestal], ps);

    _thresh_bad_bad = ps.getUntrackedParameter("thresh_bad_bad", 0.05);
}

/* virtual */ void HcalOnlineHarvesting::beginRun(
    edm::Run const& r, edm::EventSetup const& es)
{
    DQHarvester::beginRun(r,es);
    for (std::vector<DQClient*>::const_iterator it=_vsumgen.begin();
        it!=_vsumgen.end(); ++it)
        (*it)->beginRun(r,es);
}

/* virtual */ void HcalOnlineHarvesting::_dqmEndLuminosityBlock(
    DQMStore::IBooker& ib,
    DQMStore::IGetter& ig, edm::LuminosityBlock const&, 
    edm::EventSetup const&)
{
    //  DETERMINE WHICH MODULES ARE PRESENT IN DATA
    if (ig.get(_subsystem+"/"+_vnames[fRaw]+"/EventsTotal")!=nullptr)
        _vmarks[fRaw]=true;
    if (ig.get(_subsystem+"/"+_vnames[fDigi]+"/EventsTotal")!=nullptr)
        _vmarks[fDigi]=true;
    if (ig.get(_subsystem+"/"+_vnames[fTP]+"/EventsTotal")!=nullptr)
        _vmarks[fTP]=true;
    if (ig.get(_subsystem+"/"+_vnames[fReco]+"/EventsTotal")!=nullptr)
        _vmarks[fReco]=true;
    if (ig.get(_subsystem+"/"+_vnames[fPedestal]+"/EventsTotal")!=nullptr)
        _vmarks[fPedestal]=true;

    //  CREATE SUMMARY REPORT MAP FED vs LS and LOAD MODULE'S SUMMARIES
    //  NOTE: THIS STATEMENTS WILL BE EXECUTED ONLY ONCE!
    if (!_reportSummaryMap)
    {
        ig.setCurrentFolder(_subsystem+"/EventInfo");
        _reportSummaryMap = ib.book2D("reportSummaryMap", "reportSummaryMap",
            _maxLS, 1, _maxLS+1, _vFEDs.size(), 0, _vFEDs.size());
        for (uint32_t i=0; i<_vFEDs.size(); i++)
        {
            char name[5];
            sprintf(name, "%d", _vFEDs[i]);
            _reportSummaryMap->setBinLabel(i+1, name, 2);
        }
        //  set LS bit to mark Xaxis as LS axis
        _reportSummaryMap->getTH1()->SetBit(BIT(BIT_OFFSET+BIT_AXIS_LS));

        // INITIALIZE ALL THE MODULES
        for (uint32_t i=0; i<_vnames.size(); i++)
            _vcSummaryvsLS.push_back(ContainerSingle2D(_vnames[i],
                "SummaryvsLS",
                new hcaldqm::quantity::LumiSection(_maxLS),
                new hcaldqm::quantity::FEDQuantity(_vFEDs),
                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState)));

        //  LOAD ONLY THOSE MODULES THAT ARE PRESENT IN DATA
        for (uint32_t i=0; i<_vmarks.size(); i++)
        {
            if (_vmarks[i])
                _vcSummaryvsLS[i].load(ig, _subsystem);
        }

        //  Create a map of bad channels and fill
        _cKnownBadChannels_depth.initialize("RunInfo", "KnownBadChannels",
            hcaldqm::hashfunctions::fdepth,
            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),0);
        _cKnownBadChannels_depth.book(ib, _emap, _subsystem);
        for (uintCompactMap::const_iterator it=_xQuality.begin();
            it!=_xQuality.end(); ++it)
            _cKnownBadChannels_depth.fill(HcalDetId(it->first));

        ig.setCurrentFolder(_subsystem+"/EventInfo");
        _runSummary = ib.book2D("runSummary", "runSummary",
            1, 0, 1, 1, 0, 1);
    }

    int ifed=0;
    hcaldqm::flag::Flag fTotal("Status", hcaldqm::flag::fNCDAQ);
    if (_ptype != fOffline) { // hidefed2crate
        for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
            it!=_vhashFEDs.end(); ++it)
        {
            HcalElectronicsId eid(*it);
            hcaldqm::flag::Flag fSum("Status", hcaldqm::flag::fNCDAQ);
            for (uint32_t im=0; im<_vmarks.size(); im++)
                if (_vmarks[im])
                {
                    int x = _vcSummaryvsLS[im].getBinContent(eid, _currentLS);
                    hcaldqm::flag::Flag flag("Status", (hcaldqm::flag::State)x);
                    fSum+=flag;
                }
            _reportSummaryMap->setBinContent(_currentLS, ifed+1, int(fSum._state));
            ifed++;
            fTotal+=fSum;
        }
    }

    // update the Run Summary
    // ^^^TEMPORARY AT THIS POINT!
    if (fTotal._state==hcaldqm::flag::fBAD) _nBad++;
    _nTotal++;
    if (double(_nBad)/double(_nTotal)>=_thresh_bad_bad)
        _runSummary->setBinContent(1, 1, int(hcaldqm::flag::fBAD));
    else if (fTotal._state==hcaldqm::flag::fNCDAQ)
        _runSummary->setBinContent(1,1, int(hcaldqm::flag::fNCDAQ));
    else
        _runSummary->setBinContent(1,1, int(hcaldqm::flag::fGOOD));

    // HF TDC TP efficiency
    if (_vmarks[fTP]) {
        MonitorElement* meOccupancy_HF_depth = ig.get("Hcal/TPTask/OccupancyDataHF_depth/OccupancyDataHF_depth");
        MonitorElement* meOccupancyNoTDC_HF_depth = ig.get("Hcal/TPTask/OccupancyEmulHFNoTDC_depth/OccupancyEmulHFNoTDC_depth");
        MonitorElement* meOccupancy_HF_ieta = ig.get("Hcal/TPTask/OccupancyDataHF_ieta/OccupancyDataHF_ieta");
        MonitorElement* meOccupancyNoTDC_HF_ieta = ig.get("Hcal/TPTask/OccupancyEmulHFNoTDC_ieta/OccupancyEmulHFNoTDC_ieta");

        if (meOccupancy_HF_depth && meOccupancyNoTDC_HF_depth && meOccupancy_HF_ieta && meOccupancyNoTDC_HF_ieta) {
            TH2F* hOccupancy_HF_depth = meOccupancy_HF_depth->getTH2F();
            TH2F* hOccupancyNoTDC_HF_depth = meOccupancyNoTDC_HF_depth->getTH2F();
            TH1F* hOccupancy_HF_ieta = meOccupancy_HF_ieta->getTH1F();
            TH1F* hOccupancyNoTDC_HF_ieta = meOccupancyNoTDC_HF_ieta->getTH1F();

            TH2F *hEfficiency_HF_depth = (TH2F*)hOccupancy_HF_depth->Clone();
            hEfficiency_HF_depth->Divide(hOccupancyNoTDC_HF_depth);
            TH1F *hEfficiency_HF_ieta = (TH1F*)hOccupancy_HF_ieta->Clone();
            hEfficiency_HF_ieta->Divide(hOccupancyNoTDC_HF_ieta);

            ib.setCurrentFolder("Hcal/TPTask");

            MonitorElement* meEfficiency_HF_depth = ib.book2D("TDCCutEfficiency_depth", hEfficiency_HF_depth);
            meEfficiency_HF_depth->setEfficiencyFlag();
            MonitorElement* meEfficiency_HF_ieta = ib.book1D("TDCCutEfficiency_ieta", hEfficiency_HF_ieta);
            meEfficiency_HF_ieta->setEfficiencyFlag();  

            delete hEfficiency_HF_depth;
            delete hEfficiency_HF_ieta;
        }   
    }
}

/*
 *  NO END JOB PROCESSING FOR ONLINE!
 */
/* virtual */ void HcalOnlineHarvesting::_dqmEndJob(DQMStore::IBooker& ib,
    DQMStore::IGetter& ig)
{}

DEFINE_FWK_MODULE(HcalOnlineHarvesting);