egammaisolation::EgammaTowerExtractor Class Reference

#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::Track &track) const override
reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &c) 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::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
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

enum egammaisolation::EgammaTowerExtractor::HcalDepth

Enumerator
AllDepths
Undefined
Depth1
Depth2

Definition at line 36 of file EgammaTowerExtractor.h.

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

Constructor & Destructor Documentation

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

inline

Definition at line 39 of file EgammaTowerExtractor.h.

: EgammaTowerExtractor(par, iC) {}

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

inline

Definition at line 40 of file EgammaTowerExtractor.h.

References AllDepths, Depth1, Depth2, depth_, Exception, extRadius2_, and ~EgammaTowerExtractor().

        : 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;
      }
    }

EgammaTowerExtractor::~EgammaTowerExtractor ( )

override

Definition at line 30 of file EgammaTowerExtractor.cc.

Referenced by EgammaTowerExtractor().

{}

Member Function Documentation

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.

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

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

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.

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

                                                                                    {
  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;
}

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.

                                                                                                                 {
    }

Member Data Documentation

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

private

Definition at line 74 of file EgammaTowerExtractor.h.

int egammaisolation::EgammaTowerExtractor::depth_

private

Definition at line 75 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

double egammaisolation::EgammaTowerExtractor::etLow_

private

Definition at line 72 of file EgammaTowerExtractor.h.

double egammaisolation::EgammaTowerExtractor::extRadius2_

private

Definition at line 70 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

double egammaisolation::EgammaTowerExtractor::intRadius_

private

Definition at line 71 of file EgammaTowerExtractor.h.