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

Go to the documentation of this file.

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

#include "IO.h"
#include "Otile.h"
#include "Odata.h"
#include "Date.h"
#include "InjEct.h"
#include "ffl.h"

using namespace std;


class Omicron {

 public:
  
  Omicron(const string aOptionFile);

  virtual ~Omicron(void);
  bool InitSegments(Segments *aSeg, const double aTimeOffset=0.0);

  bool MakeDirectories(const int aId = 0);

  bool NewChunk(void);

  bool NewChannel(void);

  bool LoadData(double **aDataVector, int *aSize);

  int Condition(const int aInVectSize, double *aInVect);
  
  bool Project(void);
  
  bool WriteOutput(void);
  
  void PrintStatusInfo(void);
  

  Segments* GetTriggerSegments(TH1D *aThr=NULL, const double aInfValue=1e20);
  
  inline vector <string> GetChannels(void){return fChannels;};

  inline int GetChunkDuration(void){return fChunkDuration;};

  inline int GetSegmentDuration(void){return fSegmentDuration;};

  inline int GetOverlapDuration(void){return fOverlapDuration;};

  inline int GetSampleFrequency(void){return fSampleFrequency;};
 
  void PrintMessage(const string aMessage);

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

  inline int GetVerbosity(void){ return fVerbosity; };

  inline string GetVersion(void){ return "v2r1"; };

 private:

  // PROCESS STATUS
  bool status_OK;               
  time_t timer;                 
  time_t timer_start;           
  struct tm * ptm;              
  int chanindex;                
  double timeoffset;            

  // OPTIONS
  void ReadOptions(void);       
  string fOptionFile;           
  vector <string> fOptionName;  
  vector <string> fOptionType;  
  string fFflFile;              
  vector <string> fChannels;    
  int fSampleFrequency;         
  string fMaindir;              
  int fVerbosity;               
  string fOutProducts;          
  string fOutFormat;            
  string fOutStyle;             
  int fChunkDuration;           
  int fSegmentDuration;         
  int fOverlapDuration;         
  vector <double> fFreqRange;   
  vector <double> fQRange;      
  double fMismatchMax;          
  double fSNRThreshold_trigger; 
  double fSNRThreshold_map;     
  int fNTriggerMax;             
  string fClusterAlgo;          
  double fcldt;                 
  int fTileDown;                
  vector <int> fWindows;        
  int fsnrscale;                
  vector <string> fInjChan;     
  vector <double> fInjFact;     
  string fInjFilePat;           

  // PROCESS MONITORING
  Segments *inSegments;         
  Segments **outSegments;       
  int chunk_ctr;                
  int *chan_ctr;                
  int *chan_data_ctr;           
  int *chan_cond_ctr;           
  int *chan_proj_ctr;           
  int *chan_write_ctr;          
  double *chan_mapsnrmax;       
  vector <int> mapcenter;       
  vector <int> chunkstart;      
  vector <int> chunkstop;       

  // COMPONENTS
  Odata *dataseq;               
  Spectrum *spectrum;           
  ffl *FFL;                     
  Otile *tile;                  
  MakeTriggers **triggers;      
  InjEct **inject;              

  // DATA VECTORS
  double *ChunkVect;            
  double *SegVect;              
  
  // CONDITIONING
  bool Whiten(double **aDataRe, 
          double **aDataIm);
  double* GetTukeyWindow(const int aSize, const int aFractionSize); 
  double *TukeyWindow;          
  fft *offt;                    
  double **dataRe;              
  double **dataIm;              

  // OUTPUT
  string maindir;               
  vector <string> outdir;       
  void SaveAPSD(const string type="PSD");
  void SaveTS(void);            
  void MakeHtml(void);          

  // MISC
  GwollumPlot *GPlot;           
  void PrintASCIIlogo(void);    
  static const string colorcode[17];
  string GetColorCode(const double aSNRratio);


  ClassDef(Omicron,0)  
};

#endif

---