/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r2/src/Inject/InjRec.h

Go to the documentation of this file.

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

#include "TObject.h"
#include "EventMap.h"
#include "InjRea.h"

using namespace std;

class InjRec : public InjRea, public EventMap {
  friend class InjCoi;

 public:
  
  InjRec(const string aTriggerFilePattern, const string aInjFilePattern, const int aVerbose=0);

  virtual ~InjRec(void);
  int Recover(const int aInjIndex=-1, const int aClusterTagThr=0);

  void Draw(const int aPlotIndex, const string aOptions="", const int aPadIndex=0);
  
  void PrintMissedInjection(const double aAmplitudeThr=0.0);

  inline int GetNValidInjections(void){
    return ValidInjTree->GetEntries();
  }

  inline int GetNRecInjections(void){
    return (int)rec_injindex.size();
  }

  inline vector <int> GetRecInjections(void){
    return rec_injindex;
  }

  inline int GetNMissedInjections(void){
    return (int)mis_injindex.size();
  }
  
  inline double GetInjectedAmplitude(void){
    double tlocal=GetLocalTime(inj_ra,inj_dec,inj_time);
    double fplus, fcross;
    GetDetAMResponse(fplus, fcross, inj_ra, inj_dec, inj_psi, GreenwichMeanSiderealTime(tlocal));
    return fplus*GetSineGaussh0plus()+fcross*GetSineGaussh0cross();
  };


 private:
  
  TTree *ValidInjTree;       
  int injindex;              

  vector <int> rec_injindex; 
  vector <int> rec_clindex;  
  vector <int> mis_injindex; 

  TH1D *h_inj_amp;           
  TH1D *h_inj_f0;            

  TH1D *h1_rec;              
  TH2D *h2_rec;              
  TGraph *g_rec;             

  // make diagnostic plots
  void make_effamp(void);
  void make_efff0(void);
  void make_ampamp(void);
  void make_f0f0(void);

  ClassDef(InjRec,0)
};

#endif

---