#include <EgammaHcalExtractor.h>

Inheritance diagram for egammaisolation::EgammaHcalExtractor:

Public Member Functions
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const override

reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &c) const override

	EgammaHcalExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

	EgammaHcalExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &iC)

void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks) override

	~EgammaHcalExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const
	make single IsoDeposit based on trackRef as input More...

virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const
	make single IsoDeposit based on a CandidateBaseRef as input More...

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackBaseRef &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &track) const

virtual std::vector< reco::IsoDeposit >	deposits (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::CandidateBaseRef &track) const

virtual void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Attributes
double	etLow_

double	extRadius_

edm::EDGetTokenT< HBHERecHitCollection >	hcalRecHitProducerToken_

double	intRadius_

Detailed Description

Definition at line 36 of file EgammaHcalExtractor.h.

Constructor & Destructor Documentation

egammaisolation::EgammaHcalExtractor::EgammaHcalExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &&	iC
	)

inline

Definition at line 38 of file EgammaHcalExtractor.h.

References ~EgammaHcalExtractor().

38 : EgammaHcalExtractor(par, iC) {}

egammaisolation::EgammaHcalExtractor::EgammaHcalExtractor

EgammaHcalExtractor(const edm::ParameterSet &par, edm::ConsumesCollector &&iC)

Definition: EgammaHcalExtractor.h:38

EgammaHcalExtractor::EgammaHcalExtractor	(	const edm::ParameterSet &	par,
		edm::ConsumesCollector &	iC
	)

Definition at line 35 of file EgammaHcalExtractor.cc.

     : extRadius_(par.getParameter<double>("extRadius")),
       intRadius_(par.getParameter<double>("intRadius")),
       etLow_(par.getParameter<double>("etMin")),
       hcalRecHitProducerToken_(iC.consumes<HBHERecHitCollection>(par.getParameter<edm::InputTag>("hcalRecHits"))) {}

EgammaHcalExtractor::~EgammaHcalExtractor ( )

override

Definition at line 41 of file EgammaHcalExtractor.cc.

Referenced by EgammaHcalExtractor().

41 {}

Member Function Documentation

reco::IsoDeposit egammaisolation::EgammaHcalExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Track &	track
	)		const

inlineoverridevirtual

make single IsoDeposit based on track as input purely virtual: have to implement in concrete implementations

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 45 of file EgammaHcalExtractor.h.

References HltBtagPostValidation_cff::c, ev, Exception, and Skims_PA_cff::name.

Referenced by deposit().

                                                                     {
       throw cms::Exception("Configuration Error")
           << "This extractor " << (typeid(this).name()) << " is not made for tracks";
     }

reco::IsoDeposit EgammaHcalExtractor::deposit	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::Candidate &	track
	)		const

overridevirtual

make single IsoDeposit based on a candidate as input purely virtual: have to implement in concrete implementations

Reimplemented from reco::isodeposit::IsoDepositExtractor.

Definition at line 43 of file EgammaHcalExtractor.cc.

References deposit(), CaloRecHit::detid(), CaloRecHit::energy(), PV3DBase< T, PVType, FrameType >::eta(), etLow_, JetChargeProducer_cfi::exp, extRadius_, edm::EventSetup::get(), reco::Candidate::get(), edm::Event::getByToken(), CaloGeometry::getPosition(), DetId::Hcal, hcalRecHitProducerToken_, mps_fire::i, intRadius_, PV3DBase< T, PVType, FrameType >::phi(), point, edm::ESHandle< T >::product(), SimDataFormats::CaloAnalysis::sc, CaloDualConeSelector< T >::selectCallback(), and funct::sin().

                                                                                    {
   //Get MetaRecHit collection
   edm::Handle<HBHERecHitCollection> hcalRecHitHandle;
   iEvent.getByToken(hcalRecHitProducerToken_, hcalRecHitHandle);
 
   //Get Calo Geometry
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* caloGeom = pG.product();
   CaloDualConeSelector<HBHERecHit> coneSel(intRadius_, extRadius_, caloGeom, DetId::Hcal);
 
   //Take the SC position
   reco::SuperClusterRef sc = emObject.get<reco::SuperClusterRef>();
   math::XYZPoint caloPosition = sc->position();
   GlobalPoint point(caloPosition.x(), caloPosition.y(), caloPosition.z());
   // needed: coneSel.select(eta,phi,hits) is not the same!
 
   Direction candDir(caloPosition.eta(), caloPosition.phi());
   reco::IsoDeposit deposit(candDir);
   deposit.setVeto(reco::IsoDeposit::Veto(candDir, intRadius_));
   double sinTheta = sin(2 * atan(exp(-sc->eta())));
   deposit.addCandEnergy(sc->energy() * sinTheta);
 
   //Compute the HCAL energy behind ECAL
   coneSel.selectCallback(point, *hcalRecHitHandle, [&](const HBHERecHit& i) {
     const GlobalPoint& hcalHit_position = caloGeom->getPosition(i.detid());
     double hcalHit_eta = hcalHit_position.eta();
     double hcalHit_Et = i.energy() * sin(2 * atan(exp(-hcalHit_eta)));
     if (hcalHit_Et > etLow_) {
       deposit.addDeposit(Direction(hcalHit_eta, hcalHit_position.phi()), hcalHit_Et);
     }
   });
 
   return deposit;
 }

void egammaisolation::EgammaHcalExtractor::fillVetos	(	const edm::Event &	ev,
		const edm::EventSetup &	evSetup,
		const reco::TrackCollection &	tracks
	)

inlineoverridevirtual

fill vetoes: to exclude deposits at IsoDeposit creation stage check concrete extractors if it's no-op !

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 43 of file EgammaHcalExtractor.h.

43 {

44 }

Member Data Documentation

double egammaisolation::EgammaHcalExtractor::etLow_

private

Definition at line 58 of file EgammaHcalExtractor.h.

Referenced by deposit().

double egammaisolation::EgammaHcalExtractor::extRadius_

private

Definition at line 56 of file EgammaHcalExtractor.h.

Referenced by deposit().

edm::EDGetTokenT<HBHERecHitCollection> egammaisolation::EgammaHcalExtractor::hcalRecHitProducerToken_

private

Definition at line 60 of file EgammaHcalExtractor.h.

Referenced by deposit().

double egammaisolation::EgammaHcalExtractor::intRadius_

private

Definition at line 57 of file EgammaHcalExtractor.h.

Referenced by deposit().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation