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
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Track &track) const
virtual reco::IsoDeposit	deposit (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::Candidate &c) const
	EgammaTowerExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &&iC)
	EgammaTowerExtractor (const edm::ParameterSet &par, edm::ConsumesCollector &iC)
virtual void	fillVetos (const edm::Event &ev, const edm::EventSetup &evSetup, const reco::TrackCollection &tracks)
virtual	~EgammaTowerExtractor ()
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 35 of file EgammaTowerExtractor.h.

Member Enumeration Documentation

enum egammaisolation::EgammaTowerExtractor::HcalDepth

Enumerator
AllDepths
Undefined
Depth1
Depth2

Definition at line 38 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 42 of file EgammaTowerExtractor.h.

                                                                                         :
           EgammaTowerExtractor(par, iC) {}

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

inline

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

virtual

Definition at line 31 of file EgammaTowerExtractor.cc.

Referenced by EgammaTowerExtractor().

{}

Member Function Documentation

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

inlinevirtual

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

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 63 of file EgammaTowerExtractor.h.

References EnergyCorrector::c, ev, Exception, and dataset::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

virtual

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

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

inlinevirtual

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

Implements reco::isodeposit::IsoDepositExtractor.

Definition at line 61 of file EgammaTowerExtractor.h.

                                                                     { }

Member Data Documentation

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

private

Definition at line 77 of file EgammaTowerExtractor.h.

int egammaisolation::EgammaTowerExtractor::depth_

private

Definition at line 78 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

double egammaisolation::EgammaTowerExtractor::etLow_

private

Definition at line 74 of file EgammaTowerExtractor.h.

double egammaisolation::EgammaTowerExtractor::extRadius2_

private

Definition at line 72 of file EgammaTowerExtractor.h.

Referenced by EgammaTowerExtractor().

double egammaisolation::EgammaTowerExtractor::intRadius_

private

Definition at line 73 of file EgammaTowerExtractor.h.