egammaisolation::EgammaHcalExtractor Class Reference

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::ESGetToken< CaloGeometry, CaloGeometryRecord >	geometryToken_
edm::EDGetTokenT < HBHERecHitCollection >	hcalRecHitProducerToken_
double	intRadius_

Detailed Description

Definition at line 41 of file EgammaHcalExtractor.cc.

Constructor & Destructor Documentation

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

inline

Definition at line 43 of file EgammaHcalExtractor.cc.

: EgammaHcalExtractor(par, iC) {}

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

Definition at line 79 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"))),
      geometryToken_(iC.esConsumes()) {}

EgammaHcalExtractor::~EgammaHcalExtractor ( )

override

Definition at line 86 of file EgammaHcalExtractor.cc.

{}

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 50 of file EgammaHcalExtractor.cc.

References Exception, and mergeVDriftHistosByStation::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 88 of file EgammaHcalExtractor.cc.

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

                                                                                   {
  //Get MetaRecHit collection
  auto const& hcalRecHits = iEvent.get(hcalRecHitProducerToken_);

  //Get Calo Geometry
  const CaloGeometry* caloGeom = &iSetup.getData(geometryToken_);
  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, hcalRecHits, [&](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 48 of file EgammaHcalExtractor.cc.

                                                                                                                 {
    }

Member Data Documentation

double egammaisolation::EgammaHcalExtractor::etLow_

private

Definition at line 63 of file EgammaHcalExtractor.cc.

Referenced by deposit().

double egammaisolation::EgammaHcalExtractor::extRadius_

private

Definition at line 61 of file EgammaHcalExtractor.cc.

Referenced by deposit().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> egammaisolation::EgammaHcalExtractor::geometryToken_

private

Definition at line 66 of file EgammaHcalExtractor.cc.

Referenced by deposit().

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

private

Definition at line 65 of file EgammaHcalExtractor.cc.

Referenced by deposit().

double egammaisolation::EgammaHcalExtractor::intRadius_

private

Definition at line 62 of file EgammaHcalExtractor.cc.

Referenced by deposit().