/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/Omicron/v2r2/src/Oinject.h

Go to the documentation of this file.

//  Author : florent robinet (LAL - Orsay): robinet@lal.in2p3.fr
#ifndef __Oinject__
#define __Oinject__
#include "Spectrum.h"
#include "TRandom3.h"
#include "TMath.h"

#define sqrt2pi 2.5066282746310002416123552393401042

using namespace std;

class Oinject{

 public:
  
  Oinject(const double aDuration);

  virtual ~Oinject(void);
  void MakeWaveform(void);

  inline double GetWaveform(const int aIndex, const int aSamplingFrequency){
    return
      amp*
      Wg/sigma_t/sqrt2pi*exp(-(-duration/2.0+(double)aIndex/(double)aSamplingFrequency-tau)*(-duration/2.0+(double)aIndex/(double)aSamplingFrequency-tau)/2.0/sigma_t/sigma_t)* // Gaussian window
      TMath::Cos(2.0*TMath::Pi()*phi*(-duration/2.0+(double)aIndex/(double)aSamplingFrequency)+phase);// sine
  };

  double GetTrueSNR(Spectrum *aSpec1, Spectrum *aSpec2);

  inline void SetTimeRange(const double aTimeMin, const double aTimeMax){
    taumin=aTimeMin; taumax=aTimeMax;
  };

  inline void SetFrequencyRange(const double aFreqMin, const double aFreqMax){
    phimin=aFreqMin; phimax=aFreqMax;
  };

  inline void SetAmplitudeRange(const double aAmpMin, const double aAmpMax){
    ampmin=aAmpMin; ampmax=aAmpMax;
  };

  inline void SetQRange(const double aQMin, const double aQMax){
    Qmin=aQMin; Qmax=aQMax;
  };

  inline double GetTime(void){ return tau; };

  inline double GetTimeMin(void){ return taumin; };

  inline double GetTimeMax(void){ return taumax; };

  inline double GetFrequency(void){ return phi; };

  inline double GetFrequencyMin(void){ return phimin; };

  inline double GetFrequencyMax(void){ return phimax; };

  inline double GetQ(void){ return Q; };

  inline double GetQMin(void){ return Qmin; };

  inline double GetQMax(void){ return Qmax; };

  inline double GetAmplitude(void){ return amp; };

  inline double GetAmplitudeMin(void){ return ampmin; };

  inline double GetAmplitudeMax(void){ return ampmax; };

  inline double GetPhase(void){ return phase; };

  inline double GetSigmat(void){ return sigma_t; };

  inline double GetSigmaf(void){ return sigma_f; };

 private:
  
  double duration;                        
  TRandom3 *randgen;                   

  double Wg;                           
  double sigma_t;                      
  double sigma_f;                      

  void GenerateParameters(void);       
  double tau;                          
  double taumin;                       
  double taumax;                       
  double phi;                          
  double phimin;                       
  double phimax;                       
  double Q;                            
  double Qmin;                         
  double Qmax;                         
  double amp;                          
  double ampmin;                       
  double ampmax;                       
  double snr;                          
  double phase;                        

  ClassDef(Oinject,0)  
};


#endif

---