PedestalTask.cc

Go to the documentation of this file.

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

using namespace hcaldqm;
PedestalTask::PedestalTask(edm::ParameterSet const& ps):
    DQTask(ps)
{
    //  Containers
    _cPedestalMeans_SubDet.initialize(_name+"/Means/SubDet", 
        "PedestalMeans",
        mapper::fSubDet, 
        new axis::ValueAxis(axis::fXaxis, axis::fADC_15),
        new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
    _cPedestalRMSs_SubDet.initialize(_name+"/PedestalRMSs/SubDet", 
        "PedestalRMSs",
        mapper::fSubDet, 
        new axis::ValueAxis(axis::fXaxis, axis::fADC_5),
        new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
    _cPedestalMeans_depth.initialize(_name+"/PedestalMeans/depth", 
        "PedestalMeans",
        mapper::fdepth, 
        new axis::CoordinateAxis(axis::fXaxis, axis::fieta), 
        new axis::CoordinateAxis(axis::fYaxis, axis::fiphi));
    _cPedestalRMSs_depth.initialize(_name+"/PedestalRMSs/depth", "PedestalRMSs",
        mapper::fdepth, 
        new axis::CoordinateAxis(axis::fXaxis, axis::fieta), 
        new axis::CoordinateAxis(axis::fYaxis, axis::fiphi));

    //  tags
    _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE",
        edm::InputTag("hcalDigis"));
    _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO",
        edm::InputTag("hcalDigis"));
    _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF",
        edm::InputTag("hcalDigis"));
    _tagTrigger = ps.getUntrackedParameter<edm::InputTag>("tagTrigger",
        edm::InputTag("tbunpacker"));
    _tokHBHE = consumes<HBHEDigiCollection>(_tagHBHE);
    _tokHO = consumes<HODigiCollection>(_tagHO);
    _tokHF = consumes<HFDigiCollection>(_tagHF);
    _tokTrigger = consumes<HcalTBTriggerData>(_tagTrigger);
}

/* virtual */ void PedestalTask::bookHistograms(DQMStore::IBooker &ib,
    edm::Run const& r, edm::EventSetup const& es)
{
    DQTask::bookHistograms(ib, r, es);
    _cPedestalMeans_SubDet.book(ib, _subsystem);
    _cPedestalRMSs_SubDet.book(ib, _subsystem);
    _cPedestalMeans_depth.book(ib, _subsystem);
    _cPedestalRMSs_depth.book(ib, _subsystem);
}

/* virtual */ void PedestalTask::_resetMonitors(UpdateFreq uf)
{
    DQTask::_resetMonitors(uf);
}

/* virtual */ void PedestalTask::_dump()
{
    _cPedestalMeans_SubDet.reset();
    _cPedestalRMSs_SubDet.reset();
    _cPedestalMeans_depth.reset();
    _cPedestalRMSs_depth.reset();
    _cPedestals.dump(&_cPedestalMeans_SubDet, true);
    _cPedestals.dump(&_cPedestalRMSs_SubDet, false);
    _cPedestals.dump(&_cPedestalMeans_depth, true);
    _cPedestals.dump(&_cPedestalRMSs_depth, false);
}

/* virtual */ void PedestalTask::_process(edm::Event const& e,
    edm::EventSetup const& es)
{
    edm::Handle<HBHEDigiCollection>     chbhe;
    edm::Handle<HODigiCollection>       cho;
    edm::Handle<HFDigiCollection>       chf;

    if (!e.getByToken(_tokHBHE, chbhe))
        _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
            + _tagHBHE.label() + " " + _tagHBHE.instance());
    if (!e.getByToken(_tokHO, cho))
        _logger.dqmthrow("Collection HODigiCollection isn't available"
            + _tagHO.label() + " " + _tagHO.instance());
    if (!e.getByToken(_tokHF, chf))
        _logger.dqmthrow("Collection HFDigiCollection isn't available"
            + _tagHF.label() + " " + _tagHF.instance());

    for (HBHEDigiCollection::const_iterator it=chbhe->begin();
        it!=chbhe->end(); ++it)
    {
        const HBHEDataFrame digi = (const HBHEDataFrame)(*it);
        int digiSizeToUse = floor(digi.size()/constants::CAPS_NUM)*
            constants::CAPS_NUM-1;
        for (int i=0; i<digiSizeToUse; i++)
            _cPedestals.fill(digi.id(), it->sample(i).adc());
    }
    for (HODigiCollection::const_iterator it=cho->begin();
        it!=cho->end(); ++it)
    {
        const HODataFrame digi = (const HODataFrame)(*it);
        int digiSizeToUse = floor(digi.size()/constants::CAPS_NUM)*
            constants::CAPS_NUM-1;
        for (int i=0; i<digiSizeToUse; i++)
            _cPedestals.fill(digi.id(), it->sample(i).adc());
    }
    for (HFDigiCollection::const_iterator it=chf->begin();
        it!=chf->end(); ++it)
    {
        const HFDataFrame digi = (const HFDataFrame)(*it);
        int digiSizeToUse = floor(digi.size()/constants::CAPS_NUM)*
            constants::CAPS_NUM-1;
        for (int i=0; i<digiSizeToUse; i++)
            _cPedestals.fill(digi.id(), it->sample(i).adc());
    }

    if (_ptype==fOnline && _evsTotal>0 && 
        _evsTotal%constants::CALIBEVENTS_MIN==0)
        this->_dump();
}

/* virtual */ bool PedestalTask::_isApplicable(edm::Event const& e)
{
    if (_ptype==fOnline)
    {
        //  online-global
        return this->_getCalibType(e)==hc_Pedestal;
    }
    else 
    {
        //  local
        edm::Handle<HcalTBTriggerData>  ctrigger;
        if (!e.getByToken(_tokTrigger, ctrigger))
            _logger.dqmthrow("Collection HcalTBTriggerData isn't available"
                + _tagTrigger.label() + " " + _tagTrigger.instance());
        return ctrigger->wasSpillIgnorantPedestalTrigger();
    }

    return false;
}

DEFINE_FWK_MODULE(PedestalTask);