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

Go to the documentation of this file.

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

#include "TriggerBuffer.h"
#include "Spectrum.h"
#include "Omap.h"

// eq 5.95 with alpha=2
#define BIASFACT2 1.1140649371721838001292326225666329264640808105469

using namespace std;

class Oqplane: public Omap {

 public:
  friend class Otile;  

 private:
  
  // constructor - destructor
  Oqplane(const double aQ, 
      const int aSampleFrequency, 
      const int aTimeRange, 
      const double aFrequencyMin, 
      const double aFrequencyMax, 
      const double aMismatchStep);
  virtual ~Oqplane(void);

  void PrintParameters(void);
  int ProjectData(fft *aDataFft, const double aPadding=0.0);
  void FillMap(const string aContentType="snr");
  bool SaveTriggers(TriggerBuffer *aTriggers, const double aT0, Segments* aSeg);

  // GETS
  inline double GetQ(void){ return Q; };
  inline double GetSNRThr(void){ return SNRThr; };
  inline double GetTileSNR2(const int aTimeTileIndex, const int aBandIndex){
    return TMath::Max(bandFFT[aBandIndex]->GetNorm2_t(aTimeTileIndex)-2.0,0.0);
  };
  inline double GetTileAmplitude(const int aTimeTileIndex, const int aBandIndex){
    return sqrt(GetTileSNR2(aTimeTileIndex,aBandIndex))*bandNoiseAmplitude[aBandIndex];
  };

  // SETS
  bool SetPower(Spectrum *aSpec1, Spectrum *aSpec2);
  inline void SetSNRThr(const double aSNRThr){ SNRThr=aSNRThr; };


  // INTERNAL
  double GetMeanEnergy(const int aBandIndex, const double aPadding);
  double GetA1(void);
 
  // Q-PLANE
  double Q;                         
  double QPrime;                    
  double PlaneNormalization;        

  // TRIGGER SELECTION
  double SNRThr;                    
    
  // FREQUENCY BANDS
  int *bandWindowSize;              
  double **bandWindow;              
  double *bandNoiseAmplitude;       
  fft **bandFFT;                    
    
  ClassDef(Oqplane,0)  
};


#endif

---