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

Go to the documentation of this file.

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

#include "CUtils.h"
#include "TMath.h"
#include "TH2D.h"

using namespace std;

class Omap: public TH2D {

 public:
  friend class Otile;   
  friend class Oqplane; 

 private:
  
  // constructor - destructor
  Omap();
  virtual ~Omap(void);

  void SetBins(const double aQ,
           const double aFrequencyMin, const double aFrequencyMax,
           const int aTimeRange, const double aMismatchStep);

  inline double GetTimeRange(void){
    return GetXaxis()->GetBinUpEdge(GetNbinsX())-GetXaxis()->GetBinLowEdge(1); 
  };
  inline double GetTimeMin(void){ 
    return GetXaxis()->GetBinLowEdge(1); 
  };
  inline double GetTimeMax(void){ 
    return GetXaxis()->GetBinUpEdge(GetNbinsX()); 
  };
  inline double GetFrequencyMin(void){ 
    return GetYaxis()->GetBinLowEdge(1); 
  };
  inline double GetFrequencyMax(void){ 
    return GetYaxis()->GetBinUpEdge(GetNbinsY()); 
  };
  inline int GetNBands(void){ 
    return GetNbinsY();
  };
  inline long int GetNTiles(void){ 
    return Ntiles;
  };
  inline int GetBandIndex(const double aFrequency){
    return GetYaxis()->FindBin(aFrequency)-1;
  };
  inline double GetBandFrequency(const int aBandIndex){
    return GetYaxis()->GetBinCenterLog(aBandIndex+1);
  };
  inline double GetBandStart(const int aBandIndex){ 
    return GetYaxis()->GetBinLowEdge(aBandIndex+1);
  };
  inline double GetBandEnd(const int aBandIndex){ 
    return GetYaxis()->GetBinUpEdge(aBandIndex+1);
  };
  inline double GetBandWidth(const int aBandIndex){ 
    return GetYaxis()->GetBinWidth(aBandIndex+1);
  };
  inline double GetTileDuration(const int aBandIndex){
    return GetXaxis()->GetBinWidth(1)*bandMultiple[aBandIndex];
  };
  inline int GetBandNtiles(const int aBandIndex){ 
    return GetNbinsX()/bandMultiple[aBandIndex];
  };
  inline double GetTileContent(const int aTimeTileIndex, const int aBandIndex){    
    return tilecontent[aBandIndex][aTimeTileIndex];
  };
  inline double GetTilePhase(const int aTimeTileIndex, const int aBandIndex){    
    return tilephase[aBandIndex][aTimeTileIndex];
  };
  inline double GetTileTag(const int aTimeTileIndex, const int aBandIndex){
    return tiletag[aBandIndex][aTimeTileIndex];
  };
  inline double GetTileTimeStart(const int aTimeTileIndex, const int aBandIndex){
    return GetXaxis()->GetBinLowEdge(aTimeTileIndex*bandMultiple[aBandIndex]+1);
  };
  inline double GetTileTimeEnd(const int aTimeTileIndex, const int aBandIndex){
    return GetXaxis()->GetBinUpEdge((aTimeTileIndex+1)*bandMultiple[aBandIndex]);
  };
  inline double GetTileTime(const int aTimeTileIndex, const int aBandIndex){
    return (GetTileTimeStart(aTimeTileIndex,aBandIndex) + GetTileTimeEnd(aTimeTileIndex,aBandIndex)) / 2.0;
  };
  inline int GetTimeTileIndex(const int aBandIndex, const double aTime){
    return (int)floor((aTime-GetTimeMin())/GetTileDuration(aBandIndex));
  };

  // SETS
  inline void SetTileContent(const int aTimeTileIndex, const int aBandIndex, const double aContent, const double aPhase=-100.0, const bool aTag=true){
    tilecontent[aBandIndex][aTimeTileIndex]=aContent;
    tilephase[aBandIndex][aTimeTileIndex]=aPhase;
    tiletag[aBandIndex][aTimeTileIndex]=aTag;
  };
  inline void SetTileTag(const int aTimeTileIndex, const int aBandIndex, const bool aTag){ tiletag[aBandIndex][aTimeTileIndex]=aTag; };

  // MAPS
  void MakeMapContent(void);
  void MakeMapPhase(void);
  void MakeMapDisplay(void);

  
  long int Ntiles;                  
  int *bandMultiple;                
  double **tilecontent;             
  double **tilephase;               
  bool   **tiletag;                 

  ClassDef(Omap,0)  
};


#endif

---