Streams.h

Go to the documentation of this file.

 
#ifndef __Streams__
#define __Streams__
 
#include <TH2.h>
#include "TimeDelay.h"
#include "Sample.h"
#include "DetConst.h"
 
using namespace std;
 
class Streams: public Sample{
  
 public:
  
  Streams(const string aName, const unsigned int aVerbose=0);
  
  virtual ~Streams(void);
  void SetName(const string aName);
 
  inline string GetName(void){ return Name; };
 
  inline string GetNamePrefix(void){ return NamePrefix; };
 
  inline string GetNameSuffix(void){ return NameSuffix; };
 
  inline string GetNameSuffixUnderScore(void){ return NameSuffixUnderScore; };
 
  inline string GetNameConv(void){ return NamePrefix+"-"+NameSuffixUnderScore; };
 
  inline string GetDirectory(const string aOutDir="."){
    return aOutDir+"/"+Name;
  };
    
  inline string GetTriggerFileName(const unsigned int aGpsStart, const unsigned int aDuration, const string aFileFormat="root", const string aProcessName="PROC", const string aOutDir="."){
    gwl_ss<<aOutDir<<"/"<<GetNameConv()<<"_"<<aProcessName<<"-"<<aGpsStart<<"-"<<aDuration<<"."<<aFileFormat;
    string ss_ret = gwl_ss.str();
    gwl_ss.clear(); gwl_ss.str("");
    return ss_ret;
  };
    
  inline unsigned int GetDetectorIndex(void){ return DetIndex; };
 
  inline string GetDetectorPrefix(void){ return DET_PREFIX[DetIndex]; };
 
  void GetDetectorAMResponse(double &aFplus, double &aFcross,
                             const double aRa, const double aDec,
                             const double aPsi, const double aGmst);
  
  inline double GetLocalTime(const double aRa, const double aDec, const double aGeocentricTime){
    return aGeocentricTime+TimeDelayFromEarthCenter(DetLoc, aRa, aDec, aGeocentricTime);
  };
 
  bool GetSkyRingPoint(double &aThetaS, double &aPhiS, const double aGmst, const double aDeltaT, const double aOmegaS, Streams *aDprime);
 
  inline bool GetSkyRingPointRaDec(double &aDecS, double &aRaS, const double aGmst, const double aDeltaT, const double aOmegaS, Streams *aDprime){
    if(GetSkyRingPoint(aDecS, aRaS, aGmst, aDeltaT, aOmegaS, aDprime)==false) return false;
    aDecS = TMath::Pi()/2.0 - aDecS;
    return true;
  };
  
  TH2D* GetSkyAntennaFactor(const unsigned int aGpsTime);
 
protected:
  
  string Name;                 
  
private:
  
  unsigned int verbose;        
  string NamePrefix;           
  string NameSuffix;           
  string NameSuffixUnderScore; 
  unsigned int DetIndex;       
  double Response[3][3];       
  double DetLoc[3];            
 
  void MakeDetector(void);
 
  void ComputeDetectorResponse(void);
 
  ClassDef(Streams,0)
};
 
#endif