/afs/in2p3.fr/throng/virgo/virgoDev/GWOLLUM/v1r3/src/Inject/Inject.h

Go to the documentation of this file.

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

#include "Chain.h"
#include "TFile.h"
#include "TTree.h"
#include "TGraph.h"
#include "TimeDelay.h"
#include "Segments.h"
#include "Network.h"
#include "InjRea.h"

// this defines the size of the injection vector.
// this must be big enough to contain the full injection waveform...
// ... but not too big not to slow down the process
#define INJDURATIONMAX 10 // [s]

using namespace std;

class Inject {
  
 public:
  
  Inject(Network *aNet, const int aSampling, const int aVerbose=0); // constructor
  virtual ~Inject();        // destructor

  // INJECTION SET
  bool SetInjectionSet(const string aInjRootFilePattern);
  inline string GetInjectionSet(void){ return InjRootFilePattern; };
  bool InjectSet(const int aDetIndex, const int aDataSize, double *aData, const double aTimeStart);

  // INJECT
  bool InjectSignal(const int aDetIndex, const int aDataSize, double *aData, const double aTimeStart);
  bool SaveInjectionSignal(const string aRootFileName, const string aOption="recreate");
  bool SaveDetectorSignal(const string aRootFileName, const string aOption="recreate");

  // SIGNAL GENERATION
  bool SetSource(const double aTime, const double aRa, const double aDec, const double aPsi=0, const double aEcc=0);
  bool GenerateGauss(const double aAmplitude, const double aSigma);
  bool GenerateSineGauss(const double aAmplitude, const double aSigma, const double aFrequency);
  
 protected:
  
  bool status;              // status of the Inject object
  int Verbose;              // verbosity level

  // PARAMETERS
  Network *Net;             // network - DO NOT DELETE
  int fSampling;            // sampling frequency
  int DataSize;             // data vector size

  // INJECTION SET
  InjRea *InjSet;           // injection set
  string InjRootFilePattern;// root file name

  // INJECTION
  int InjSize;              // vector size for injection waveform
  double *hplus;            // injection waveform (+)
  double *hcross;           // injection waveform (x)
  double *hdet;             // injection waveform projected on the detector
  double *hset;             // set of injection waveform projected on the detector
  double inj_time;          // Geocentric time of injection
  double inj_ra;            // right-ascension
  double inj_dec;           // declination
  double inj_psi;           // psi
  double inj_ecc;           // eccentricity

  bool ProjectSignal(const int aDetIndex, const double aTimeStart);// fill hdet[]

  // WINDOW  
  double *window;
  bool CreateTukeyWindow(const int size);


  ClassDef(Inject,0)
};

#endif

---