Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-12/MdcHoughFinder/HoughZsFit.h

Go to the documentation of this file.

#ifndef  HoughZsFit_H
#define HoughZsFit_H
#include "MdcHoughFinder/HoughStereo.h"
class HoughZsFit
{
  public:
    HoughZsFit(vector<HoughRecHit>* recHitCol);
    ~HoughZsFit();
//  typedef std::pair<double,double> Point;
    struct Line{
      double _k;
      double _b;
      double _chi;
      //std::vector<HoughRecHit> _recCol;
      std::vector<HoughRecHit> _pointCol;
      std::vector<int> _ambig;
      std::vector<HoughRecHit*> _clusterCol;
    };
    void leastFit(Line &linefit, int N);
    int leastLine(int ,double x[],double y[],double &,double &,double& );
    void  doit();
    void  initPoint();
    double getTanl()  const {return _tanl;}
    double getZ0()    const {return _z0;}
    double getPro()    const {return _pro_correct;}
    void sortHit();
    void print();
    static int m_debug;
    static int m_recMethod;
    static double m_factor;
    static int m_nPoint3D;
    int cgemfit();
    void hough();
    static void setmethod( int m){m_recMethod = m;}
    static vector<double> m_cut_deltaZ;
 
  private:
    vector<HoughRecHit*> _recStereoHit; 
//  Point**  _vecPoint;
    int    _hitSize;
    double _k;
    double _b;
    double _chi;
//  int    _ambig;
    double _z0;
    double _tanl;
    double _pro_correct;
    vector<HoughRecHit*> _cgemCluster[3];
 
};
bool compare_zsfit(const HoughZsFit::Line& a,const HoughZsFit::Line& b);
bool layer_in_track(const HoughRecHit* hita,const HoughRecHit* hitb);
#endif