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
// $Id: SiStripGainCosmicCalculator.h,v 1.2 2007/09/13 20:02:52 dkcira Exp $
#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>

class SiStripGainCosmicCalculator : public ConditionDBWriter<SiStripApvGain> {
public:
  explicit SiStripGainCosmicCalculator(const edm::ParameterSet&);
  ~SiStripGainCosmicCalculator();
private:
  void algoAnalyze(const edm::Event &, const edm::EventSetup &);
  void algoBeginJob(const edm::EventSetup&);
  void algoEndJob();
  SiStripApvGain * getNewObject();
private:
  std::pair<double,double> getPeakOfLandau( TH1F * inputHisto );
  double moduleWidth(const uint32_t detid, const edm::EventSetup* iSetup);
  double moduleThickness(const uint32_t detid, const edm::EventSetup* iSetup);
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;
  const edm::EventSetup * eventSetupCopy_;
  unsigned int MinNrEntries;
  double MaxChi2OverNDF;
  bool outputHistogramsInRootFile;
  TString outputFileName ;
  bool printdebug_;
};
#endif