PixelThresholdClusterizer.h

Go to the documentation of this file.

#ifndef RecoLocalTracker_SiPixelClusterizer_PixelThresholdClusterizer_H
#define RecoLocalTracker_SiPixelClusterizer_PixelThresholdClusterizer_H

//-----------------------------------------------------------------------
//-----------------------------------------------------------------------

// Base class, defines SiPixelDigi and SiPixelCluster.  The latter includes
// Pixel, PixelPos and Shift as inner classes.
//
#include "DataFormats/Common/interface/DetSetVector.h"
#include "PixelClusterizerBase.h"

// The private pixel buffer
#include "SiPixelArrayBuffer.h"

// Parameter Set:
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include <vector>


class dso_hidden PixelThresholdClusterizer final : public PixelClusterizerBase {
 public:

  PixelThresholdClusterizer(edm::ParameterSet const& conf);
  ~PixelThresholdClusterizer();

  // Full I/O in DetSet
  void clusterizeDetUnit( const edm::DetSet<PixelDigi> & input, 
                  const PixelGeomDetUnit * pixDet,
                  const std::vector<short>& badChannels,
                  edmNew::DetSetVector<SiPixelCluster>::FastFiller& output) { clusterizeDetUnitT(input, pixDet, badChannels, output); }
  void clusterizeDetUnit( const edmNew::DetSet<SiPixelCluster> & input,
                          const PixelGeomDetUnit * pixDet,
                          const std::vector<short>& badChannels,
                          edmNew::DetSetVector<SiPixelCluster>::FastFiller& output) { clusterizeDetUnitT(input, pixDet, badChannels, output); }

 private:

  template<typename T>
  void clusterizeDetUnitT( const T & input,
                           const PixelGeomDetUnit * pixDet,
                           const std::vector<short>& badChannels,
                           edmNew::DetSetVector<SiPixelCluster>::FastFiller& output);

  SiPixelArrayBuffer               theBuffer;         // internal nrow * ncol matrix
  bool                             bufferAlreadySet;  // status of the buffer array
  std::vector<SiPixelCluster::PixelPos>  theSeeds;          // cached seed pixels
  std::vector<SiPixelCluster>            theClusters;       // resulting clusters  
  
  float thePixelThresholdInNoiseUnits;    // Pixel threshold in units of noise
  float theSeedThresholdInNoiseUnits;     // Pixel cluster seed in units of noise
  float theClusterThresholdInNoiseUnits;  // Cluster threshold in units of noise

  const int   thePixelThreshold;  // Pixel threshold in electrons
  const int   theSeedThreshold;   // Seed threshold in electrons
  const float theClusterThreshold;  // Cluster threshold in electrons
  const int   theConversionFactor;  // adc to electron conversion factor
  const int   theOffset;            // adc to electron conversion offset

  const int   theStackADC_;          // The maximum ADC count for the stack layers
  const int   theFirstStack_;        // The index of the first stack layer
  const double theElectronPerADCGain_;  //  ADC to electrons conversion

  int  theNumOfRows;
  int  theNumOfCols;
  uint32_t detid_;
  bool dead_flag;
  const bool doMissCalibrate; // Use calibration or not
  const bool doSplitClusters;
  bool setup(const PixelGeomDetUnit * pixDet);
  void copy_to_buffer( DigiIterator begin, DigiIterator end );   
  void copy_to_buffer( ClusterIterator begin, ClusterIterator end );
  void clear_buffer( DigiIterator begin, DigiIterator end );
  void clear_buffer( ClusterIterator begin, ClusterIterator end );
  SiPixelCluster make_cluster( const SiPixelCluster::PixelPos& pix, edmNew::DetSetVector<SiPixelCluster>::FastFiller& output);
  // Calibrate the ADC charge to electrons 
  int calibrate(int adc, int col, int row);

};

#endif