/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r1/src/Streams/Spectrum.h

Go to the documentation of this file.

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

#include "CUtils.h"
#include "FFT.h"
#include "TMath.h"
#include "TGraph.h"
#include "TH1.h"
#include <gsl/gsl_errno.h>
#include <gsl/gsl_spline.h>

using namespace std;

class Spectrum{

 public:

  Spectrum(const int aSamplingFrequency, 
       const int aNpoints,
       const int aPaddingDuration, 
       const int aVerbosity=0);

  virtual ~Spectrum(void);
  bool LoadData(const int aDataSize, double *aData);

  double GetPower(const double aFrequency);
  
  TGraph* GetPSD(void);

  TGraph* GetASD(void);

  inline bool GetStatus(void){ return status_OK; };

 private:

  bool status_OK;              
  int fVerbosity;              
   
  int fSamplingFrequency;      
  int PaddingDuration;         
  int PaddingSize;             

  int DataSize;                
  double *Data;                
  double *Datain;              
  double *HannWindow;          
  int fPsdSize;                
  double *PSD;                 
  double *f_PSD;               
  fft *fft_psd;                

  gsl_interp_accel *acc;       
  gsl_spline *interp_psd;      
  
  double* GetHannWindow(const int aSize);
  bool ComputePSD(void);

  ClassDef(Spectrum,0)  
};

#endif

---