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

inline

Definition at line 20 of file EGHcalRecHitSelector.h.

References addDetIds().

20 {}

Member Function Documentation

template<typename CollType >

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

Definition at line 47 of file EGHcalRecHitSelector.h.

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

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

 {
   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 26 of file EGHcalRecHitSelector.cc.

References particleFlow_cfi::dEta.

Referenced by addDetIds(), and setup().

 {  
   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 37 of file EGHcalRecHitSelector.cc.

References particleFlow_cfi::dPhi.

Referenced by addDetIds(), and setup().

 {
   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 45 of file EGHcalRecHitSelector.cc.

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

Referenced by addDetIds(), and setup().

 {
   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 15 of file EGHcalRecHitSelector.cc.

References edm::ParameterSetDescription::add().

Referenced by EgammaIsoHcalDetIdCollectionProducer::fillDescriptions(), and setup().

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