#include <EgammaTowerExtractor.h>

Inheritance diagram for egammaisolation::EgammaTowerExtractor:

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

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

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

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

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

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

	~EgammaTowerExtractor () override

Public Member Functions inherited from reco::isodeposit::IsoDepositExtractor
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 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 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 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 void	initEvent (const edm::Event &, const edm::EventSetup &)
	perform some constly initialization More...

virtual	~IsoDepositExtractor ()
	Destructor. More...

Private Attributes
edm::EDGetTokenT< CaloTowerCollection >	caloTowerToken

int	depth_

double	etLow_

double	extRadius2_

double	intRadius_

Detailed Description

Definition at line 34 of file EgammaTowerExtractor.h.

Member Enumeration Documentation

◆ HcalDepth

enum egammaisolation::EgammaTowerExtractor::HcalDepth

Enumerator
AllDepths
Undefined
Depth1
Depth2

Definition at line 36 of file EgammaTowerExtractor.h.

36 { AllDepths = -1, Undefined = 0, Depth1 = 1, Depth2 = 2 };

Constructor & Destructor Documentation

◆ EgammaTowerExtractor() [1/2]

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

inline

Definition at line 39 of file EgammaTowerExtractor.h.

39 : EgammaTowerExtractor(par, iC) {}

◆ EgammaTowerExtractor() [2/2]

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

inline

Definition at line 40 of file EgammaTowerExtractor.h.

         : extRadius2_(par.getParameter<double>("extRadius")),
           intRadius_(par.getParameter<double>("intRadius")),
           etLow_(par.getParameter<double>("etMin")),
           caloTowerToken(iC.consumes<CaloTowerCollection>(par.getParameter<edm::InputTag>("caloTowers"))),
           depth_(par.getParameter<int>("hcalDepth")) {
       extRadius2_ *= extRadius2_;
       //lets just check we have a valid depth
       //should we throw an exception or just warn and then fail gracefully later?
       if (depth_ != AllDepths && depth_ != Depth1 && depth_ != Depth2) {
         throw cms::Exception("Configuration Error")
             << "hcalDepth passed to EgammaTowerExtractor is invalid " << std::endl;
       }
     }

References AllDepths, Depth1, Depth2, depth_, Exception, and extRadius2_.

◆ ~EgammaTowerExtractor()

EgammaTowerExtractor::~EgammaTowerExtractor ( )

override

Definition at line 30 of file EgammaTowerExtractor.cc.

30 {}

Member Function Documentation

◆ deposit() [1/2]

reco::IsoDeposit EgammaTowerExtractor::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 32 of file EgammaTowerExtractor.cc.

                                                                                     {
   edm::Handle<CaloTowerCollection> towercollectionH;
   iEvent.getByToken(caloTowerToken, towercollectionH);
  
   //Take the SC position
   reco::SuperClusterRef sc = emObject.get<reco::SuperClusterRef>();
   math::XYZPoint caloPosition = sc->position();
  
   Direction candDir(caloPosition.eta(), caloPosition.phi());
   reco::IsoDeposit deposit(candDir);
   deposit.setVeto(reco::IsoDeposit::Veto(candDir, intRadius_));
   deposit.addCandEnergy(sc->energy() * sin(2 * atan(exp(-sc->eta()))));
  
   //loop over tracks
   for (CaloTowerCollection::const_iterator trItr = towercollectionH->begin(), trEnd = towercollectionH->end();
        trItr != trEnd;
        ++trItr) {
     double depEt = 0;
     //the hcal can be seperated into different depths
     //currently it is setup to check that the depth is valid in constructor
     //if the depth is not valid it fails gracefully
     //small bug fix, hadEnergyHeInnerLater returns zero for towers which are only depth 1
     //but we want Depth1 isolation to include these so we have to manually check for this
     if (depth_ == AllDepths)
       depEt = trItr->hadEt();
     else if (depth_ == Depth1)
       depEt = trItr->ietaAbs() < 18 || trItr->ietaAbs() > 29
                   ? trItr->hadEt()
                   : trItr->hadEnergyHeInnerLayer() * sin(trItr->p4().theta());
     else if (depth_ == Depth2)
       depEt = trItr->hadEnergyHeOuterLayer() * sin(trItr->p4().theta());
  
     if (depEt < etLow_)
       continue;
  
     Direction towerDir(trItr->eta(), trItr->phi());
     double dR2 = candDir.deltaR2(towerDir);
  
     if (dR2 < extRadius2_) {
       deposit.addDeposit(towerDir, depEt);
     }
  
   }  //end loop over tracks
  
   return deposit;
 }

References edm::SortedCollection< T, SORT >::begin(), reco::isodeposit::Direction::deltaR2(), edm::SortedCollection< T, SORT >::end(), JetChargeProducer_cfi::exp, reco::Candidate::get(), iEvent, and funct::sin().

◆ deposit() [2/2]

reco::IsoDeposit egammaisolation::EgammaTowerExtractor::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 59 of file EgammaTowerExtractor.h.

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

References Exception, and Skims_PA_cff::name.

◆ fillVetos()

void egammaisolation::EgammaTowerExtractor::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 57 of file EgammaTowerExtractor.h.

57 {

58 }

Member Data Documentation

◆ caloTowerToken

edm::EDGetTokenT<CaloTowerCollection> egammaisolation::EgammaTowerExtractor::caloTowerToken

private

Definition at line 74 of file EgammaTowerExtractor.h.

◆ depth_

int egammaisolation::EgammaTowerExtractor::depth_

private

Definition at line 75 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

◆ etLow_

double egammaisolation::EgammaTowerExtractor::etLow_

private

Definition at line 72 of file EgammaTowerExtractor.h.

◆ extRadius2_

double egammaisolation::EgammaTowerExtractor::extRadius2_

private

Definition at line 70 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

◆ intRadius_

double egammaisolation::EgammaTowerExtractor::intRadius_

private

Definition at line 71 of file EgammaTowerExtractor.h.

Public Types

Public Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ HcalDepth

Constructor & Destructor Documentation

◆ EgammaTowerExtractor() [1/2]

◆ EgammaTowerExtractor() [2/2]

◆ ~EgammaTowerExtractor()

Member Function Documentation

◆ deposit() [1/2]

◆ deposit() [2/2]

◆ fillVetos()

Member Data Documentation

◆ caloTowerToken

◆ depth_

◆ etLow_

◆ extRadius2_

◆ intRadius_