MinL3Algorithm.h

Go to the documentation of this file.

#ifndef MinL3Algorithm_H
#define MinL3Algorithm_H

#include <vector>
#include <iostream>


class MinL3Algorithm
{
public:
  MinL3Algorithm(float kweight_ = 0., int squareMode_ = 5, int mineta_ = 1, int maxeta_ = 85, int minphi_ = 1, int maxphi_ = 20);

  ~MinL3Algorithm();

  std::vector<float> iterate(const std::vector<std::vector<float> >& eventMatrix, const std::vector<int>& VmaxCeta, const std::vector<int>& VmaxCphi, const std::vector<float>& energyVector, const int& nIter, const bool& normalizeFlag = false);


  void addEvent(const std::vector<float>& eventSquare, const int& maxCeta, const int& maxCphi, const float& energy);

  std::vector<float> recalibrateEvent(const std::vector<float>& eventSquare, const int& maxCeta, const int& maxCphi, const std::vector<float>& recalibrateVector); 

  std::vector<float> getSolution(bool resetsolution=true);

  void resetSolution(); 

  int indexSqr2Reg(const int& sqrIndex, const int& maxCeta, const int& maxCphi);


private:

  float kweight;
  int squareMode;
  int mineta, maxeta, minphi, maxphi;
  int countEvents;
  int Nchannels, Nxtals;
  std::vector<float> wsum;
  std::vector<float> Ewsum;

};

#endif // MinL3Algorithm_H