WeightedPoints3D.h

#ifndef WeightedPoints3D_h
#define WeightedPoints3D_h


#include <stdio.h>
#include <iostream>
#include <iomanip>
#include <string>
#include <sstream>
#include <cstdlib>
#include <math.h>
#include<vector>


class WeightedPoints3D {

public:

  WeightedPoints3D(int nhits, double* a, double* x, double* y, double* z);
  WeightedPoints3D(const std::vector<double> &cog,
                   const std::vector<double> &cov, const std::vector<double> &mayor_axis_error =  std::vector<double>() , 
                   int npnt = 0 ,double wgtsum =0.0 , double wgt2sum=0.0 , double wgt4sum=0.0 );

  ~WeightedPoints3D();

  void setErrors(double *ex, double* ey, double *ez);



  int getNumberOfPoints();

  double getTotalWeight();

  double getTotalSquaredWeight();

  double getTotalQuarticWeight();

  double* getCentreOfGravity();
 
  double* getCentreOfGravityErrors();

  inline double* getCenterOfGravity() { return getCentreOfGravity() ; }
  inline double* getCenterOfGravityErrors() { return getCentreOfGravityErrors() ; }
  
  double* getEigenVal();

  double* getEigenValErrors();

  double* getEigenVecCartesian();

  double* getEigenVecPolar();
 
  double* getEigenVecCartesianErrors();
  double* getEigenVecPolarErrors();

  double getWidth();

 
  inline double radius() { return _radius; }


private:

  int _nPoints;

  double* _Wgt;
  double* _x;
  double* _y;
  double* _z;
  double* _ex;
  double* _ey;
  double* _ez;
  double* _xl;
  double* _xt;
  double* _t;
  double* _s;

  int   _ifNotPointsGiven;
  int   _ifNotCOG;
  int   _ifNotCOGErrors;
  double _sumWgt;
  double _wgt_sqr_sum ;
  double _wgt_4_sum ;
  double _radius;
  double _xgr;
  double _ygr;
  double _zgr;
  double _COG[3];
  double _COGCov[3][3];
  double _COGCovCov[6][6];


  int   _ifNotWidth;
  double _Width;


  int   _ifNotEigenSolved;
  int   _ifNotEigenToPolarDone ;
  int   _ifNotCovErrors;
  int   _ifNotVecErrorsPolarPropagated;
  int   _ifNotVecErrorsCartesianPropagated;
  int   _ifNotValErrorsPropagated;
  int   _ifNot_dVec_dCov;
  int   _ifNot_fourthmom;
  int   _last_evec_to_error_propagate;
  double _EigenVal[3];
  double _EigenValCov[3][3];
  double _EigenVec[3][3];
  double _EigenVecAngle[2][3];
  double _EigenVecCov[9][9];
  double _theta_phi_cov[2][2][3];
  double _xyz_cov[3][3][3];
  double _var_lam[3];
  double _nfact1;
  double _nfact2;
  double _nfact3;
  double _nfact4;
  double _fourth_mom[3][3][3][3];
  double dang[6][2][3]; 



  void  findCOG();
  void  findCOGErrors();
  double findDistance(int i);
  void  solveEigenValEq();
  void  findCovErrors();
  void  findWidth();
  void  propagateVecErrorsPolar();
  void  propagateVecErrorsCartesian();
  void  propagateValErrors();
  void  find_dVec_dCov();

};


#endif