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_
bool	debug_
std::vector< double >	dEdXWeights_
edm::InputTag	eeRecHitInputTag_
edm::InputTag	fhRecHitInputTag_
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_, debug_, dEdXWeights_, eeRecHitInputTag_, fhRecHitInputTag_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), 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") ),
        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_);
    pcaHelper_.setdEdXWeights(dEdXWeights_);
    debug_ = iConfig.getUntrackedParameter<bool>("debug", false);
}

HGCalEgammaIDHelper::~HGCalEgammaIDHelper ( )

inline

Definition at line 35 of file HGCalEgammaIDHelper.h.

References computeHGCAL(), eventInit(), iEvent, and TCMET_cfi::radius.

{}

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 37 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().

Referenced by ~HGCalEgammaIDHelper().

                                                                                 {
    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 22 of file HGCalEgammaIDHelper.cc.

References hgcal::EGammaPCAHelper::fillHitMap(), edm::Event::getByToken(), hgcal::RecHitTools::getEventSetup(), isoHelper_, pcaHelper_, recHitsBH_, recHitsEE_, recHitsFH_, recHitTools_, HGCalIsoCalculator::setRecHits(), HGCalIsoCalculator::setRecHitTools(), and hgcal::EGammaPCAHelper::setRecHitTools().

Referenced by ~HGCalEgammaIDHelper().

                                                                                      {
    edm::Handle<HGCRecHitCollection> recHitHandleEE;
    iEvent.getByToken(recHitsEE_, recHitHandleEE);
    edm::Handle<HGCRecHitCollection> recHitHandleFH;
    iEvent.getByToken(recHitsFH_, recHitHandleFH);
    edm::Handle<HGCRecHitCollection> recHitHandleBH;
    iEvent.getByToken(recHitsBH_, recHitHandleBH);

    recHitTools_.getEventSetup(iSetup);
    pcaHelper_.setRecHitTools(&recHitTools_);
    isoHelper_.setRecHitTools(&recHitTools_);
    pcaHelper_.fillHitMap(*recHitHandleEE,*recHitHandleFH,*recHitHandleBH);
    isoHelper_.setRecHits(recHitHandleEE, recHitHandleFH, recHitHandleBH);
}

float HGCalEgammaIDHelper::getIsolationRing ( unsigned int  ring ) const

inline

Definition at line 64 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().

bool HGCalEgammaIDHelper::debug_

private

Definition at line 83 of file HGCalEgammaIDHelper.h.

Referenced by computeHGCAL(), and HGCalEgammaIDHelper().

std::vector<double> HGCalEgammaIDHelper::dEdXWeights_

private

Definition at line 76 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().

HGCalIsoCalculator HGCalEgammaIDHelper::isoHelper_

private

Definition at line 78 of file HGCalEgammaIDHelper.h.

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

hgcal::EGammaPCAHelper HGCalEgammaIDHelper::pcaHelper_

private

Definition at line 77 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 81 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

edm::EDGetTokenT<HGCRecHitCollection> HGCalEgammaIDHelper::recHitsEE_

private

Definition at line 79 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

edm::EDGetTokenT<HGCRecHitCollection> HGCalEgammaIDHelper::recHitsFH_

private

Definition at line 80 of file HGCalEgammaIDHelper.h.

Referenced by eventInit(), and HGCalEgammaIDHelper().

hgcal::RecHitTools HGCalEgammaIDHelper::recHitTools_

private

Definition at line 82 of file HGCalEgammaIDHelper.h.

Referenced by eventInit().