DQHarvester.cc

Go to the documentation of this file.

#include "DQM/HcalCommon/interface/DQHarvester.h"
#include "FWCore/Framework/interface/Run.h"

namespace hcaldqm
{
    using namespace constants;

    DQHarvester::DQHarvester(edm::ParameterSet const& ps) :
        DQModule(ps)
    {}

    void DQHarvester::beginRun(edm::Run const& r,
        edm::EventSetup const& es)
    {
        if (_ptype==fLocal)
            if (r.runAuxiliary().run()==1)
                return;

        //      TEMPORARY FIX
        if (_ptype != fOffline) { // hidefed2crate
            _vhashFEDs.clear();
            _vcdaqEids.clear();
        }

        //  - get the Hcal Electronics Map
        //  - collect all the FED numbers and FED's rawIds
        edm::ESHandle<HcalDbService> dbs;
        es.get<HcalDbRecord>().get(dbs);
        _emap = dbs->getHcalMapping();

        if (_ptype != fOffline) { // hidefed2crate
            _vFEDs = utilities::getFEDList(_emap);
            for (std::vector<int>::const_iterator it=_vFEDs.begin();
                it!=_vFEDs.end(); ++it)
            {
                //
                //  FIXME
                //  until there exists a map of FED2Crate and Crate2FED,
                //  all the unknown Crates will be mapped to 0...
                //
                if (*it==0)
                {
                    _vhashFEDs.push_back(HcalElectronicsId(
                        0, SLOT_uTCA_MIN, FIBER_uTCA_MIN1,
                        FIBERCH_MIN, false).rawId());
                    continue;
                }
        
                if (*it>FED_VME_MAX)
                {
                    std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
                    _vhashFEDs.push_back(HcalElectronicsId(
                        cspair.first, cspair.second, FIBER_uTCA_MIN1,
                        FIBERCH_MIN, false).rawId());
                }
                else
                    _vhashFEDs.push_back(HcalElectronicsId(FIBERCH_MIN,
                        FIBER_VME_MIN, SPIGOT_MIN, (*it)-FED_VME_MIN).rawId());
            }

            //  get the FEDs registered at cDAQ
                        if (auto runInfoRec = es.tryToGet<RunInfoRcd>())
            {
                edm::ESHandle<RunInfo> ri;
                                runInfoRec->get(ri);
                std::vector<int> vfeds= ri->m_fed_in;
                for (std::vector<int>::const_iterator it=vfeds.begin();
                    it!=vfeds.end(); ++it)
                {
                    if (*it>=constants::FED_VME_MIN && *it<=FED_VME_MAX)
                        _vcdaqEids.push_back(HcalElectronicsId(
                            constants::FIBERCH_MIN,
                            constants::FIBER_VME_MIN, SPIGOT_MIN,
                            (*it)-FED_VME_MIN).rawId());
                    else if (*it>=constants::FED_uTCA_MIN && 
                        *it<=FEDNumbering::MAXHCALuTCAFEDID)
                    {
                        std::pair<uint16_t, uint16_t> cspair = utilities::fed2crate(*it);
                        _vcdaqEids.push_back(HcalElectronicsId(
                            cspair.first, cspair.second,
                            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
                    }
                }
            }
        } else {
            _vCrates = utilities::getCrateList(_emap);
            std::map<int, uint32_t> crateHashMap = utilities::getCrateHashMap(_emap);
            for (auto& it_crate : _vCrates) {
                _vhashCrates.push_back(crateHashMap[it_crate]);
            }
        }

        // Initialize channel quality masks, but do not load (changed for 10_4_X, moving to LS granularity)
        _xQuality.initialize(hashfunctions::fDChannel);
    }

    void DQHarvester::dqmBeginLuminosityBlock(
        DQMStore::IBooker& ib, DQMStore::IGetter& ig,
        edm::LuminosityBlock const& lb, edm::EventSetup const& es)
    {
        //  get the Hcal Channels Quality for channels that are not 0
        edm::ESHandle<HcalChannelQuality> hcq;
        es.get<HcalChannelQualityRcd>().get("withTopo", hcq);
        const HcalChannelQuality *cq = hcq.product();
        std::vector<DetId> detids = cq->getAllChannels();
        for (std::vector<DetId>::const_iterator it=detids.begin();
            it!=detids.end(); ++it)
        {
             if (HcalGenericDetId(*it).genericSubdet()==
                HcalGenericDetId::HcalGenUnknown)
                 continue;

             if (HcalGenericDetId(*it).isHcalDetId())
             {
                 HcalDetId did(*it);
                  uint32_t mask = (cq->getValues(did))->getValue();
                   if (mask!=0)
                       _xQuality.push(did, mask);
             }
        }
    }

    void DQHarvester::dqmEndLuminosityBlock(
        DQMStore::IBooker& ib, DQMStore::IGetter& ig,
        edm::LuminosityBlock const& lb, edm::EventSetup const& es)
    {
        //  should be the same - just in case!
        _currentLS = lb.luminosityBlock();
        _totalLS++;
        _dqmEndLuminosityBlock(ib, ig, lb, es);
    }
    void DQHarvester::dqmEndJob(DQMStore::IBooker& ib, 
        DQMStore::IGetter& ig)
    {
        _dqmEndJob(ib, ig);
    }
}