HLTHcalTowerNoiseCleaner.h

Go to the documentation of this file.

#ifndef HLTHcalTowerNoiseCleaner_h
#define HLTHcalTowerNoiseCleaner_h

#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"   
#include "RecoMET/METAlgorithms/interface/HcalNoiseAlgo.h"

namespace edm {
   class ConfigurationDescriptions;
}

class HLTHcalTowerNoiseCleaner : public edm::EDProducer {
  
 public:
  explicit HLTHcalTowerNoiseCleaner(const edm::ParameterSet&);
  ~HLTHcalTowerNoiseCleaner();
  static void fillDescriptions(edm::ConfigurationDescriptions & descriptions);
  virtual void produce(edm::Event&, const edm::EventSetup&);

 private:
  // parameters
  edm::InputTag HcalNoiseRBXCollectionTag_;
  edm::InputTag TowerCollectionTag_;
  int severity_;
  int maxNumRBXs_;
  int numRBXsToConsider_;
  bool needEMFCoincidence_;
  double minRBXEnergy_;
  double minRatio_;
  double maxRatio_;
  int minHPDHits_;
  int minRBXHits_;
  int minHPDNoOtherHits_;
  int minZeros_;
  double minHighEHitTime_;
  double maxHighEHitTime_;
  double maxRBXEMF_;

  // imported from the RecoMET/METProducers/python/hcalnoiseinfoproducer_cfi
  double minRecHitE_, minLowHitE_, minHighHitE_;

  double TS4TS5EnergyThreshold_;
  std::vector<std::pair<double, double> > TS4TS5UpperCut_;
  std::vector<std::pair<double, double> > TS4TS5LowerCut_;  

  // helper function to compare noise data energies
  struct noisedatacomp {
    inline bool operator() ( const CommonHcalNoiseRBXData& t1, const CommonHcalNoiseRBXData& t2) {
      return t1.energy()>t2.energy();
    }
  };
  typedef std::set<CommonHcalNoiseRBXData, noisedatacomp> noisedataset_t;
};

#endif //HLTHcalTowerNoiseCleaner_h