RecoSummaryTask.cc

Go to the documentation of this file.

#include "DQM/EcalMonitorTasks/interface/RecoSummaryTask.h"

#include "DQM/EcalCommon/interface/EcalDQMCommonUtils.h"
#include "DataFormats/EcalRawData/interface/EcalDCCHeaderBlock.h"
#include "DataFormats/EcalDetId/interface/EcalDetIdCollections.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"

#include "RecoEcal/EgammaCoreTools/interface/EcalTools.h"
#include "RecoEcal/EgammaCoreTools/interface/EcalClusterTools.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

namespace ecaldqm {
  RecoSummaryTask::RecoSummaryTask()
      : DQWorkerTask(), rechitThresholdEB_(0.), rechitThresholdEE_(0.), ebHits_(nullptr), eeHits_(nullptr) {}

  void RecoSummaryTask::setParams(edm::ParameterSet const& _params) {
    rechitThresholdEB_ = _params.getUntrackedParameter<double>("rechitThresholdEB");
    rechitThresholdEE_ = _params.getUntrackedParameter<double>("rechitThresholdEE");
  }

  void RecoSummaryTask::addDependencies(DependencySet& _dependencies) {
    _dependencies.push_back(Dependency(kEBBasicCluster, kEBRecHit));
    _dependencies.push_back(Dependency(kEEBasicCluster, kEERecHit));
  }

  bool RecoSummaryTask::filterRunType(short const* _runType) {
    for (unsigned iFED(0); iFED != ecaldqm::nDCC; iFED++) {
      if (_runType[iFED] == EcalDCCHeaderBlock::COSMIC || _runType[iFED] == EcalDCCHeaderBlock::MTCC ||
          _runType[iFED] == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
          _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_GLOBAL || _runType[iFED] == EcalDCCHeaderBlock::COSMICS_LOCAL ||
          _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_LOCAL)
        return true;
    }

    return false;
  }

  void RecoSummaryTask::endEvent(edm::Event const&, edm::EventSetup const&) {
    ebHits_ = nullptr;
    eeHits_ = nullptr;
  }

  void RecoSummaryTask::runOnRecHits(EcalRecHitCollection const& _hits, Collections _collection) {
    bool isBarrel(_collection == kEBRecHit);

    MESet& meEnergyMax(MEs_.at("EnergyMax"));
    MESet& meChi2(MEs_.at("Chi2"));
    MESet& meTime(MEs_.at("Time"));
    MESet* meSwissCross(isBarrel ? &MEs_.at("SwissCross") : nullptr);
    MESet& meRecoFlag(MEs_.at("RecoFlagAll"));

    double maxE[2] = {-1., -1};
    int subdet(isBarrel ? EcalBarrel : EcalEndcap);

    for (EcalRecHitCollection::const_iterator hitItr(_hits.begin()); hitItr != _hits.end(); ++hitItr) {
      meRecoFlag.fill(subdet, hitItr->recoFlag());
      float energy(hitItr->energy());

      int signedSubdet;
      float rechitThreshold;

      if (isBarrel) {
        signedSubdet = EcalBarrel;
        rechitThreshold = rechitThresholdEB_;

        if (energy > 3.) {
          EBDetId ebId(hitItr->id());
          if (ebId.ieta() != 85)
            meSwissCross->fill(EcalTools::swissCross(ebId, _hits, 0.));
        }

        if (energy > maxE[0])
          maxE[0] = energy;
      } else {
        rechitThreshold = rechitThresholdEE_;

        EEDetId eeId(hitItr->id());
        if (eeId.zside() < 0) {
          signedSubdet = -EcalEndcap;
          if (energy > maxE[0])
            maxE[0] = energy;
        } else {
          signedSubdet = EcalEndcap;
          if (energy > maxE[1])
            maxE[1] = energy;
        }
      }

      if (energy > rechitThreshold) {
        meChi2.fill(signedSubdet, hitItr->chi2());
        meTime.fill(signedSubdet, hitItr->time());
      }
    }

    if (isBarrel) {
      meEnergyMax.fill(EcalBarrel, maxE[0]);

      ebHits_ = &_hits;
    } else {
      meEnergyMax.fill(-EcalEndcap, maxE[0]);
      meEnergyMax.fill(EcalEndcap, maxE[1]);

      eeHits_ = &_hits;
    }
  }

  void RecoSummaryTask::runOnReducedRecHits(EcalRecHitCollection const& _hits, Collections _collections) {
    MESet& meRecoFlag(MEs_.at("RecoFlagReduced"));

    int subdet(_collections == kEBReducedRecHit ? EcalBarrel : EcalEndcap);

    for (EcalRecHitCollection::const_iterator hitItr(_hits.begin()); hitItr != _hits.end(); ++hitItr)
      meRecoFlag.fill(subdet, hitItr->recoFlag());
  }

  void RecoSummaryTask::runOnBasicClusters(edm::View<reco::CaloCluster> const& _bcs, Collections _collection) {
    bool isBarrel(_collection == kEBBasicCluster);

    MESet& meRecoFlag(MEs_.at("RecoFlagBasicCluster"));

    EcalRecHitCollection const* hitCol(isBarrel ? ebHits_ : eeHits_);
    int subdet(isBarrel ? EcalBarrel : EcalEndcap);

    for (edm::View<reco::CaloCluster>::const_iterator bcItr(_bcs.begin()); bcItr != _bcs.end(); ++bcItr) {
      if (bcItr->caloID().detectors() != 0) {
        if (isBarrel && !bcItr->caloID().detector(reco::CaloID::DET_ECAL_BARREL))
          continue;
        if (!isBarrel && !bcItr->caloID().detector(reco::CaloID::DET_ECAL_ENDCAP))
          continue;
      }

      std::vector<std::pair<DetId, float> > const& haf(bcItr->hitsAndFractions());
      for (unsigned iH(0); iH != haf.size(); ++iH) {
        if (isBarrel && haf[iH].first.subdetId() != EcalBarrel)
          continue;
        if (!isBarrel && haf[iH].first.subdetId() != EcalEndcap)
          continue;
        EcalRecHitCollection::const_iterator hItr(hitCol->find(haf[iH].first));
        if (hItr == hitCol->end())
          continue;
        meRecoFlag.fill(subdet, hItr->recoFlag());
      }
    }
  }

  DEFINE_ECALDQM_WORKER(RecoSummaryTask);
}  // namespace ecaldqm