EdbMath.h

Go to the documentation of this file.

#ifndef ROOT_EdbMath
#define ROOT_EdbMath
 
//                                                                      //
// EdbMath                                                              //
//                                                                      //
// Collection of matematical algorithns                                 //
//                                                                      //
#include "TObject.h"
#include "TArrayI.h"
#include "TArrayD.h"
#include "TH2D.h"
 
//______________________________________________________________________________
class EdbMath {
 
 public:
  EdbMath(){}  
  virtual ~EdbMath(){}
 
  static double Magnitude3( float Point1[3], float Point2[3] );
  static double DistancePointLine3( float Point[3], 
                    float LineStart[3], float LineEnd[3], 
                    bool inside );
 
  static double Angle3( float tx1, float ty1, float tx2, float ty2);
 
  static bool LineLineIntersect( float p1[3], float p2[3], float p3[3], float p4[3],
                 float pa[3], float pb[3],
                 double &mua, double &mub );
 
  static bool LinFitDiag( int n, float *x, float *y, float *e, 
              float p[2], float d[2][2], float *chi2);
  static bool LinFitCOV( int n, float *x, float *y, double *c, 
             float *p, float *d, float *chi2);
  static void LFITW( float *X, float *Y, float *W, int L, int KEY, float &A, float &B, float &E );
 
  static int LFIT3( float *X, float *Y, float *Z, float *W, int L, 
            float &X0, float &Y0, float &Z0, float &TX, float &TY, float &EX, float &EY );
  
  static int ArrStat( int n, float *x, float par[4] );
 
  ClassDef(EdbMath,2)  // general matematical algorithms
};
 
//______________________________________________________________________________
class TIndex2 : public TArrayD {
 private:
 
 public:
  TIndex2() : TArrayD() {}
  TIndex2(int size) : TArrayD(size) {}
  virtual ~TIndex2(){}
  
  static Double_t BuildValue( Int_t major, Int_t minor ) {return major+minor*1e-9;}
  void   SetMinor(int i, Int_t minor) { (*this)[i] = BuildValue( Major(i), minor ); }
 
  Int_t  Major(int i) { return (Int_t)((*this)[i]); }
  Int_t  Minor(int i) { return (Int_t)(((*this)[i]- (Double_t)Major(i)+1e-10)*1e+9 ); }
  void   BuildIndex( int n, double *w );
  Int_t  Find(Int_t major, Int_t minor);
  Int_t  FindIndex(Int_t major);
  Int_t  FindIndexArr(Int_t major, TArrayI &mina );
  void   Print();
 
  ClassDef(TIndex2,1)  // simple index class for fast search in linear cases
};
 
//______________________________________________________________________________
class EdbFilter2D {
 
  private:
    Int_t   eKX,eKY;
    Double_t *eKernel; 
 
    public:
      EdbFilter2D(){eKernel=0;eKX=0;eKY=0;}
      virtual ~EdbFilter2D(){}
    
     void SetKernel(int k);
     void Smooth0(TH2D &h);
     void Smooth(TH2D &h);
     void Print();
    
     ClassDef(EdbFilter2D,1)  // 2D histogram filtering (i.e. smoothing)
};
 
#endif /* ROOT_EdbMath */