QIE11Task.cc

Go to the documentation of this file.

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

using namespace hcaldqm;
using namespace hcaldqm::constants;
QIE11Task::QIE11Task(edm::ParameterSet const& ps):
    DQTask(ps)
{
    
    //  tags
    _tagQIE11 = ps.getUntrackedParameter<edm::InputTag>("tagQIE11",
        edm::InputTag("hcalDigis"));
    _tokQIE11 = consumes<QIE11DigiCollection>(_tagQIE11);

    _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN",
                              edm::InputTag("hcalDigis"));
    _tokuMN = consumes<HcalUMNioDigi>(_taguMN);

    //  cuts
    _cut = ps.getUntrackedParameter<double>("cut", 50.0);
    _ped = ps.getUntrackedParameter<int>("ped", 4);
    _laserType = ps.getUntrackedParameter<int32_t>("laserType", -1);
    _eventType = ps.getUntrackedParameter<int32_t>("eventType", -1);
}
/* virtual */ void QIE11Task::bookHistograms(DQMStore::IBooker &ib,
    edm::Run const& r, edm::EventSetup const& es)
{
    if (_ptype==fLocal)
        if (r.runAuxiliary().run()==1)
            return;

    DQTask::bookHistograms(ib, r, es);

    //  GET WHAT YOU NEED
    edm::ESHandle<HcalDbService> dbs;
    es.get<HcalDbRecord>().get(dbs);
    _emap = dbs->getHcalMapping();
    std::vector<uint32_t> vhashC34;
    vhashC34.push_back(HcalElectronicsId(34, 11,
        FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashC34.push_back(HcalElectronicsId(34, 12,
        FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    _filter_C34.initialize(filter::fPreserver, hcaldqm::hashfunctions::fCrateSlot,
        vhashC34);

    std::vector<std::pair<int, int> > timingChannels;
    timingChannels.push_back(std::pair<int, int>(28, 63));
    timingChannels.push_back(std::pair<int, int>(28, 65));
    timingChannels.push_back(std::pair<int, int>(20, 63));
    timingChannels.push_back(std::pair<int, int>(20, 65));
    for (int iChan = 0; iChan < 4; ++iChan) {
        std::vector<uint32_t> vhashTimingChannel;
        for (int depth = 1; depth <= 7; ++depth) {
            vhashTimingChannel.push_back(HcalDetId(HcalEndcap, timingChannels[iChan].first, timingChannels[iChan].second, depth));
        }
        _filter_timingChannels[iChan].initialize(filter::fPreserver, hcaldqm::hashfunctions::fDChannel, vhashTimingChannel);
    }

    //  INITIALIZE what you need
    
    // EChannel plots, online+local only
    if (_ptype != fOffline) {
        unsigned int itr = 0;
        for (unsigned int crate = 34; crate <= 34; ++crate) {
            for (unsigned int slot = 11; slot <= 12; ++slot) {
                std::vector<uint32_t> vhashSlot;
                vhashSlot.push_back(HcalElectronicsId(crate, slot, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
                _filter_slot[itr].initialize(filter::fPreserver, hashfunctions::fCrateSlot, vhashSlot);
                _cShapeCut_EChannel[itr].initialize(_name,
                    "ShapeCut", hcaldqm::hashfunctions::fEChannel,
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000));
                _cLETDCvsTS_EChannel[itr].initialize(_name, 
                    "LETDCvsTS", hcaldqm::hashfunctions::fEChannel, 
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true), 0);
                _cLETDCTime_EChannel[itr].initialize(_name,
                    "LETDCTime", hcaldqm::hashfunctions::fEChannel,
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
                for (unsigned int j=0; j<10; j++) {
                    _cLETDCvsADC_EChannel[j][itr].initialize(_name,
                        "LETDCvsADC", hcaldqm::hashfunctions::fEChannel,
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
                    _cADC_EChannel[j][itr].initialize(_name,
                        "ADC", hcaldqm::hashfunctions::fEChannel,
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
                    _cLETDC_EChannel[j][itr].initialize(_name,
                        "LETDC", hcaldqm::hashfunctions::fEChannel,
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
                }
                ++itr;
            }
        }
    }
    _cShapeCut.initialize(_name,
        "ShapeCut", 
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000));
    _cLETDCvsADC.initialize(_name, "LETDCvsADC",
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
    _cLETDCTimevsADC.initialize(_name, "LETDCTimevsADC",
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
    _cLETDC.initialize(_name, "LETDC",
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);
    _cADC.initialize(_name, "ADC",
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),0);

    if (_ptype != fOffline) {
        unsigned int itr = 0;
        std::map<std::pair<unsigned int, unsigned int>, unsigned int> itr_map;
        for(unsigned int crate = 34; crate <= 34; ++crate) {
            for(unsigned int slot=11; slot<=12; ++slot) {
                char aux[100];
                sprintf(aux, "/Crate%d_Slot%d", crate, slot);
                _cShapeCut_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
                _cLETDCvsTS_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
                _cLETDCTime_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
                for (unsigned int j=0; j<10; j++) {
                    char aux2[100];
                    sprintf(aux2, "/Crate%d_Slot%d/TS%d", crate, slot, j);
                    _cLETDCvsADC_EChannel[j][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux2);
                    _cLETDC_EChannel[j][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux2);
                    _cADC_EChannel[j][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux2);
                }
                itr_map[std::make_pair(crate, slot)] = itr;
                ++itr;
            }
        }
    }
    _cShapeCut.book(ib, _subsystem);
    _cLETDCvsADC.book(ib, _subsystem);
    _cLETDCTimevsADC.book(ib, _subsystem);
    _cLETDC.book(ib, _subsystem);
    _cADC.book(ib, _subsystem);

    _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap, _filter_C34);
}

/* virtual */ void QIE11Task::endLuminosityBlock(edm::LuminosityBlock const& lb,
    edm::EventSetup const& es)
{
    
    //  finish
    DQTask::endLuminosityBlock(lb, es);
}

/* virtual */ void QIE11Task::_process(edm::Event const& e, 
    edm::EventSetup const&)
{
    edm::Handle<QIE11DigiCollection> cqie11;
    if (!e.getByToken(_tokQIE11, cqie11))
        return;

    for (uint32_t i=0; i<cqie11->size(); i++)
    {
        QIE11DataFrame frame = static_cast<QIE11DataFrame>((*cqie11)[i]);
        DetId did = frame.detid();
        if (HcalDetId(frame.detid()).subdet() != HcalEndcap) {
            continue;
        }

        HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
        if (!eid.rawId()) {
          continue;
        }
        int fakecrate = -1;
        if (eid.crateId() == 34) fakecrate = 0;
        int index = fakecrate * 12 + (eid.slot() - 10) - 1;

        //  compute the signal, ped subracted
        CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, frame);
//      double q = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(frame,
//          constants::adc2fC[_ped], 0, frame.samples()-1);

        //  iterate thru all TS and fill
        for (int j=0; j<frame.samples(); j++)
        {
            double q_pedsub = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, frame, j);
            if (_ptype != fOffline) {
                if (index == 0 || index == 1) {
                    //  shapes are after the cut
                    _cShapeCut_EChannel[index].fill(eid, j, q_pedsub);  
                    _cLETDCvsTS_EChannel[index].fill(eid, j, frame[j].tdc());   

                    //  w/o a cut
                    _cLETDCvsADC_EChannel[j][index].fill(eid, frame[j].adc(), 
                        frame[j].tdc());
                    _cLETDC_EChannel[j][index].fill(eid, frame[j].tdc());
                    if (frame[j].tdc() < 50) {
                        // Each TDC count is 0.5 ns. 
                        // tdc == 62 or 63 means value was below or above threshold for whole time slice. 
                        _cLETDCTime_EChannel[index].fill(eid, j*25. + (frame[j].tdc() / 2.));
                    }
                    _cADC_EChannel[j][index].fill(eid, frame[j].adc());
                }
            }
            _cShapeCut.fill(eid, j, q_pedsub);

            _cLETDCvsADC.fill(frame[j].adc(), frame[j].tdc());
            if (frame[j].tdc() < 50) {
                _cLETDCTimevsADC.fill(frame[j].adc(), j*25. + (frame[j].tdc() / 2.));
            }
            _cLETDC.fill(eid, frame[j].tdc());

            _cADC.fill(eid, frame[j].adc());

        }
    }
}


/* virtual */ bool QIE11Task::_isApplicable(edm::Event const& e)
{
  if (_ptype!=fOnline || (_laserType < 0 && _eventType < 0))
    return true;
  else
    {
      //      fOnline mode
      edm::Handle<HcalUMNioDigi> cumn;
      if (!e.getByToken(_tokuMN, cumn))
    return false;

      //      event type check first
      int eventType = cumn->eventType();
      if (eventType==_eventType)
    return true;

      //      check if this analysis task is of the right laser type
      int laserType = cumn->valueUserWord(0);
      if (laserType==_laserType)
    return true;
    }

  return false;
}


/* virtual */ void QIE11Task::_resetMonitors(hcaldqm::UpdateFreq)
{
}

DEFINE_FWK_MODULE(QIE11Task);