/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r1/src/Inject/InjRea.h

Go to the documentation of this file.

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

#include "Chain.h"
#include "TH1.h"
#include "TMath.h"

using namespace std;

class InjRea {

  friend class InjRec;

 public:
  
  InjRea(const string aPattern, const int aVerbose=0);

  virtual ~InjRea();
  inline bool LoadInjection(const int aInjIndex){if(InjTree->GetEntry(aInjIndex)>0){ InjIndex=aInjIndex; return true; } return false; };

  double GetSineGaussh0plus(void);

  double GetSineGaussh0cross(void);

  inline double GetInjectionTime(void){ return inj_time; };

  inline double GetInjectionRa(void){ return inj_ra; };

  inline double GetInjectionDec(void){ return inj_dec; };

  inline double GetInjectionEccentricity(void){ return inj_ecc; };

  inline double GetInjectionPolarization(void){ return inj_psi; };

  inline double GetInjectionAmplitude(void){ return inj_amp; };

  inline double GetInjectionAmplitudeMin(void){ return inj_ampmin; };

  inline double GetInjectionAmplitudeMax(void){ return inj_ampmax; };
  
  inline double GetInjectionSigma(void){ return inj_sigma; };

  inline double GetInjectionFrequency(void){ return inj_f0; };

  inline double GetInjectionFrequencyMin(void){ return inj_f0min; };

  inline double GetInjectionFrequencyMax(void){ return inj_f0max; };

  inline bool GetInjectionTag(void){ return inj_tag[InjIndex]; };

  inline void SetInjectionTag(const bool aNewTag){ inj_tag[InjIndex]=aNewTag; };

  int GetFirstInjectionIndex(const double aTimeStart, const double aTimeEnd);
  
  int GetLastInjectionIndex(const double aTimeStart, const double aTimeEnd);
  
  inline bool GetStatus(void){ return status_OK; };
  
  inline int GetNInjections(void){ return InjTree->GetEntries(); };
  
  inline string GetInjectionType(void){ return (string)InjTree->GetName(); };
 
  inline double GetInjectionSigmaMax(void){ return inj_sigmamax; };

  TH1D* GetInjectionParamDist(const string aParamName, const int aNbins=1, const string aBinType="UNIFORM");

 
 protected:
  
  bool status_OK;            
  int Verbose;               

  //inline int GetTreeNumber(void){ return InjTree->GetTreeNumber(); };
  inline Long64_t* GetTreeOffset(void){ return InjTree->GetTreeOffset(); };

  double inj_ra,             
    inj_dec,                 
    inj_psi,                 
    inj_time,                
    inj_amp,                 
    inj_f0,                  
    inj_sigma,               
    inj_ecc;                 

  double inj_sigmamin;       
  double inj_sigmamax;       
  double inj_ampmin;         
  double inj_ampmax;         
  double inj_f0min;          
  double inj_f0max;          
 
 private:

  Chain *InjTree;            
  bool *inj_tag;             
  int InjIndex;              
  double inj_timefirst;      
  double inj_timelast;       
 
  ClassDef(InjRea,0)
};

#endif

---