DQM/HcalTasks/plugins/LEDTask.cc - master - githubmirror/cmssw

DQM/HcalTasks/plugins/LEDTask.cc
957 строк43 KB
Deniz Sunar Cerci
[DQM] Add LS-normalized HCAL LED trend plots
01 апр 2026, 20:37
01 апр 2026, 20:373f2a43d

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

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

namespace {
  constexpr double kTiny = 1e-9;

  class FineTimingDiffNs : public hcaldqm::quantity::ValueQuantity {
  public:
    FineTimingDiffNs() : ValueQuantity(hcaldqm::quantity::fTimingDiff_ns) {}

    FineTimingDiffNs* makeCopy() override { return new FineTimingDiffNs(); }
    int nbins() override { return 100; }
    double min() override { return -50.; }
    double max() override { return 50.; }
  };
}  // namespace

LEDTask::LEDTask(edm::ParameterSet const& ps)
    : DQTask(ps), hcalDbServiceToken_(esConsumes<HcalDbService, HcalDbRecord, edm::Transition::BeginRun>()) {
  _nevents = ps.getUntrackedParameter<int>("nevents", 2000);
  //	tags
  _tagQIE11 = ps.getUntrackedParameter<edm::InputTag>("tagHE", edm::InputTag("hcalDigis"));
  _tagHO = ps.getUntrackedParameter<edm::InputTag>("tagHO", edm::InputTag("hcalDigis"));
  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hcalDigis"));
  _tagTrigger = ps.getUntrackedParameter<edm::InputTag>("tagTrigger", edm::InputTag("tbunpacker"));
  _taguMN = ps.getUntrackedParameter<edm::InputTag>("taguMN", edm::InputTag("hcalDigis"));
  _tokQIE11 = consumes<QIE11DigiCollection>(_tagQIE11);
  _tokHO = consumes<HODigiCollection>(_tagHO);
  _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
  _tokTrigger = consumes<HcalTBTriggerData>(_tagTrigger);
  _tokuMN = consumes<HcalUMNioDigi>(_taguMN);

  //	constants
  _lowHBHE = ps.getUntrackedParameter<double>("lowHBHE", 20);
  _lowHO = ps.getUntrackedParameter<double>("lowHO", 20);
  _lowHF = ps.getUntrackedParameter<double>("lowHF", 20);
}

/* virtual */ void LEDTask::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);

  edm::ESHandle<HcalDbService> dbService = es.getHandle(hcalDbServiceToken_);
  _emap = dbService->getHcalMapping();

  // Book LED calibration channels from emap
  std::vector<HcalElectronicsId> eids = _emap->allElectronicsId();
  for (unsigned i = 0; i < eids.size(); i++) {
    HcalElectronicsId eid = eids[i];
    DetId id = _emap->lookup(eid);
    if (HcalGenericDetId(id.rawId()).isHcalCalibDetId()) {
      HcalCalibDetId calibId(id);
      if (calibId.calibFlavor() == HcalCalibDetId::CalibrationBox) {
        auto cUch = calibId.cboxChannel();
        bool isLED(false);
        HcalSubdetector this_subdet = HcalEmpty;

        switch (calibId.hcalSubdet()) {
          case HcalBarrel:
            this_subdet = HcalBarrel;
            if (cUch == 0 || cUch == 1) {
              isLED = true;
            }
            break;
          case HcalEndcap:
            this_subdet = HcalEndcap;
            if (cUch == 0 || cUch == 1) {
              isLED = true;
            }
            break;
          case HcalOuter:
            this_subdet = HcalOuter;
            isLED = true;
            break;
          case HcalForward:
            this_subdet = HcalForward;
            if (cUch == 0 || cUch == 8) {
              isLED = true;
            }
            break;
          default:
            this_subdet = HcalEmpty;
            break;
        }

        if (isLED) {
          _ledCalibrationChannels[this_subdet].push_back(HcalDetId(id.rawId()));
        }
      }
    }
  }

  std::vector<uint32_t> vhashVME;
  std::vector<uint32_t> vhashuTCA;
  std::vector<uint32_t> vhashC36;
  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, hcaldqm::hashfunctions::fElectronics, vhashVME);
  _filter_uTCA.initialize(filter::fFilter, hcaldqm::hashfunctions::fElectronics, vhashuTCA);

  //	INITIALIZE
  _cSignalMean_Subdet.initialize(_name,
                                 "SignalMean",
                                 hcaldqm::hashfunctions::fSubdet,
                                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
                                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                 0);
  _cSignalRMS_Subdet.initialize(_name,
                                "SignalRMS",
                                hcaldqm::hashfunctions::fSubdet,
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_3000),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                0);
  _cTimingMean_Subdet.initialize(_name,
                                 "TimingMean",
                                 hcaldqm::hashfunctions::fSubdet,
                                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                                 new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                 0);
  _cTimingRMS_Subdet.initialize(_name,
                                "TimingRMS",
                                hcaldqm::hashfunctions::fSubdet,
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                0);

  if (_ptype == fLocal) {  // hidefed2crate
    _cSignalMean_FEDuTCA.initialize(_name,
                                    "SignalMean",
                                    hcaldqm::hashfunctions::fFED,
                                    new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                    new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),
                                    0);
    _cSignalRMS_FEDuTCA.initialize(_name,
                                   "SignalRMS",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_10000, true),
                                   0);
    _cTimingMean_FEDuTCA.initialize(_name,
                                    "TimingMean",
                                    hcaldqm::hashfunctions::fFED,
                                    new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                    new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                                    0);
    _cTimingRMS_FEDuTCA.initialize(_name,
                                   "TimingRMS",
                                   hcaldqm::hashfunctions::fFED,
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                   new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fFiberuTCAFiberCh),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                                   0);

    _cShapeCut_FEDSlot.initialize(_name,
                                  "Shape",
                                  hcaldqm::hashfunctions::fFEDSlot,
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000),
                                  0);
  }

  _cSignalMean_depth.initialize(_name,
                                "SignalMean",
                                hcaldqm::hashfunctions::fdepth,
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000, true),
                                0);
  _cSignalRMS_depth.initialize(_name,
                               "SignalRMS",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_10000, true),
                               0);
  _cTimingMean_depth.initialize(_name,
                                "TimingMean",
                                hcaldqm::hashfunctions::fdepth,
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                                0);
  _cTimingRMS_depth.initialize(_name,
                               "TimingRMS",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                               0);
  _cChargeCorr_depth.initialize(_name,
                                "ChargeCorrected",
                                hcaldqm::hashfunctions::fdepth,
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000),
                                0);
  _cCWTCorr_depth.initialize(_name,
                             "CWTCorrected",
                             hcaldqm::hashfunctions::fdepth,
                             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                             new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                             0);
  _cChargeNormLS2_depth.initialize(_name,
                                   "ChargeNormLS2",
                                   hcaldqm::hashfunctions::fdepth,
                                   new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                   new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2),
                                   0);
  _cCWTNormLS2_depth.initialize(_name,
                                "CWTNormLS2",
                                hcaldqm::hashfunctions::fdepth,
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTimingDiff_ns),
                                0);
  _cChargeAbs_depth.initialize(_name,
                               "ChargeAbs",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::ffC_generic_400000),
                               0);
  _cCWTAbs_depth.initialize(_name,
                            "CWTAbs",
                            hcaldqm::hashfunctions::fdepth,
                            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                            new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS200),
                            0);

  _cMissing_depth.initialize(_name,
                             "Missing",
                             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 (_ptype != fOffline) {  // hidefed2crate
    _cMissing_FEDuTCA.initialize(_name,
                                 "Missing",
                                 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);
  }
  _cChargeNormvsLS_SubdetPM.initialize(_name,
                                       "ChargeNormvsLS",
                                       hcaldqm::hashfunctions::fSubdetPM,
                                       new hcaldqm::quantity::LumiSection(_maxLS),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fRatio_0to2),
                                       new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                       0);
  _cCWTNormvsLS_SubdetPM.initialize(_name,
                                    "CWTNormvsLS",
                                    hcaldqm::hashfunctions::fSubdetPM,
                                    new hcaldqm::quantity::LumiSection(_maxLS),
                                    new FineTimingDiffNs(),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                    0);

  // Plots for LED in global
  if (_ptype == fOnline || _ptype == fLocal) {
    _cADCvsTS_SubdetPM.initialize(_name,
                                  "ADCvsTS",
                                  hcaldqm::hashfunctions::fSubdetPM,
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                  0);
  }
  if (_ptype == fOnline) {
    _cLowSignal_CrateSlot.initialize(_name,
                                     "LowSignal",
                                     new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fCrateuTCA),
                                     new hcaldqm::quantity::ElectronicsQuantity(hcaldqm::quantity::fSlotuTCA),
                                     new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                     0);
    _cSumQ_SubdetPM.initialize(_name,
                               "SumQ",
                               hcaldqm::hashfunctions::fSubdetPM,
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                               0);
    _cTDCTime_SubdetPM.initialize(_name,
                                  "TDCTime",
                                  hcaldqm::hashfunctions::fSubdetPM,
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true),
                                  0);
    _cTDCTime_depth.initialize(_name,
                               "TDCTime",
                               hcaldqm::hashfunctions::fdepth,
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                               new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
                               0);
    _LED_ADCvsTS_Subdet.initialize(_name + "/CU_LED",
                                   "CU_ADCvsTS",
                                   hcaldqm::hashfunctions::fSubdet,
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                                   new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                   0);
  } else if (_ptype == fLocal) {
    _LED_ADCvsEvn_Subdet.initialize(_name + "CU_LED",
                                    "CU_ADCvsEvn",
                                    hcaldqm::hashfunctions::fSubdet,
                                    new hcaldqm::quantity::EventNumber(_nevents),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fADC_256_4),
                                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN),
                                    0);
  }

  //	initialize compact containers
  _xSignalSum.initialize(hcaldqm::hashfunctions::fDChannel);
  _xSignalSum2.initialize(hcaldqm::hashfunctions::fDChannel);
  _xTimingSum.initialize(hcaldqm::hashfunctions::fDChannel);
  _xTimingSum2.initialize(hcaldqm::hashfunctions::fDChannel);
  _xEntries.initialize(hcaldqm::hashfunctions::fDChannel);
  _xPedestalChargeSumLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xPedestalChargeEntriesLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xPedestalChargePrevLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xChargeRefLS2.initialize(hcaldqm::hashfunctions::fDChannel);
  _xCWTRefLS2.initialize(hcaldqm::hashfunctions::fDChannel);
  _xChargeRefLS2Entries.initialize(hcaldqm::hashfunctions::fDChannel);
  _xCWTRefLS2Entries.initialize(hcaldqm::hashfunctions::fDChannel);
  _xChargeNormSumLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xCWTNormSumLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xChargeNormEntriesLS.initialize(hcaldqm::hashfunctions::fDChannel);
  _xCWTNormEntriesLS.initialize(hcaldqm::hashfunctions::fDChannel);

  //	BOOK
  _cSignalMean_Subdet.book(ib, _emap, _subsystem);
  _cSignalRMS_Subdet.book(ib, _emap, _subsystem);
  _cTimingMean_Subdet.book(ib, _emap, _subsystem);
  _cTimingRMS_Subdet.book(ib, _emap, _subsystem);

  _cSignalMean_depth.book(ib, _emap, _subsystem);
  _cSignalRMS_depth.book(ib, _emap, _subsystem);
  _cTimingMean_depth.book(ib, _emap, _subsystem);
  _cTimingRMS_depth.book(ib, _emap, _subsystem);
  _cChargeCorr_depth.book(ib, _emap, _subsystem);
  _cCWTCorr_depth.book(ib, _emap, _subsystem);
  _cChargeNormLS2_depth.book(ib, _emap, _subsystem);
  _cCWTNormLS2_depth.book(ib, _emap, _subsystem);
  _cChargeAbs_depth.book(ib, _emap, _subsystem);
  _cCWTAbs_depth.book(ib, _emap, _subsystem);
  _cChargeNormvsLS_SubdetPM.book(ib, _emap, _subsystem);
  _cCWTNormvsLS_SubdetPM.book(ib, _emap, _subsystem);

  _cMissing_depth.book(ib, _emap, _subsystem);
  if (_ptype == fLocal) {  // hidefed2crate
    _cSignalMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cSignalRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cTimingMean_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cTimingRMS_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
    _cShapeCut_FEDSlot.book(ib, _emap, _subsystem);
    _cMissing_FEDuTCA.book(ib, _emap, _filter_VME, _subsystem);
  }
  if (_ptype == fOnline || _ptype == fLocal) {
    _cADCvsTS_SubdetPM.book(ib, _emap, _subsystem);
  }
  if (_ptype == fOnline) {
    _cLowSignal_CrateSlot.book(ib, _subsystem);
    _cSumQ_SubdetPM.book(ib, _emap, _subsystem);
    _cTDCTime_SubdetPM.book(ib, _emap, _subsystem);
    _cTDCTime_depth.book(ib, _emap, _subsystem);
  }

  if (_ptype == fOnline) {
    _LED_ADCvsTS_Subdet.book(ib, _emap, _subsystem);
  } else if (_ptype == fLocal) {
    _LED_ADCvsEvn_Subdet.book(ib, _emap, _subsystem);
  }

  _xSignalSum.book(_emap);
  _xSignalSum2.book(_emap);
  _xEntries.book(_emap);
  _xTimingSum.book(_emap);
  _xTimingSum2.book(_emap);
  _xPedestalChargeSumLS.book(_emap);
  _xPedestalChargeEntriesLS.book(_emap);
  _xPedestalChargePrevLS.book(_emap);
  _xChargeRefLS2.book(_emap);
  _xCWTRefLS2.book(_emap);
  _xChargeRefLS2Entries.book(_emap);
  _xCWTRefLS2Entries.book(_emap);
  _xChargeNormSumLS.book(_emap);
  _xCWTNormSumLS.book(_emap);
  _xChargeNormEntriesLS.book(_emap);
  _xCWTNormEntriesLS.book(_emap);

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

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

/* virtual */ void LEDTask::_dump() {
  _cSignalMean_Subdet.reset();
  _cSignalRMS_Subdet.reset();
  _cTimingMean_Subdet.reset();
  _cTimingRMS_Subdet.reset();
  _cSignalMean_depth.reset();
  _cSignalRMS_depth.reset();
  _cTimingMean_depth.reset();
  _cTimingRMS_depth.reset();

  if (_ptype == fLocal) {  // hidefed2crate
    _cSignalMean_FEDuTCA.reset();
    _cSignalRMS_FEDuTCA.reset();
    _cTimingMean_FEDuTCA.reset();
    _cTimingRMS_FEDuTCA.reset();
  }

  std::vector<HcalGenericDetId> dids = _emap->allPrecisionId();
  for (std::vector<HcalGenericDetId>::const_iterator it = dids.begin(); it != dids.end(); ++it) {
    if (!it->isHcalDetId())
      continue;
    HcalDetId did = HcalDetId(it->rawId());
    HcalElectronicsId eid(_ehashmap.lookup(did));
    int n = _xEntries.get(did);
    double msig = _xSignalSum.get(did) / n;
    double mtim = _xTimingSum.get(did) / n;
    double rsig = sqrt(_xSignalSum2.get(did) / n - msig * msig);
    double rtim = sqrt(_xTimingSum2.get(did) / n - mtim * mtim);

    //	channels missing or low signal
    if (n == 0) {
      _cMissing_depth.fill(did);
      if (_ptype == fLocal) {  // hidefed2crate
        if (!eid.isVMEid())
          _cMissing_FEDuTCA.fill(eid);
      }
      continue;
    }
    _cSignalMean_Subdet.fill(did, msig);
    _cSignalMean_depth.fill(did, msig);
    _cSignalRMS_Subdet.fill(did, rsig);
    _cSignalRMS_depth.fill(did, rsig);
    _cTimingMean_Subdet.fill(did, mtim);
    _cTimingMean_depth.fill(did, mtim);
    _cTimingRMS_Subdet.fill(did, rtim);
    _cTimingRMS_depth.fill(did, rtim);
    if (_ptype == fLocal) {  // hidefed2crate
      if (!eid.isVMEid()) {
        _cSignalMean_FEDuTCA.fill(eid, msig);
        _cSignalRMS_FEDuTCA.fill(eid, rsig);
        _cTimingMean_FEDuTCA.fill(eid, mtim);
        _cTimingRMS_FEDuTCA.fill(eid, rtim);
      }
    }
  }
}

/* virtual */ void LEDTask::globalEndLuminosityBlock(edm::LuminosityBlock const& lb, edm::EventSetup const& es) {
  auto lumiCache = luminosityBlockCache(lb.index());
  _currentLS = lumiCache->currentLS;

  if (_emap) {
    std::vector<HcalGenericDetId> dids = _emap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = dids.begin(); it != dids.end(); ++it) {
      if (!it->isHcalDetId())
        continue;
      HcalDetId did(it->rawId());
      int nChargeNorm = _xChargeNormEntriesLS.get(did);
      if (nChargeNorm > 0) {
        _cChargeNormvsLS_SubdetPM.fill(did, _currentLS, _xChargeNormSumLS.get(did) / nChargeNorm);
      }
      int nCWTNorm = _xCWTNormEntriesLS.get(did);
      if (nCWTNorm > 0) {
        _cCWTNormvsLS_SubdetPM.fill(did, _currentLS, _xCWTNormSumLS.get(did) / nCWTNorm);
      }
      int nPed = _xPedestalChargeEntriesLS.get(did);
      if (nPed > 0) {
        _xPedestalChargePrevLS.set(did, _xPedestalChargeSumLS.get(did) / nPed);
      }
    }
  }
  _xChargeNormSumLS.reset();
  _xCWTNormSumLS.reset();
  _xChargeNormEntriesLS.reset();
  _xCWTNormEntriesLS.reset();
  _xPedestalChargeSumLS.reset();
  _xPedestalChargeEntriesLS.reset();

  DQTask::globalEndLuminosityBlock(lb, es);
}

/* virtual */ void LEDTask::_process(edm::Event const& e, edm::EventSetup const& es) {
  edm::Handle<HODigiCollection> c_ho;
  edm::Handle<QIE10DigiCollection> c_QIE10;
  edm::Handle<QIE11DigiCollection> c_QIE11;

  if (!e.getByToken(_tokHO, c_ho))
    _logger.dqmthrow("Collection HODigiCollection isn't available " + _tagHO.label() + " " + _tagHO.instance());
  if (!e.getByToken(_tokQIE10, c_QIE10))
    _logger.dqmthrow("Collection QIE10DigiCollection isn't available " + _tagQIE10.label() + " " +
                     _tagQIE10.instance());
  if (!e.getByToken(_tokQIE11, c_QIE11))
    _logger.dqmthrow("Collection QIE11DigiCollection isn't available " + _tagQIE11.label() + " " +
                     _tagQIE11.instance());

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

  for (QIE11DigiCollection::const_iterator it = c_QIE11->begin(); it != c_QIE11->end(); ++it) {
    const QIE11DataFrame digi = static_cast<const QIE11DataFrame>(*it);
    HcalDetId const& did = digi.detid();
    HcalCalibDetId hcdid(digi.id());
    if ((did.subdet() != HcalBarrel) && (did.subdet() != HcalEndcap)) {
      // LED monitoring from calibration channels
      if (did.subdet() == HcalOther) {
        HcalOtherDetId hodid(digi.detid());
        if (hodid.subdet() == HcalCalibration) {
          if (std::find(_ledCalibrationChannels[HcalEndcap].begin(), _ledCalibrationChannels[HcalEndcap].end(), did) !=
              _ledCalibrationChannels[HcalEndcap].end()) {
            for (int i = 0; i < digi.samples(); i++) {
              if (_ptype == fOnline) {
                _LED_ADCvsTS_Subdet.fill(HcalDetId(HcalEndcap, 16, 1, 1), i, digi[i].adc());
              } else if (_ptype == fLocal) {
                _LED_ADCvsEvn_Subdet.fill(
                    HcalDetId(HcalEndcap, 16, 1, 1), e.eventAuxiliary().id().event(), digi[i].adc());
              }
            }
          } else if (std::find(_ledCalibrationChannels[HcalBarrel].begin(),
                               _ledCalibrationChannels[HcalBarrel].end(),
                               did) != _ledCalibrationChannels[HcalBarrel].end()) {
            if (hcdid.rawId() == constants::HBLasMon.rawId())
              continue;
            for (int i = 0; i < digi.samples(); i++) {
              if (_ptype == fOnline) {
                _LED_ADCvsTS_Subdet.fill(HcalDetId(HcalBarrel, 1, 1, 1), i, digi[i].adc());
              } else if (_ptype == fLocal) {
                _LED_ADCvsEvn_Subdet.fill(
                    HcalDetId(HcalBarrel, 1, 1, 1), e.eventAuxiliary().id().event(), digi[i].adc());
              }
            }
          }
        }
      }
      continue;
    }
    uint32_t rawid = _ehashmap.lookup(did);
    if (!rawid) {
      std::string unknown_id_string = "Detid " + std::to_string(int(did)) + ", ieta " + std::to_string(did.ieta());
      unknown_id_string += ", iphi " + std::to_string(did.iphi()) + ", depth " + std::to_string(did.depth());
      unknown_id_string += ", is not in emap. Skipping.";
      _logger.warn(unknown_id_string);
      continue;
    }
    HcalElectronicsId const& eid(rawid);

    CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE11DataFrame>(_dbService, did, digi);
    //double sumQ = hcaldqm::utilities::sumQ_v10<QIE11DataFrame>(digi, 2.5, 0, digi.samples()-1);
    double sumQ = hcaldqm::utilities::sumQDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
    int const nSamples = digi.samples();
    if (nSamples <= 0)
      continue;
    double rawCharge = 0;
    double pedestalChargeEvent = 0;
    for (int i = 0; i < nSamples; ++i) {
      rawCharge += digi_fC[i];
    }
    int nPedestalTS = 1;
    for (int i = 0; i < nPedestalTS; ++i)
      pedestalChargeEvent += digi_fC[i];
    pedestalChargeEvent = (nPedestalTS > 0 ? pedestalChargeEvent * nSamples / nPedestalTS : 0);
    _xPedestalChargeSumLS.get(did) += pedestalChargeEvent;
    _xPedestalChargeEntriesLS.get(did)++;

    if (sumQ >= _lowHBHE) {
      //double aveTS = hcaldqm::utilities::aveTS_v10<QIE11DataFrame>(digi, 2.5, 0,digi.samples()-1);
      double aveTS = hcaldqm::utilities::aveTSDB<QIE11DataFrame>(_dbService, digi_fC, did, digi, 0, digi.size() - 1);

      _xSignalSum.get(did) += sumQ;
      _xSignalSum2.get(did) += sumQ * sumQ;
      _xTimingSum.get(did) += aveTS;
      _xTimingSum2.get(did) += aveTS * aveTS;
      _xEntries.get(did)++;

      double pedestalPrevLS = _xPedestalChargePrevLS.get(did);
      if (pedestalPrevLS <= 0)
        pedestalPrevLS = pedestalChargeEvent;
      double pedestalPerTS = pedestalPrevLS / nSamples;
      double chargeAbs = rawCharge - pedestalPrevLS;
      double cwtDenom = 0;
      double cwtNumer = 0;
      for (int i = 0; i < nSamples; ++i) {
        double const qCorr = digi_fC[i] - pedestalPerTS;
        cwtDenom += qCorr;
        cwtNumer += i * qCorr;
      }
      double cwtCorr = (cwtDenom > 0 ? cwtNumer / cwtDenom : constants::GARBAGE_VALUE);

      _cChargeCorr_depth.fill(did, chargeAbs);
      _cChargeAbs_depth.fill(did, chargeAbs);
      if (cwtCorr != constants::GARBAGE_VALUE) {
        _cCWTCorr_depth.fill(did, cwtCorr);
        _cCWTAbs_depth.fill(did, cwtCorr);
      }

      if (_currentLS == 2) {
        int nChargeRef = _xChargeRefLS2Entries.get(did);
        double chargeRef = _xChargeRefLS2.get(did);
        _xChargeRefLS2.set(did, (chargeRef * nChargeRef + chargeAbs) / (nChargeRef + 1));
        _xChargeRefLS2Entries.get(did)++;
        if (cwtCorr != constants::GARBAGE_VALUE) {
          int nCWTRef = _xCWTRefLS2Entries.get(did);
          double cwtRef = _xCWTRefLS2.get(did);
          _xCWTRefLS2.set(did, (cwtRef * nCWTRef + cwtCorr) / (nCWTRef + 1));
          _xCWTRefLS2Entries.get(did)++;
        }
      }

      if (_xChargeRefLS2Entries.get(did) > 0) {
        double chargeRef = _xChargeRefLS2.get(did);
        if (std::abs(chargeRef) > kTiny) {
          double chargeNorm = chargeAbs / chargeRef;
          _cChargeNormLS2_depth.fill(did, chargeNorm);
          _xChargeNormSumLS.get(did) += chargeNorm;
          _xChargeNormEntriesLS.get(did)++;
        }
      }
      if (cwtCorr != constants::GARBAGE_VALUE && _xCWTRefLS2Entries.get(did) > 0) {
        double cwtNormNS = (cwtCorr - _xCWTRefLS2.get(did)) * 25.;
        _cCWTNormLS2_depth.fill(did, cwtNormNS);
        _xCWTNormSumLS.get(did) += cwtNormNS;
        _xCWTNormEntriesLS.get(did)++;
      }

      if (_ptype == fLocal) {  // hidefed2crate
        for (int i = 0; i < digi.samples(); i++) {
          //_cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC()-2.5);
          _cShapeCut_FEDSlot.fill(
              eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE11DataFrame>(_dbService, digi_fC, did, digi, i));
        }
      }
      if (_ptype == fOnline || _ptype == fLocal) {
        for (int iTS = 0; iTS < digi.samples(); ++iTS) {
          _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc());
        }
      }
      if (_ptype == fOnline) {
        for (int iTS = 0; iTS < digi.samples(); ++iTS) {
          if (digi[iTS].tdc() < 50) {
            double time = iTS * 25. + (digi[iTS].tdc() / 2.);
            _cTDCTime_SubdetPM.fill(did, time);
            _cTDCTime_depth.fill(did, time);
          }
        }
        _cSumQ_SubdetPM.fill(did, sumQ);

        // Low signal in SOI
        short soi = -1;
        for (int i = 0; i < digi.samples(); i++) {
          if (digi[i].soi()) {
            soi = i;
            break;
          }
        }
        if (digi[soi].adc() < 30) {
          _cLowSignal_CrateSlot.fill(eid);
        }
      }
    }
  }
  for (HODigiCollection::const_iterator it = c_ho->begin(); it != c_ho->end(); ++it) {
    const HODataFrame digi = (const HODataFrame)(*it);
    HcalDetId did = digi.id();
    if (did.subdet() != HcalOuter) {
      // LED monitoring from calibration channels
      if (did.subdet() == HcalOther) {
        HcalOtherDetId hodid(did);
        if (hodid.subdet() == HcalCalibration) {
          if (std::find(_ledCalibrationChannels[HcalOuter].begin(), _ledCalibrationChannels[HcalOuter].end(), did) !=
              _ledCalibrationChannels[HcalOuter].end()) {
            for (int i = 0; i < digi.size(); i++) {
              if (_ptype == fOnline) {
                _LED_ADCvsTS_Subdet.fill(HcalDetId(HcalOuter, 1, 1, 4), i, digi[i].adc());
              } else if (_ptype == fLocal) {
                _LED_ADCvsEvn_Subdet.fill(
                    HcalDetId(HcalOuter, 1, 1, 4), e.eventAuxiliary().id().event(), digi[i].adc());
              }
            }
          }
        }
      }
      continue;
    }
    HcalElectronicsId eid = digi.elecId();
    //double sumQ = hcaldqm::utilities::sumQ<HODataFrame>(digi, 8.5, 0, digi.size()-1);
    CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<HODataFrame>(_dbService, did, digi);
    double sumQ = hcaldqm::utilities::sumQDB<HODataFrame>(_dbService, digi_fC, did, digi, 0, digi.size() - 1);
    int const nSamples = digi.size();
    if (nSamples <= 0)
      continue;
    double rawCharge = 0;
    double pedestalChargeEvent = 0;
    for (int i = 0; i < nSamples; ++i) {
      rawCharge += digi_fC[i];
    }
    int nPedestalTS = 1;
    for (int i = 0; i < nPedestalTS; ++i)
      pedestalChargeEvent += digi_fC[i];
    pedestalChargeEvent = (nPedestalTS > 0 ? pedestalChargeEvent * nSamples / nPedestalTS : 0);
    _xPedestalChargeSumLS.get(did) += pedestalChargeEvent;
    _xPedestalChargeEntriesLS.get(did)++;

    if (sumQ >= _lowHO) {
      //double aveTS = hcaldqm::utilities::aveTS<HODataFrame>(digi, 8.5, 0, digi.size()-1);
      double aveTS = hcaldqm::utilities::aveTSDB<HODataFrame>(_dbService, digi_fC, did, digi, 0, digi.size() - 1);

      _xSignalSum.get(did) += sumQ;
      _xSignalSum2.get(did) += sumQ * sumQ;
      _xTimingSum.get(did) += aveTS;
      _xTimingSum2.get(did) += aveTS * aveTS;
      _xEntries.get(did)++;

      double pedestalPrevLS = _xPedestalChargePrevLS.get(did);
      if (pedestalPrevLS <= 0)
        pedestalPrevLS = pedestalChargeEvent;
      double pedestalPerTS = pedestalPrevLS / nSamples;
      double chargeAbs = rawCharge - pedestalPrevLS;
      double cwtDenom = 0;
      double cwtNumer = 0;
      for (int i = 0; i < nSamples; ++i) {
        double const qCorr = digi_fC[i] - pedestalPerTS;
        cwtDenom += qCorr;
        cwtNumer += i * qCorr;
      }
      double cwtCorr = (cwtDenom > 0 ? cwtNumer / cwtDenom : constants::GARBAGE_VALUE);

      _cChargeCorr_depth.fill(did, chargeAbs);
      _cChargeAbs_depth.fill(did, chargeAbs);
      if (cwtCorr != constants::GARBAGE_VALUE) {
        _cCWTCorr_depth.fill(did, cwtCorr);
        _cCWTAbs_depth.fill(did, cwtCorr);
      }

      if (_currentLS == 2) {
        int nChargeRef = _xChargeRefLS2Entries.get(did);
        double chargeRef = _xChargeRefLS2.get(did);
        _xChargeRefLS2.set(did, (chargeRef * nChargeRef + chargeAbs) / (nChargeRef + 1));
        _xChargeRefLS2Entries.get(did)++;
        if (cwtCorr != constants::GARBAGE_VALUE) {
          int nCWTRef = _xCWTRefLS2Entries.get(did);
          double cwtRef = _xCWTRefLS2.get(did);
          _xCWTRefLS2.set(did, (cwtRef * nCWTRef + cwtCorr) / (nCWTRef + 1));
          _xCWTRefLS2Entries.get(did)++;
        }
      }

      if (_xChargeRefLS2Entries.get(did) > 0) {
        double chargeRef = _xChargeRefLS2.get(did);
        if (std::abs(chargeRef) > kTiny) {
          double chargeNorm = chargeAbs / chargeRef;
          _cChargeNormLS2_depth.fill(did, chargeNorm);
          _xChargeNormSumLS.get(did) += chargeNorm;
          _xChargeNormEntriesLS.get(did)++;
        }
      }
      if (cwtCorr != constants::GARBAGE_VALUE && _xCWTRefLS2Entries.get(did) > 0) {
        double cwtNormNS = (cwtCorr - _xCWTRefLS2.get(did)) * 25.;
        _cCWTNormLS2_depth.fill(did, cwtNormNS);
        _xCWTNormSumLS.get(did) += cwtNormNS;
        _xCWTNormEntriesLS.get(did)++;
      }

      if (_ptype == fLocal) {  // hidefed2crate
        for (int i = 0; i < digi.size(); i++) {
          //_cShapeCut_FEDSlot.fill(eid, i, digi.sample(i).nominal_fC()-8.5);
          _cShapeCut_FEDSlot.fill(
              eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<HODataFrame>(_dbService, digi_fC, did, digi, i));
        }
      }
      if (_ptype == fOnline || _ptype == fLocal) {
        for (int iTS = 0; iTS < digi.size(); ++iTS) {
          _cADCvsTS_SubdetPM.fill(did, iTS, digi.sample(iTS).adc());
        }
      }
      if (_ptype == fOnline) {
        _cSumQ_SubdetPM.fill(did, sumQ);
      }
    }
  }

  for (QIE10DigiCollection::const_iterator it = c_QIE10->begin(); it != c_QIE10->end(); ++it) {
    const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
    HcalDetId did = digi.detid();
    if (did.subdet() != HcalForward) {
      // LED monitoring from calibration channels
      if (did.subdet() == HcalOther) {
        HcalOtherDetId hodid(digi.detid());
        if (hodid.subdet() == HcalCalibration) {
          if (std::find(_ledCalibrationChannels[HcalForward].begin(),
                        _ledCalibrationChannels[HcalForward].end(),
                        did) != _ledCalibrationChannels[HcalForward].end()) {
            for (int i = 0; i < digi.samples(); i++) {
              if (_ptype == fOnline) {
                _LED_ADCvsTS_Subdet.fill(HcalDetId(HcalForward, 29, 1, 1), i, digi[i].adc());
              } else if (_ptype == fLocal) {
                _LED_ADCvsEvn_Subdet.fill(
                    HcalDetId(HcalForward, 29, 1, 1), e.eventAuxiliary().id().event(), digi[i].adc());
              }
            }
          }
        }
      }
      continue;
    }
    HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
    //double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
    CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, digi);
    double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.samples() - 1);
    int const nSamples = digi.samples();
    if (nSamples <= 0)
      continue;
    double rawCharge = 0;
    double pedestalChargeEvent = 0;
    for (int i = 0; i < nSamples; ++i) {
      rawCharge += digi_fC[i];
    }
    int nPedestalTS = 1;
    for (int i = 0; i < nPedestalTS; ++i)
      pedestalChargeEvent += digi_fC[i];
    pedestalChargeEvent = (nPedestalTS > 0 ? pedestalChargeEvent * nSamples / nPedestalTS : 0);
    _xPedestalChargeSumLS.get(did) += pedestalChargeEvent;
    _xPedestalChargeEntriesLS.get(did)++;

    if (sumQ >= _lowHF) {
      //double aveTS = hcaldqm::utilities::aveTS_v10<QIE10DataFrame>(digi, 2.5, 0, digi.samples()-1);
      double aveTS = hcaldqm::utilities::aveTSDB<QIE10DataFrame>(_dbService, digi_fC, did, digi, 0, digi.size() - 1);

      _xSignalSum.get(did) += sumQ;
      _xSignalSum2.get(did) += sumQ * sumQ;
      _xTimingSum.get(did) += aveTS;
      _xTimingSum2.get(did) += aveTS * aveTS;
      _xEntries.get(did)++;

      double pedestalPrevLS = _xPedestalChargePrevLS.get(did);
      if (pedestalPrevLS <= 0)
        pedestalPrevLS = pedestalChargeEvent;
      double pedestalPerTS = pedestalPrevLS / nSamples;
      double chargeAbs = rawCharge - pedestalPrevLS;
      double cwtDenom = 0;
      double cwtNumer = 0;
      for (int i = 0; i < nSamples; ++i) {
        double const qCorr = digi_fC[i] - pedestalPerTS;
        cwtDenom += qCorr;
        cwtNumer += i * qCorr;
      }
      double cwtCorr = (cwtDenom > 0 ? cwtNumer / cwtDenom : constants::GARBAGE_VALUE);

      _cChargeCorr_depth.fill(did, chargeAbs);
      _cChargeAbs_depth.fill(did, chargeAbs);
      if (cwtCorr != constants::GARBAGE_VALUE) {
        _cCWTCorr_depth.fill(did, cwtCorr);
        _cCWTAbs_depth.fill(did, cwtCorr);
      }

      if (_currentLS == 2) {
        int nChargeRef = _xChargeRefLS2Entries.get(did);
        double chargeRef = _xChargeRefLS2.get(did);
        _xChargeRefLS2.set(did, (chargeRef * nChargeRef + chargeAbs) / (nChargeRef + 1));
        _xChargeRefLS2Entries.get(did)++;
        if (cwtCorr != constants::GARBAGE_VALUE) {
          int nCWTRef = _xCWTRefLS2Entries.get(did);
          double cwtRef = _xCWTRefLS2.get(did);
          _xCWTRefLS2.set(did, (cwtRef * nCWTRef + cwtCorr) / (nCWTRef + 1));
          _xCWTRefLS2Entries.get(did)++;
        }
      }

      if (_xChargeRefLS2Entries.get(did) > 0) {
        double chargeRef = _xChargeRefLS2.get(did);
        if (std::abs(chargeRef) > kTiny) {
          double chargeNorm = chargeAbs / chargeRef;
          _cChargeNormLS2_depth.fill(did, chargeNorm);
          _xChargeNormSumLS.get(did) += chargeNorm;
          _xChargeNormEntriesLS.get(did)++;
        }
      }
      if (cwtCorr != constants::GARBAGE_VALUE && _xCWTRefLS2Entries.get(did) > 0) {
        double cwtNormNS = (cwtCorr - _xCWTRefLS2.get(did)) * 25.;
        _cCWTNormLS2_depth.fill(did, cwtNormNS);
        _xCWTNormSumLS.get(did) += cwtNormNS;
        _xCWTNormEntriesLS.get(did)++;
      }

      if (_ptype == fLocal) {  // hidefed2crate
        for (int i = 0; i < digi.samples(); ++i) {
          // Note: this used to be digi.sample(i).nominal_fC() - 2.5, but this branch doesn't exist in QIE10DataFrame.
          // Instead, use lookup table.
          //_cShapeCut_FEDSlot.fill(eid, i, constants::adc2fC[digi[i].adc()]);
          _cShapeCut_FEDSlot.fill(
              eid, i, hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, digi, i));
        }
      }
      if (_ptype == fOnline || _ptype == fLocal) {
        for (int iTS = 0; iTS < digi.samples(); ++iTS) {
          _cADCvsTS_SubdetPM.fill(did, iTS, digi[iTS].adc());
        }
      }
      if (_ptype == fOnline) {
        for (int iTS = 0; iTS < digi.samples(); ++iTS) {
          if (digi[iTS].le_tdc() < 50) {
            double time = iTS * 25. + (digi[iTS].le_tdc() / 2.);
            _cTDCTime_SubdetPM.fill(did, time);
            _cTDCTime_depth.fill(did, time);
          }
        }
        _cSumQ_SubdetPM.fill(did, sumQ);
      }
    }
  }

  if (_ptype == fOnline && _evsTotal > 0 && _evsTotal % constants::CALIBEVENTS_MIN == 0)
    this->_dump();
}

/* virtual */ bool LEDTask::_isApplicable(edm::Event const& e) {
  if (_ptype != fOnline) {
    //	local
    edm::Handle<HcalTBTriggerData> ctrigger;
    if (!e.getByToken(_tokTrigger, ctrigger))
      _logger.dqmthrow("Collection HcalTBTriggerData isn't available " + _tagTrigger.label() + " " +
                       _tagTrigger.instance());
    return ctrigger->wasLEDTrigger();
  } else {
    //	fOnline mode
    edm::Handle<HcalUMNioDigi> cumn;
    if (!e.getByToken(_tokuMN, cumn)) {
      return false;
    }

    return (cumn->eventType() == constants::EVENTTYPE_LED);
  }

  return false;
}

DEFINE_FWK_MODULE(LEDTask);