/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r2/src/Streams/Sample.h

Go to the documentation of this file.

//  Author : florent robinet (LAL - Orsay): robinet@lal.in2p3.fr
#ifndef __Sample__
#define __Sample__

#include "CUtils.h"
#include "TMath.h"

#define NORDER 20
#define NORDERHP 12

using namespace std;

class ZPGFilter {

 public:
  friend class Sample;  
  friend class IIRFilter;

 private:

  ZPGFilter(const unsigned int numZeros, const unsigned int numPoles);
  virtual ~ZPGFilter(void);

  bool WToZ(void);     
  
  unsigned int nzeros; 
  unsigned int npoles; 
  double *zeros[2];    
  double *poles[2];    
  double gain[2];      
};


class IIRFilter {

 public:
  friend class Sample;           

 private:

  IIRFilter(ZPGFilter *zpgFilter);
  virtual ~IIRFilter(void);       

  unsigned int nzeros;            
  unsigned int npoles;            
  double *zeros[2];               
  double *poles[2];               
  double gain[2];                 
  unsigned int numDirect;         
  unsigned int numRecurs;         
  double *directCoef;             
  double *recursCoef;             
  double *history;                
  double *history_init;           

};


class Sample {

 public:

  Sample(const int aVerbosity=0);

  virtual ~Sample(void);
  void SetDCRemoval(const bool aDCremove=true){ fRemoveDC=aDCremove; };

  bool SetWindows(const double aTukeyFraction);

  bool SetFrequencies(const int aNativeFrequency, const int aWorkingFrequency, const double aHighPassFrequency=0.0);

  bool SetNativeFrequency(const int aNativeFrequency);

  bool SetWorkingFrequency(const int aWorkingFrequency);

  bool SetHighPassFrequency(const double aHighPassFrequency);

  bool Transform(const int invectorsize, double *invector, const int outvectorsize, double *outvector);

  inline int GetNativeFrequency(void){ return fNativeFrequency; };

  inline int GetWorkingFrequency(void){ return fWorkingFrequency; };

  inline double GetHighPassFrequency(void){ return fHighPassFrequency; };

  inline double GetTukeyFraction(void){ return fTukeyFraction; };


 protected:

  bool fRemoveDC;              
  double fTukeyFraction;       
  int fWorkingFrequency;       
  int fNativeFrequency;        
  int fVerbosity;              
  double fHighPassFrequency;   

 private:

  // RESAMPLING
  bool MakeSamplingFilters(void);  
  bool Resample(int invectorsize, double *invector, double *outvector);
  bool presample;                  
  int fResampleFactor;             
  int fIntermediateFrequency;      
  IIRFilter *iirFilter[NORDER];    
  bool PreSample(int &datavectorsize, double *datavector); 
  
  
  // HIGHPASS
  bool MakeHighPassFilters(void);  
  bool HighPass(const int invectorsize, double *invector);
  IIRFilter *iirFilterHP[NORDERHP];

  // FILTERING
  bool ApplyFilter(const int datavectorsize, double *datavector, IIRFilter *filter); 

  // WINDOWING
  bool ApplyTukeyWindow(const int datavectorsize, double *datavector);

  // DC REMOVAL
  void RemoveDC(const int datavectorsize, double *datavector);
  
  ClassDef(Sample,0)  

};

#endif

---