NoiseRates.h

Go to the documentation of this file.

#ifndef _ANOMALOUSNOISE_ANALYZERS_NOISERATES_H_
#define _ANOMALOUSNOISE_ANALYZERS_NOISERATES_H_

//
// NoiseRates.h
//
//    description: Makes plots to calculate the anomalous noise rates
//
//    author: J.P. Chou, Brown
//
//

// system include files
#include <memory>
#include <string>
#include <vector>

// user include files
#include "DQMServices/Core/interface/DQMEDAnalyzer.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DQMServices/Core/interface/DQMStore.h"

#include "FWCore/ServiceRegistry/interface/Service.h"

#include "DataFormats/METReco/interface/HcalNoiseRBX.h"

// Hcal Hoise Summary
#include "DataFormats/METReco/interface/HcalNoiseSummary.h"

//
// class declaration
//

class NoiseRates : public DQMEDAnalyzer {
public:
  explicit NoiseRates(const edm::ParameterSet &);
  ~NoiseRates() override;

  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;

private:
  void analyze(const edm::Event &, const edm::EventSetup &) override;

  std::string outputFile_;

  // parameters
  edm::InputTag rbxCollName_;  // label for the rbx collection
  edm::EDGetTokenT<reco::HcalNoiseRBXCollection> tok_rbx_;
  double minRBXEnergy_;  // RBX energy threshold
  double minHitEnergy_;  // RecHit energy threshold

  bool useAllHistos_;

  // Hcal Noise Summary Parameters
  edm::EDGetTokenT<HcalNoiseSummary> noisetoken_;

  MonitorElement *hLumiBlockCount_;
  MonitorElement *hRBXEnergy_;
  MonitorElement *hRBXEnergyType1_;
  MonitorElement *hRBXEnergyType2_;
  MonitorElement *hRBXEnergyType3_;
  MonitorElement *hRBXNHits_;

  // Hcal Noise Summary Plots

  MonitorElement *nNNumChannels_;
  MonitorElement *nNSumE_;
  MonitorElement *nNSumEt_;

  MonitorElement *sNNumChannels_;
  MonitorElement *sNSumE_;
  MonitorElement *sNSumEt_;

  MonitorElement *iNNumChannels_;
  MonitorElement *iNSumE_;
  MonitorElement *iNSumEt_;

  MonitorElement *hNoise_maxZeros_;
  MonitorElement *hNoise_maxHPDHits_;
  MonitorElement *hNoise_maxHPDNoOtherHits_;

  // count lumi segments
  std::map<int, int> lumiCountMap_;
};

#endif