HcalRechitIsoCalculator.h

Go to the documentation of this file.

#ifndef HiEgammaAlgos_HcalRechitIsoCalculator_h
#define HiEgammaAlgos_HcalRechitIsoCalculator_h

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "FWCore/Utilities/interface/InputTag.h"

#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
#include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"

#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"

#include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

class HcalRechitIsoCalculator {
public:
  HcalRechitIsoCalculator(const CaloGeometry *geometry,
                          const edm::Handle<HBHERecHitCollection> hbhe,
                          const edm::Handle<HFRecHitCollection> hfLabel,
                          const edm::Handle<HORecHitCollection> hoLabel);

  double getHcalRechitIso(const reco::SuperClusterRef clus,
                          const double i,
                          const double threshold,
                          const double innerR = 0.0);
  double getBkgSubHcalRechitIso(const reco::SuperClusterRef clus,
                                const double i,
                                const double threshold,
                                const double innerR = 0.0);

private:
  const HBHERecHitCollection *fHBHERecHits_;
  const HORecHitCollection *fHORecHits_;
  const HFRecHitCollection *fHFRecHits_;
  const CaloGeometry *geometry_;
};

#endif