/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r2/src/Triggers/EventMap.h

Go to the documentation of this file.

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

#include "GwollumPlot.h"
#include "ReadTriggers.h"
#include "TriggerSelect.h"
#include "TH2.h"

using namespace std;

class EventMap: public GwollumPlot, public ReadTriggers {
  
 public:
  
  EventMap(const int aNmaps, const string aPattern, const string aDirectory="", const int aVerbose=0);

  virtual ~EventMap(void);
  bool BuildMap(const int aMapIndex, const double aTime);

  bool BuildClusterMap(const int aMapIndex, const int aCl, const double aTimeCenter=0.0);

  bool PrintMap(const int aMapIndex);

  bool PrintMapNext(const double aTimeFraction=0.25);
  
  bool PrintMapLive(const double aDuration, const string aOutFile="none");

  inline void SetMapAmplitudeContent(void){ MapContentType=true; ResetAll(); };

  inline void SetMapSNRContent(void){ MapContentType=false; ResetAll(); };

  bool SetMapTimeRange(const int aMapIndex, const double aTimeRange);
 
  bool SetMapFrequencyRange(const int aMapIndex, const double aFreqMin, const double aFreqMax);
 
  bool SetMapQRange(const int aMapIndex, const double aQMin, const double aQMax);
 
  bool SetMapResolution(const int aMapIndex, const int aTresolution, const int aFresolution);
 
  double GetMapLoudestSNR(const int aMapIndex);

  double GetMapLoudestTime(const int aMapIndex);

  double GetMapLoudestFrequency(const int aMapIndex);

  bool NormalizeMap(const int aMapIndex);

  bool NormalizeMapLoudest(const int aMapIndex);

  double GetMapMatch(double &aTimeShift, double &aMatchSumw2, const int aMapIndex1, const int aMapIndex2, const double aTimeJitter=0.0);

  double GetMapMatch(double &aTimeShift, double &aMatchSumw2, const int aMapIndex, EventMap *aExtMap, const int aExtMapIndex, const double aTimeJitter=0.0);

  inline int GetMapTimeResolution(const int aMapIndex){
    if(aMapIndex>=0&&aMapIndex<Nmaps) return Map[aMapIndex]->GetNbinsX();
    return -1;
  }

  inline int GetMapFrequencyResolution(const int aMapIndex){
    if(aMapIndex>=0&&aMapIndex<Nmaps) return Map[aMapIndex]->GetNbinsY();
    return -1;
  }
    
  inline int GetMapTimeRange(const int aMapIndex){
    if(aMapIndex>=0&&aMapIndex<Nmaps) return MapParam[aMapIndex]->GetTimeRange();
    return -1;
  }

  inline int GetMapFrequencyMin(const int aMapIndex){
    if(aMapIndex>=0&&aMapIndex<Nmaps) return MapParam[aMapIndex]->GetFrequencyMin();
    return -1;
  }

  inline int GetMapFrequencyMax(const int aMapIndex){
    if(aMapIndex>=0&&aMapIndex<Nmaps) return MapParam[aMapIndex]->GetFrequencyMax();
    return -1;
  }

 private:

  void ResetMap(const int aMapIndex);// delete map and rebuild it with new parameters
  void ResetMapContent(const int aMapIndex);// reset map content
  inline void ResetAll(void){
    for(int m=0; m<Nmaps; m++) ResetMap(m);
  };// reset all maps

  int Nmaps;                
  TH2D **Map;               
  double ***MapSNR;         
  double **MapNoiseAmp;     
  TriggerSelect **MapParam; 
  double *MapCentralTime;   
  int *MapLoudestBin;       
  double **MapLoudestSNR_bf;
  int MapPrinted;           
  bool MapContentType;      
  int *MapBinXmin,          
    *MapBinXmax,            
    *MapBinYmin,            
    *MapBinYmax;            

  double GetMatch(double &aTimeShift,
          double &aSumw2,
          EventMap *aEMap1, const int aMapIndex1,
          EventMap *aEMap2, const int aMapIndex2,
          const double aTimeJitter);

  ClassDef(EventMap,0)  
};

#endif

---