EgammaTowerIsolation Class Reference

#include <EgammaTowerIsolation.h>

Public Types
enum	HcalDepth { AllDepths =-1, Undefined =0, Depth1 =1, Depth2 =2 }

Public Member Functions
	EgammaTowerIsolation (float extRadiusI, float intRadiusI, float etLow, signed int depth, const CaloTowerCollection *towers)
double	getTowerESum (const reco::Candidate *cand, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
double	getTowerESum (reco::SuperCluster const *sc, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
double	getTowerEtSum (const reco::Candidate *cand, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
double	getTowerEtSum (reco::SuperCluster const *sc, const std::vector< CaloTowerDetId > *detIdToExclude=0) const

Private Member Functions
double	getSum (bool et, reco::SuperCluster const &sc, const std::vector< CaloTowerDetId > *detIdToExclude) const

Private Attributes
signed int	depth_
float	extRadius
float	intRadius

Static Private Attributes
static thread_local uint32_t	id15 =0
static thread_local EgammaTowerIsolationNew< 1 > *	newAlgo =nullptr
static thread_local const CaloTowerCollection *	oldTowers =nullptr

Detailed Description

Definition at line 208 of file EgammaTowerIsolation.h.

Member Enumeration Documentation

enum EgammaTowerIsolation::HcalDepth

Enumerator
AllDepths
Undefined
Depth1
Depth2

Definition at line 211 of file EgammaTowerIsolation.h.

{AllDepths=-1,Undefined=0,Depth1=1,Depth2=2};

Constructor & Destructor Documentation

EgammaTowerIsolation::EgammaTowerIsolation	(	float	extRadiusI,
		float	intRadiusI,
		float	etLow,
		signed int	depth,
		const CaloTowerCollection *	towers
	)

Definition at line 31 of file EgammaTowerIsolation.cc.

References extRadius, id15, intRadius, newAlgo, EgammaTowerIsolationNew< NC >::nt, oldTowers, and edm::SortedCollection< T, SORT >::size().

                                                            :  
  depth_(depth), 
  extRadius(extRadiusI),
  intRadius(intRadiusI)
{
  assert(0==etLow);

  // extremely poor in quality  (test of performance)
  if (newAlgo==nullptr ||  towers!=oldTowers || towers->size()!=newAlgo->nt || (towers->size()>15 && (*towers)[15].id()!=id15)) {
    delete newAlgo;
    newAlgo = new EgammaTowerIsolationNew<1>(&extRadius,&intRadius,*towers);
    oldTowers=towers;
    id15 = towers->size()>15 ? (*towers)[15].id() : 0;
  }
}

Member Function Documentation

double EgammaTowerIsolation::getSum ( bool  et, reco::SuperCluster const &  sc, const std::vector< CaloTowerDetId > *  detIdToExclude  ) const

private

Definition at line 52 of file EgammaTowerIsolation.cc.

References AllDepths, EgammaTowerIsolationNew< NC >::compute(), Depth1, Depth2, depth_, extRadius, EgammaTowerIsolationNew< NC >::h2, EgammaTowerIsolationNew< NC >::he, intRadius, newAlgo, and EgammaTowerIsolationNew< NC >::setRadius().

Referenced by getTowerESum(), and getTowerEtSum().

                                                                                                                                 {

  if (0!=detIdToExclude) assert(0==intRadius);

  // hack
  newAlgo->setRadius(&extRadius,&intRadius);

  EgammaTowerIsolationNew<1>::Sum sum;
  newAlgo->compute(et, sum, sc, 
          (detIdToExclude==0) ? nullptr : &((*detIdToExclude).front()),
          (detIdToExclude==0) ? nullptr : (&(*detIdToExclude).back())+1
          );
  
  switch(depth_){
  case AllDepths: return detIdToExclude==0 ? sum.he[0] : sum.heBC[0]; 
  case Depth1: return detIdToExclude==0 ? sum.he[0]-sum.h2[0] : sum.heBC[0]-sum.h2BC[0]; 
  case Depth2:return detIdToExclude==0 ? sum.h2[0] : sum.h2BC[0]; 
  default: return 0;
  }
  return 0;
}

double EgammaTowerIsolation::getTowerESum ( const reco::Candidate *  cand, const std::vector< CaloTowerDetId > *  detIdToExclude = 0  ) const

inline

Definition at line 224 of file EgammaTowerIsolation.h.

References reco::Candidate::get(), and getSum().

Referenced by ConvertedPhotonProducer::buildCollections(), ConversionTrackCandidateProducer::buildCollections(), PhotonProducer::fillPhotonCollection(), GEDPhotonProducer::fillPhotonCollection(), ElectronHcalHelper::hcalESumDepth1(), and ElectronHcalHelper::hcalESumDepth2().

                                                                                                            {
    reco::SuperCluster const & sc =  *cand->get<reco::SuperClusterRef>().get();
    return getSum(false,sc,detIdToExclude);
  }

double EgammaTowerIsolation::getTowerESum ( reco::SuperCluster const *  sc, const std::vector< CaloTowerDetId > *  detIdToExclude = 0  ) const

inline

Definition at line 231 of file EgammaTowerIsolation.h.

References getSum().

                                                                                                              {
    return getSum(false,*sc,detIdToExclude);
  }

double EgammaTowerIsolation::getTowerEtSum ( const reco::Candidate *  cand, const std::vector< CaloTowerDetId > *  detIdToExclude = 0  ) const

inline

Definition at line 220 of file EgammaTowerIsolation.h.

References reco::Candidate::get(), and getSum().

Referenced by PhotonIsolationCalculator::calculateHcalTowerIso(), GsfElectronAlgo::createElectron(), and EgammaHLTBcHcalIsolationProducersRegional::produce().

                                                                                                             {
    reco::SuperCluster const & sc =  *cand->get<reco::SuperClusterRef>().get();
    return getSum(true,sc,detIdToExclude);
  }

double EgammaTowerIsolation::getTowerEtSum ( reco::SuperCluster const *  sc, const std::vector< CaloTowerDetId > *  detIdToExclude = 0  ) const

inline

Definition at line 228 of file EgammaTowerIsolation.h.

References getSum().

                                                                                                               {
    return getSum(true,*sc,detIdToExclude);
  }

Member Data Documentation

signed int EgammaTowerIsolation::depth_

private

Definition at line 242 of file EgammaTowerIsolation.h.

Referenced by getSum().

float EgammaTowerIsolation::extRadius

private

Definition at line 243 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation(), and getSum().

thread_local uint32_t EgammaTowerIsolation::id15 =0

staticprivate

Definition at line 241 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation().

float EgammaTowerIsolation::intRadius

private

Definition at line 244 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation(), and getSum().

thread_local EgammaTowerIsolationNew< 1 > * EgammaTowerIsolation::newAlgo =nullptr

staticprivate

Definition at line 239 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation(), and getSum().

thread_local const CaloTowerCollection * EgammaTowerIsolation::oldTowers =nullptr

staticprivate

Definition at line 240 of file EgammaTowerIsolation.h.

Referenced by EgammaTowerIsolation().