HLTHcalMETNoiseCleaner.h

Go to the documentation of this file.

#ifndef HLTHcalMETNoiseCleaner_h
#define HLTHcalMETNoiseCleaner_h
 
#include "FWCore/Framework/interface/EDFilter.h"
#include "RecoMET/METAlgorithms/interface/HcalNoiseAlgo.h"
#include "DataFormats/METReco/interface/HcalNoiseRBX.h"
#include "DataFormats/METReco/interface/CaloMET.h"
#include "DataFormats/METReco/interface/CaloMETCollection.h"
 
namespace edm {
  class ConfigurationDescriptions;
}
 
class HLTHcalMETNoiseCleaner : public edm::EDFilter {
public:
  explicit HLTHcalMETNoiseCleaner(const edm::ParameterSet&);
  ~HLTHcalMETNoiseCleaner() override;
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
  bool filter(edm::Event&, const edm::EventSetup&) override;
 
private:
  edm::EDGetTokenT<reco::CaloMETCollection> m_theCaloMetToken;
  edm::EDGetTokenT<reco::HcalNoiseRBXCollection> m_theHcalNoiseToken;
  // parameters
  edm::InputTag HcalNoiseRBXCollectionTag_;
  edm::InputTag CaloMetCollectionTag_;
  double CaloMetCut_;
  int severity_;
  int maxNumRBXs_;
  int numRBXsToConsider_;
  bool accept2NoiseRBXEvents_;
  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_, minR45HitE_;
 
  double TS4TS5EnergyThreshold_;
  std::vector<std::pair<double, double> > TS4TS5UpperCut_;
  std::vector<std::pair<double, double> > TS4TS5LowerCut_;
 
  reco::CaloMET BuildCaloMet(float sumet, float pt, float phi);
 
  // helper function to compare noise data energies
  struct noisedatacomp {
    inline bool operator()(const CommonHcalNoiseRBXData& t1, const CommonHcalNoiseRBXData& t2) const {
      return t1.energy() > t2.energy();
    }
  };
  typedef std::set<CommonHcalNoiseRBXData, noisedatacomp> noisedataset_t;
};
 
#endif  //HLTHcalMETNoiseCleaner_h