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