QVector.hh

Go to the documentation of this file.

 
#ifndef __QVECTOR_HH_
#define __QVECTOR_HH_
 
#include <gsl/gsl_vector.h>
#include <iostream>
#include <string>
#include "QObject.hh"
#include <vector>
 
class TGraph;
Q_BEGIN_NAMESPACE
 
//FIXME CMT: Mult and Div definition has changed
 
class QMatrix;
class QVectorI;
 
class QVector : public QObject {
public:
  QVector();
  explicit QVector(const UInt_t size);
  QVector(const QVector& orig);
  virtual ~QVector();
  inline UInt_t Size() const {return fSize;}
  void Initialize(const double val = 0);
  double& operator[](const UInt_t i);
  double& operator()(const UInt_t i){return operator[](i);}
  const double &operator[] (const UInt_t i) const;
 
  const double &operator() (const UInt_t i) const { return operator[](i);}
  const QVector& operator = (const QVector& orig);
  const QVector& operator = (const QVectorI& orig);
  
  const QVector& operator = (const std::vector<double>& orig);
  QVector operator-() const;
  //was:                const QVector& Mult(const QVector& other);
  QVector Mult(const QVector& other) const;
  //was:                const QVector& Div(const QVector& other);
  QVector Div(const QVector& other) const;
  const QVector& operator*=(const double t);
  const QVector & operator/=(const double v);
  const QVector & operator+=(const double v);
  const QVector & operator-=(const double v);
 
 
  const QVector & operator+=(const QVector&v);
  const QVector & operator-=(const QVector&v);
  const QVector & operator*=(const QVector &v);
  const QVector & operator/=(const QVector &v);
 
  double operator*(const QVector& other) const;
  QVector operator+(const QVector& other) const;
  QVector operator-(const QVector& other) const;
  QVector operator+(double t) const;
  QVector operator-(double t) const;
  QVector operator*(double t) const;
  QVector operator/(double t) const;
  QMatrix T() const;
 
  QMatrix Reshape(UInt_t n,UInt_t m);
  const QMatrix Reshape(UInt_t n,UInt_t m) const;
  
  const QVector& Rebin(const int rebin,bool sum=false);
 
  const QVector& Shift(const int nstep);
  const QVector& Reverse();
 
  const QVector& ShiftReal(const double fstep);
  const QVector& Differentiate();
  
  QVector Derivative2P(const double delta = 1.) const;
  QVector Derivative2PF(const double delta = 1.) const;
  
  QVector Derivative3P(const double delta = 1.) const;
  
  QVector Derivative3PB(const double delta = 1.) const;
  QVector Derivative3PF(const double delta = 1.) const;
  
  QVector Derivative5P(const double delta = 1.) const;
  
  QVector Derivative5PB(const double delta = 1.) const;
  
  QVector DerivativeEulero(const double delta = 1.) const;
  QVector Abs() const;
 
  const QVector& Sort(const bool asc=true);
  const QVector& Integrate(const double constTerm = 0);
  void Resize(const UInt_t newsize);
  void Append(const double val);
  double Norma() const;
  double Sum(UInt_t nelem, UInt_t first = 0) const; 
  double Sum() const { return fSize*GetMean(); }
  
  double GetMean(const UInt_t nelem,const UInt_t first = 0) const {
    return Sum(nelem,first)/nelem;
  }
  double GetMean() const { 
    return GetMean(fSize);
  }
  double GetRMS(const UInt_t nelem,const UInt_t first = 0) const; 
  double GetRMS() const { 
    return GetRMS(fSize);
  }
  double GetMedian() const;
  double GetMedianAbsoluteDeviation() const;
  double GetMax() const 
  { SetMathVector(); return gsl_vector_max(fMathVec);  }
  double GetMin() const 
  { SetMathVector(); return gsl_vector_min(fMathVec);  }
  int GetMaxIndex(const UInt_t nelem,const UInt_t first) const;
  int GetMinIndex(const UInt_t nelem,const UInt_t first) const;
  UInt_t GetMaxIndex() const 
  {SetMathVector(); return gsl_vector_max_index(fMathVec);}
  UInt_t GetMinIndex() const {SetMathVector(); return gsl_vector_min_index(fMathVec);}
  double* GetArray() const;
  std::vector<double> GetArrayVector()const{return std::vector<double>(fData,fData + Size());}
  const double* GetConstArray() const { return fData;}
  void SetArray(const double *orig,const UInt_t size);
  virtual void Clear() {Resize(0);}
  std::string sqlString() const;
  
  virtual void Dump(std::ostream& o) const;
  
  void Draw(Option_t* option=""); // *MENU*
  
  TGraph* GetGraph(double samplingFrequency = 1., double scale = 1.) const;
 
  QVector GetSubVector(UInt_t size,UInt_t start=0,UInt_t stride=1);
 
  const QVector GetSubVector(UInt_t size,UInt_t start=0,UInt_t stride=1) const;
 
protected: 
  QVector(const QVector& orig,const bool dataOwner);
  
  QVector(double* vec,const UInt_t size,const UInt_t stride);
 
  const bool fDataOwner; 
 
  UInt_t fSize; 
  
  UInt_t fAllocSize; 
 
  UInt_t fStride;
  
  Double_t *fData; //[fSize]
 
  mutable gsl_vector* fMathVec; 
  
private:
  void SetMathVector() const;
  static double * ArrayAlloc(const UInt_t size);
  static void ArrayCopy(const double* orig,const UInt_t origStride, double* dest,const UInt_t destStride,const UInt_t size);
  static void ArrayFree(double* array);
  
  friend class QVectorView;
  friend class QVectorConstView;
  friend class QVectorC;
  friend class QMatrix;
  friend class QMatrixC;
  friend class QMatrixCArray;
  friend class QRealComplexFFTW3;
  friend class QFFTWRealComplexPlan;
  QObjectDef(QVector,4);
};
 
 
Q_END_NAMESPACE
 
// operations other than member functions
Diana::QVector operator+(double t, const Diana::QVector&v);
Diana::QVector operator-(double t, const Diana::QVector&v);
Diana::QVector operator*(double t, const Diana::QVector&v);
Diana::QVector operator/(double t, const Diana::QVector&v);
Diana::QVector abs(const Diana::QVector &v);
Diana::QVector sin(const Diana::QVector &v);
Diana::QVector cos(const Diana::QVector &v);
Diana::QVector tan(const Diana::QVector &v);
Diana::QVector exp(const Diana::QVector &v);
Diana::QVector log(const Diana::QVector &v);
Diana::QVector log10(const Diana::QVector &v);
Diana::QVector pow(const Diana::QVector &v,double p);
Diana::QVector pow(double p,const Diana::QVector &v);
Diana::QVector pow(const Diana::QVector &u,const Diana::QVector &v);
double max(const Diana::QVector &v);
double min(const Diana::QVector &v);
 
 
 
#endif