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

Go to the documentation of this file.

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

#include "Segments.h"
#include "TimeDelay.h"
#include "ReadCoinc.h"
#include "InjRea.h"
#include "TH2.h"

using namespace std;

class InjCoi : public ReadCoinc, public InjRea{
  
 public:
  
  InjCoi(const string aCoincFilePattern, const string aInjFilePattern, Network *aNet, const int aVerbose=0);// constructor for single detector
  virtual ~InjCoi(void);        // destructor

  // RECOVER
  bool Recover(void);
  inline void SetTimeWindow(const double aTimeWindow){ TimeWindow=fabs(aTimeWindow); };
  inline void SetFrequencyRatioMax(const double aFrequencyRatioMax){ FrequencyRatioMax=fabs(aFrequencyRatioMax); };

  // RANGE
  bool SetAmplitudeRange(const double aMin, const double aMax);
  bool SetSigmaRange(const double aMin, const double aMax);
  bool SetFrequencyRange(const double aMin, const double aMax);
  bool SetEccentricityRange(const double aMin, const double aMax);
  bool SetPolarizationRange(const double aMin, const double aMax);
  bool SetRecoveryType(const int aRecType);

  bool Reset(void);// reset all
  bool SaveOutput(const string outFileName);// save results

 protected:
  
  bool status;                // class status

  // RECOVERY
  vector<int> sel_inj;        // indexes of selected injections
  vector<int> rec_inj;        // indexes of recovered injections
  vector<int> mis_inj;        // indexes of missed injections
  vector<int> rec_c;          // indexes of coinc
  double TimeWindow;          // time window
  double FrequencyRatioMax;   // frequency ratio
  int RecType;                // recovery type: 0=exact, 1=near

  // RANGES
  double ampMin, ampMax;      // amplitude range
  double sigmaMin, sigmaMax;  // sigma range
  double f0Min, f0Max;        // frequency range
  double eccMin, eccMax;      // eccentricity range
  double psiMin, psiMax;      // polarization range

  // Monitors
  bool MakeMonitors(void);
  TH1D* Nsel_amp;// number of selected inj in amplitude bins
  TH1D* Nrec_amp;// number of recovered inj in amplitude bins
  TH1D* Nmis_amp;// number of missed inj in amplitude bins
  TH1D* Nsel_freq;// number of selected inj in frequency bins
  TH1D* Nrec_freq;// number of recovered inj in frequency bins
  TH1D* Nmis_freq;// number of missed inj in frequency bins
  TH1D* Eff_amp;// recovery efficiency in amplitude bins
  TH1D* Eff_freq;// recovery efficiency in frequency bins
  TH1D** time_rec;// time reconstruction by stream
  TH2D** freq_rec;// frequency reconstruction by stream

  ClassDef(InjCoi,0)
};

#endif

---