SiStripCorrelateBadStripAndNoise.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripCorrelateBadStripAndNoise
// Class:      SiStripCorrelateBadStripAndNoise
//
//
// Original Author:  Domenico GIORDANO
//         Created:  Mon Aug 10 10:42:04 CEST 2009
//
//

// system include files
#include <memory>

// user include files
#include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

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

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

#include "CalibTracker/Records/interface/SiStripDependentRecords.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
#include "CommonTools/TrackerMap/interface/TrackerMap.h"

#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"

//
// class decleration
//
class TrackerTopology;
class SiStripCorrelateBadStripAndNoise : public edm::EDAnalyzer {
public:
  explicit SiStripCorrelateBadStripAndNoise(const edm::ParameterSet &);
  ~SiStripCorrelateBadStripAndNoise() override;

private:
  void beginRun(const edm::Run &run, const edm::EventSetup &es) override;
  void analyze(const edm::Event &, const edm::EventSetup &) override{};
  void endJob() override;

  void DoAnalysis(const edm::EventSetup &);
  void getHistos(const uint32_t &detid, const TrackerTopology *tTopo, std::vector<TH2F *> &histos);
  TH2F *getHisto(const long unsigned int &index);

  unsigned long long getNoiseCache(const edm::EventSetup &eSetup) {
    return eSetup.get<SiStripNoisesRcd>().cacheIdentifier();
  }
  unsigned long long getQualityCache(const edm::EventSetup &eSetup) {
    return eSetup.get<SiStripQualityRcd>().cacheIdentifier();
  }

  void iterateOnDets(const TrackerTopology *tTopo);
  void iterateOnBadStrips(const uint32_t &detid, const TrackerTopology *tTopo, SiStripQuality::Range &sqrange);
  void correlateWithNoise(const uint32_t &detid,
                          const TrackerTopology *tTopo,
                          const uint32_t &firstStrip,
                          const uint32_t &range);
  float getMeanNoise(const SiStripNoises::Range &noiseRange, const uint32_t &first, const uint32_t &range);

  // ----------member data ---------------------------

  SiStripDetInfoFileReader *fr;
  edm::ESHandle<SiStripQuality> qualityHandle_;
  edm::ESHandle<SiStripNoises> noiseHandle_;

  TFile *file;
  std::vector<TH2F *> vTH2;

  TrackerMap *tkmap;

  unsigned long long cacheID_quality;
  unsigned long long cacheID_noise;
};