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

Go to the documentation of this file.

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

#include "GwollumPlot.h"
#include "IO.h"
#include "ffl.h"
#include "Streams.h"
#include "Spectrum.h"
#include "Sample.h"
#include "Triggers.h"
#include "EventMap.h"
#include "TMath.h"
#include "Otile.h"
#include "Odata.h"
#include "Date.h"

using namespace std;


class Omicron {

 public:
  
  Omicron(const string aOptionFile);

  virtual ~Omicron(void);
  bool Process(Segments *aSeg);

  bool Scan(const double aTimeCenter);

  bool ScanTriggers(const double aTimeCenter);

  bool InitSegments(Segments *aSeg);

  bool MakeDirectories(const double aGPS = 0);

  bool NewChunk(void);

  int ConditionVector(const int aChNumber, const int aInVectSize, double *aInVect);
  
  int MakeTriggers(const int aChNumber);
  
  bool WriteTriggers(const int aChNumber);
  
  bool MakeMaps(const int aChNumber, const double aTimeCenter);
  
  bool WriteMaps(void);

  bool ScanReport(const string aScanDir="");


  //Segments* GetOnlineSegments(const int aChNumber, TH1D *aThr, const double aPadding=0.0, const double aInfValue=1e20);
  
  inline vector <string> GetChannelList(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 PrintStatusInfo(void);
  
  inline bool GetStatus(void){ return status_OK; };

 private:

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

  // INPUT OPTIONS
  bool ReadOptions(void);       
  string fOptionFile;           
  vector <string> fOptionName;  
  vector <string> fOptionType;  
  string fMaindir;              
  vector <string> fOutdir;      
  vector <string> fChannels;    
  vector <string> fInjChan;     
  vector <double> fInjFact;     
  string fFflFile;              
  string fTrigDir;              
  int fSampleFrequency;         
  vector <double> fFreqRange;   
  vector <double> fQRange;      
  int fChunkDuration;           
  int fSegmentDuration;         
  int fOverlapDuration;         
  double fMismatchMax;          
  vector <int> fWindows;        
  int fWindowMax;               
  string ffftplan;              
  double fSNRThreshold;         
  int fNtriggerMax;             
  vector <string> fClusterAlgo; 
  double fcldt;                 
  int fVerbosity;               
  string fOutFormat;            
  string fStyle;                
  bool writepsd;                
  bool writetimeseries;         
  string fWriteMode;            

  // PROCESS MONITORING
  Segments *inSegments;         
  Segments **outSegments;       
  int *chunk_ctr;               
  int *cor_chunk_ctr;           
  int *cor_data_ctr;            
  int *max_chunk_ctr;           

  // COMPONENTS
  Odata *dataseq;               
  Sample **sample;              
  Streams **streams;            
  Spectrum *spectrum;           
  ffl *FFL;                     
  Otile *tile;                  
  Triggers **triggers;          

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

  // SCANS
  int MapChNumber;              
  TH2D **Qmap;                  
  TH2D **Qmap_full;             
  double Qmap_center;           
  int *loudest_qmap;            
  GwollumPlot *GPlot;           
  string fScandir;              

  // MISC
  void SaveAPSD(const int c, const string type="PSD");    
  void SaveTS(const int c, double tcenter=0);
  bool *first_save;             
  void PrintASCIIlogo(void);    

  ClassDef(Omicron,0)  
};

#endif

---