Oinject.h

Go to the documentation of this file.

 
#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 unsigned int aIndex, const unsigned int aSamplingFrequency){
    return
      amp* // signal amplitude
      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;      
 
  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;        
 
  void GenerateParameters(void);
 
 
  ClassDef(Oinject,0)  
};
 
 
#endif