/afs/in2p3.fr/throng/virgo/virgoDev/Omicron/v1r3/src/Oqplane.h

Go to the documentation of this file.

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

#include "CUtils.h"
#include "Triggers.h"
#include "Spectrum.h"
#include "FFT.h"
#include "TFile.h"
#include "TGraph.h"
#include "TH2.h"


#define NFROWMAX 1000


using namespace std;

class FreqRow {

 public:
  friend class Oqplane;  
  
 private:
  
  FreqRow(const int aTimeRange, 
      const int aTimePad, 
      const int aSampleFrequency, 
      const double aF, 
      const double aQ, 
      const double aMismatchStep, 
      const double aSNRThreshold);
  virtual ~FreqRow(void);

  double* GetSNRs(double *aDataRe, double *aDataIm);// Get SNRs in FD
  bool GetTriggers(Triggers *aTriggers, double *aDataRe, double *aDataIm, const int aStartTime);
  inline void SetPower(const double aPower){ fPower=aPower; };
  
  // PARAMETERS
  double fF;                
  double fQ;                
  int fTimeRange;           
  int fTimePad;             
  int fSampleFrequency;     
  double fMismatchStep;     
  double fSNRThreshold;     

  // DERIVED PARAMETERS
  double fBandWidth;        
  double fDuration;         

  // F-ROW
  vector <double> fTime;    
  int fNumberOfTiles;       
  bool *ValidIndices;       
  double fPower;            

  // Q-TRANSFORM
  vector <double> fWindow;  
  vector <double> fWindowFrequency; 
  vector <int> fDataIndices;
  int fZeroPadLength;       
  fft *offt;                
};


class Oqplane {

 public:
  friend class Otile;  

 private:

  Oqplane(const double aQ, 
      const int aSampleFrequency, 
      const int aTimeRange, 
      const int aTimePad, 
      const double aFrequencyMin, 
      const double aFrequencyMax, 
      const double aMismatchStep, 
      const double aSNRThreshold);
  virtual ~Oqplane(void);
    
  bool GetTriggers(Triggers *aTriggers, double *aDataRe, double *aDataIm, const int aTimeStart);
  bool SetPowerSpectrum(Spectrum *aSpec);
  TH2D* GetMap(double *aDataRe, double *aDataIm, const double time_offset=0.0, const bool printamplitude=false);
  
  // STATUS
  bool status_OK;                  

  // PARAMETERS
  double fQ;
  int fTimeRange;                  
  int fTimePad;                    
  int fSampleFrequency;            
  double fFrequencyMin,            
    fFrequencyMax;                 
  double fMismatchStep;            
  double fSNRThreshold;            

  // DERIVED PARAMETERS
  double fQPrime;                  

  // Q-PLANE
  TH2D *hplane;                    
  int fNumberOfRows;               
  double fPlaneNormalization;      
  vector <double> fFreq;           
  FreqRow *freqrow[NFROWMAX];      
  int fNumberOfTiles;              
  
  bool CheckParameters(void);      
  void GetPlaneNormalization(void);
  void BuildTiles(void);

  ClassDef(Oqplane,0)  
};


#endif

---