RFSignal.h

#ifndef RFSIGNAL_H
#define RFSIGNAL_H
#include "TGraph.h"
#include "FFTWComplex.h"

class RFFilter;


class RFSignal : public TGraph {
public:
  RFSignal(); 
  RFSignal(TGraph *grWave, Int_t mvNs=0); 
  RFSignal(RFSignal *rfWave); 
  RFSignal(Int_t numPoints,const Double_t *tVals, const Double_t *vVals, Int_t mvNs=0);
  RFSignal(Int_t numFreqs, const  Double_t *freqVals, const FFTWComplex *complexNums, Int_t mvNs=0);
  ~RFSignal(); 
  
  TGraph *getFreqMagGraph();
  Double_t *getFreqs();
  Double_t *getMags();
  Double_t *getPhases();
  void setFreqs(Int_t nfreq, Double_t *freqs);
  void setMagsPhases(Double_t *mags, Double_t *phases);
  FFTWComplex *getComplexNums();
  Int_t getNumFreqs();
  void addToSignal(RFSignal *grSignal);
  void applyFilter(RFFilter *theFilter);
  void updateTimeDomain();//because we could operate on freq domain once it has already been extracted!
  
 private:  
  void fillFreqStuff();
  void extractFromComplex();
  void reExtractFromComplex();
  Int_t fGotFreqs;
  Int_t fNumFreqs;
  Int_t fMvNs;
  Double_t *fFreqs;  //[fNumFreqs]
  Double_t *fMags; //[fNumFreqs]
  Double_t *fPhases; //[fNumFreqs]
  FFTWComplex *fComplexNums; //[fNumFreqs]
  
ClassDef(RFSignal,2) 

};

#endif // RFSIGNAL_H