HouseholderDecomposition.h

Go to the documentation of this file.

#ifndef HouseholderDecomposition_H
#define HouseholderDecomposition_H
 
#include <vector>
#include <iostream>
 
class HouseholderDecomposition {
public:
  HouseholderDecomposition(int squareMode_ = 5, int mineta_ = 1, int maxeta_ = 85, int minphi_ = 1, int maxphi_ = 20);
 
  ~HouseholderDecomposition();
 
  std::vector<float> runRegional(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 int& regLength = 5);
 
  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);
 
  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>& energyVectorOrig);
 
  std::vector<float> recalibrateEvent(const std::vector<float>& eventSquare,
                                      const int& maxCeta,
                                      const int& maxCphi,
                                      const std::vector<float>& recalibrateVector);
 
  int indexSqr2Reg(const int& sqrIndex, const int& maxCeta, const int& maxCphi);
 
private:
  bool decompose();
 
  void solve(std::vector<float>& y);
 
  std::vector<std::vector<float> > unzipMatrix(const std::vector<std::vector<float> >& eventMatrix,
                                               const std::vector<int>& VmaxCeta,
                                               const std::vector<int>& VmaxCphi);
 
  void makeRegions(const int& regLength);
 
  int squareMode, countEvents;
  int mineta, maxeta, minphi, maxphi, Neta, Nphi;
  int Nchannels, Nxtals, Nevents;
  std::vector<std::vector<float> > eventMatrixOrig;
  std::vector<std::vector<float> > eventMatrixProc;
  std::vector<float> energyVectorProc;
  std::vector<float> alpha;
  std::vector<int> pivot;
 
  std::vector<int> regMinPhi, regMaxPhi, regMinEta, regMaxEta;
  float sigmaReplacement;
};
 
#endif  // HouseholderDecomposition_H