HGCalEgammaIDHelper Class Reference

#include <HGCalEgammaIDHelper.h>

Public Member Functions
const math::XYZVector &	axis () const
const math::XYZPoint &	barycenter () const
float	clusterDepthCompatibility (const hgcal::LongDeps &ld, float &measDepth, float &expDepth, float &expSigma)
void	computeHGCAL (const reco::Photon &thePhoton, float radius)
void	computeHGCAL (const reco::GsfElectron &theElectron, float radius)
const TVectorD &	eigenValues () const
hgcal::LongDeps	energyPerLayer (float radius, bool withHalo=true)
void	eventInit (const edm::Event &iEvent, const edm::EventSetup &iSetup)
float	getIsolationRing (unsigned int ring) const
	HGCalEgammaIDHelper ()
	HGCalEgammaIDHelper (const edm::ParameterSet &, edm::ConsumesCollector &&iC)
const hgcal::EGammaPCAHelper *	pcaHelper () const
void	printHits (float radius) const
double	sigmaEE () const
double	sigmaPP () const
const TVectorD &	sigmas () const
double	sigmaUU () const
double	sigmaVV () const
	~HGCalEgammaIDHelper ()

Private Attributes
edm::InputTag	bhRecHitInputTag_
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometry_
bool	debug_
std::vector< double >	dEdXWeights_
edm::InputTag	eeRecHitInputTag_
edm::InputTag	fhRecHitInputTag_
edm::EDGetTokenT < std::unordered_map< DetId, const HGCRecHit * > >	hitMap_
edm::InputTag	hitMapInputTag_
HGCalIsoCalculator	isoHelper_
hgcal::EGammaPCAHelper	pcaHelper_
edm::EDGetTokenT < HGCRecHitCollection >	recHitsBH_
edm::EDGetTokenT < HGCRecHitCollection >	recHitsEE_
edm::EDGetTokenT < HGCRecHitCollection >	recHitsFH_
hgcal::RecHitTools	recHitTools_

Detailed Description

Definition at line 31 of file HGCalEgammaIDHelper.h.

Constructor & Destructor Documentation

HGCalEgammaIDHelper::HGCalEgammaIDHelper ( )

inline

Definition at line 33 of file HGCalEgammaIDHelper.h.

{}

HGCalEgammaIDHelper::HGCalEgammaIDHelper	(	const edm::ParameterSet &	iConfig,
		edm::ConsumesCollector &&	iC
	)

Definition at line 5 of file HGCalEgammaIDHelper.cc.

References bhRecHitInputTag_, caloGeometry_, debug_, dEdXWeights_, eeRecHitInputTag_, fhRecHitInputTag_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hitMap_, hitMapInputTag_, isoHelper_, pcaHelper_, recHitsBH_, recHitsEE_, recHitsFH_, hgcal::EGammaPCAHelper::setdEdXWeights(), HGCalIsoCalculator::setDeltaR(), HGCalIsoCalculator::setMinDeltaR(), and HGCalIsoCalculator::setNRings().

    : eeRecHitInputTag_(iConfig.getParameter<edm::InputTag>("EERecHits")),
      fhRecHitInputTag_(iConfig.getParameter<edm::InputTag>("FHRecHits")),
      bhRecHitInputTag_(iConfig.getParameter<edm::InputTag>("BHRecHits")),
      hitMapInputTag_(iConfig.getParameter<edm::InputTag>("hitMapTag")),
      dEdXWeights_(iConfig.getParameter<std::vector<double>>("dEdXWeights")) {
  isoHelper_.setDeltaR(iConfig.getParameter<double>("isoDeltaR"));
  isoHelper_.setNRings(iConfig.getParameter<unsigned int>("isoNRings"));
  isoHelper_.setMinDeltaR(iConfig.getParameter<double>("isoDeltaRmin"));

  recHitsEE_ = iC.consumes<HGCRecHitCollection>(eeRecHitInputTag_);
  recHitsFH_ = iC.consumes<HGCRecHitCollection>(fhRecHitInputTag_);
  recHitsBH_ = iC.consumes<HGCRecHitCollection>(bhRecHitInputTag_);
  hitMap_ = iC.consumes<std::unordered_map<DetId, const HGCRecHit*>>(hitMapInputTag_);
  caloGeometry_ = iC.esConsumes();
  pcaHelper_.setdEdXWeights(dEdXWeights_);
  debug_ = iConfig.getUntrackedParameter<bool>("debug", false);
}

HGCalEgammaIDHelper::~HGCalEgammaIDHelper ( )

inline

Definition at line 35 of file HGCalEgammaIDHelper.h.

{}

Member Function Documentation

const math::XYZVector& HGCalEgammaIDHelper::axis ( ) const

inline

Definition at line 52 of file HGCalEgammaIDHelper.h.

References hgcal::EGammaPCAHelper::axis(), and pcaHelper_.

{ return pcaHelper_.axis(); }

const math::XYZPoint& HGCalEgammaIDHelper::barycenter ( ) const

inline

Definition at line 51 of file HGCalEgammaIDHelper.h.

References hgcal::EGammaPCAHelper::barycenter(), and pcaHelper_.

{ return pcaHelper_.barycenter(); }

float HGCalEgammaIDHelper::clusterDepthCompatibility ( const hgcal::LongDeps &  ld, float &  measDepth, float &  expDepth, float &  expSigma  )

inline

Definition at line 60 of file HGCalEgammaIDHelper.h.

References hgcal::EGammaPCAHelper::clusterDepthCompatibility(), and pcaHelper_.

                                                                                                               {
    return pcaHelper_.clusterDepthCompatibility(ld, measDepth, expDepth, expSigma);
  }

void HGCalEgammaIDHelper::computeHGCAL	(	const reco::Photon &	thePhoton,
		float	radius
	)

Definition at line 36 of file HGCalEgammaIDHelper.cc.

References hgcal::EGammaPCAHelper::clear(), hgcal::EGammaPCAHelper::computePCA(), hgcal::EGammaPCAHelper::computeShowerWidth(), gather_cfg::cout, debug_, reco::Photon::isEB(), isoHelper_, pcaHelper_, hgcal::EGammaPCAHelper::pcaInitialComputation(), HGCalIsoCalculator::produceHGCalIso(), hgcal::EGammaPCAHelper::storeRecHits(), and reco::Photon::superCluster().

                                                                                {
  if (thePhoton.isEB()) {
    if (debug_)
      std::cout << "The photon is in the barrel" << std::endl;
    pcaHelper_.clear();
    return;
  }

  pcaHelper_.storeRecHits(*thePhoton.superCluster()->seed());
  if (debug_)
    std::cout << " Stored the hits belonging to the photon superCluster seed " << std::endl;

  // initial computation, no radius cut, but halo hits not taken
  if (debug_)
    std::cout << " Calling PCA initial computation" << std::endl;
  pcaHelper_.pcaInitialComputation();
  // first computation within cylinder, halo hits included
  pcaHelper_.computePCA(radius);
  // second computation within cylinder, halo hits included
  pcaHelper_.computePCA(radius);
  pcaHelper_.computeShowerWidth(radius);

  // isolation
  isoHelper_.produceHGCalIso(thePhoton.superCluster()->seed());
}

void HGCalEgammaIDHelper::computeHGCAL	(	const reco::GsfElectron &	theElectron,
		float	radius
	)

Definition at line 62 of file HGCalEgammaIDHelper.cc.

References hgcal::EGammaPCAHelper::clear(), hgcal::EGammaPCAHelper::computePCA(), hgcal::EGammaPCAHelper::computeShowerWidth(), gather_cfg::cout, debug_, reco::GsfElectron::electronCluster(), reco::GsfElectron::isEB(), isoHelper_, pcaHelper_, hgcal::EGammaPCAHelper::pcaInitialComputation(), HGCalIsoCalculator::produceHGCalIso(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                                                       {
  if (theElectron.isEB()) {
    if (debug_)
      std::cout << "The electron is in the barrel" << std::endl;
    pcaHelper_.clear();
    return;
  }

  pcaHelper_.storeRecHits(*theElectron.electronCluster());
  if (debug_)
    std::cout << " Stored the hits belonging to the electronCluster " << std::endl;

  // initial computation, no radius cut, but halo hits not taken
  if (debug_)
    std::cout << " Calling PCA initial computation" << std::endl;
  pcaHelper_.pcaInitialComputation();
  // first computation within cylinder, halo hits included
  pcaHelper_.computePCA(radius);
  // second computation within cylinder, halo hits included
  pcaHelper_.computePCA(radius);
  pcaHelper_.computeShowerWidth(radius);
  isoHelper_.produceHGCalIso(theElectron.electronCluster());
}

const TVectorD& HGCalEgammaIDHelper::eigenValues ( ) const

inline

Definition at line 49 of file HGCalEgammaIDHelper.h.

References hgcal::EGammaPCAHelper::eigenValues(), and pcaHelper_.

{ return pcaHelper_.eigenValues(); }

hgcal::LongDeps HGCalEgammaIDHelper::energyPerLayer ( float  radius, bool  withHalo = true  )

inline

Definition at line 55 of file HGCalEgammaIDHelper.h.

References hgcal::EGammaPCAHelper::energyPerLayer(), and pcaHelper_.

                                                                   {
    return pcaHelper_.energyPerLayer(radius, withHalo);
  }

void HGCalEgammaIDHelper::eventInit	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 24 of file HGCalEgammaIDHelper.cc.

References caloGeometry_, edm::Event::get(), edm::EventSetup::getData(), edm::Event::getHandle(), hitMap_, isoHelper_, pcaHelper_, recHitsBH_, recHitsEE_, recHitsFH_, recHitTools_, hgcal::RecHitTools::setGeometry(), hgcal::EGammaPCAHelper::setHitMap(), HGCalIsoCalculator::setRecHits(), hgcal::EGammaPCAHelper::setRecHitTools(), and HGCalIsoCalculator::setRecHitTools().

                                                                                       {
  auto recHitHandleEE = iEvent.getHandle(recHitsEE_);
  auto recHitHandleFH = iEvent.getHandle(recHitsFH_);
  auto recHitHandleBH = iEvent.getHandle(recHitsBH_);

  recHitTools_.setGeometry(iSetup.getData(caloGeometry_));
  pcaHelper_.setRecHitTools(&recHitTools_);
  isoHelper_.setRecHitTools(&recHitTools_);
  pcaHelper_.setHitMap(&iEvent.get(hitMap_));
  isoHelper_.setRecHits(recHitHandleEE, recHitHandleFH, recHitHandleBH);
}

float HGCalEgammaIDHelper::getIsolationRing ( unsigned int  ring ) const

inline

Definition at line 65 of file HGCalEgammaIDHelper.h.

References HGCalIsoCalculator::getIso(), and isoHelper_.

{ return isoHelper_.getIso(ring); };

const hgcal::EGammaPCAHelper* HGCalEgammaIDHelper::pcaHelper ( ) const

inline

Definition at line 69 of file HGCalEgammaIDHelper.h.

References pcaHelper_.

{ return &pcaHelper_; }

void HGCalEgammaIDHelper::printHits ( float  radius ) const

inline

Definition at line 68 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::printHits().

Referenced by ntuplePrintersDiff.TrackingParticlePrinter::__call__(), ntuplePrintersDiff.SeedPrinter::diff(), ntuplePrintersDiff.TrackPrinter::diff(), ntuplePrintersDiff.TrackingParticlePrinter::diff(), ntuplePrintersDiff.SeedPrinter::printSeed(), and ntuplePrintersDiff.TrackPrinter::printTrack().

{ pcaHelper_.printHits(radius); }

double HGCalEgammaIDHelper::sigmaEE ( ) const

inline

Definition at line 47 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::sigmaEE().

{ return pcaHelper_.sigmaEE(); }

double HGCalEgammaIDHelper::sigmaPP ( ) const

inline

Definition at line 48 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::sigmaPP().

{ return pcaHelper_.sigmaPP(); }

const TVectorD& HGCalEgammaIDHelper::sigmas ( ) const

inline

Definition at line 50 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::sigmas().

{ return pcaHelper_.sigmas(); }

double HGCalEgammaIDHelper::sigmaUU ( ) const

inline

Definition at line 45 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::sigmaUU().

{ return pcaHelper_.sigmaUU(); }

double HGCalEgammaIDHelper::sigmaVV ( ) const

inline

Definition at line 46 of file HGCalEgammaIDHelper.h.

References pcaHelper_, and hgcal::EGammaPCAHelper::sigmaVV().

{ return pcaHelper_.sigmaVV(); }

Member Data Documentation

edm::InputTag HGCalEgammaIDHelper::bhRecHitInputTag_

private

Definition at line 74 of file HGCalEgammaIDHelper.h.

Referenced by HGCalEgammaIDHelper().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HGCalEgammaIDHelper::caloGeometry_

private

Definition at line 84 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

bool HGCalEgammaIDHelper::debug_

private

Definition at line 86 of file HGCalEgammaIDHelper.h.

Referenced by computeHGCAL(), and HGCalEgammaIDHelper().

std::vector<double> HGCalEgammaIDHelper::dEdXWeights_

private

Definition at line 77 of file HGCalEgammaIDHelper.h.

Referenced by HGCalEgammaIDHelper().

edm::InputTag HGCalEgammaIDHelper::eeRecHitInputTag_

private

Definition at line 72 of file HGCalEgammaIDHelper.h.

Referenced by HGCalEgammaIDHelper().

edm::InputTag HGCalEgammaIDHelper::fhRecHitInputTag_

private

Definition at line 73 of file HGCalEgammaIDHelper.h.

Referenced by HGCalEgammaIDHelper().

edm::EDGetTokenT<std::unordered_map<DetId, const HGCRecHit*> > HGCalEgammaIDHelper::hitMap_

private

Definition at line 83 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

edm::InputTag HGCalEgammaIDHelper::hitMapInputTag_

private

Definition at line 75 of file HGCalEgammaIDHelper.h.

Referenced by HGCalEgammaIDHelper().

HGCalIsoCalculator HGCalEgammaIDHelper::isoHelper_

private

Definition at line 79 of file HGCalEgammaIDHelper.h.

Referenced by computeHGCAL(), eventInit(), getIsolationRing(), and HGCalEgammaIDHelper().

hgcal::EGammaPCAHelper HGCalEgammaIDHelper::pcaHelper_

private

Definition at line 78 of file HGCalEgammaIDHelper.h.

Referenced by axis(), barycenter(), clusterDepthCompatibility(), computeHGCAL(), eigenValues(), energyPerLayer(), eventInit(), HGCalEgammaIDHelper(), pcaHelper(), printHits(), sigmaEE(), sigmaPP(), sigmas(), sigmaUU(), and sigmaVV().

edm::EDGetTokenT<HGCRecHitCollection> HGCalEgammaIDHelper::recHitsBH_

private

Definition at line 82 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

edm::EDGetTokenT<HGCRecHitCollection> HGCalEgammaIDHelper::recHitsEE_

private

Definition at line 80 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

edm::EDGetTokenT<HGCRecHitCollection> HGCalEgammaIDHelper::recHitsFH_

private

Definition at line 81 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

hgcal::RecHitTools HGCalEgammaIDHelper::recHitTools_

private

Definition at line 85 of file HGCalEgammaIDHelper.h.

Referenced by eventInit().