HBHEIsolatedNoiseReflagger.h

Go to the documentation of this file.

#ifndef __HBHE_ISOLATED_NOISE_REFLAGGER_H__
#define __HBHE_ISOLATED_NOISE_REFLAGGER_H__
 
/*
Description: "Reflags" HB/HE hits based on their ECAL, HCAL, and tracking isolation.
 
Original Author: John Paul Chou (Brown University)
                 Thursday, September 2, 2010
*/
#include <memory>
 
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"
#include "RecoLocalCalo/HcalRecAlgos/interface/HBHEIsolatedNoiseAlgos.h"
 
class EcalChannelStatusRcd;
class HcalChannelQuality;
class HcalChannelQualityRcd;
class HcalSeverityLevelComputer;
class HcalSeverityLevelComputerRcd;
class EcalSeverityLevelAlgo;
class EcalSeverityLevelAlgoRcd;
class CaloTowerConstituentsMap;
class CaloGeometryRecord;
class HcalFrontEndMap;
class HcalFrontEndMapRcd;
 
class HBHEIsolatedNoiseReflagger : public edm::stream::EDProducer<> {
public:
  explicit HBHEIsolatedNoiseReflagger(const edm::ParameterSet&);
  ~HBHEIsolatedNoiseReflagger() override;
 
private:
  void produce(edm::Event&, const edm::EventSetup&) override;
 
  void DumpHBHEHitMap(std::vector<HBHEHitMap>& i) const;
 
  // parameters
  edm::EDGetTokenT<HBHERecHitCollection> tok_hbhe_;
  edm::EDGetTokenT<EcalRecHitCollection> tok_EB_;
  edm::EDGetTokenT<EcalRecHitCollection> tok_EE_;
  edm::EDGetTokenT<std::vector<reco::TrackExtrapolation> > tok_trackExt_;
  const HcalFrontEndMap* hfemap;
 
  double LooseHcalIsol_;
  double LooseEcalIsol_;
  double LooseTrackIsol_;
  double TightHcalIsol_;
  double TightEcalIsol_;
  double TightTrackIsol_;
 
  double LooseRBXEne1_, LooseRBXEne2_;
  int LooseRBXHits1_, LooseRBXHits2_;
  double TightRBXEne1_, TightRBXEne2_;
  int TightRBXHits1_, TightRBXHits2_;
  double LooseHPDEne1_, LooseHPDEne2_;
  int LooseHPDHits1_, LooseHPDHits2_;
  double TightHPDEne1_, TightHPDEne2_;
  int TightHPDHits1_, TightHPDHits2_;
  double LooseDiHitEne_;
  double TightDiHitEne_;
  double LooseMonoHitEne_;
  double TightMonoHitEne_;
 
  double RBXEneThreshold_;
 
  bool debug_;
 
  // object validator
  ObjectValidator objvalidator_;
 
  // ES tokens
  edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> ecalChStatusToken_;
  edm::ESGetToken<HcalChannelQuality, HcalChannelQualityRcd> hcalChStatusToken_;
  edm::ESGetToken<HcalSeverityLevelComputer, HcalSeverityLevelComputerRcd> hcalSevToken_;
  edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> ecalSevToken_;
  edm::ESGetToken<CaloTowerConstituentsMap, CaloGeometryRecord> ctcmToken_;
  edm::ESGetToken<HcalFrontEndMap, HcalFrontEndMapRcd> hfemapToken_;
  edm::ESGetToken<CaloGeometry, CaloGeometryRecord> geoToken_;
};
 
#endif