#include <EGHcalRecHitSelector.h>

Public Member Functions
template<typename CollType >
void	addDetIds (const reco::SuperCluster &superClus, const HBHERecHitCollection &recHits, CollType &detIdsToStore) const

	EGHcalRecHitSelector (const edm::ParameterSet &config, edm::ConsumesCollector)

void	setup (const edm::EventSetup &iSetup)

	~EGHcalRecHitSelector ()

Static Public Member Functions
static edm::ParameterSetDescription	makePSetDescription ()

Private Member Functions
float	getMinEnergyHCAL_ (HcalDetId id) const

Static Private Member Functions
static int	calDIEta (int iEta1, int iEta2)

static int	calDIPhi (int iPhi1, int iPhi2)

Private Attributes
int	maxDIEta_

int	maxDIPhi_

float	minEnergyHB_

float	minEnergyHEDefault_

float	minEnergyHEDepth1_

edm::ESHandle < CaloTowerConstituentsMap >	towerMap_

edm::ESGetToken < CaloTowerConstituentsMap, CaloGeometryRecord >	towerMapToken_

Detailed Description

Definition at line 19 of file EGHcalRecHitSelector.h.

Constructor & Destructor Documentation

EGHcalRecHitSelector::EGHcalRecHitSelector	(	const edm::ParameterSet &	config,
		edm::ConsumesCollector	cc
	)

explicit

Definition at line 5 of file EGHcalRecHitSelector.cc.

     : maxDIEta_(config.getParameter<int>("maxDIEta")),
       maxDIPhi_(config.getParameter<int>("maxDIPhi")),
       minEnergyHB_(config.getParameter<double>("minEnergyHB")),
       minEnergyHEDepth1_(config.getParameter<double>("minEnergyHEDepth1")),
       minEnergyHEDefault_(config.getParameter<double>("minEnergyHEDefault")),
       towerMapToken_(cc.esConsumes<CaloTowerConstituentsMap, CaloGeometryRecord, edm::Transition::BeginRun>()) {}

EGHcalRecHitSelector::~EGHcalRecHitSelector ( )

inline

Definition at line 22 of file EGHcalRecHitSelector.h.

22 {}

Member Function Documentation

template<typename CollType >

void EGHcalRecHitSelector::addDetIds	(	const reco::SuperCluster &	superClus,
		const HBHERecHitCollection &	recHits,
		CollType &	detIdsToStore
	)		const

Definition at line 49 of file EGHcalRecHitSelector.h.

References funct::abs(), calDIEta(), calDIPhi(), DetId::det(), DetId::Ecal, DetId::Forward, getMinEnergyHCAL_(), maxDIEta_, maxDIPhi_, reco::SuperCluster::seed(), ecaldqm::towerId(), and towerMap_.

Referenced by EgammaIsoHcalDetIdCollectionProducer::produce().

                                                                     {
   DetId seedId = superClus.seed()->seed();
   if (seedId.det() != DetId::Ecal && seedId.det() != DetId::Forward) {
     edm::LogError("EgammaIsoHcalDetIdCollectionProducerError")
         << "Somehow the supercluster has a seed which is not ECAL, something is badly wrong";
     return;
   }
   //so we are using CaloTowers to get the iEta/iPhi of the HCAL rec hit behind the seed cluster, this might go funny on tower 28 but shouldnt matter there
 
   if (seedId.det() == DetId::Forward)
     return;
 
   CaloTowerDetId towerId(towerMap_->towerOf(seedId));
   int seedHcalIEta = towerId.ieta();
   int seedHcalIPhi = towerId.iphi();
 
   for (auto& recHit : recHits) {
     int dIEtaAbs = std::abs(calDIEta(seedHcalIEta, recHit.id().ieta()));
     int dIPhiAbs = std::abs(calDIPhi(seedHcalIPhi, recHit.id().iphi()));
 
     if (dIEtaAbs <= maxDIEta_ && dIPhiAbs <= maxDIPhi_ && recHit.energy() > getMinEnergyHCAL_(recHit.id()))
       detIdsToStore.insert(detIdsToStore.end(), recHit.id().rawId());
   }
 }

int EGHcalRecHitSelector::calDIEta	(	int	iEta1,
		int	iEta2
	)

staticprivate

Definition at line 23 of file EGHcalRecHitSelector.cc.

Referenced by addDetIds().

                                                        {
   int dEta = iEta1 - iEta2;
   if (iEta1 * iEta2 < 0) {  //-ve to +ve transistion and no crystal at zero
     if (dEta < 0)
       dEta++;
     else
       dEta--;
   }
   return dEta;
 }

int EGHcalRecHitSelector::calDIPhi	(	int	iPhi1,
		int	iPhi2
	)

staticprivate

Definition at line 34 of file EGHcalRecHitSelector.cc.

Referenced by addDetIds().

                                                        {
   int dPhi = iPhi1 - iPhi2;
   if (dPhi > 72 / 2)
     dPhi -= 72;
   else if (dPhi < -72 / 2)
     dPhi += 72;
   return dPhi;
 }

float EGHcalRecHitSelector::getMinEnergyHCAL_ ( HcalDetId id ) const

private

Definition at line 43 of file EGHcalRecHitSelector.cc.

References HLT_FULL_cff::depth, HcalBarrel, HcalEndcap, SiStripPI::max, minEnergyHB_, minEnergyHEDefault_, and minEnergyHEDepth1_.

Referenced by addDetIds().

                                                                 {
   if (id.subdetId() == HcalBarrel)
     return minEnergyHB_;
   else if (id.subdetId() == HcalEndcap) {
     if (id.depth() == 1)
       return minEnergyHEDepth1_;
     else
       return minEnergyHEDefault_;
   } else
     return std::numeric_limits<float>::max();
 }

edm::ParameterSetDescription EGHcalRecHitSelector::makePSetDescription ( )

static

Definition at line 13 of file EGHcalRecHitSelector.cc.

References edm::ParameterSetDescription::add(), and submitPVResolutionJobs::desc.

Referenced by EgammaIsoHcalDetIdCollectionProducer::fillDescriptions().

                                                                      {
   edm::ParameterSetDescription desc;
   desc.add<int>("maxDIEta", 5);
   desc.add<int>("maxDIPhi", 5);
   desc.add<double>("minEnergyHB", 0.8);
   desc.add<double>("minEnergyHEDepth1", 0.1);
   desc.add<double>("minEnergyHEDefault", 0.2);
   return desc;
 }