SiStripCorrelateNoise.h

Go to the documentation of this file.

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

// system include files
#include <memory>

// user include files
#include "CondFormats/SiStripObjects/interface/SiStripApvGain.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 SiStripCorrelateNoise : public edm::EDAnalyzer {
public:
  explicit SiStripCorrelateNoise(const edm::ParameterSet &);
  ~SiStripCorrelateNoise() 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 DoPlots();
  void DoAnalysis(const edm::EventSetup &, const SiStripNoises &, SiStripNoises &);
  void getHistos(const uint32_t &detid, const TrackerTopology *tTopo, std::vector<TH1F *> &histos);
  TH1F *getHisto(const long unsigned int &index);

  unsigned long long getNoiseCache(const edm::EventSetup &eSetup) {
    return eSetup.get<SiStripNoisesRcd>().cacheIdentifier();
  }
  unsigned long long getGainCache(const edm::EventSetup &eSetup) {
    return eSetup.get<SiStripApvGainRcd>().cacheIdentifier();
  }
  void checkGainCache(const edm::EventSetup &es);
  float getMeanNoise(const SiStripNoises::Range &noiseRange, const uint32_t &first, const uint32_t &range);

  float getGainRatio(const uint32_t &detid, const uint16_t &apv);

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

  struct Data {
    uint32_t detid;
    std::vector<float> values;
  };

  SiStripDetInfoFileReader *fr;
  edm::ESHandle<SiStripApvGain> gainHandle_;
  edm::ESHandle<SiStripNoises> noiseHandle_;

  uint32_t theRun;
  SiStripNoises *refNoise;

  SiStripApvGain *oldGain, *newGain;
  bool equalGain;

  TFile *file;
  std::vector<TH1F *> vTH1;

  TrackerMap *tkmap;

  unsigned long long cacheID_noise;
  unsigned long long cacheID_gain;
};