DigiComparisonTask.cc

Go to the documentation of this file.

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

using namespace hcaldqm;
using namespace hcaldqm::constants;

DigiComparisonTask::DigiComparisonTask(edm::ParameterSet const& ps):
    DQTask(ps)
{
    //  tags and tokens
    _tagHBHE1 = ps.getUntrackedParameter<edm::InputTag>("tagHBHE1",
        edm::InputTag("hcalDigis"));
    _tagHBHE2 = ps.getUntrackedParameter<edm::InputTag>("tagHBHE2",
        edm::InputTag("vmeDigis"));
    _tokHBHE1 = consumes<HBHEDigiCollection>(_tagHBHE1);
    _tokHBHE2 = consumes<HBHEDigiCollection>(_tagHBHE2);
}
/* virtual */ void DigiComparisonTask::bookHistograms(DQMStore::IBooker &ib,
    edm::Run const& r, edm::EventSetup const& es)
{
    DQTask::bookHistograms(ib, r, es);

    //  GET WHAT YOU NEED
    edm::ESHandle<HcalDbService> dbs;
    es.get<HcalDbRecord>().get(dbs);
    edm::ESHandle<HcalElectronicsMap> item;
    es.get<HcalElectronicsMapRcd>().get("full", item);
    _emap = item.product();
    std::vector<int> vFEDs = utilities::getFEDList(_emap);
    std::vector<int> vFEDsVME = utilities::getFEDVMEList(_emap);
    std::vector<int> vFEDsuTCA = utilities::getFEDuTCAList(_emap);
    std::vector<uint32_t> vhashVME;
    std::vector<uint32_t> vhashuTCA;
    vhashVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
        constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
    vhashuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
        FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    _filter_VME.initialize(filter::fFilter, hashfunctions::fElectronics,
        vhashVME);
    _filter_uTCA.initialize(filter::fFilter, hashfunctions::fElectronics,
        vhashuTCA);
    
    //  INITIALIZE
    for  (unsigned int i=0; i<10; i++)
    {
        _cADC_Subdet[i].initialize(_name, "ADC",
            hashfunctions::fSubdet,
            new quantity::ValueQuantity(quantity::fADCCorr_128),
            new quantity::ValueQuantity(quantity::fADCCorr_128),
            new quantity::ValueQuantity(quantity::fN, true));
    }
    _cADCall_Subdet.initialize(_name, "ADC",
        hashfunctions::fSubdet,
        new quantity::ValueQuantity(quantity::fADCCorr_128),
        new quantity::ValueQuantity(quantity::fADCCorr_128),
        new quantity::ValueQuantity(quantity::fN, true));
    _cADCMsnuTCA_Subdet.initialize(_name, "ADCMsnuTCA",
        hashfunctions::fSubdet,
        new quantity::ValueQuantity(quantity::fADC_128),
        new quantity::ValueQuantity(quantity::fN, true));
    _cADCMsnVME_Subdet.initialize(_name, "ADCMsnVME",
        hashfunctions::fSubdet,
        new quantity::ValueQuantity(quantity::fADC_128),
        new quantity::ValueQuantity(quantity::fN, true));
    _cMsm_depth.initialize(_name, "Mismatched",
        hashfunctions::fdepth,
        new quantity::DetectorQuantity(quantity::fieta),
        new quantity::DetectorQuantity(quantity::fiphi),
        new quantity::ValueQuantity(quantity::fN));
    _cMsm_FEDVME.initialize(_name, "Mismatched",
        hashfunctions::fFED,
        new quantity::ElectronicsQuantity(quantity::fSpigot),
        new quantity::ElectronicsQuantity(quantity::fFiberVMEFiberCh),
        new quantity::ValueQuantity(quantity::fN));
    _cMsm_FEDuTCA.initialize(_name, "Mismatched",
        hashfunctions::fFED,
        new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
        new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
        new quantity::ValueQuantity(quantity::fN));
    _cMsnVME_depth.initialize(_name, "Missing",
        hashfunctions::fdepth,
        new quantity::DetectorQuantity(quantity::fieta),
        new quantity::DetectorQuantity(quantity::fiphi),
        new quantity::ValueQuantity(quantity::fN));
    _cMsnuTCA_depth.initialize(_name, "Missing",
        hashfunctions::fdepth,
        new quantity::DetectorQuantity(quantity::fieta),
        new quantity::DetectorQuantity(quantity::fiphi),
        new quantity::ValueQuantity(quantity::fN));
    _cMsn_FEDVME.initialize(_name, "Missing",
        hashfunctions::fFED,
        new quantity::ElectronicsQuantity(quantity::fSpigot),
        new quantity::ElectronicsQuantity(quantity::fFiberVMEFiberCh),
        new quantity::ValueQuantity(quantity::fN));
    _cMsn_FEDuTCA.initialize(_name, "Missing",
        hashfunctions::fFED,
        new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
        new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
        new quantity::ValueQuantity(quantity::fN));

    //  BOOK
    char aux[20];
    for (unsigned int i=0; i<10; i++)
    {
        sprintf(aux, "TS%d", i);
        _cADC_Subdet[i].book(ib, _emap, _subsystem, aux);
    }
    _cADCall_Subdet.book(ib, _emap, _subsystem);
    _cADCMsnVME_Subdet.book(ib, _emap, _subsystem);
    _cADCMsnuTCA_Subdet.book(ib, _emap, _subsystem);
    _cMsm_depth.book(ib, _emap, _subsystem);
    _cMsnVME_depth.book(ib, _emap, _subsystem, std::string("VME"));
    _cMsnuTCA_depth.book(ib, _emap, _subsystem, std::string("uTCA"));
    _cMsm_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
    _cMsn_FEDVME.book(ib, _emap, _filter_uTCA, _subsystem);
    _cMsm_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cMsn_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);

    _ehashmapuTCA.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap,
        _filter_VME);
    _ehashmapVME.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap,
        _filter_uTCA);
}

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

/* virtual */ void DigiComparisonTask::_process(edm::Event const& e,
    edm::EventSetup const& es)
{
    edm::Handle<HBHEDigiCollection> chbhe1;
    edm::Handle<HBHEDigiCollection> chbhe2;
    
    if (!e.getByToken(_tokHBHE1, chbhe1))
        _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
            + _tagHBHE1.label() + " " + _tagHBHE1.instance());
    if (!e.getByToken(_tokHBHE2, chbhe2))
        _logger.dqmthrow("Collection HBHEDigiCollection isn't available"
            + _tagHBHE2.label() + " " + _tagHBHE2.instance());

    //  assume that coll1 is primary(uTCA) and coll2 is secondary(VME)
    //  uTCA is X and VME is Y axis
    for (HBHEDigiCollection::const_iterator it1=chbhe1->begin();
        it1!=chbhe1->end(); ++it1)
    {
        //  iterate thru the utca collection
        //  get the same detid digi from vme collection
        //  if missing - fill vme missing
        //  else correlate
        HcalDetId did = it1->id();
        HcalElectronicsId eid1 = it1->elecId();
        HBHEDigiCollection::const_iterator it2 = chbhe2->find(did);

        //  get the eid for vme by did
        HcalElectronicsId eid2 = HcalElectronicsId(_ehashmapVME.lookup(did));
        if (it2==chbhe2->end())
        {
            //  fill the depth plot
            _cMsnVME_depth.fill(did);
            _cMsn_FEDVME.fill(eid2);
            for (int i=0; i<it1->size(); i++)
            {
                _cADCMsnVME_Subdet.fill(did, it1->sample(i).adc());
                _cADCall_Subdet.fill(did, it1->sample(i).adc(), -2);
                _cADC_Subdet[i].fill(did, it1->sample(i).adc(), -2);
            }
        }
        else
            for (int i=0; i<it1->size(); i++)
            {
                _cADCall_Subdet.fill(did, double(it1->sample(i).adc()),
                    double(it2->sample(i).adc()));
                _cADC_Subdet[i].fill(did, double(it1->sample(i).adc()),
                    double(it2->sample(i).adc()));
                if (it1->sample(i).adc()!=it2->sample(i).adc())
                {
                    //  fill depth, uTCA and VME as well for which guys
                    //  mismatches happen
                    _cMsm_depth.fill(did);
                    _cMsm_FEDVME.fill(eid2);
                    _cMsm_FEDuTCA.fill(eid1);
                }
            }
    }
    for (HBHEDigiCollection::const_iterator it2=chbhe2->begin();
        it2!=chbhe2->end(); ++it2)
    {
        //  itearte thru VME 
        //  find utca digi by detid
        //  check if present or missing
        HcalDetId did = it2->id();
        HBHEDigiCollection::const_iterator it1 = chbhe1->find(did);
        if (it1==chbhe1->end())
        {
            HcalElectronicsId eid1 = HcalElectronicsId(_ehashmapuTCA.lookup(did));
            _cMsn_FEDuTCA.fill(eid1);
            for (int i=0; i<it2->size(); i++)
            {
                _cADCMsnuTCA_Subdet.fill(did, it2->sample(i).adc());
                _cADCall_Subdet.fill(did, -2, it2->sample(i).adc());
                _cADC_Subdet[i].fill(did, -2, it2->sample(i).adc());
            }
        }
    }
}

/* virtual */ void DigiComparisonTask::endLuminosityBlock(
    edm::LuminosityBlock const& lb, edm::EventSetup const& es)
{
    //  in the end always
    DQTask::endLuminosityBlock(lb, es);
}

DEFINE_FWK_MODULE(DigiComparisonTask);