SiStripGainCosmicCalculator.h

Go to the documentation of this file.

#ifndef CalibTracker_SiStripChannelGain_SiStripGainCosmicCalculator_h
#define CalibTracker_SiStripChannelGain_SiStripGainCosmicCalculator_h
// -*- C++ -*-
// Package:    SiStripApvGainCalculator
// Class:      SiStripApvGainCalculator
// Original Author:  Dorian Kcira, Pierre Rodeghiero
//         Created:  Mon Nov 20 10:04:31 CET 2006
#include "CommonTools/ConditionDBWriter/interface/ConditionDBWriter.h"
#include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
#include <vector>
#include "TF1.h"
#include "TH1F.h"
#include "TObjArray.h"
#include "TFile.h"
#include "TString.h"
#include <fstream>
#include <sstream>
#include <memory>

class TrackerTopology;
class SiStripDetCabling;
class TrackerGeometry;

class SiStripGainCosmicCalculator : public ConditionDBWriter<SiStripApvGain> {
public:
  explicit SiStripGainCosmicCalculator(const edm::ParameterSet&);
  ~SiStripGainCosmicCalculator() override;

private:
  void algoAnalyze(const edm::Event&, const edm::EventSetup&) override;
  void algoBeginJob(const edm::EventSetup&) override;
  void algoEndJob() override;
  std::unique_ptr<SiStripApvGain> getNewObject() override;

private:
  std::pair<double, double> getPeakOfLandau(TH1F* inputHisto);
  double moduleWidth(const uint32_t detid);
  double moduleThickness(const uint32_t detid);

private:
  std::string TrackProducer;
  std::string TrackLabel;
  //
  TObjArray* HlistAPVPairs;
  TObjArray* HlistOtherHistos;
  uint32_t total_nr_of_events;
  double ExpectedChargeDeposition;
  std::map<uint32_t, double> thickness_map;  // map of detector id to respective thickness
  std::vector<uint32_t> SelectedDetIds;
  std::vector<uint32_t> detModulesToBeExcluded;
  SiStripDetCabling const* siStripDetCabling = nullptr;
  TrackerGeometry const* tkGeom = nullptr;
  unsigned int MinNrEntries;
  double MaxChi2OverNDF;
  bool outputHistogramsInRootFile;
  TString outputFileName;
  bool printdebug_;
  const TrackerTopology* tTopo;
};
#endif