RecHitTask.cc

Go to the documentation of this file.

#include "DQM/HcalTasks/interface/RecHitTask.h"
#include <cmath>

using namespace hcaldqm;
using namespace hcaldqm::constants;
using namespace hcaldqm::filter;

RecHitTask::RecHitTask(edm::ParameterSet const& ps)
    : DQTask(ps), hcalDbServiceToken_(esConsumes<HcalDbService, HcalDbRecord, edm::Transition::BeginRun>()) {
  _tagHBHE = ps.getUntrackedParameter<edm::InputTag>("tagHBHE", edm::InputTag("hbhereco"));
  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO", edm::InputTag("horeco"));
  _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hfreco"));
  _tagPreHF = ps.getUntrackedParameter<edm::InputTag>("tagPreHF", edm::InputTag(""));
  _hfPreRecHitsAvailable = ps.getUntrackedParameter<bool>("hfPreRecHitsAvailable", false);

  _tokHBHE = consumes<HBHERecHitCollection>(_tagHBHE);
  _tokHO = consumes<HORecHitCollection>(_tagHO);
  _tokHF = consumes<HFRecHitCollection>(_tagHF);
  _tokPreHF = consumes<HFPreRecHitCollection>(_tagPreHF);

  _cutE_HBHE = ps.getUntrackedParameter<double>("cutE_HBHE", 5);
  _cutE_HO = ps.getUntrackedParameter<double>("cutE_HO", 5);
  _cutE_HF = ps.getUntrackedParameter<double>("cutE_HF", 5);
  _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);

  //    order must be the same as in RecoFlag enum
  _vflags.resize(nRecoFlag);
  _vflags[fUni] = flag::Flag("UniSlotHF");
  _vflags[fTCDS] = flag::Flag("TCDS");
  _vflags[fUnknownIds] = flag::Flag("UnknownIds");
}

/* virtual */ void RecHitTask::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.getHandle(hcalDbServiceToken_);
  _emap = dbs->getHcalMapping();

  std::vector<uint32_t> vVME;
  std::vector<uint32_t> vuTCA;
  vVME.push_back(
      HcalElectronicsId(constants::FIBERCH_MIN, constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
  vuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
  _filter_VME.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vVME);
  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vuTCA);
  std::vector<uint32_t> vhashHF;
  vhashHF.push_back(
      hcaldqm::hashfunctions::hash_did[hcaldqm::hashfunctions::fSubdet](HcalDetId(HcalForward, 29, 1, 1)));

  _filter_HF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fSubdet, vhashHF);

  //    INITIALIZE FIRST
  //    Energy
  _cEnergy_Subdet.initialize(_name,
                             "Energy",
                             hcaldqm::hashfunctions::fSubdet,
                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy),
                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                             0);
  _cEnergy_depth.initialize(_name,
                            "Energy",
                            hcaldqm::hashfunctions::fdepth,
                            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
                            0);

  //    Timing
  _cTimingCut_SubdetPM.initialize(_name,
                                  "TimingCut",
                                  hcaldqm::hashfunctions::fSubdetPM,
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                  0);
  _cTimingvsEnergy_SubdetPM.initialize(_name,
                                       "TimingvsEnergy",
                                       hcaldqm::hashfunctions::fSubdetPM,
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy, true),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                       0);

  _cTimingCut_HBHEPartition.initialize(_name,
                                       "TimingCut",
                                       hcaldqm::hashfunctions::fHBHEPartition,
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                       0);
  _cTimingCut_depth.initialize(_name,
                               "TimingCut",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                               0);
  _cTimingCutvsLS_SubdetPM.initialize(_name,
                                      "TimingCutvsLS",
                                      hcaldqm::hashfunctions::fSubdetPM,
                                      new hcaldqm::quantity::LumiSection(_maxLS),
                                      new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                      0);

  //    Occupancy
  _cOccupancy_depth.initialize(_name,
                               "Occupancy",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                               0);

  _cOccupancyvsLS_Subdet.initialize(_name,
                                    "OccupancyvsLS",
                                    hcaldqm::hashfunctions::fSubdet,
                                    new hcaldqm::quantity::LumiSection(_maxLS),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
                                    0);

  _cOccupancyCut_depth.initialize(_name,
                                  "OccupancyCut",
                                  hcaldqm::hashfunctions::fdepth,
                                  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                  new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                  0);

  if (_hfPreRecHitsAvailable) {
    _cDAAsymmetryVsCharge_SubdetPM.initialize(
        _name,
        "ChargeVsAsymmetry",
        hcaldqm::hashfunctions::fSubdetPM,
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
        0);
    _cDAAsymmetryMean_cut_depth.initialize(_name,
                                           "AsymmetryMean",
                                           hcaldqm::hashfunctions::fdepth,
                                           new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                           new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
                                           0);
    _cDAAsymmetry_cut_SubdetPM.initialize(_name,
                                          "Asymmetry",
                                          hcaldqm::hashfunctions::fSubdetPM,
                                          new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fDualAnodeAsymmetry),
                                          new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                          0);
  }

  //    INITIALIZE HISTOGRAMS to be used only in Online
  if (_ptype == fOnline) {
    _cEnergyvsieta_Subdet.initialize(_name,
                                     "Energyvsieta",
                                     hcaldqm::hashfunctions::fSubdet,
                                     new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                     new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
                                     0);
    _cEnergyvsiphi_SubdetPM.initialize(_name,
                                       "Energyvsiphi",
                                       hcaldqm::hashfunctions::fSubdetPM,
                                       new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
                                       0);
    _cEnergyvsLS_SubdetPM.initialize(_name,
                                     "EnergyvsLS",
                                     hcaldqm::hashfunctions::fSubdetPM,
                                     new hcaldqm::quantity::LumiSection(_maxLS),
                                     new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
                                     0);
    _cEnergyvsBX_SubdetPM.initialize(_name,
                                     "EnergyvsBX",
                                     hcaldqm::hashfunctions::fSubdetPM,
                                     new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
                                     new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fEnergy_1TeV),
                                     0);
    _cTimingCutvsieta_Subdet.initialize(_name,
                                        "TimingCutvsieta",
                                        hcaldqm::hashfunctions::fSubdet,
                                        new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                        0);
    _cTimingCutvsiphi_SubdetPM.initialize(_name,
                                          "TimingCutvsiphi",
                                          hcaldqm::hashfunctions::fSubdetPM,
                                          new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                          new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                          0);
    _cTimingCutvsBX_SubdetPM.initialize(_name,
                                        "TimingCutvsBX",
                                        hcaldqm::hashfunctions::fSubdetPM,
                                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
                                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                        0);
    _cOccupancyvsiphi_SubdetPM.initialize(_name,
                                          "Occupancyvsiphi",
                                          hcaldqm::hashfunctions::fSubdetPM,
                                          new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                          new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                          0);
    _cOccupancyvsieta_Subdet.initialize(_name,
                                        "Occupancyvsieta",
                                        hcaldqm::hashfunctions::fSubdet,
                                        new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                        new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                        0);
    _cOccupancyCutvsiphi_SubdetPM.initialize(_name,
                                             "OccupancyCutvsiphi",
                                             hcaldqm::hashfunctions::fSubdetPM,
                                             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                             0);
    _cOccupancyCutvsieta_Subdet.initialize(_name,
                                           "OccupancyCutvsieta",
                                           hcaldqm::hashfunctions::fSubdet,
                                           new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                           0);
    _cOccupancyCutvsBX_Subdet.initialize(_name,
                                         "OccupancyCutvsBX",
                                         hcaldqm::hashfunctions::fSubdet,
                                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fBX),
                                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                         0);
    _cOccupancyCutvsiphivsLS_SubdetPM.initialize(_name,
                                                 "OccupancyCutvsiphivsLS",
                                                 hcaldqm::hashfunctions::fSubdetPM,
                                                 new hcaldqm::quantity::LumiSection(_maxLS),
                                                 new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                                 0);
    _cOccupancyCutvsLS_Subdet.initialize(_name,
                                         "OccupancyCutvsLS",
                                         hcaldqm::hashfunctions::fSubdet,
                                         new hcaldqm::quantity::LumiSection(_maxLS),
                                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN_to8000),
                                         0);
  }

  // FED-based plots
  if (_ptype != fOffline) {  // hidefed2crate
    std::vector<int> vFEDs = hcaldqm::utilities::getFEDList(_emap);
    std::vector<int> vFEDsVME = hcaldqm::utilities::getFEDVMEList(_emap);
    std::vector<int> vFEDsuTCA = hcaldqm::utilities::getFEDuTCAList(_emap);

    //  push the rawIds of each fed into the vector
    for (std::vector<int>::const_iterator it = vFEDsVME.begin(); it != vFEDsVME.end(); ++it)
      _vhashFEDs.push_back(HcalElectronicsId(FIBERCH_MIN, FIBER_VME_MIN, SPIGOT_MIN, (*it) - FED_VME_MIN).rawId());
    for (std::vector<int>::const_iterator it = vFEDsuTCA.begin(); it != vFEDsuTCA.end(); ++it) {
      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());
    }

    _cTimingCut_FEDuTCA.initialize(_name,
                                   "TimingCut",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                   0);
    _cTimingCutvsLS_FED.initialize(_name,
                                   "TimingCutvsLS",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::LumiSection(_maxLS),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                   0);
    _cTimingCut_ElectronicsuTCA.initialize(_name,
                                           "TimingCut",
                                           hcaldqm::hashfunctions::fElectronics,
                                           new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
                                           new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_ns),
                                           0);

    _cOccupancy_FEDuTCA.initialize(_name,
                                   "Occupancy",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                   0);
    _cOccupancy_ElectronicsuTCA.initialize(_name,
                                           "Occupancy",
                                           hcaldqm::hashfunctions::fElectronics,
                                           new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
                                           new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                           0);

    _cOccupancyCut_FEDuTCA.initialize(_name,
                                      "OccupancyCut",
                                      hcaldqm::hashfunctions::fFED,
                                      new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                      new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                      new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                      0);
    _cOccupancyCut_ElectronicsuTCA.initialize(_name,
                                              "OccupancyCut",
                                              hcaldqm::hashfunctions::fElectronics,
                                              new hcaldqm::quantity::FEDQuantity(vFEDsuTCA),
                                              new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                              new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                              0);
    if (_ptype == fOnline) {
      _cSummaryvsLS_FED.initialize(_name,
                                   "SummaryvsLS",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::LumiSection(_maxLS),
                                   new hcaldqm::quantity::FlagQuantity(_vflags),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
                                   0);
      _cSummaryvsLS.initialize(_name,
                               "SummaryvsLS",
                               new hcaldqm::quantity::LumiSection(_maxLS),
                               new hcaldqm::quantity::FEDQuantity(vFEDs),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fState),
                               0);

      _xUniHF.initialize(hcaldqm::hashfunctions::fFEDSlot);
      _xUni.initialize(hcaldqm::hashfunctions::fFED);
    }
  }

  //    BOOK HISTOGRAMS
  char cutstr[200];
  sprintf(cutstr, "_EHBHE%dHO%dHF%d", int(_cutE_HBHE), int(_cutE_HO), int(_cutE_HF));
  char cutstr2[200];
  sprintf(cutstr2, "_EHF%d", int(_cutE_HF));

  //    Energy
  _cEnergy_Subdet.book(ib, _emap, _subsystem);
  _cEnergy_depth.book(ib, _emap, _subsystem);

  //    Timing
  _cTimingCut_SubdetPM.book(ib, _emap, _subsystem);
  _cTimingvsEnergy_SubdetPM.book(ib, _emap, _subsystem);
  if (_ptype != fOffline) {  // hidefed2crate
    _cTimingCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cTimingCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
  }
  _cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
  _cTimingCut_depth.book(ib, _emap, _subsystem);
  _cTimingCutvsLS_SubdetPM.book(ib, _emap, _subsystem);
  if (_ptype != fOffline) {  // hidefed2crate
    _cTimingCutvsLS_FED.book(ib, _emap, _subsystem);
  }

  //    Occupancy
  _cOccupancy_depth.book(ib, _emap, _subsystem);
  if (_ptype != fOffline) {  // hidefed2crate
    _cOccupancy_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cOccupancy_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
  }
  _cOccupancyvsLS_Subdet.book(ib, _emap, _subsystem);
  _cOccupancyCut_depth.book(ib, _emap, _subsystem);
  if (_ptype != fOffline) {  // hidefed2crate
    _cOccupancyCut_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cOccupancyCut_ElectronicsuTCA.book(ib, _emap, _filter_VME, _subsystem);
  }

  if (_hfPreRecHitsAvailable) {
    _cDAAsymmetryVsCharge_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
    _cDAAsymmetryMean_cut_depth.book(ib, _emap, _filter_HF, _subsystem);
    _cDAAsymmetry_cut_SubdetPM.book(ib, _emap, _filter_HF, _subsystem);
  }

  //    BOOK HISTOGRAMS to be used only in Online
  if (_ptype == fOnline) {
    _cEnergyvsLS_SubdetPM.book(ib, _emap, _subsystem);
    _cEnergyvsieta_Subdet.book(ib, _emap, _subsystem);
    _cEnergyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
    _cEnergyvsBX_SubdetPM.book(ib, _emap, _subsystem);
    _cTimingCutvsieta_Subdet.book(ib, _emap, _subsystem);
    _cTimingCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
    _cTimingCutvsBX_SubdetPM.book(ib, _emap, _subsystem);
    _cOccupancyvsiphi_SubdetPM.book(ib, _emap, _subsystem);
    _cOccupancyvsieta_Subdet.book(ib, _emap, _subsystem);
    _cOccupancyCutvsiphi_SubdetPM.book(ib, _emap, _subsystem);
    _cOccupancyCutvsieta_Subdet.book(ib, _emap, _subsystem);
    _cOccupancyCutvsBX_Subdet.book(ib, _emap, _subsystem);
    _cOccupancyCutvsiphivsLS_SubdetPM.book(ib, _emap, _subsystem);
    _cOccupancyCutvsLS_Subdet.book(ib, _emap, _subsystem);
    _cSummaryvsLS_FED.book(ib, _emap, _subsystem);
    _cSummaryvsLS.book(ib, _subsystem);

    std::vector<uint32_t> vhashFEDHF;
    vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN + 6, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    HashFilter filter_FEDHF;
    filter_FEDHF.initialize(filter::fPreserver, hcaldqm::hashfunctions::fFED, vhashFEDHF);

    _gids = _emap->allPrecisionId();
    if (_ptype != fOffline) {  // hidefed2crate
      _xUniHF.book(_emap, filter_FEDHF);
      _xUni.book(_emap);
    }
  }

  //    initialize hash map
  _ehashmap.initialize(_emap, hcaldqm::electronicsmap::fD2EHashMap);

  //    book some mes...
  ib.setCurrentFolder(_subsystem + "/" + _name);
  auto scope = DQMStore::IBooker::UseLumiScope(ib);
  meUnknownIds1LS = ib.book1DD("UnknownIds", "UnknownIds", 1, 0, 1);
  _unknownIdsPresent = false;
}

/* virtual */ void RecHitTask::_resetMonitors(hcaldqm::UpdateFreq uf) {
  switch (uf) {
    case hcaldqm::f1LS:
      _unknownIdsPresent = false;
      break;
    default:
      break;
  }

  DQTask::_resetMonitors(uf);
}

/* virtual */ void RecHitTask::_process(edm::Event const& e, edm::EventSetup const&) {
  edm::Handle<HBHERecHitCollection> chbhe;
  edm::Handle<HORecHitCollection> cho;
  edm::Handle<HFRecHitCollection> chf;

  if (!(e.getByToken(_tokHBHE, chbhe)))
    _logger.dqmthrow("Collection HBHERecHitCollection not available " + _tagHBHE.label() + " " + _tagHBHE.instance());
  if (!(e.getByToken(_tokHO, cho)))
    _logger.dqmthrow("Collection HORecHitCollection not available " + _tagHO.label() + " " + _tagHO.instance());
  if (!(e.getByToken(_tokHF, chf)))
    _logger.dqmthrow("Collection HFRecHitCollection not available " + _tagHF.label() + " " + _tagHF.instance());

  edm::Handle<HFPreRecHitCollection> cprehf;
  if (_hfPreRecHitsAvailable) {
    if (!(e.getByToken(_tokPreHF, cprehf)))
      _logger.dqmthrow("Collection HFPreRecHitCollection not available " + _tagPreHF.label() + " " +
                       _tagPreHF.instance());
  }

  //    extract some info per event
  int bx = e.bunchCrossing();

  auto lumiCache = luminosityBlockCache(e.getLuminosityBlock().index());
  _currentLS = lumiCache->currentLS;

  //  To fill histograms outside of the loop, you need to determine if there were
  //  any valid det ids first
  uint32_t rawidHBValid = 0;
  uint32_t rawidHEValid = 0;

  double ehbm = 0;
  double ehbp = 0;
  double ehem = 0;
  double ehep = 0;
  int nChsHB = 0;
  int nChsHE = 0;
  int nChsHBCut = 0;
  int nChsHECut = 0;
  for (HBHERecHitCollection::const_iterator it = chbhe->begin(); it != chbhe->end(); ++it) {
    double energy = it->energy();
    double timing = it->time();

    //  Explicit check on the DetIds present in the Collection
    HcalDetId did = it->id();
    uint32_t rawid = _ehashmap.lookup(did);
    /*
         * Needs to be removed as DetIds that belong to the HEP17 after combination
         * are not present in the emap
         * Removed until further notice!
         *
         */
    //if (rawid==0)
    //{meUnknownIds1LS->Fill(1); _unknownIdsPresent=true;continue;}

    HcalElectronicsId const& eid(rawid);
    if (did.subdet() == HcalBarrel)
      rawidHBValid = did.rawId();
    else if (did.subdet() == HcalEndcap)
      rawidHEValid = did.rawId();

    _cEnergy_Subdet.fill(did, energy);
    _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
    _cOccupancy_depth.fill(did);
    did.subdet() == HcalBarrel ? did.ieta() > 0 ? ehbp += energy : ehbm += energy
    : did.ieta() > 0           ? ehep += energy
                               : ehem += energy;

    //  ONLINE ONLY!
    if (_ptype == fOnline) {
      _cOccupancyvsiphi_SubdetPM.fill(did);
      _cOccupancyvsieta_Subdet.fill(did);
    }
    //  ^^^ONLINE ONLY!
    //
    if (_ptype != fOffline) {  // hidefed2crate
      //    Also, for these electronics plots, require that the channel was found in the emap.
      if (rawid != 0) {
        if (!eid.isVMEid()) {
          _cOccupancy_FEDuTCA.fill(eid);
          _cOccupancy_ElectronicsuTCA.fill(eid);
        }
      }
    }

    if (energy > _cutE_HBHE) {
      //    ONLINE ONLY!
      if (_ptype == fOnline) {
        _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
        _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
        _cEnergyvsieta_Subdet.fill(did, energy);
        _cEnergyvsiphi_SubdetPM.fill(did, energy);
        _cTimingCutvsieta_Subdet.fill(did, timing);
        _cTimingCutvsiphi_SubdetPM.fill(did, timing);
        _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
        _cOccupancyCutvsiphi_SubdetPM.fill(did);
        _cOccupancyCutvsieta_Subdet.fill(did);
        _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
      }
      //    ^^^ONLINE ONLY!
      _cEnergy_depth.fill(did, energy);
      _cTimingCut_SubdetPM.fill(did, timing);
      _cTimingCut_HBHEPartition.fill(did, timing);
      _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);

      //    ONLINE
      if (_ptype == fOnline) {
        if (rawid != 0) {
          _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
        }
        _cTimingCut_depth.fill(did, timing);
        //  ^^^ONLINE
      } else {
        if (_ptype != fOffline) {  // hidefed2crate
          if (rawid != 0) {
            _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
          }
        }
        _cTimingCut_depth.fill(did, timing);
      }
      _cOccupancyCut_depth.fill(did);
      if (_ptype != fOffline) {  // hidefed2crate
        if (rawid != 0) {
          if (eid.isUTCAid()) {
            if (_ptype == fOnline) {
              //    time constraints are explicit!
              _cTimingCut_FEDuTCA.fill(eid, timing);
              _cTimingCut_ElectronicsuTCA.fill(eid, timing);
            } else {
              _cTimingCut_FEDuTCA.fill(eid, timing);
              _cTimingCut_ElectronicsuTCA.fill(eid, timing);
            }
            _cOccupancyCut_FEDuTCA.fill(eid);
            _cOccupancyCut_ElectronicsuTCA.fill(eid);
          }
        }
      }
      did.subdet() == HcalBarrel ? nChsHBCut++ : nChsHECut++;
    }
    did.subdet() == HcalBarrel ? nChsHB++ : nChsHE++;
  }

  if (rawidHBValid != 0 && rawidHEValid != 0) {
    _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, nChsHB);
    _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, nChsHE);

    //  ONLINE ONLY!
    if (_ptype == fOnline) {
      _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHBValid), bx, nChsHBCut);
      _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidHEValid), bx, nChsHECut);
      _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHBValid), _currentLS, nChsHBCut);
      _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidHEValid), _currentLS, nChsHECut);
    }
    //  ^^^ONLINE ONLY!
  }

  //  To fill histograms outside of the loop, you need to determine if there were
  //  any valid det ids first
  uint32_t rawidValid = 0;

  int nChsHO = 0;
  int nChsHOCut = 0;
  double ehop = 0;
  double ehom = 0;
  for (HORecHitCollection::const_iterator it = cho->begin(); it != cho->end(); ++it) {
    double energy = it->energy();
    double timing = it->time();

    //  Explicit check on the DetIds present in the Collection
    HcalDetId did = it->id();
    uint32_t rawid = _ehashmap.lookup(did);
    if (rawid == 0) {
      meUnknownIds1LS->Fill(1);
      _unknownIdsPresent = true;
      continue;
    }
    HcalElectronicsId const& eid(rawid);
    if (did.subdet() == HcalOuter)
      rawidValid = did.rawId();

    _cEnergy_Subdet.fill(did, energy);
    _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
    _cOccupancy_depth.fill(did);
    did.ieta() > 0 ? ehop += energy : ehom += energy;

    //  IMPORTANT: ONLINE ONLY!
    if (_ptype == fOnline) {
      _cOccupancyvsiphi_SubdetPM.fill(did);
      _cOccupancyvsieta_Subdet.fill(did);
    }
    //  ONLINE ONLY!

    if (_ptype != fOffline) {  // hidefed2crate
      if (!eid.isVMEid()) {
        _cOccupancy_FEDuTCA.fill(eid);
        _cOccupancy_ElectronicsuTCA.fill(eid);
      }
    }

    if (energy > _cutE_HO) {
      //    ONLINE ONLY!
      if (_ptype == fOnline) {
        _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
        _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
        _cEnergyvsieta_Subdet.fill(did, energy);
        _cEnergyvsiphi_SubdetPM.fill(did, energy);
        _cTimingCutvsieta_Subdet.fill(did, timing);
        _cTimingCutvsiphi_SubdetPM.fill(did, timing);
        _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
        _cOccupancyCutvsiphi_SubdetPM.fill(did);
        _cOccupancyCutvsieta_Subdet.fill(did);
        _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
      }
      //    ^^^ONLINE ONLY!

      _cEnergy_depth.fill(did, energy);
      _cTimingCut_SubdetPM.fill(did, timing);
      _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
      if (_ptype != fOffline) {  // hidefed2crate
        _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
      }
      _cOccupancyCut_depth.fill(did);
      _cTimingCut_depth.fill(did, timing);
      if (_ptype != fOffline) {  // hidefed2crate
        if (!eid.isVMEid()) {
          _cTimingCut_FEDuTCA.fill(eid, timing);
          _cTimingCut_ElectronicsuTCA.fill(eid, timing);
          _cOccupancyCut_FEDuTCA.fill(eid);
          _cOccupancyCut_ElectronicsuTCA.fill(eid);
        }
      }
      nChsHOCut++;
    }
    nChsHO++;
  }

  if (rawidValid != 0) {
    _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHO);
    //  ONLINE ONLY!
    if (_ptype == fOnline) {
      _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, nChsHOCut);
      _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHOCut);
    }
    //  ^^^ONLINE ONLY!
  }

  //reset
  rawidValid = 0;

  int nChsHF = 0;
  int nChsHFCut = 0;
  double ehfp = 0;
  double ehfm = 0;
  for (HFRecHitCollection::const_iterator it = chf->begin(); it != chf->end(); ++it) {
    double energy = it->energy();
    double timing = it->time();

    //  Explicit check on the DetIds present in the Collection
    HcalDetId did = it->id();
    uint32_t rawid = _ehashmap.lookup(did);
    if (rawid == 0) {
      meUnknownIds1LS->Fill(1);
      _unknownIdsPresent = true;
      continue;
    }
    HcalElectronicsId const& eid(rawid);
    if (did.subdet() == HcalForward)
      rawidValid = did.rawId();

    _cEnergy_Subdet.fill(did, energy);
    _cTimingvsEnergy_SubdetPM.fill(did, energy, timing);
    _cOccupancy_depth.fill(did);
    did.ieta() > 0 ? ehfp += energy : ehfm += energy;

    //  IMPORTANT:
    //  only for Online Processing
    //
    if (_ptype == fOnline) {
      _cOccupancyvsiphi_SubdetPM.fill(did);
      _cOccupancyvsieta_Subdet.fill(did);
    }
    //  ONLINE ONLY!

    if (_ptype != fOffline) {  // hidefed2crate
      if (!eid.isVMEid()) {
        _cOccupancy_FEDuTCA.fill(eid);
        _cOccupancy_ElectronicsuTCA.fill(eid);
      }
    }

    if (energy > _cutE_HF) {
      //    ONLINE ONLY!
      if (_ptype == fOnline) {
        _cEnergyvsLS_SubdetPM.fill(did, _currentLS, energy);
        _cEnergyvsBX_SubdetPM.fill(did, bx, energy);
        _cEnergyvsieta_Subdet.fill(did, energy);
        _cEnergyvsiphi_SubdetPM.fill(did, energy);
        _cTimingCutvsieta_Subdet.fill(did, timing);
        _cTimingCutvsiphi_SubdetPM.fill(did, timing);
        _cTimingCutvsBX_SubdetPM.fill(did, bx, timing);
        _cOccupancyCutvsiphi_SubdetPM.fill(did);
        _cOccupancyCutvsieta_Subdet.fill(did);
        _cOccupancyCutvsiphivsLS_SubdetPM.fill(did, _currentLS);
        if (_ptype != fOffline) {  // hidefed2crate
          _xUniHF.get(eid)++;
        }
      }
      //    ^^^ONLINE ONLY!
      _cEnergy_depth.fill(did, energy);
      _cTimingCut_SubdetPM.fill(did, timing);
      _cTimingCutvsLS_SubdetPM.fill(did, _currentLS, timing);
      if (_ptype != fOffline) {  // hidefed2crate
        _cTimingCutvsLS_FED.fill(eid, _currentLS, timing);
      }
      _cOccupancyCut_depth.fill(did);
      _cTimingCut_depth.fill(did, timing);
      if (_ptype != fOffline) {  // hidefed2crate
        if (!eid.isVMEid()) {
          _cTimingCut_FEDuTCA.fill(eid, timing);
          _cTimingCut_ElectronicsuTCA.fill(eid, timing);
          _cOccupancyCut_FEDuTCA.fill(eid);
          _cOccupancyCut_ElectronicsuTCA.fill(eid);
        }
      }
      nChsHFCut++;
    }
    nChsHF++;
  }

  if (rawidValid != 0) {
    _cOccupancyvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHF);
    //  ONLINE ONLY!
    if (_ptype == fOnline) {
      _cOccupancyCutvsBX_Subdet.fill(HcalDetId(rawidValid), bx, nChsHFCut);
      _cOccupancyCutvsLS_Subdet.fill(HcalDetId(rawidValid), _currentLS, nChsHFCut);
    }
    //  ^^^ONLINE ONLY!
  }

  // Loop over HFPreRecHits to get charge and charge asymmetry
  if (_hfPreRecHitsAvailable) {
    for (HFPreRecHitCollection::const_iterator it = cprehf->begin(); it != cprehf->end(); ++it) {
      HcalDetId did = it->id();
      if (_filter_HF.filter(did)) {
        continue;
      }
      std::pair<float, bool> chargeAsymmetry = it->chargeAsymmetry(0.);
      std::pair<float, bool> chargeAsymmetryCut = it->chargeAsymmetry(20.);

      if (chargeAsymmetry.second) {
        _cDAAsymmetryVsCharge_SubdetPM.fill(did, chargeAsymmetry.first, it->charge());
      }
      if (chargeAsymmetryCut.second) {
        _cDAAsymmetryMean_cut_depth.fill(did, chargeAsymmetryCut.first);
        _cDAAsymmetry_cut_SubdetPM.fill(did, chargeAsymmetryCut.first);
      }
    }
  }
}

std::shared_ptr<hcaldqm::Cache> RecHitTask::globalBeginLuminosityBlock(edm::LuminosityBlock const& lb,
                                                                       edm::EventSetup const& es) const {
  return DQTask::globalBeginLuminosityBlock(lb, es);
}

/* virtual */ void RecHitTask::globalEndLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
  if (_ptype != fOnline)
    return;

  auto lumiCache = luminosityBlockCache(lb.index());
  _currentLS = lumiCache->currentLS;

  //
  //    GENERATE STATUS ONLY FOR ONLINE
  //
  //    for (std::vector<HcalGenericDetId>::const_iterator it=gids.begin();
  //        it!=gids.end(); ++it)
  //    {}

  for (uintCompactMap::const_iterator it = _xUniHF.begin(); it != _xUniHF.end(); ++it) {
    uint32_t hash1 = it->first;
    HcalElectronicsId eid1(hash1);
    double x1 = it->second;

    for (uintCompactMap::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)++;
    }
  }

  if (_ptype != fOffline) {  // hidefed2crate
    for (std::vector<uint32_t>::const_iterator it = _vhashFEDs.begin(); it != _vhashFEDs.end(); ++it) {
      flag::Flag fSum("RECO");
      HcalElectronicsId eid = HcalElectronicsId(*it);

      std::vector<uint32_t>::const_iterator cit = std::find(_vcdaqEids.begin(), _vcdaqEids.end(), *it);
      if (cit == _vcdaqEids.end()) {
        //  not @cDAQ
        for (uint32_t iflag = 0; iflag < _vflags.size(); iflag++)
          _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(flag::fNCDAQ));
        _cSummaryvsLS.setBinContent(eid, _currentLS, int(flag::fNCDAQ));
        continue;
      }

      //    FED is @cDAQ
      if (hcaldqm::utilities::isFEDHF(eid) && (_runkeyVal == 0 || _runkeyVal == 4)) {
        if (_xUni.get(eid) > 0)
          _vflags[fUni]._state = flag::fPROBLEMATIC;
        else
          _vflags[fUni]._state = flag::fGOOD;
      }

      if (_unknownIdsPresent)
        _vflags[fUnknownIds]._state = flag::fBAD;
      else
        _vflags[fUnknownIds]._state = flag::fGOOD;

      int iflag = 0;
      for (std::vector<flag::Flag>::iterator ft = _vflags.begin(); ft != _vflags.end(); ++ft) {
        _cSummaryvsLS_FED.setBinContent(eid, _currentLS, int(iflag), int(ft->_state));
        fSum += (*ft);
        iflag++;

        //  reset after using
        ft->reset();
      }
      _cSummaryvsLS.setBinContent(eid, _currentLS, fSum._state);
    }
    _xUniHF.reset();
    _xUni.reset();
  }

  //    in the end always do the DQTask::endLumi
  DQTask::globalEndLuminosityBlock(lb, es);
}

DEFINE_FWK_MODULE(RecHitTask);