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 "FWCore/Framework/interface/ESHandle.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&,  edm::ESHandle<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