/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/Omicron/v2r2/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 "Oinject.h"
#include "Date.h"
#include "InjEct.h"
#include "ffl.h"
#include "TRandom.h"
#include "TRandom3.h"
#include "TriggerBuffer.h"

using namespace std;


class Omicron {

 public:
  
  Omicron(const string aOptionFile);

  virtual ~Omicron(void);
  bool InitSegments(Segments *aInSeg, Segments *aOutSeg=NULL);

  bool MakeDirectories(const int aId = 0);

  bool NewChunk(void);

  bool DefineNewChunk(const int aTimeStart, const int aTimeEnd, const bool aResetPSDBuffer=false);

  bool NewChannel(void);

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

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

  string WriteTriggers(const bool aLVDirConvention = false);

  bool WriteOutput(void);

  void PrintStatusInfo(void);
  
  Segments* GetTriggerSegments(TH1D *aThr=NULL, const double aInfValue=1e20);
  
  vector <string> GetChannels(void);

  inline int GetChunkDuration(void){return tile->GetTimeRange();};

  inline int GetOverlapDuration(void){return tile->GetOverlapDuration();};

  inline int GetSampleFrequency(void){return triggers[0]->GetWorkingFrequency();};
 
  inline int GetNTriggers(void){ return triggers[chanindex]->GetNtriggers(); };

  
  inline double GetTriggerBufferLiveTime(void){ return triggers[chanindex]->GetBufferLiveTime(); };

  inline void ResetTriggerBuffer(void){ return triggers[chanindex]->ResetBuffer(); };

  inline void ResetTriggerBufferAfter(const double aTime){ return triggers[chanindex]->ResetAfter(aTime); };

  inline void ResetPSDBuffer(void) { spectrum1[chanindex]->Reset(); };

  void PrintMessage(const string aMessage);

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

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

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

 private:

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

  // OPTIONS
  void ReadOptions(void);       
  string fOptionFile;           
  int fVerbosity;               
  string fMaindir;              
  string fOutFormat;            
  string fOutProducts;          
  bool fNoLogo;                 
  vector <int> fWindows;        
  string fClusterAlgo;          
  string fftplan;               
  double fratemax;              
  vector <string> fInjChan;     
  vector <double> fInjFact;     
  int fsginj;                   

  // 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;          
  int *trig_ctr;                
  double *chan_mapsnrmax;       
  vector <int> chunkcenter;     
  vector <string> chunktfile;   
  ofstream oinjfile;            
  
  // COMPONENTS
  int nchannels;                
  GwollumPlot *GPlot;           
  vector <string> outdir;       
  Spectrum **spectrum1;         
  Spectrum **spectrum2;         
  Spectrum *spectrumw;          
  ffl *FFL;                     
  ffl *FFL_inject;              
  fft *offt;                    
  Otile *tile;                  
  TriggerBuffer **triggers;     
  Oinject *oinj;                
  InjEct **inject;              

  // DATA VECTORS
  double *ChunkVect;            
    
  // CONDITIONING & WHITENING
  void Whiten(Spectrum *aSpec); 
  double* GetTukeyWindow(const int aSize, const int aFractionSize); 
  double *TukeyWindow;          
 
  // OUTPUT
  string maindir;               
  void SaveAPSD(const string aType);
  void SaveWPSD(void);          
  void SaveTS(const bool aWhite=false); 
  void SaveSG(void);            
  void SaveSpectral(void);      
  void MakeHtml(void);          

  // MISC
  void PrintASCIIlogo(void);    
  static const string colorcode[17];
  string GetColorCode(const double aSNRratio);
  bool IsFlat(const int aInVectSize, double *aInVect);

  ClassDef(Omicron,0)  
};

#endif

---