/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/SimCalorimetry/EcalTrigPrimAlgos/src/EcalFenixPeakFinder.cc

Go to the documentation of this file.

#include <SimCalorimetry/EcalTrigPrimAlgos/interface/EcalFenixPeakFinder.h>
  
EcalFenixPeakFinder::EcalFenixPeakFinder()
  :inputsAlreadyIn_(0){ 
}

  
EcalFenixPeakFinder::~EcalFenixPeakFinder(){
}

int EcalFenixPeakFinder::setInput(int input)
{
  if(inputsAlreadyIn_<3)
    {
      buffer_[inputsAlreadyIn_]=input;
      inputsAlreadyIn_++;
    }
  else
    {
      for(int i=0; i<2; i++) buffer_[i]=buffer_[i+1];
      buffer_[2]=input;
    }
  return 1;
}
  
int EcalFenixPeakFinder::process()
{

  if(inputsAlreadyIn_<3) return 0;
  if(buffer_[1]>buffer_[0] && buffer_[1]>buffer_[2]) return 1;
  else return 0;
    
}

std::vector<int> EcalFenixPeakFinder::process(std::vector<int> &filtout, std::vector<int> & output)
{
  
  // FIXME: 3
  inputsAlreadyIn_=0;
  for (unsigned int i =0;i<3;i++) buffer_[i]=0;

  //  std::vector<int> output;
  
  // attention, we have to shift by one, because the peak is found one too late
  for (unsigned int i =0;i<filtout.size();i++){
    
    setInput(filtout[i]);
    if (i>0) {
      //      int outone = process();
      //      output.push_back(outone);
      output[i-1]= process();
    }
  }
  //  output.resize(filtout.size());

  return output;
}