RxCalculator Class Reference

#include <RxCalculator.h>

List of all members.

Public Member Functions
double	getCRx (const reco::SuperClusterRef clus, double i, double threshold)
double	getRFx (const reco::SuperClusterRef clus, double i, double threshold)
double	getROx (const reco::SuperClusterRef clus, double i, double threshold)
double	getRx (const reco::SuperClusterRef clus, double i, double threshold)
	RxCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, edm::InputTag hbheLabel, edm::InputTag hfLabel, edm::InputTag hoLabel)
Private Attributes
const HBHERecHitCollection *	fHBHERecHits_
const HFRecHitCollection *	fHFRecHits_
const HORecHitCollection *	fHORecHits_
const CaloGeometry *	geometry_

---

Detailed Description

Definition at line 19 of file RxCalculator.h.

---

Constructor & Destructor Documentation

RxCalculator::RxCalculator	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		edm::InputTag	hbheLabel,
		edm::InputTag	hfLabel,
		edm::InputTag	hoLabel
	)

Definition at line 20 of file RxCalculator.cc.

References edm::EventSetup::get(), edm::Event::getByLabel(), edm::ESHandle< T >::product(), and edm::Handle< T >::product().

{
   Handle<HFRecHitCollection> hfhandle;
   iEvent.getByLabel(hfLabel, hfhandle);
   fHFRecHits_ = hfhandle.product();

   Handle<HORecHitCollection> hohandle;
   iEvent.getByLabel(hoLabel, hohandle);
   fHORecHits_ = hohandle.product();

   Handle<HBHERecHitCollection> hehbhandle;
   iEvent.getByLabel(hbheLabel, hehbhandle);
   fHBHERecHits_ = hehbhandle.product();

   ESHandle<CaloGeometry> geometryHandle;
   iSetup.get<CaloGeometryRecord>().get(geometryHandle);
   geometry_ = geometryHandle.product();

} 

---

Member Function Documentation

double RxCalculator::getCRx	(	const reco::SuperClusterRef	clus,
		double	i,
		double	threshold
	)

Definition at line 156 of file RxCalculator.cc.

References cond::rpcobgas::detid, dPhi(), CaloRecHit::energy(), ExpressReco_HICollisions_FallBack::et, eta(), PV3DBase< T, PVType, FrameType >::eta(), HBHERecHit::id(), getHLTprescales::index, PV3DBase< T, PVType, FrameType >::phi(), phi, PI, dt_offlineAnalysis_common_cff::reco, and ExpressReco_HICollisions_FallBack::x.

Referenced by HiEgammaIsolationProducer::produce().

{
   using namespace edm;
   using namespace reco;


   if(!fHBHERecHits_) {       
      LogError("RxCalculator") << "Error! Can't get HBHERecHits for event.";
      return -100;
   }

   double SClusterEta = cluster->eta();
   double SClusterPhi = cluster->phi();
   double TotalEt = 0;

   for(size_t index = 0; index < fHBHERecHits_->size(); index++) {
      const HBHERecHit &rechit = (*fHBHERecHits_)[index];
      const DetId &detid = rechit.id();
      const GlobalPoint& hitpoint = geometry_->getPosition(detid);
      double eta = hitpoint.eta();
      double phi = hitpoint.phi();
      double dEta = fabs(eta-SClusterEta);
      double dPhi = fabs(phi-SClusterPhi);
      while (dPhi>2*PI) dPhi-=2*PI;
      if (dPhi>PI) dPhi=2*PI-dPhi;

      if (dEta<x*0.1) {
         double et = rechit.energy()/cosh(eta);
         if (et<threshold) et=0;
         TotalEt += et;
      }
   }

   double Rx = getRx(cluster,x,threshold);
   double CRx = Rx - TotalEt / 40.0 * x;

   return CRx;
}

double RxCalculator::getRFx	(	const reco::SuperClusterRef	clus,
		double	i,
		double	threshold
	)

Definition at line 116 of file RxCalculator.cc.

References cond::rpcobgas::detid, dPhi(), CaloRecHit::energy(), ExpressReco_HICollisions_FallBack::et, eta(), PV3DBase< T, PVType, FrameType >::eta(), HFRecHit::id(), getHLTprescales::index, PV3DBase< T, PVType, FrameType >::phi(), phi, PI, dt_offlineAnalysis_common_cff::reco, and mathSSE::sqrt().

{
   using namespace edm;
   using namespace reco;

   if(!fHFRecHits_) {       
      LogError("RxCalculator") << "Error! Can't get HFRecHits for event.";
      return -100;
   }

   double SClusterEta = cluster->eta();
   double SClusterPhi = cluster->phi();
   double TotalEt = 0;

   for(size_t index = 0; index < fHFRecHits_->size(); index++) {
      const HFRecHit &rechit = (*fHFRecHits_)[index];
      const DetId &detid = rechit.id();
      const GlobalPoint& hitpoint = geometry_->getPosition(detid);
      double eta = hitpoint.eta();
      double phi = hitpoint.phi();
      double dEta = fabs(eta-SClusterEta);
      double dPhi = fabs(phi-SClusterPhi);
      while (dPhi>2*PI) dPhi-=2*PI;
      if (dPhi>PI) dPhi=2*PI-dPhi;


      double dR = sqrt(dEta * dEta + dPhi * dPhi);
      if (dR<x*0.1) {
         double et = rechit.energy()/cosh(eta);
         if (et<threshold) et=0;
         TotalEt += et;
      }
   }



   return TotalEt;
}

double RxCalculator::getROx	(	const reco::SuperClusterRef	clus,
		double	i,
		double	threshold
	)

Definition at line 81 of file RxCalculator.cc.

References cond::rpcobgas::detid, dPhi(), CaloRecHit::energy(), ExpressReco_HICollisions_FallBack::et, eta(), PV3DBase< T, PVType, FrameType >::eta(), HORecHit::id(), getHLTprescales::index, PV3DBase< T, PVType, FrameType >::phi(), phi, PI, dt_offlineAnalysis_common_cff::reco, and mathSSE::sqrt().

{
   using namespace edm;
   using namespace reco;

   if(!fHORecHits_) {       
      LogError("RxCalculator") << "Error! Can't get HORecHits for event.";
      return -100;
   }

   double SClusterEta = cluster->eta();
   double SClusterPhi = cluster->phi();
   double TotalEt = 0;

   for(size_t index = 0; index < fHORecHits_->size(); index++) {
      const HORecHit &rechit = (*fHORecHits_)[index];
      const DetId &detid = rechit.id();
      const GlobalPoint& hitpoint = geometry_->getPosition(detid);
      double eta = hitpoint.eta();
      double phi = hitpoint.phi();
      double dEta = fabs(eta-SClusterEta);
      double dPhi = fabs(phi-SClusterPhi);
      while (dPhi>2*PI) dPhi-=2*PI;
      if (dPhi>PI) dPhi=2*PI-dPhi;

      double dR = sqrt(dEta * dEta + dPhi * dPhi);
      if (dR<x*0.1) {
         double et = rechit.energy()/cosh(eta);
         if (et<threshold) et=0;
         TotalEt += et;
      }
   }
   return TotalEt;
}

double RxCalculator::getRx	(	const reco::SuperClusterRef	clus,
		double	i,
		double	threshold
	)

Definition at line 41 of file RxCalculator.cc.

References cond::rpcobgas::detid, dPhi(), CaloRecHit::energy(), ExpressReco_HICollisions_FallBack::et, eta(), PV3DBase< T, PVType, FrameType >::eta(), HBHERecHit::id(), getHLTprescales::index, PV3DBase< T, PVType, FrameType >::phi(), phi, PI, dt_offlineAnalysis_common_cff::reco, and mathSSE::sqrt().

Referenced by HiEgammaIsolationProducer::produce().

{
   using namespace edm;
   using namespace reco;


   if(!fHBHERecHits_) {       
      LogError("RxCalculator") << "Error! Can't get HBHERecHits for event.";
      return -100;
   }

   double SClusterEta = cluster->eta();
   double SClusterPhi = cluster->phi();
   double TotalEt = 0;

   for(size_t index = 0; index < fHBHERecHits_->size(); index++) {
      const HBHERecHit &rechit = (*fHBHERecHits_)[index];
      const DetId &detid = rechit.id();
      const GlobalPoint& hitpoint = geometry_->getPosition(detid);
      double eta = hitpoint.eta();
      double phi = hitpoint.phi();
      double dEta = fabs(eta-SClusterEta);
      double dPhi = fabs(phi-SClusterPhi);
      while (dPhi>2*PI) dPhi-=2*PI;
      if (dPhi>PI) dPhi=2*PI-dPhi;

      if (dPhi>PI) dPhi=2*PI-dPhi;

      double dR = sqrt(dEta * dEta + dPhi * dPhi);

      if (dR<x*0.1) {
         double et = rechit.energy()/cosh(eta);
         if (et<threshold) et=0;
         TotalEt += et;
      }
   }

   return TotalEt;
}

---

Member Data Documentation

const HBHERecHitCollection* RxCalculator::fHBHERecHits_ [private]

Definition at line 32 of file RxCalculator.h.

const HFRecHitCollection* RxCalculator::fHFRecHits_ [private]

Definition at line 34 of file RxCalculator.h.

const HORecHitCollection* RxCalculator::fHORecHits_ [private]

Definition at line 33 of file RxCalculator.h.

const CaloGeometry* RxCalculator::geometry_ [private]

Definition at line 35 of file RxCalculator.h.