SiStripHotStripAlgorithmFromClusterOccupancy.h

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    SiStripQuality
// Class:      SiStripHotStripAlgorithmFromClusterOccupancy
//
//
// Original Author:  Domenico GIORDANO
//         Created:  Wed Oct  3 12:11:10 CEST 2007
//
//
 
#ifndef CalibTracker_SiStripQuality_SiStripHotStripAlgorithmFromClusterOccupancy_H
#define CalibTracker_SiStripQuality_SiStripHotStripAlgorithmFromClusterOccupancy_H
 
// system include files
#include <memory>
#include <vector>
#include <sstream>
#include <iostream>
 
#include "TMath.h"
#include "TTree.h"
#include "TFile.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
#include "Geometry/CommonTopologies/interface/StripTopology.h"
#include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
 
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "CalibTracker/SiStripQuality/interface/SiStripQualityHistos.h"
#include "CalibTracker/SiStripCommon/interface/SiStripDetInfoFileReader.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "CalibTracker/Records/interface/SiStripQualityRcd.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "CommonTools/TrackerMap/interface/TrackerMap.h"
 
class SiStripQuality;
class TrackerTopology;
 
class SiStripHotStripAlgorithmFromClusterOccupancy {
public:
  typedef SiStrip::QualityHistosMap HistoMap;
 
  SiStripHotStripAlgorithmFromClusterOccupancy(const edm::ParameterSet&, const TrackerTopology*);
 
  virtual ~SiStripHotStripAlgorithmFromClusterOccupancy();
 
  void setProbabilityThreshold(long double prob) { prob_ = prob; }
  void setMinNumEntries(unsigned short m) { MinNumEntries_ = m; }
  void setMinNumEntriesPerStrip(unsigned short m) { MinNumEntriesPerStrip_ = m; }
  void setOccupancyThreshold(long double occupancy) {
    occupancy_ = occupancy;
    minNevents_ = occupancy_ * Nevents_;
  }
  void setNumberOfEvents(double Nevents);
  void setOutputFileName(std::string OutputFileName, bool WriteOutputFile) {
    OutFileName_ = OutputFileName;
    WriteOutputFile_ = WriteOutputFile;
  }
  void setTrackerGeometry(const TrackerGeometry* tkgeom) { TkGeom = tkgeom; }
  void extractBadStrips(SiStripQuality*, HistoMap&, const SiStripQuality*);
 
private:
  // unsigned long long m_cacheID_;
  //std::string dataLabel_;
  //edm::FileInPath fp_;
  //SiStripDetInfoFileReader* reader;
 
  struct pHisto {
    pHisto() : _NEntries(0), _NEmptyBins(0), _SubdetId(0){};
    TH1F* _th1f;
    int _NEntries;
    int _NEmptyBins;
    int _SubdetId;
  };
 
  void iterativeSearch(pHisto&, std::vector<unsigned int>&, int);
  void evaluatePoissonian(std::vector<long double>&, long double& meanVal);
 
  long double prob_;
  long double ratio_;
  unsigned short MinNumEntries_;
  unsigned short MinNumEntriesPerStrip_;
  double Nevents_;
  double minNevents_;
  long double occupancy_;
  std::string OutFileName_;
  bool WriteOutputFile_;
  const TrackerGeometry* TkGeom;
  const TrackerTopology* tTopo;
 
  SiStripQuality* pQuality;
 
  TFile* f;
  TTree* striptree;
  bool UseInputDB_;
  int detrawid;
  int subdetid;
  int layer_ring;
  int disc;
  int isback;
  int isexternalstring;
  int iszminusside;
  int rodstringpetal;
  int isstereo;
  int module_position;
  int number_strips;
  int strip_number;
  int apv_channel;
 
  float global_position_x;
  float global_position_y;
  float global_position_z;
 
  int isHot;
  int hotStripsPerAPV;
  int hotStripsPerModule;
  double stripOccupancy;
  int stripHits;
  double medianAPVHits;
  double avgAPVHits;
  double poissonProb;
 
  int ishot[768];
  int hotstripsperapv[6];
  int hotstripspermodule;
  double stripoccupancy[768];
  int striphits[768];
  double poissonprob[768];
  double medianapvhits[6];
  double avgapvhits[6];
 
  std::stringstream ss;
};
#endif