#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)

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_

Detailed Description

Definition at line 17 of file EGHcalRecHitSelector.h.

Constructor & Destructor Documentation

◆ EGHcalRecHitSelector()

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

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")) {}

◆ ~EGHcalRecHitSelector()

EGHcalRecHitSelector::~EGHcalRecHitSelector ( )

inline

Definition at line 20 of file EGHcalRecHitSelector.h.

20 {}

Member Function Documentation

◆ addDetIds()

template<typename CollType >

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

Definition at line 46 of file EGHcalRecHitSelector.h.

                                                                     {
   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());
   }
 }

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

Referenced by ReducedEGProducer::linkHcalHits(), and EgammaIsoHcalDetIdCollectionProducer::produce().

◆ calDIEta()

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

staticprivate

Definition at line 22 of file EGHcalRecHitSelector.cc.

                                                        {
   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;
 }

References HLT_2018_cff::dEta.

Referenced by addDetIds().

◆ calDIPhi()

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

staticprivate

Definition at line 33 of file EGHcalRecHitSelector.cc.

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

References HLT_2018_cff::dPhi.

Referenced by addDetIds().

◆ getMinEnergyHCAL_()

float EGHcalRecHitSelector::getMinEnergyHCAL_ ( HcalDetId id ) const

private

Definition at line 42 of file EGHcalRecHitSelector.cc.

                                                                 {
   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();
 }

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

Referenced by addDetIds().

◆ makePSetDescription()

edm::ParameterSetDescription EGHcalRecHitSelector::makePSetDescription ( )

static

Definition at line 12 of file EGHcalRecHitSelector.cc.

                                                                      {
   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;
 }