HiEgammaAlgos/interface/HcalRechitIsoCalculator.h

0001 #ifndef HiEgammaAlgos_HcalRechitIsoCalculator_h
0002 #define HiEgammaAlgos_HcalRechitIsoCalculator_h
0003
0004 #include "FWCore/Framework/interface/Event.h"
0005 #include "FWCore/Framework/interface/EventSetup.h"
0006 #include "FWCore/Framework/interface/ESHandle.h"
0007
0008 #include "FWCore/Utilities/interface/InputTag.h"
0009
0010 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
0011 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
0012
0013 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
0014 #include "Geometry/Records/interface/IdealGeometryRecord.h"
0015
0016 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
0017 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
0018
0019 class HcalRechitIsoCalculator {
0020 public:
0021   HcalRechitIsoCalculator(const CaloGeometry *geometry,
0022                           const edm::Handle<HBHERecHitCollection> hbhe,
0023                           const edm::Handle<HFRecHitCollection> hfLabel,
0024                           const edm::Handle<HORecHitCollection> hoLabel);
0025
0026   /// Return the hcal rechit energy in a cone around the SC
0027   double getHcalRechitIso(const reco::SuperClusterRef clus,
0028                           const double i,
0029                           const double threshold,
0030                           const double innerR = 0.0);
0031   /// Return the background-subtracted hcal rechit energy in a cone around the SC
0032   double getBkgSubHcalRechitIso(const reco::SuperClusterRef clus,
0033                                 const double i,
0034                                 const double threshold,
0035                                 const double innerR = 0.0);
0036
0037 private:
0038   const HBHERecHitCollection *fHBHERecHits_;
0039   const HORecHitCollection *fHORecHits_;
0040   const HFRecHitCollection *fHFRecHits_;
0041   const CaloGeometry *geometry_;
0042 };
0043
0044 #endif