RecoRunSummary.cc

Go to the documentation of this file.

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

namespace hcaldqm {
  using namespace constants;

  RecoRunSummary::RecoRunSummary(std::string const& name, std::string const& taskname, edm::ParameterSet const& ps)
      : DQClient(name, taskname, ps) {
    _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
    _thresh_tcds = ps.getUntrackedParameter<double>("thresh_tcds", 1.5);
  }

  /* virtual */ void RecoRunSummary::beginRun(edm::Run const& r, edm::EventSetup const& es) {
    DQClient::beginRun(r, es);
  }

  /*
     *
     */
  /* virtual */ void RecoRunSummary::endLuminosityBlock(DQMStore::IBooker& ib,
                                                        DQMStore::IGetter& ig,
                                                        edm::LuminosityBlock const& lb,
                                                        edm::EventSetup const& es) {
    DQClient::endLuminosityBlock(ib, ig, lb, es);
  }

  /*
     *
     */
  /* virtual */ std::vector<flag::Flag> RecoRunSummary::endJob(DQMStore::IBooker& ib, DQMStore::IGetter& ig) {
    if (_ptype != fOffline)
      return std::vector<flag::Flag>();

    // FILTERS, some useful vectors, hash maps
    std::vector<uint32_t> vhashCrateHF;
    vhashCrateHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashCrateHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashCrateHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    filter::HashFilter filter_CrateHF;
    filter_CrateHF.initialize(filter::fPreserver, hashfunctions::fCrate,
                              vhashCrateHF);  // preserve only HF crates
    electronicsmap::ElectronicsMap ehashmap;
    ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
    bool tcdsshift = false;
    std::vector<flag::Flag> vflags;
    vflags.resize(nRecoFlag);
    vflags[fUniSlotHF] = flag::Flag("UniSlotHF");
    vflags[fTCDS] = flag::Flag("TCDS");

    //  INITIALIZE
    Container2D cOccupancy_depth, cOccupancyCut_depth;
    ContainerSingle2D cSummary;
    Container1D cTimingCut_HBHEPartition;
    ContainerXXX<double> xUniHF, xUni;
    xUni.initialize(hashfunctions::fCrate);
    xUniHF.initialize(hashfunctions::fCrateSlot);
    cOccupancy_depth.initialize(_taskname,
                                "Occupancy",
                                hashfunctions::fdepth,
                                new quantity::DetectorQuantity(quantity::fieta),
                                new quantity::DetectorQuantity(quantity::fiphi),
                                new quantity::ValueQuantity(quantity::fN),
                                0);
    cOccupancyCut_depth.initialize(_taskname,
                                   "OccupancyCut",
                                   hashfunctions::fdepth,
                                   new quantity::DetectorQuantity(quantity::fieta),
                                   new quantity::DetectorQuantity(quantity::fiphi),
                                   new quantity::ValueQuantity(quantity::fN),
                                   0);
    cTimingCut_HBHEPartition.initialize(_taskname,
                                        "TimingCut",
                                        hashfunctions::fHBHEPartition,
                                        new quantity::ValueQuantity(quantity::fTiming_ns),
                                        new quantity::ValueQuantity(quantity::fN),
                                        0);

    cSummary.initialize(_name,
                        "Summary",
                        new quantity::CrateQuantity(_emap),
                        new quantity::FlagQuantity(vflags),
                        new quantity::ValueQuantity(quantity::fState),
                        0);

    //  BOOK
    xUniHF.book(_emap, filter_CrateHF);

    //  LOAD
    cOccupancy_depth.load(ig, _emap, _subsystem);
    cOccupancyCut_depth.load(ig, _emap, _subsystem);
    cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
    cSummary.book(ib, _subsystem);

    //  iterate over all channels
    std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = gids.begin(); it != gids.end(); ++it) {
      if (!it->isHcalDetId())
        continue;

      HcalDetId did(it->rawId());
      HcalElectronicsId eid = HcalElectronicsId(ehashmap.lookup(did));

      if (did.subdet() == HcalForward)
        xUniHF.get(eid) += cOccupancyCut_depth.getBinContent(did);
    }

    //  iphi/slot HF non uniformity
    for (doubleCompactMap::const_iterator it = xUniHF.begin(); it != xUniHF.end(); ++it) {
      uint32_t hash1 = it->first;
      HcalElectronicsId eid1(hash1);
      double x1 = it->second;
      for (doubleCompactMap::const_iterator jt = xUniHF.begin(); jt != xUniHF.end(); ++jt) {
        if (jt == it)
          continue;

        double x2 = jt->second;
        if (x2 == 0)
          continue;
        if (x1 / x2 < _thresh_unihf)
          xUni.get(eid1)++;
      }
    }

    //  TCDS shift
    double a = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 5, 1));
    double b = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 30, 1));
    double c = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 55, 1));
    double dab = fabs(a - b);
    double dac = fabs(a - c);
    double dbc = fabs(b - c);
    if (dab >= _thresh_tcds || dac >= _thresh_tcds || dbc >= _thresh_tcds)
      tcdsshift = true;

    //  summary flags
    std::vector<flag::Flag> sumflags;
    int icrate = 0;
    for (std::vector<uint32_t>::const_iterator it = _vhashCrates.begin(); it != _vhashCrates.end(); ++it) {
      flag::Flag fSum("RECO");
      HcalElectronicsId eid(*it);
      HcalDetId did = HcalDetId(_emap->lookup(eid));

      //    registered @cDAQ
      if (did.subdet() == HcalBarrel || did.subdet() == HcalEndcap) {
        if (tcdsshift)
          vflags[fTCDS]._state = flag::fBAD;
        else
          vflags[fTCDS]._state = flag::fGOOD;
      }
      if (did.subdet() == HcalForward) {
        if (xUni.get(eid) > 0)
          vflags[fUniSlotHF]._state = flag::fBAD;
        else
          vflags[fUniSlotHF]._state = flag::fGOOD;
      }

      //    combine
      int iflag = 0;
      for (std::vector<flag::Flag>::iterator ft = vflags.begin(); ft != vflags.end(); ++ft) {
        cSummary.setBinContent(eid, iflag, ft->_state);
        fSum += (*ft);
        iflag++;
        ft->reset();
      }
      sumflags.push_back(fSum);
      icrate++;
    }

    return sumflags;
  }
}  // namespace hcaldqm