ClusterFillMap.cc

Go to the documentation of this file.

#include "Calibration/EcalCalibAlgos/interface/ClusterFillMap.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "FWCore/Utilities/interface/isFinite.h"
 
ClusterFillMap::ClusterFillMap(int WindowX,
                               int WindowY,
                               const std::map<int, int>& xtalReg,
                               double minE,
                               double maxE,
                               const std::map<int, int>& IndexReg,
                               EcalIntercalibConstantMap* barrelMap,
                               EcalIntercalibConstantMap* endcapMap)
    : VFillMap(WindowX, WindowY, xtalReg, minE, maxE, IndexReg, barrelMap, endcapMap)
 
{}
 
ClusterFillMap::~ClusterFillMap() {}
 
void ClusterFillMap::fillMap(const std::vector<std::pair<DetId, float> >& v1,
                             const DetId Max,
                             const EcalRecHitCollection* barrelHitsCollection,
                             const EcalRecHitCollection* endcapHitsCollection,
                             std::map<int, double>& xtlMap,
                             double& pSubtract) {
  for (std::vector<std::pair<DetId, float> >::const_iterator idsIt = v1.begin(); idsIt != v1.end(); ++idsIt) {
    int RegionNumber = m_xtalRegionId[Max.rawId()];
    EcalRecHitCollection::const_iterator itrechit;
    double dummy = 0.;
    if (idsIt->first.subdetId() == EcalBarrel) {
      itrechit = barrelHitsCollection->find(idsIt->first);
      dummy = itrechit->energy();
      dummy *= (*m_barrelMap)[idsIt->first];
    } else if (idsIt->first.subdetId() == EcalEndcap) {
      itrechit = endcapHitsCollection->find(idsIt->first);
      dummy = itrechit->energy();
      dummy *= (*m_endcapMap)[idsIt->first];
    }
    if (edm::isNotFinite(dummy)) {
      dummy = 0;
    }
    if (dummy < m_minEnergyPerCrystal)
      continue;  //return 1;
    if (dummy > m_maxEnergyPerCrystal) {
      continue;
    }
    int ID = idsIt->first.rawId();
    if (m_xtalRegionId[ID] == RegionNumber)
      xtlMap[m_IndexInRegion[ID]] += dummy;
    else
      pSubtract += dummy;
  }
}