CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Private Member Functions | Private Attributes
EgammaHcalIsolation Class Reference

#include <EgammaHcalIsolation.h>

Public Types

using arrayHB = std::array< double, 4 >
 
using arrayHE = std::array< double, 7 >
 
enum  InclusionRule : int { InclusionRule::withinConeAroundCluster = 0, InclusionRule::isBehindClusterSeed = 1 }
 

Public Member Functions

 EgammaHcalIsolation (InclusionRule extIncRule, double extRadius, InclusionRule intIncRule, double intRadius, const arrayHB &eThresHB, const arrayHB &etThresHB, int maxSeverityHB, const arrayHE &eThresHE, const arrayHE &etThresHE, int maxSeverityHE, const HBHERecHitCollection &mhbhe, edm::ESHandle< CaloGeometry > caloGeometry, edm::ESHandle< HcalTopology > hcalTopology, edm::ESHandle< HcalChannelQuality > hcalChStatus, edm::ESHandle< HcalSeverityLevelComputer > hcalSevLvlComputer, edm::ESHandle< CaloTowerConstituentsMap > towerMap)
 
 EgammaHcalIsolation (InclusionRule extIncRule, double extRadius, InclusionRule intIncRule, double intRadius, const arrayHB &eThresHB, const arrayHB &etThresHB, int maxSeverityHB, const arrayHE &eThresHE, const arrayHE &etThresHE, int maxSeverityHE, const HBHERecHitCollection &mhbhe, const CaloGeometry &caloGeometry, const HcalTopology &hcalTopology, const HcalChannelQuality &hcalChStatus, const HcalSeverityLevelComputer &hcalSevLvlComputer, const CaloTowerConstituentsMap &towerMap)
 
double getHcalESum (const reco::Candidate *c, int depth) const
 
double getHcalESum (const reco::SuperCluster *sc, int depth) const
 
double getHcalESum (const math::XYZPoint &p, int depth) const
 
double getHcalESum (const GlobalPoint &pclu, int depth) const
 
double getHcalESumBc (const reco::Candidate *c, int depth) const
 
double getHcalESumBc (const reco::SuperCluster *sc, int depth) const
 
double getHcalESumBc (const math::XYZPoint &p, int depth, int ieta, int iphi, int include_or_exclude) const
 
double getHcalESumBc (const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude) const
 
double getHcalEtSum (const reco::Candidate *c, int depth) const
 
double getHcalEtSum (const reco::SuperCluster *sc, int depth) const
 
double getHcalEtSum (const math::XYZPoint &p, int depth) const
 
double getHcalEtSum (const GlobalPoint &pclu, int depth) const
 
double getHcalEtSumBc (const reco::Candidate *c, int depth) const
 
double getHcalEtSumBc (const reco::SuperCluster *sc, int depth) const
 
double getHcalEtSumBc (const math::XYZPoint &p, int depth, int ieta, int iphi, int include_or_exclude) const
 
double getHcalEtSumBc (const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude) const
 

Private Member Functions

double getHcalSum (const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
 
double goodHitEnergy (float pcluEta, float pcluPhi, const HBHERecHit &hit, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
 

Private Attributes

const CaloGeometry & caloGeometry_
 
arrayHB eThresHB_
 
arrayHE eThresHE_
 
arrayHB etThresHB_
 
arrayHE etThresHE_
 
InclusionRule extIncRule_
 
double extRadius_
 
const HcalChannelQuality & hcalChStatus_
 
const HcalSeverityLevelComputer & hcalSevLvlComputer_
 
const HcalTopology & hcalTopology_
 
InclusionRule intIncRule_
 
double intRadius_
 
int maxSeverityHB_
 
int maxSeverityHE_
 
const HBHERecHitCollection & mhbhe_
 
const CaloTowerConstituentsMap & towerMap_
 

Detailed Description

Definition at line 38 of file EgammaHcalIsolation.h.

Member Typedef Documentation

◆ arrayHB

using EgammaHcalIsolation::arrayHB = std::array<double, 4>

Definition at line 41 of file EgammaHcalIsolation.h.

◆ arrayHE

using EgammaHcalIsolation::arrayHE = std::array<double, 7>

Definition at line 42 of file EgammaHcalIsolation.h.

Member Enumeration Documentation

◆ InclusionRule

enum EgammaHcalIsolation::InclusionRule : int
strong
Enumerator
withinConeAroundCluster 
isBehindClusterSeed 

Definition at line 40 of file EgammaHcalIsolation.h.

40 : int { withinConeAroundCluster = 0, isBehindClusterSeed = 1 };

Constructor & Destructor Documentation

◆ EgammaHcalIsolation() [1/2]

EgammaHcalIsolation::EgammaHcalIsolation ( InclusionRule  extIncRule,
double  extRadius,
InclusionRule  intIncRule,
double  intRadius,
const arrayHB &  eThresHB,
const arrayHB &  etThresHB,
int  maxSeverityHB,
const arrayHE &  eThresHE,
const arrayHE &  etThresHE,
int  maxSeverityHE,
const HBHERecHitCollection &  mhbhe,
edm::ESHandle< CaloGeometry >  caloGeometry,
edm::ESHandle< HcalTopology >  hcalTopology,
edm::ESHandle< HcalChannelQuality >  hcalChStatus,
edm::ESHandle< HcalSeverityLevelComputer >  hcalSevLvlComputer,
edm::ESHandle< CaloTowerConstituentsMap >  towerMap 
)

◆ EgammaHcalIsolation() [2/2]

EgammaHcalIsolation::EgammaHcalIsolation ( InclusionRule  extIncRule,
double  extRadius,
InclusionRule  intIncRule,
double  intRadius,
const arrayHB &  eThresHB,
const arrayHB &  etThresHB,
int  maxSeverityHB,
const arrayHE &  eThresHE,
const arrayHE &  etThresHE,
int  maxSeverityHE,
const HBHERecHitCollection &  mhbhe,
const CaloGeometry &  caloGeometry,
const HcalTopology &  hcalTopology,
const HcalChannelQuality &  hcalChStatus,
const HcalSeverityLevelComputer &  hcalSevLvlComputer,
const CaloTowerConstituentsMap &  towerMap 
)

Member Function Documentation

◆ getHcalESum() [1/4]

double EgammaHcalIsolation::getHcalESum ( const reco::Candidate *  c,
int  depth 
) const
inline

Definition at line 79 of file EgammaHcalIsolation.h.

References HltBtagPostValidation_cff::c, LEDCalibrationChannels::depth, and edm::Ref< C, T, F >::get().

Referenced by getHcalESum(), and EgammaHLTHcalVarProducerFromRecHit::produce().

79  {
80  return getHcalESum(c->get<reco::SuperClusterRef>().get(), depth);
81  }
EgammaHcalIsolation::getHcalESum
double getHcalESum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:79
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:232

◆ getHcalESum() [2/4]

double EgammaHcalIsolation::getHcalESum ( const reco::SuperCluster *  sc,
int  depth 
) const
inline

Definition at line 85 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalESum(), and reco::CaloCluster::position().

Referenced by getHcalESum().

85 { return getHcalESum(sc->position(), depth); }
reco::CaloCluster::position
const math::XYZPoint & position() const
cluster centroid position
Definition: CaloCluster.h:154
EgammaHcalIsolation::getHcalESum
double getHcalESum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:79

◆ getHcalESum() [3/4]

double EgammaHcalIsolation::getHcalESum ( const math::XYZPoint &  p,
int  depth 
) const
inline

Definition at line 87 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalESum(), and AlCaHLTBitMon_ParallelJobs::p.

87  {
88  return getHcalESum(GlobalPoint(p.x(), p.y(), p.z()), depth);
89  }
EgammaHcalIsolation::getHcalESum
double getHcalESum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:79
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ getHcalESum() [4/4]

double EgammaHcalIsolation::getHcalESum ( const GlobalPoint &  pclu,
int  depth 
) const
inline

Definition at line 93 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalSum(), and scaleToE().

93 { return getHcalSum(pclu, depth, 0, 0, 0, &scaleToE); }
EgammaHcalIsolation::getHcalSum
double getHcalSum(const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
Definition: EgammaHcalIsolation.cc:180
scaleToE
double scaleToE(const double &eta)
Definition: EgammaHcalIsolation.cc:21

◆ getHcalESumBc() [1/4]

double EgammaHcalIsolation::getHcalESumBc ( const reco::Candidate *  c,
int  depth 
) const
inline

Definition at line 96 of file EgammaHcalIsolation.h.

References HltBtagPostValidation_cff::c, LEDCalibrationChannels::depth, and edm::Ref< C, T, F >::get().

Referenced by getHcalESumBc(), and EgammaHLTHcalVarProducerFromRecHit::produce().

96  {
97  return getHcalESumBc(c->get<reco::SuperClusterRef>().get(), depth);
98  }
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
EgammaHcalIsolation::getHcalESumBc
double getHcalESumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:96
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:232

◆ getHcalESumBc() [2/4]

double EgammaHcalIsolation::getHcalESumBc ( const reco::SuperCluster *  sc,
int  depth 
) const
inline

Definition at line 102 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, extIncRule_, getHcalESumBc(), intIncRule_, isBehindClusterSeed, reco::CaloCluster::position(), reco::SuperCluster::seed(), l1tHGCalTowerProducer_cfi::tower, towerMap_, egamma::towerOf(), and withinConeAroundCluster.

102  {
103  const auto tower = egamma::towerOf(*(sc->seed()), towerMap_);
104 
105  if (extIncRule_ == InclusionRule::isBehindClusterSeed and intIncRule_ == InclusionRule::withinConeAroundCluster)
106  return getHcalESumBc(sc->position(), depth, tower.ieta(), tower.iphi(), -1);
107  else if (extIncRule_ == InclusionRule::withinConeAroundCluster and
108  intIncRule_ == InclusionRule::isBehindClusterSeed)
109  return getHcalESumBc(sc->position(), depth, tower.ieta(), tower.iphi(), 1);
110 
111  return getHcalESumBc(sc->position(), depth, tower.ieta(), tower.iphi(), 0);
112  }
reco::CaloCluster::position
const math::XYZPoint & position() const
cluster centroid position
Definition: CaloCluster.h:154
EgammaHcalIsolation::getHcalESumBc
double getHcalESumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:96
egamma::towerOf
CaloTowerDetId towerOf(reco::CaloCluster const &cluster, CaloTowerConstituentsMap const &towerMap)
Definition: EgammaHadTower.cc:10
EgammaHcalIsolation::towerMap_
const CaloTowerConstituentsMap & towerMap_
Definition: EgammaHcalIsolation.h:172
reco::SuperCluster::seed
const CaloClusterPtr & seed() const
seed BasicCluster
Definition: SuperCluster.h:79
l1tHGCalTowerProducer_cfi.tower
tower
Definition: l1tHGCalTowerProducer_cfi.py:4

◆ getHcalESumBc() [3/4]

double EgammaHcalIsolation::getHcalESumBc ( const math::XYZPoint &  p,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude 
) const
inline

Definition at line 124 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalESumBc(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and AlCaHLTBitMon_ParallelJobs::p.

124  {
125  return getHcalESumBc(GlobalPoint(p.x(), p.y(), p.z()), depth, ieta, iphi, include_or_exclude);
126  }
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
EgammaHcalIsolation::getHcalESumBc
double getHcalESumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:96
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ getHcalESumBc() [4/4]

double EgammaHcalIsolation::getHcalESumBc ( const GlobalPoint &  pclu,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude 
) const
inline

Definition at line 130 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalSum(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and scaleToE().

130  {
131  return getHcalSum(pclu, depth, ieta, iphi, include_or_exclude, &scaleToE);
132  }
EgammaHcalIsolation::getHcalSum
double getHcalSum(const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
Definition: EgammaHcalIsolation.cc:180
scaleToE
double scaleToE(const double &eta)
Definition: EgammaHcalIsolation.cc:21

◆ getHcalEtSum() [1/4]

double EgammaHcalIsolation::getHcalEtSum ( const reco::Candidate *  c,
int  depth 
) const
inline

Definition at line 82 of file EgammaHcalIsolation.h.

References HltBtagPostValidation_cff::c, LEDCalibrationChannels::depth, and edm::Ref< C, T, F >::get().

Referenced by GsfElectronAlgo::createElectron(), getHcalEtSum(), and EgammaHLTHcalVarProducerFromRecHit::produce().

82  {
83  return getHcalEtSum(c->get<reco::SuperClusterRef>().get(), depth);
84  }
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
EgammaHcalIsolation::getHcalEtSum
double getHcalEtSum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:82
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:232

◆ getHcalEtSum() [2/4]

double EgammaHcalIsolation::getHcalEtSum ( const reco::SuperCluster *  sc,
int  depth 
) const
inline

Definition at line 86 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalEtSum(), and reco::CaloCluster::position().

Referenced by getHcalEtSum().

86 { return getHcalEtSum(sc->position(), depth); }
reco::CaloCluster::position
const math::XYZPoint & position() const
cluster centroid position
Definition: CaloCluster.h:154
EgammaHcalIsolation::getHcalEtSum
double getHcalEtSum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:82

◆ getHcalEtSum() [3/4]

double EgammaHcalIsolation::getHcalEtSum ( const math::XYZPoint &  p,
int  depth 
) const
inline

Definition at line 90 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalEtSum(), and AlCaHLTBitMon_ParallelJobs::p.

90  {
91  return getHcalEtSum(GlobalPoint(p.x(), p.y(), p.z()), depth);
92  }
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
EgammaHcalIsolation::getHcalEtSum
double getHcalEtSum(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:82
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ getHcalEtSum() [4/4]

double EgammaHcalIsolation::getHcalEtSum ( const GlobalPoint &  pclu,
int  depth 
) const
inline

Definition at line 94 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalSum(), and scaleToEt().

94 { return getHcalSum(pclu, depth, 0, 0, 0, &scaleToEt); }
EgammaHcalIsolation::getHcalSum
double getHcalSum(const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
Definition: EgammaHcalIsolation.cc:180
scaleToEt
double scaleToEt(const double &eta)
Definition: EgammaHcalIsolation.cc:22

◆ getHcalEtSumBc() [1/4]

double EgammaHcalIsolation::getHcalEtSumBc ( const reco::Candidate *  c,
int  depth 
) const
inline

Definition at line 99 of file EgammaHcalIsolation.h.

References HltBtagPostValidation_cff::c, LEDCalibrationChannels::depth, and edm::Ref< C, T, F >::get().

Referenced by GsfElectronAlgo::createElectron(), getHcalEtSumBc(), and EgammaHLTHcalVarProducerFromRecHit::produce().

99  {
100  return getHcalEtSumBc(c->get<reco::SuperClusterRef>().get(), depth);
101  }
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
EgammaHcalIsolation::getHcalEtSumBc
double getHcalEtSumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:99
edm::Ref::get
T const * get() const
Returns C++ pointer to the item.
Definition: Ref.h:232

◆ getHcalEtSumBc() [2/4]

double EgammaHcalIsolation::getHcalEtSumBc ( const reco::SuperCluster *  sc,
int  depth 
) const
inline

Definition at line 113 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, extIncRule_, getHcalEtSumBc(), intIncRule_, isBehindClusterSeed, reco::CaloCluster::position(), reco::SuperCluster::seed(), l1tHGCalTowerProducer_cfi::tower, towerMap_, egamma::towerOf(), and withinConeAroundCluster.

113  {
114  const auto tower = egamma::towerOf(*(sc->seed()), towerMap_);
115 
116  if (extIncRule_ == InclusionRule::isBehindClusterSeed and intIncRule_ == InclusionRule::withinConeAroundCluster)
117  return getHcalEtSumBc(sc->position(), depth, tower.ieta(), tower.iphi(), -1);
118  else if (extIncRule_ == InclusionRule::withinConeAroundCluster and
119  intIncRule_ == InclusionRule::isBehindClusterSeed)
120  return getHcalEtSumBc(sc->position(), depth, tower.ieta(), tower.iphi(), 1);
121 
122  return getHcalEtSumBc(sc->position(), depth, tower.ieta(), tower.iphi(), 0);
123  }
reco::CaloCluster::position
const math::XYZPoint & position() const
cluster centroid position
Definition: CaloCluster.h:154
egamma::towerOf
CaloTowerDetId towerOf(reco::CaloCluster const &cluster, CaloTowerConstituentsMap const &towerMap)
Definition: EgammaHadTower.cc:10
EgammaHcalIsolation::getHcalEtSumBc
double getHcalEtSumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:99
EgammaHcalIsolation::towerMap_
const CaloTowerConstituentsMap & towerMap_
Definition: EgammaHcalIsolation.h:172
reco::SuperCluster::seed
const CaloClusterPtr & seed() const
seed BasicCluster
Definition: SuperCluster.h:79
l1tHGCalTowerProducer_cfi.tower
tower
Definition: l1tHGCalTowerProducer_cfi.py:4

◆ getHcalEtSumBc() [3/4]

double EgammaHcalIsolation::getHcalEtSumBc ( const math::XYZPoint &  p,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude 
) const
inline

Definition at line 127 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalEtSumBc(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and AlCaHLTBitMon_ParallelJobs::p.

127  {
128  return getHcalEtSumBc(GlobalPoint(p.x(), p.y(), p.z()), depth, ieta, iphi, include_or_exclude);
129  }
GlobalPoint
Global3DPoint GlobalPoint
Definition: GlobalPoint.h:10
EgammaHcalIsolation::getHcalEtSumBc
double getHcalEtSumBc(const reco::Candidate *c, int depth) const
Definition: EgammaHcalIsolation.h:99
AlCaHLTBitMon_ParallelJobs.p
def p
Definition: AlCaHLTBitMon_ParallelJobs.py:153

◆ getHcalEtSumBc() [4/4]

double EgammaHcalIsolation::getHcalEtSumBc ( const GlobalPoint &  pclu,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude 
) const
inline

Definition at line 133 of file EgammaHcalIsolation.h.

References LEDCalibrationChannels::depth, getHcalSum(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and scaleToEt().

133  {
134  return getHcalSum(pclu, depth, ieta, iphi, include_or_exclude, &scaleToEt);
135  }
EgammaHcalIsolation::getHcalSum
double getHcalSum(const GlobalPoint &pclu, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
Definition: EgammaHcalIsolation.cc:180
scaleToEt
double scaleToEt(const double &eta)
Definition: EgammaHcalIsolation.cc:22

◆ getHcalSum()

double EgammaHcalIsolation::getHcalSum ( const GlobalPoint &  pclu,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude,
double(*)(const double &)  scale 
) const
private

Definition at line 180 of file EgammaHcalIsolation.cc.

References LEDCalibrationChannels::depth, PV3DBase< T, PVType, FrameType >::eta(), goodHitEnergy(), LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, mhbhe_, PV3DBase< T, PVType, FrameType >::phi(), and L1EGammaClusterEmuProducer_cfi::scale.

Referenced by getHcalESum(), getHcalESumBc(), getHcalEtSum(), and getHcalEtSumBc().

185  {
186  double sum = 0.;
187  const float pcluEta = pclu.eta();
188  const float pcluPhi = pclu.phi();
189  for (const auto &hit : mhbhe_)
190  sum += goodHitEnergy(pcluEta, pcluPhi, hit, depth, ieta, iphi, include_or_exclude, scale);
191 
192  return sum;
193 }
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:66
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:73
EgammaHcalIsolation::goodHitEnergy
double goodHitEnergy(float pcluEta, float pcluPhi, const HBHERecHit &hit, int depth, int ieta, int iphi, int include_or_exclude, double(*scale)(const double &)) const
Definition: EgammaHcalIsolation.cc:122
EgammaHcalIsolation::mhbhe_
const HBHERecHitCollection & mhbhe_
Definition: EgammaHcalIsolation.h:167
L1EGammaClusterEmuProducer_cfi.scale
scale
Definition: L1EGammaClusterEmuProducer_cfi.py:15

◆ goodHitEnergy()

double EgammaHcalIsolation::goodHitEnergy ( float  pcluEta,
float  pcluPhi,
const HBHERecHit &  hit,
int  depth,
int  ieta,
int  iphi,
int  include_or_exclude,
double(*)(const double &)  scale 
) const
private

Definition at line 122 of file EgammaHcalIsolation.cc.

References caloGeometry_, HLTMuonOfflineAnalyzer_cfi::deltaR2, LEDCalibrationChannels::depth, HcalDetId::depth(), eThresHB_, eThresHE_, etThresHB_, etThresHE_, extIncRule_, extRadius_, RemoveAddSevLevel::flag, CaloGeometry::getGeometry(), HcalSeverityLevelComputer::getSeverityLevel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), HcalBarrel, hcalChStatus_, HcalEndcap, hcalSevLvlComputer_, hcalTopology_, hcalSimParameters_cfi::he, trackingPlots::hp, HcalTopology::idFront(), LEDCalibrationChannels::ieta, intIncRule_, intRadius_, LEDCalibrationChannels::iphi, maxSeverityHB_, maxSeverityHE_, or, HcalSeverityLevelComputer::recoveredRecHit(), L1EGammaClusterEmuProducer_cfi::scale, scaleToEt(), ErrorSummaryFilter_cfi::severity, and withinConeAroundCluster.

Referenced by getHcalSum().

129  {
130  const HcalDetId hid(hit.detid());
131  const int hd = hid.depth(), he = hid.ieta(), hp = hid.iphi();
132  const int h1 = hd - 1;
133 
134  if (include_or_exclude == -1 and (he != ieta or hp != iphi))
135  return 0.;
136 
137  if (include_or_exclude == 1 and (he == ieta and hp == iphi))
138  return 0.;
139 
140  if ((hid.subdet() == HcalBarrel and (hd < 1 or hd > int(eThresHB_.size()))) or
141  (hid.subdet() == HcalEndcap and (hd < 1 or hd > int(eThresHE_.size()))))
142  edm::LogWarning("EgammaHcalIsolation")
143  << " hit in subdet " << hid.subdet() << " has an unaccounted for depth of " << hd << "!!";
144 
145  const bool right_depth = (depth == 0 or hd == depth);
146  if (!right_depth)
147  return 0.;
148 
149  const bool goodHBe = hid.subdet() == HcalBarrel and hit.energy() > eThresHB_[h1];
150  const bool goodHEe = hid.subdet() == HcalEndcap and hit.energy() > eThresHE_[h1];
151  if (!(goodHBe or goodHEe))
152  return 0.;
153 
154  const auto phit = caloGeometry_.getGeometry(hit.detid())->repPos();
155  const float phitEta = phit.eta();
156 
157  if (extIncRule_ == InclusionRule::withinConeAroundCluster or intIncRule_ == InclusionRule::withinConeAroundCluster) {
158  auto const dR2 = deltaR2(pcluEta, pcluPhi, phitEta, phit.phi());
159  if ((extIncRule_ == InclusionRule::withinConeAroundCluster and dR2 > extRadius_) or
160  (intIncRule_ == InclusionRule::withinConeAroundCluster and dR2 < intRadius_))
161  return 0.;
162  }
163 
164  DetId did = hcalTopology_.idFront(hid);
165  const uint32_t flag = hit.flags();
166  const uint32_t dbflag = hcalChStatus_.getValues(did)->getValue();
167  int severity = hcalSevLvlComputer_.getSeverityLevel(did, flag, dbflag);
168  bool recovered = hcalSevLvlComputer_.recoveredRecHit(did, flag);
169 
170  const double het = hit.energy() * scaleToEt(phitEta);
171  const bool goodHB = goodHBe and (severity <= maxSeverityHB_ or recovered) and het > etThresHB_[h1];
172  const bool goodHE = goodHEe and (severity <= maxSeverityHE_ or recovered) and het > etThresHE_[h1];
173 
174  if (goodHB or goodHE)
175  return hit.energy() * scale(phitEta);
176 
177  return 0.;
178 }
scaleToEt
double scaleToEt(const double &eta)
Definition: EgammaHcalIsolation.cc:22
EgammaHcalIsolation::caloGeometry_
const CaloGeometry & caloGeometry_
Definition: EgammaHcalIsolation.h:168
HcalSeverityLevelComputer::recoveredRecHit
bool recoveredRecHit(const DetId &myid, const uint32_t &myflag) const
Definition: HcalSeverityLevelComputer.cc:363
EgammaHcalIsolation::hcalChStatus_
const HcalChannelQuality & hcalChStatus_
Definition: EgammaHcalIsolation.h:170
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
HLTMuonOfflineAnalyzer_cfi.deltaR2
deltaR2
Definition: HLTMuonOfflineAnalyzer_cfi.py:53
HcalTopology::idFront
HcalDetId idFront(const HcalDetId &id) const
Definition: HcalTopology.h:170
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
ErrorSummaryFilter_cfi.severity
severity
Definition: ErrorSummaryFilter_cfi.py:5
CaloGeometry::getGeometry
std::shared_ptr< const CaloCellGeometry > getGeometry(const DetId &id) const
Get the cell geometry of a given detector id.
Definition: CaloGeometry.cc:60
HcalChannelStatus::getValue
uint32_t getValue() const
Definition: HcalChannelStatus.h:60
hcalSimParameters_cfi.he
he
Definition: hcalSimParameters_cfi.py:79
DetId
Definition: DetId.h:17
L1EGammaClusterEmuProducer_cfi.scale
scale
Definition: L1EGammaClusterEmuProducer_cfi.py:15
HcalSeverityLevelComputer::getSeverityLevel
int getSeverityLevel(const DetId &myid, const uint32_t &myflag, const uint32_t &mystatus) const
Definition: HcalSeverityLevelComputer.cc:304
EgammaHcalIsolation::hcalTopology_
const HcalTopology & hcalTopology_
Definition: EgammaHcalIsolation.h:169
EgammaHcalIsolation::hcalSevLvlComputer_
const HcalSeverityLevelComputer & hcalSevLvlComputer_
Definition: EgammaHcalIsolation.h:171
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
HcalDetId::depth
constexpr int depth() const
get the tower depth
Definition: HcalDetId.h:164

Member Data Documentation

◆ caloGeometry_

const CaloGeometry& EgammaHcalIsolation::caloGeometry_
private

Definition at line 168 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ eThresHB_

arrayHB EgammaHcalIsolation::eThresHB_
private

Definition at line 159 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ eThresHE_

arrayHE EgammaHcalIsolation::eThresHE_
private

Definition at line 163 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ etThresHB_

arrayHB EgammaHcalIsolation::etThresHB_
private

Definition at line 160 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ etThresHE_

arrayHE EgammaHcalIsolation::etThresHE_
private

Definition at line 164 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ extIncRule_

InclusionRule EgammaHcalIsolation::extIncRule_
private

Definition at line 154 of file EgammaHcalIsolation.h.

Referenced by getHcalESumBc(), getHcalEtSumBc(), and goodHitEnergy().

◆ extRadius_

double EgammaHcalIsolation::extRadius_
private

Definition at line 155 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ hcalChStatus_

const HcalChannelQuality& EgammaHcalIsolation::hcalChStatus_
private

Definition at line 170 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ hcalSevLvlComputer_

const HcalSeverityLevelComputer& EgammaHcalIsolation::hcalSevLvlComputer_
private

Definition at line 171 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ hcalTopology_

const HcalTopology& EgammaHcalIsolation::hcalTopology_
private

Definition at line 169 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ intIncRule_

InclusionRule EgammaHcalIsolation::intIncRule_
private

Definition at line 156 of file EgammaHcalIsolation.h.

Referenced by getHcalESumBc(), getHcalEtSumBc(), and goodHitEnergy().

◆ intRadius_

double EgammaHcalIsolation::intRadius_
private

Definition at line 157 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ maxSeverityHB_

int EgammaHcalIsolation::maxSeverityHB_
private

Definition at line 161 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ maxSeverityHE_

int EgammaHcalIsolation::maxSeverityHE_
private

Definition at line 165 of file EgammaHcalIsolation.h.

Referenced by goodHitEnergy().

◆ mhbhe_

const HBHERecHitCollection& EgammaHcalIsolation::mhbhe_
private

Definition at line 167 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

◆ towerMap_

const CaloTowerConstituentsMap& EgammaHcalIsolation::towerMap_
private

Definition at line 172 of file EgammaHcalIsolation.h.

Referenced by getHcalESumBc(), and getHcalEtSumBc().

Generated for CMSSW Reference Manual by doxygen 1.8.14