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

Go to the documentation of this file.

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

#include "Chain.h"
#include "TTreeIndex.h"
#include "Segments.h"

using namespace std;

class Triggers {

 public:

  Triggers(const int aVerbose=0);

  virtual ~Triggers(void);
  bool AddTrigger(const double aTime,   const double aFrequency,
          const double aSNR,    const double aQ,
          const double aTstart, const double aTend,
          const double aFstart, const double aFend,
          const double aAmplitude, const double aPhase);

  bool SortTriggers(void);

  bool SetTriggerBranchStatus(const string aBname, const bool aStatus = true);

  inline int GetNTrig(void){ return Ttree->GetEntries(); };

  inline int GetNtrig(void){ return Ttree->GetEntries(); };

  inline int GetNTriggers(void){ return Ttree->GetEntries(); };

  inline int GetNtriggers(void){ return Ttree->GetEntries(); };

  string GetTriggerFile(const int aTr);

  inline double GetTriggerTime(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Ttime;
    return 0.0;
  };

  inline double GetTriggerFrequency(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tfreq;
    return 0.0;
  };

  double GetTriggerSNR(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tsnr;
    return 0.0;
  };

  double GetTriggerAmplitude(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tamp;
    return 0.0;
  };

  double GetTriggerQ(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tq;
    return 0.0;
  };

  double GetTriggerPhase(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tph;
    return 0.0;
  };

  double GetTriggerTimeStart(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Ttstart;
    return 0.0;
  };

  double GetTriggerTimeEnd(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Ttend;
    return 0.0;
  };

  inline double GetTriggerDuration(const int aTr){
    return GetTriggerTimeEnd(aTr)-GetTriggerTimeStart(aTr);
  };

  double GetTriggerFrequencyStart(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tfstart;
    return 0.0;
  };

  double GetTriggerFrequencyEnd(const int aTr){
    if(Ttree->GetEntry(aTr)>0) return Tfend;
    return 0.0;
  };

  double GetTriggerBandWidth(const int aTr){
    return GetTriggerFrequencyEnd(aTr)-GetTriggerFrequencyStart(aTr);
  };

  int GetTriggerIndex(const double aTime);

  Segments* GetTriggerSegments(TH1D *aThr=NULL, const double aInfValue=1e20);

  void ResetClusters(void);
  
  bool Clusterize(const int aTag=-1);

  inline void SetClusterizeDt(const double aDt){Cdelta_t=aDt; return;};

  inline void SetClusterizeSizeMin(const double aSizeMin){Csizemin=aSizeMin; return;};

  inline void SetClusterizeSNRthr(const double aSNRthr){cl_snr_thr=aSNRthr; return;};

  bool SetClusterTag(const int aCl, const int aNewTag);

  inline void SetClusterTag(const int aNewTag){
    for(int c=0; c<Ctree->GetEntries(); c++) cl_tag[c]=aNewTag;
  };

  bool SetClusterBranchStatus(const string aBname, const bool aStatus = true);

  inline int GetNClust(void){return Ctree->GetEntries();};
  inline int GetNclust(void){return Ctree->GetEntries();};

  inline int GetNClusters(void){return Ctree->GetEntries();};
  inline int GetNclusters(void){return Ctree->GetEntries();};

  inline double GetClusterizeDt(void){return Cdelta_t;};
  
  inline double GetClusterTime(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Ctime;
    return 0.0;
  };

  inline double GetClusterFrequency(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Cfreq;
    return 0.0;
  };

  inline double GetClusterSNR(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Csnr;
    return 0.0;
  };

  inline double GetClusterAmplitude(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Camp;
    return 0.0;
  };

  inline double GetClusterQ(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Cq;
    return 0.0;
  };

  inline double GetClusterPhase(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Cph;
    return -100.0;
  };

  inline double GetClusterTimeStart(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Ctstart;
    return 0.0;
  };

  inline double GetClusterTimeEnd(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Ctend;
    return 0.0;
  };

  inline double GetClusterDuration(const int aCl){
    return GetClusterTimeEnd(aCl)-GetClusterTimeStart(aCl);
  };

  inline double GetClusterFrequencyStart(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Cfstart;
    return 0.0;
  };

  inline double GetClusterFrequencyEnd(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return Cfend;
    return 0.0;
  };

  inline double GetClusterBandWidth(const int aCl){
    return GetClusterFrequencyEnd(aCl)-GetClusterFrequencyStart(aCl);
  };

  inline int GetClusterSize(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return (int)Csize;
    return 0;
  };

  inline int GetClusterFirstTrigger(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return (int)Cfirstentry;
    return -1;
  };

  int GetClusterIndex(const double aTime);

  inline int GetClusterTag(const int aCl){
    if(Ctree->GetEntry(aCl)>0) return cl_tag[aCl];
    return -100;
  }

  inline double GetClusterMeanDuration(void){ return Cmeandur; };
  
  inline double GetClusterMaxDuration(void){ return Cmaxdur; };
  
  void ResetTriggers(void);
  
  bool ReadTriggerFiles(const string aPattern, const string aDirectory);

 protected:

  // GENERAL
  int Verbose;             

  // TRIGGERS
  TTree *Ttree;            
  Chain *Ttree_read;       
  TTree *Ttree_make;       
  double Ttime,            
    Tfreq,                 
    Tq,                    
    Tsnr,                  
    Ttstart,               
    Ttend,                 
    Tfstart,               
    Tfend,                 
    Tamp,                  
    Tph;                   

  // CLUSTERS
  TTree *Ctree;            
  double Ctime,            
    Ctstart,               
    Ctend,                 
    Cfreq,                 
    Cfstart,               
    Cfend,                 
    Csnr,                  
    Camp,                  
    Cq,                    
    Cph;                   
  Long64_t Cfirstentry,    
    Csize;                 
  int *cl_tag;             
  double Cmeandur;         
  double Cmaxdur;          
  
  // READING-MODE ONLY
  Long64_t *Toffsets;      

  // BUILDING-MODE ONLY
  int Ttstart_ms;          
  
 private:

  // GENERAL
  bool readmode;           
  int randid;              
  string srandid;          

  // CLUSTERIZE
  double Cdelta_t;         
  int Csizemin;            
  double cl_snr_thr;       
  bool Clusterize_t(void); 
  
  ClassDef(Triggers,0)  
};

#endif

---