RPixDetClusterizer.h

Go to the documentation of this file.

/**********************************************************************
 *
 * Author: F.Ferro - INFN Genova
 * September 2016
 *
 **********************************************************************/
#ifndef RecoCTPPS_PixelLocal_RPixDetClusterizer_H
#define RecoCTPPS_PixelLocal_RPixDetClusterizer_H

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/Common/interface/DetSetVector.h"

#include "DataFormats/CTPPSDigi/interface/CTPPSPixelDigi.h"
#include "DataFormats/CTPPSDigi/interface/CTPPSPixelDigiCollection.h"
#include "DataFormats/CTPPSReco/interface/CTPPSPixelCluster.h"

#include "CondFormats/CTPPSReadoutObjects/interface/CTPPSPixelGainCalibrations.h"
#include "RecoCTPPS/PixelLocal/interface/CTPPSPixelGainCalibrationDBService.h"

#include "CondFormats/CTPPSReadoutObjects/interface/CTPPSPixelAnalysisMask.h"

#include <vector>
#include <set>

class RPixCalibDigi : public CTPPSPixelDigi {

public:

RPixCalibDigi(unsigned char row, unsigned char col, unsigned short adc, unsigned short ele) : CTPPSPixelDigi(row,col,adc){
    electrons_ = ele;
  }

RPixCalibDigi() : CTPPSPixelDigi(){
    electrons_ = 0;
}

  int electrons() const {
    return electrons_;
  }
  void set_electrons(int a) {
    electrons_=a;
  }

private:

  int electrons_;

};


class RPixDetClusterizer{

public:

  RPixDetClusterizer(edm::ParameterSet const& conf);

  void buildClusters(unsigned int detId, const std::vector<CTPPSPixelDigi> &digi, std::vector<CTPPSPixelCluster> &clusters, const CTPPSPixelGainCalibrations * pcalibration, const CTPPSPixelAnalysisMask*  mask);
  int calibrate(unsigned int, int, int, int ,const CTPPSPixelGainCalibrations * pcalibration);

  void make_cluster( RPixCalibDigi const &aSeed,  std::vector<CTPPSPixelCluster> &clusters );
  ~RPixDetClusterizer();


private:

  std::set<CTPPSPixelDigi> rpix_digi_set_;
  std::map<unsigned int,RPixCalibDigi> calib_rpix_digi_map_;
  const edm::ParameterSet &params_;
  int verbosity_;
  unsigned short SeedADCThreshold_;
  unsigned short ADCThreshold_;
  double ElectronADCGain_;
  int VcaltoElectronGain_;
  int VcaltoElectronOffset_;
  bool doSingleCalibration_;
  std::string CalibrationFile_;
  std::vector<RPixCalibDigi> SeedVector_;
};



class RPixTempCluster{

public:

  RPixTempCluster()
  {
    isize=0; 
    curr=0; 
  }
  ~RPixTempCluster(){}

  static constexpr unsigned short MAXSIZE = 256;
  unsigned short adc[MAXSIZE];
  uint8_t row[MAXSIZE];
  uint8_t col[MAXSIZE];
  unsigned short isize;
  unsigned short curr;

  // stack interface (unsafe ok for use below)
  unsigned short top() const { return curr;}
  void pop() { ++curr;}   
  bool empty() { return curr==isize;}

  bool addPixel(unsigned char myrow, unsigned char mycol, unsigned short const iadc) {
    if (isize==MAXSIZE) return false;
    adc[isize]=iadc;
    row[isize]=myrow;
    col[isize++]=mycol;
    return true;
  }
 

};


#endif