/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r2/src/Detchar/GlitchVeto.h

Go to the documentation of this file.

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

#include "Coinc.h"
#include "TH2.h"
#include "TGraph.h"

using namespace std;

class GlitchVeto{

 public:

  GlitchVeto(const int averbose=0);

  virtual ~GlitchVeto(void);
  bool SetMainTriggers(ReadTriggers *Amaintriggers);

  bool SetAuxTriggers(ReadTriggers *Aauxtriggers);

  bool SetSegments(Segments *aNewSeg);

  bool ProcessGlitch2Glitch(void);

  bool SetThreshold(const string aVariableName, const double aThreshold);

  bool MakeVeto(void); 
  
  bool WriteReport(const string aRootFileName="./myreport.root");

  TH1* GetMonitor(const string MonitorName);

  double GetFullDeadTime(void);

  double GetVetoEfficiency(const bool normalized=true);

  double GetVetoDeadTime(void);

  double GetVetoUsePercentage(void);

  inline Segments* GetVetoSegments(void){ return VetoSeg; };


  inline void SetFixBins(const int nfixbins=0){ nfix=nfixbins; };

 private:
  
  int fverbose;                    
  bool status_OK;                  

  // BINNING
  bool MakeBinning(const int nfixbins=0);
  double *tbins;                   
  double *fauxbins;                
  double *fmainbins;               
  double *sauxbins;                
  double *smainbins;               
  int nfix;                        
  int ntbins;                      
  int nfauxbins;                   
  int nfmainbins;                  
  int nsauxbins;                   
  int nsmainbins;                  

  // TRIGGERS
  ReadTriggers *MainTriggers;      
  ReadTriggers *AuxTriggers;       
  Coinc *coinc;                    
  vector <int> MainIndexes;        
  vector <int> AuxIndexes;         
  Segments *Seg;                   

  // PARAMETER SPACE
  bool ResetParameterSpace(void);  
  bool CreateParameterSpace(void); 
  bool FillParameterSpace(void);   

  // MAIN CHANNEL
  int nmaintrig;                   
  TH2D *h_Mstat, *h_Mstat_cum;     
  TH1D *h_mr;                      
  TGraph *g_Msnrtime;              

  // AUX CHANNEL
  TH2D *h_Astat, *h_Astat_cum,     
    *h_dead, *h_dead_cum, *h_d;    
  TH1D *h_ar;                      
  TGraph *g_snrsnr;                
  TH2D *h_ff;                      

  // VETO TUNING
  TH2D *h_nveto, *h_nveto_cum,     
    *h_use, *h_use_cum,            
    *h_up, *h_e, *h_ed;            
  
  // THRESHOLDS
  TH1D *h_thr;                     
  TH1I *h_thr_bin;                 

  // VETO PERF
  TH1D *h_Me;                      
  TH1D *h_Ause;                    
  TH1D *h_vr;                      
  TGraph *g_snrsnr_veto;           
  TGraph *g_ff_veto;               
  TH2D *h_ff_veto;                 
  TGraph *g_Msnrtime_veto;         

  // THRESHOLDS
  double GetUPthr(const double scalefactor=2);
  double UPthr;                    
  double EDthr;                    
  double Dthr;                     
  double NVthr;                    
  double NUthr;                    

  // VETO
  Segments *VetoSeg;               
  TTree *Tfom;                     
  double fom_dt, fom_eff, fom_up;
  
  ClassDef(GlitchVeto,0)  
};

#endif

---