EgammaHcalIsolation Class Reference

#include <EgammaHcalIsolation.h>

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

Public Member Functions
	EgammaHcalIsolation (double extRadius, double intRadius, double eLowB, double eLowE, double etLowB, double etLowE, edm::ESHandle< CaloGeometry > theCaloGeom, const HBHERecHitCollection &mhbhe)
double	getHcalESum (const GlobalPoint &pclu) const
double	getHcalESum (const math::XYZPoint &p) const
double	getHcalESum (const reco::Candidate *c) const
double	getHcalESum (const reco::SuperCluster *sc) const
double	getHcalESumDepth1 (const GlobalPoint &pclu) const
double	getHcalESumDepth1 (const math::XYZPoint &p) const
double	getHcalESumDepth1 (const reco::Candidate *c) const
double	getHcalESumDepth1 (const reco::SuperCluster *sc) const
double	getHcalESumDepth2 (const GlobalPoint &pclu) const
double	getHcalESumDepth2 (const math::XYZPoint &p) const
double	getHcalESumDepth2 (const reco::Candidate *c) const
double	getHcalESumDepth2 (const reco::SuperCluster *sc) const
double	getHcalEtSum (const GlobalPoint &pclu) const
double	getHcalEtSum (const math::XYZPoint &p) const
double	getHcalEtSum (const reco::Candidate *c) const
double	getHcalEtSum (const reco::SuperCluster *sc) const
double	getHcalEtSumDepth1 (const GlobalPoint &pclu) const
double	getHcalEtSumDepth1 (const math::XYZPoint &p) const
double	getHcalEtSumDepth1 (const reco::Candidate *c) const
double	getHcalEtSumDepth1 (const reco::SuperCluster *sc) const
double	getHcalEtSumDepth2 (const GlobalPoint &pclu) const
double	getHcalEtSumDepth2 (const math::XYZPoint &p) const
double	getHcalEtSumDepth2 (const reco::Candidate *c) const
double	getHcalEtSumDepth2 (const reco::SuperCluster *sc) const
	~EgammaHcalIsolation ()

Private Member Functions
double	getHcalSum (const GlobalPoint &, const HcalDepth &, double(*)(const double &)) const
bool	isDepth2 (const DetId &) const

Private Attributes
CaloDualConeSelector< HBHERecHit > *	doubleConeSel_
double	eLowB_
double	eLowE_
double	etLowB_
double	etLowE_
double	extRadius_
double	intRadius_
const HBHERecHitCollection &	mhbhe_
edm::ESHandle< CaloGeometry >	theCaloGeom_

Detailed Description

Definition at line 30 of file EgammaHcalIsolation.h.

Member Enumeration Documentation

HcalDepth

enum EgammaHcalIsolation::HcalDepth

Enumerator
AllDepths
Depth1
Depth2

Definition at line 32 of file EgammaHcalIsolation.h.

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

Constructor & Destructor Documentation

EgammaHcalIsolation()

EgammaHcalIsolation::EgammaHcalIsolation	(	double	extRadius,
		double	intRadius,
		double	eLowB,
		double	eLowE,
		double	etLowB,
		double	etLowE,
		edm::ESHandle< CaloGeometry >	theCaloGeom,
		const HBHERecHitCollection &	mhbhe
	)

Definition at line 27 of file EgammaHcalIsolation.cc.

    : extRadius_(extRadius),
      intRadius_(intRadius),
      eLowB_(eLowB),
      eLowE_(eLowE),
      etLowB_(etLowB),
      etLowE_(etLowE),
      theCaloGeom_(theCaloGeom),
      mhbhe_(mhbhe) {
  //set up the geometry and selector
  const CaloGeometry* caloGeom = theCaloGeom_.product();
  doubleConeSel_ = new CaloDualConeSelector<HBHERecHit>(intRadius_, extRadius_, caloGeom, DetId::Hcal);
}

References doubleConeSel_, extRadius_, DetId::Hcal, intRadius_, edm::ESHandle< T >::product(), and theCaloGeom_.

~EgammaHcalIsolation()

EgammaHcalIsolation::~EgammaHcalIsolation

(

)

Definition at line 48 of file EgammaHcalIsolation.cc.

{ delete doubleConeSel_; }

References doubleConeSel_.

Member Function Documentation

getHcalESum() [1/4]

double EgammaHcalIsolation::getHcalESum ( const GlobalPoint &  pclu ) const

inline

Definition at line 54 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, AllDepths, &scaleToE); }

References AllDepths, getHcalSum(), and scaleToE().

getHcalESum() [2/4]

double EgammaHcalIsolation::getHcalESum ( const math::XYZPoint &  p ) const

inline

Definition at line 52 of file EgammaHcalIsolation.h.

{ return getHcalESum(GlobalPoint(p.x(), p.y(), p.z())); }

References getHcalESum(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by getHcalESum().

getHcalESum() [3/4]

double EgammaHcalIsolation::getHcalESum ( const reco::Candidate *  c ) const

inline

Definition at line 48 of file EgammaHcalIsolation.h.

{ return getHcalESum(c->get<reco::SuperClusterRef>().get()); }

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

Referenced by getHcalESum().

getHcalESum() [4/4]

double EgammaHcalIsolation::getHcalESum ( const reco::SuperCluster *  sc ) const

inline

Definition at line 50 of file EgammaHcalIsolation.h.

{ return getHcalESum(sc->position()); }

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

Referenced by getHcalESum().

getHcalESumDepth1() [1/4]

double EgammaHcalIsolation::getHcalESumDepth1 ( const GlobalPoint &  pclu ) const

inline

Definition at line 72 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, Depth1, &scaleToE); }

References Depth1, getHcalSum(), and scaleToE().

getHcalESumDepth1() [2/4]

double EgammaHcalIsolation::getHcalESumDepth1 ( const math::XYZPoint &  p ) const

inline

Definition at line 66 of file EgammaHcalIsolation.h.

                                                        {
    return getHcalESumDepth1(GlobalPoint(p.x(), p.y(), p.z()));
  }

References getHcalESumDepth1(), and AlCaHLTBitMon_ParallelJobs::p.

getHcalESumDepth1() [3/4]

double EgammaHcalIsolation::getHcalESumDepth1 ( const reco::Candidate *  c ) const

inline

Definition at line 58 of file EgammaHcalIsolation.h.

                                                         {
    return getHcalESumDepth1(c->get<reco::SuperClusterRef>().get());
  }

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

Referenced by getHcalESumDepth1().

getHcalESumDepth1() [4/4]

double EgammaHcalIsolation::getHcalESumDepth1 ( const reco::SuperCluster *  sc ) const

inline

Definition at line 64 of file EgammaHcalIsolation.h.

{ return getHcalESumDepth1(sc->position()); }

References getHcalESumDepth1(), and reco::CaloCluster::position().

Referenced by getHcalESumDepth1().

getHcalESumDepth2() [1/4]

double EgammaHcalIsolation::getHcalESumDepth2 ( const GlobalPoint &  pclu ) const

inline

Definition at line 90 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, Depth2, &scaleToE); }

References Depth2, getHcalSum(), and scaleToE().

getHcalESumDepth2() [2/4]

double EgammaHcalIsolation::getHcalESumDepth2 ( const math::XYZPoint &  p ) const

inline

Definition at line 84 of file EgammaHcalIsolation.h.

                                                        {
    return getHcalESumDepth2(GlobalPoint(p.x(), p.y(), p.z()));
  }

References getHcalESumDepth2(), and AlCaHLTBitMon_ParallelJobs::p.

getHcalESumDepth2() [3/4]

double EgammaHcalIsolation::getHcalESumDepth2 ( const reco::Candidate *  c ) const

inline

Definition at line 76 of file EgammaHcalIsolation.h.

                                                         {
    return getHcalESumDepth2(c->get<reco::SuperClusterRef>().get());
  }

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

Referenced by getHcalESumDepth2().

getHcalESumDepth2() [4/4]

double EgammaHcalIsolation::getHcalESumDepth2 ( const reco::SuperCluster *  sc ) const

inline

Definition at line 82 of file EgammaHcalIsolation.h.

{ return getHcalESumDepth2(sc->position()); }

References getHcalESumDepth2(), and reco::CaloCluster::position().

Referenced by getHcalESumDepth2().

getHcalEtSum() [1/4]

double EgammaHcalIsolation::getHcalEtSum ( const GlobalPoint &  pclu ) const

inline

Definition at line 55 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, AllDepths, &scaleToEt); }

References AllDepths, getHcalSum(), and scaleToEt().

getHcalEtSum() [2/4]

double EgammaHcalIsolation::getHcalEtSum ( const math::XYZPoint &  p ) const

inline

Definition at line 53 of file EgammaHcalIsolation.h.

{ return getHcalEtSum(GlobalPoint(p.x(), p.y(), p.z())); }

References getHcalEtSum(), and AlCaHLTBitMon_ParallelJobs::p.

Referenced by getHcalEtSum().

getHcalEtSum() [3/4]

double EgammaHcalIsolation::getHcalEtSum ( const reco::Candidate *  c ) const

inline

Definition at line 49 of file EgammaHcalIsolation.h.

{ return getHcalEtSum(c->get<reco::SuperClusterRef>().get()); }

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

Referenced by getHcalEtSum().

getHcalEtSum() [4/4]

double EgammaHcalIsolation::getHcalEtSum ( const reco::SuperCluster *  sc ) const

inline

Definition at line 51 of file EgammaHcalIsolation.h.

{ return getHcalEtSum(sc->position()); }

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

Referenced by getHcalEtSum().

getHcalEtSumDepth1() [1/4]

double EgammaHcalIsolation::getHcalEtSumDepth1 ( const GlobalPoint &  pclu ) const

inline

Definition at line 73 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, Depth1, &scaleToEt); }

References Depth1, getHcalSum(), and scaleToEt().

getHcalEtSumDepth1() [2/4]

double EgammaHcalIsolation::getHcalEtSumDepth1 ( const math::XYZPoint &  p ) const

inline

Definition at line 69 of file EgammaHcalIsolation.h.

                                                         {
    return getHcalEtSumDepth1(GlobalPoint(p.x(), p.y(), p.z()));
  }

References getHcalEtSumDepth1(), and AlCaHLTBitMon_ParallelJobs::p.

getHcalEtSumDepth1() [3/4]

double EgammaHcalIsolation::getHcalEtSumDepth1 ( const reco::Candidate *  c ) const

inline

Definition at line 61 of file EgammaHcalIsolation.h.

                                                          {
    return getHcalEtSumDepth1(c->get<reco::SuperClusterRef>().get());
  }

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

Referenced by getHcalEtSumDepth1().

getHcalEtSumDepth1() [4/4]

double EgammaHcalIsolation::getHcalEtSumDepth1 ( const reco::SuperCluster *  sc ) const

inline

Definition at line 65 of file EgammaHcalIsolation.h.

{ return getHcalEtSumDepth1(sc->position()); }

References getHcalEtSumDepth1(), and reco::CaloCluster::position().

Referenced by getHcalEtSumDepth1().

getHcalEtSumDepth2() [1/4]

double EgammaHcalIsolation::getHcalEtSumDepth2 ( const GlobalPoint &  pclu ) const

inline

Definition at line 91 of file EgammaHcalIsolation.h.

{ return getHcalSum(pclu, Depth2, &scaleToEt); }

References Depth2, getHcalSum(), and scaleToEt().

getHcalEtSumDepth2() [2/4]

double EgammaHcalIsolation::getHcalEtSumDepth2 ( const math::XYZPoint &  p ) const

inline

Definition at line 87 of file EgammaHcalIsolation.h.

                                                         {
    return getHcalEtSumDepth2(GlobalPoint(p.x(), p.y(), p.z()));
  }

References getHcalEtSumDepth2(), and AlCaHLTBitMon_ParallelJobs::p.

getHcalEtSumDepth2() [3/4]

double EgammaHcalIsolation::getHcalEtSumDepth2 ( const reco::Candidate *  c ) const

inline

Definition at line 79 of file EgammaHcalIsolation.h.

                                                          {
    return getHcalEtSumDepth2(c->get<reco::SuperClusterRef>().get());
  }

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

Referenced by getHcalEtSumDepth2().

getHcalEtSumDepth2() [4/4]

double EgammaHcalIsolation::getHcalEtSumDepth2 ( const reco::SuperCluster *  sc ) const

inline

Definition at line 83 of file EgammaHcalIsolation.h.

{ return getHcalEtSumDepth2(sc->position()); }

References getHcalEtSumDepth2(), and reco::CaloCluster::position().

Referenced by getHcalEtSumDepth2().

getHcalSum()

double EgammaHcalIsolation::getHcalSum ( const GlobalPoint &  pclu, const HcalDepth &  depth, double(*)(const double &)  scale  ) const

private

Definition at line 50 of file EgammaHcalIsolation.cc.

                                                                             {
  double sum = 0.;
  if (!mhbhe_.empty()) {
    //Compute the HCAL energy behind ECAL
    doubleConeSel_->selectCallback(pclu, mhbhe_, [this, &sum, &depth, &scale](const HBHERecHit& i) {
      double eta = theCaloGeom_.product()->getPosition(i.detid()).eta();
      HcalDetId hcalDetId(i.detid());
      if (hcalDetId.subdet() == HcalBarrel &&         //Is it in the barrel?
          i.energy() > eLowB_ &&                      //Does it pass the min energy?
          i.energy() * scaleToEt(eta) > etLowB_ &&    //Does it pass the min et?
          (depth == AllDepths || depth == Depth1)) {  //Are we asking for the first depth?
        sum += i.energy() * scale(eta);
      }
      if (hcalDetId.subdet() == HcalEndcap &&       //Is it in the endcap?
          i.energy() > eLowE_ &&                    //Does it pass the min energy?
          i.energy() * scaleToEt(eta) > etLowE_) {  //Does it pass the min et?
        switch (depth) {                            //Which depth?
          case AllDepths:
            sum += i.energy() * scale(eta);
            break;
          case Depth1:
            sum += (isDepth2(i.detid())) ? 0 : i.energy() * scale(eta);
            break;
          case Depth2:
            sum += (isDepth2(i.detid())) ? i.energy() * scale(eta) : 0;
            break;
        }
      }
    });
  }
 
  return sum;
}

References AllDepths, LEDCalibrationChannels::depth, Depth1, Depth2, doubleConeSel_, eLowB_, eLowE_, edm::SortedCollection< T, SORT >::empty(), PVValHelper::eta, etLowB_, etLowE_, CaloGeometry::getPosition(), HcalBarrel, HcalEndcap, mps_fire::i, isDepth2(), mhbhe_, edm::ESHandle< T >::product(), L1EGammaCrystalsEmulatorProducer_cfi::scale, scaleToEt(), CaloDualConeSelector< T >::selectCallback(), and theCaloGeom_.

Referenced by getHcalESum(), getHcalESumDepth1(), getHcalESumDepth2(), getHcalEtSum(), getHcalEtSumDepth1(), and getHcalEtSumDepth2().

isDepth2()

bool EgammaHcalIsolation::isDepth2 ( const DetId &  detId ) const

private

Definition at line 86 of file EgammaHcalIsolation.cc.

                                                           {
  if ((HcalDetId(detId).depth() == 2 && HcalDetId(detId).ietaAbs() >= 18 && HcalDetId(detId).ietaAbs() < 27) ||
      (HcalDetId(detId).depth() == 3 && HcalDetId(detId).ietaAbs() == 27)) {
    return true;
 
  } else {
    return false;
  }
}

References LEDCalibrationChannels::depth.

Referenced by getHcalSum().

Member Data Documentation

doubleConeSel_

CaloDualConeSelector<HBHERecHit>* EgammaHcalIsolation::doubleConeSel_

private

Definition at line 107 of file EgammaHcalIsolation.h.

Referenced by EgammaHcalIsolation(), getHcalSum(), and ~EgammaHcalIsolation().

eLowB_

double EgammaHcalIsolation::eLowB_

private

Definition at line 99 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

eLowE_

double EgammaHcalIsolation::eLowE_

private

Definition at line 100 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

etLowB_

double EgammaHcalIsolation::etLowB_

private

Definition at line 101 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

etLowE_

double EgammaHcalIsolation::etLowE_

private

Definition at line 102 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

extRadius_

double EgammaHcalIsolation::extRadius_

private

Definition at line 97 of file EgammaHcalIsolation.h.

Referenced by EgammaHcalIsolation().

intRadius_

double EgammaHcalIsolation::intRadius_

private

Definition at line 98 of file EgammaHcalIsolation.h.

Referenced by EgammaHcalIsolation().

mhbhe_

const HBHERecHitCollection& EgammaHcalIsolation::mhbhe_

private

Definition at line 105 of file EgammaHcalIsolation.h.

Referenced by getHcalSum().

theCaloGeom_

edm::ESHandle<CaloGeometry> EgammaHcalIsolation::theCaloGeom_

private

Definition at line 104 of file EgammaHcalIsolation.h.

Referenced by EgammaHcalIsolation(), and getHcalSum().