EgammaHadTower Class Reference

#include <EgammaHadTower.h>

Public Types
enum	HoeMode { SingleTower = 0, TowersBehindCluster = 1 }

Public Member Functions
	EgammaHadTower (const edm::EventSetup &es, HoeMode mode=SingleTower)
double	getDepth1HcalESum (const std::vector< CaloTowerDetId > &towers, CaloTowerCollection const &) const
double	getDepth2HcalESum (const std::vector< CaloTowerDetId > &towers, CaloTowerCollection const &) const
bool	hasActiveHcal (const std::vector< CaloTowerDetId > &towers) const
CaloTowerDetId	towerOf (const reco::CaloCluster &cluster) const
std::vector< CaloTowerDetId >	towersOf (const reco::SuperCluster &sc) const

Private Attributes
const HcalChannelQuality *	hcalQuality_
const HcalTopology *	hcalTopology_
HoeMode	mode_
unsigned int	nMaxClusters_
const CaloTowerConstituentsMap *	towerMap_

Detailed Description

Definition at line 17 of file EgammaHadTower.h.

Member Enumeration Documentation

HoeMode

enum EgammaHadTower::HoeMode

Enumerator
SingleTower
TowersBehindCluster

Definition at line 19 of file EgammaHadTower.h.

{ SingleTower = 0, TowersBehindCluster = 1 };

Constructor & Destructor Documentation

EgammaHadTower()

EgammaHadTower::EgammaHadTower	(	const edm::EventSetup &	es,
		HoeMode	mode = SingleTower
	)

Definition at line 20 of file EgammaHadTower.cc.

                                                                    : mode_(mode) {
  edm::ESHandle<CaloTowerConstituentsMap> ctmaph;
  es.get<CaloGeometryRecord>().get(ctmaph);
  towerMap_ = &(*ctmaph);
  nMaxClusters_ = 4;
 
  edm::ESHandle<HcalChannelQuality> hQuality;
  es.get<HcalChannelQualityRcd>().get("withTopo", hQuality);
  hcalQuality_ = hQuality.product();
 
  edm::ESHandle<HcalTopology> hcalTopology;
  es.get<HcalRecNumberingRecord>().get(hcalTopology);
  hcalTopology_ = hcalTopology.product();
}

References edm::EventSetup::get(), get, hcalQuality_, hcalTopology_, nMaxClusters_, edm::ESHandle< T >::product(), and towerMap_.

Member Function Documentation

getDepth1HcalESum()

double EgammaHadTower::getDepth1HcalESum	(	const std::vector< CaloTowerDetId > &	towers,
		CaloTowerCollection const &	towerCollection
	)		const

Definition at line 94 of file EgammaHadTower.cc.

                                                                                           {
  double esum = 0.;
  for (auto const& tower : towerCollection) {
    if (std::find(towers.begin(), towers.end(), tower.id()) != towers.end()) {
      esum += tower.ietaAbs() < 18 || tower.ietaAbs() > 29 ? tower.hadEnergy() : tower.hadEnergyHeInnerLayer();
    }
  }
  return esum;
}

References spr::find(), hgcalTowerProducer_cfi::tower, HLT_2018_cff::towerCollection, and HLT_2018_cff::towers.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth1BehindClusters(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

getDepth2HcalESum()

double EgammaHadTower::getDepth2HcalESum	(	const std::vector< CaloTowerDetId > &	towers,
		CaloTowerCollection const &	towerCollection
	)		const

Definition at line 105 of file EgammaHadTower.cc.

                                                                                           {
  double esum = 0.;
  for (auto const& tower : towerCollection) {
    if (std::find(towers.begin(), towers.end(), tower.id()) != towers.end()) {
      esum += tower.hadEnergyHeOuterLayer();
    }
  }
  return esum;
}

References spr::find(), hgcalTowerProducer_cfi::tower, HLT_2018_cff::towerCollection, and HLT_2018_cff::towers.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth2BehindClusters(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

hasActiveHcal()

bool EgammaHadTower::hasActiveHcal

(

const std::vector< CaloTowerDetId > &

towers

)

const

Definition at line 116 of file EgammaHadTower.cc.

                                                                                {
  bool active = false;
  int statusMask = ((1 << HcalChannelStatus::HcalCellOff) | (1 << HcalChannelStatus::HcalCellMask) |
                    (1 << HcalChannelStatus::HcalCellDead));
#ifdef EDM_ML_DEBUG
  std::cout << "DEBUG: hasActiveHcal called with " << towers.size() << " detids. First tower detid ieta "
            << towers.front().ieta() << " iphi " << towers.front().iphi() << std::endl;
#endif
  for (auto towerid : towers) {
    unsigned int ngood = 0, nbad = 0;
    for (DetId id : towerMap_->constituentsOf(towerid)) {
      if (id.det() != DetId::Hcal) {
        continue;
      }
      HcalDetId hid(id);
      if (hid.subdet() != HcalBarrel && hid.subdet() != HcalEndcap)
        continue;
#ifdef EDM_ML_DEBUG
      std::cout << "EgammaHadTower DetId " << std::hex << id.rawId() << "  hid.rawId  " << hid.rawId() << std::dec
                << "   sub " << hid.subdet() << "   ieta " << hid.ieta() << "   iphi " << hid.iphi() << "   depth "
                << hid.depth() << std::endl;
#endif
      // Sunanda's fix for 2017 Plan1
      // and removed protection
      int status =
          hcalQuality_->getValues((DetId)(hcalTopology_->idFront(HcalDetId(id))), /*throwOnFail=*/true)->getValue();
 
#ifdef EDM_ML_DEBUG
      std::cout << "channels status = " << std::hex << status << std::dec << "  int value = " << status << std::endl;
#endif
 
      if (status & statusMask) {
#ifdef EDM_ML_DEBUG
        std::cout << "          BAD!" << std::endl;
#endif
        nbad++;
      } else {
        ngood++;
      }
    }
#ifdef EDM_ML_DEBUG
    std::cout << "    overall ngood " << ngood << " nbad " << nbad << "\n";
#endif
    if (nbad == 0 || (ngood > 0 && nbad < ngood)) {
      active = true;
    }
  }
  return active;
}

References CaloTowerConstituentsMap::constituentsOf(), gather_cfg::cout, TauDecayModes::dec, HcalDetId::depth(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalChannelStatus::HcalCellDead, HcalChannelStatus::HcalCellMask, HcalChannelStatus::HcalCellOff, HcalEndcap, hcalQuality_, hcalTopology_, HcalTopology::idFront(), HcalDetId::ieta(), HcalDetId::iphi(), DetId::rawId(), mps_update::status, HcalDetId::subdet(), towerMap_, and HLT_2018_cff::towers.

Referenced by GEDPhotonProducer::fillPhotonCollection(), and ElectronHcalHelper::hasActiveHcal().

towerOf()

CaloTowerDetId EgammaHadTower::towerOf

(

const reco::CaloCluster &

cluster

)

const

Definition at line 35 of file EgammaHadTower.cc.

                                                                           {
  DetId detid = cluster.seed();
  if (detid.det() != DetId::Ecal) {
    // Basic clusters of hybrid super-cluster do not have the seed set; take the first DetId instead
    // Should be checked . The single Tower Mode should be favoured until fixed
    detid = cluster.hitsAndFractions()[0].first;
    if (detid.det() != DetId::Ecal) {
      CaloTowerDetId tower;
      return tower;
    }
  }
  CaloTowerDetId id(towerMap_->towerOf(detid));
  return id;
}

References DetId::det(), DetId::Ecal, reco::CaloCluster::hitsAndFractions(), triggerObjects_cff::id, reco::CaloCluster::seed(), hgcalTowerProducer_cfi::tower, towerMap_, and CaloTowerConstituentsMap::towerOf().

Referenced by towersOf().

towersOf()

std::vector< CaloTowerDetId > EgammaHadTower::towersOf

(

const reco::SuperCluster &

sc

)

const

Definition at line 50 of file EgammaHadTower.cc.

                                                                                   {
  std::vector<CaloTowerDetId> towers;
  std::vector<reco::CaloClusterPtr> orderedClusters;
 
  // in this mode, check only the tower behind the seed
  if (mode_ == SingleTower) {
    towers.push_back(towerOf(*sc.seed()));
  }
 
  // in this mode check the towers behind each basic cluster
  if (mode_ == TowersBehindCluster) {
    // Loop on the basic clusters
    for (auto it = sc.clustersBegin(); it != sc.clustersEnd(); ++it) {
      orderedClusters.push_back(*it);
    }
    std::sort(orderedClusters.begin(), orderedClusters.end(), [](auto& c1, auto& c2) { return (*c1 > *c2); });
    unsigned nclusters = orderedClusters.size();
    for (unsigned iclus = 0; iclus < nclusters && iclus < nMaxClusters_; ++iclus) {
      // Get the tower
      CaloTowerDetId id = towerOf(*(orderedClusters[iclus]));
#ifdef EDM_ML_DEBUG
      std::cout << "CaloTowerId " << id << std::endl;
#endif
      if (std::find(towers.begin(), towers.end(), id) == towers.end()) {
        towers.push_back(id);
      }
    }
  }
  //  if(towers.size() > 4) {
  //    std::cout << " NTOWERS " << towers.size() << " ";
  //    for(unsigned i=0; i<towers.size() ; ++i) {
  //      std::cout << towers[i] << " ";
  //    }
  //    std::cout <<  std::endl;
  //    for ( unsigned iclus=0 ; iclus < orderedClusters.size(); ++iclus) {
  //      // Get the tower
  //      CaloTowerDetId id = towerOf(*(orderedClusters[iclus]));
  //      std::cout << " Pos " << orderedClusters[iclus]->position() << " " << orderedClusters[iclus]->energy() << " " << id ;
  //    }
  //    std::cout << std::endl;
  //  }
  return towers;
}

References alignmentValidation::c1, reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), gather_cfg::cout, spr::find(), triggerObjects_cff::id, mode_, nMaxClusters_, reco::SuperCluster::seed(), SingleTower, towerOf(), HLT_2018_cff::towers, and TowersBehindCluster.

Referenced by PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hasActiveHcal(), ElectronHcalHelper::hcalTowersBehindClusters(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

Member Data Documentation

hcalQuality_

const HcalChannelQuality* EgammaHadTower::hcalQuality_

private

Definition at line 32 of file EgammaHadTower.h.

Referenced by EgammaHadTower(), and hasActiveHcal().

hcalTopology_

const HcalTopology* EgammaHadTower::hcalTopology_

private

Definition at line 33 of file EgammaHadTower.h.

Referenced by EgammaHadTower(), and hasActiveHcal().

mode_

HoeMode EgammaHadTower::mode_

private

Definition at line 30 of file EgammaHadTower.h.

Referenced by towersOf().

nMaxClusters_

unsigned int EgammaHadTower::nMaxClusters_

private

Definition at line 31 of file EgammaHadTower.h.

Referenced by EgammaHadTower(), and towersOf().

towerMap_

const CaloTowerConstituentsMap* EgammaHadTower::towerMap_

private

Definition at line 29 of file EgammaHadTower.h.

Referenced by EgammaHadTower(), hasActiveHcal(), and towerOf().