Coinc2.h

Go to the documentation of this file.

 
#ifndef __Coinc2__
#define __Coinc2__
 
#include "GwollumPlot.h"
#include "TriggerSelect.h"
 
using namespace std;
 
#define COINC2_NT 2
 
class Coinc2: public GwollumPlot {
 
 public:
 
  Coinc2(const unsigned int aVerbose=0);
 
  virtual ~Coinc2(void);
  bool SetTriggers(ReadTriggers *aTrigger0, ReadTriggers *aTrigger1,
                   const double aTimeOffset0=0.0, const double aTimeOffset1=0.0,
                   Segments *aValidSegments=NULL);
  
  unsigned int MakeCoinc(void);
 
  void PrintCoinc(void);
 
  void PrintCoincNot(const unsigned int aSampleIndex);
 
  inline double GetCoincTimeWindowStart(const unsigned int aSampleIndex, const unsigned int aClusterIndex){
    return TMath::Max(triggers[aSampleIndex]->GetClusterTime(aClusterIndex)+toffset[aSampleIndex]-coinc_dt/2.0,
                      triggers[aSampleIndex]->GetClusterTimeStart(aClusterIndex)+toffset[aSampleIndex]);
  };
 
  inline double GetCoincTimeWindowEnd(const unsigned int aSampleIndex, const unsigned int aClusterIndex){
    return TMath::Min(triggers[aSampleIndex]->GetClusterTime(aClusterIndex)+toffset[aSampleIndex]+coinc_dt/2.0,
                      triggers[aSampleIndex]->GetClusterTimeEnd(aClusterIndex)+toffset[aSampleIndex]);
  };
 
  inline void SetCoincDeltat(const double aDeltat){ coinc_dt=aDeltat; };
 
  inline unsigned int GetCoincN(void){ return CoC[0].size(); };
 
  inline unsigned int GetClusterIndex(const unsigned int aSampleIndex, const unsigned int aCoincIndex){ 
    return CoC[aSampleIndex%COINC2_NT][aCoincIndex];
  };
 
  inline void SetCoincTag(const unsigned int aCoincIndex, const bool aNewTag){ 
    CoTag[aCoincIndex]=aNewTag;
  };
 
  inline void SetCoincTag(const bool aNewTag){ 
    for(unsigned int c=0; c<CoC[0].size(); c++) CoTag[c]=aNewTag;
  };
 
  inline bool GetCoincTag(const unsigned int aCoincIndex){ return CoTag[aCoincIndex]; };
 
  inline unsigned int GetActiveClusterN(const unsigned int aSampleIndex){ 
    return nactive[aSampleIndex%COINC2_NT];
  };
 
  inline double GetTimeOffset(const unsigned int aSampleIndex){ 
    return toffset[aSampleIndex%COINC2_NT];
  };
 
  inline double GetCoincDeltat(void){ return coinc_dt; };
 
  inline double GetCoincLiveTime(void){ return CoSeg->GetLiveTime(); };
 
  inline Segments* GetCoincSegments(void){
    return CoSeg;
  };
 
  inline Segments* GetCoincSegmentsCopy(void){
    return new Segments(CoSeg->GetStarts(), CoSeg->GetEnds());
  };
 
  inline void SetSnrMin(const unsigned int aSampleIndex, const double aSnrMin){ 
    snrmin[aSampleIndex%COINC2_NT]=aSnrMin;
    return;
  };
 
  inline double GetSnrMin(const unsigned int aSampleIndex){ 
    return snrmin[aSampleIndex%COINC2_NT];
  };
 
  inline void SetFrequencyRange(const unsigned int aSampleIndex, const double aFrequencyMin, const double aFrequencyMax){ 
    freqmin[aSampleIndex%COINC2_NT]=aFrequencyMin;
    freqmax[aSampleIndex%COINC2_NT]=aFrequencyMax;
    return;
  };
 
  inline double GetFrequencyMin(const unsigned int aSampleIndex){ 
    return freqmin[aSampleIndex%COINC2_NT];
  };
 
  inline double GetFrequencyMax(const unsigned int aSampleIndex){ 
    return freqmax[aSampleIndex%COINC2_NT];
  };
 
  bool MakeComparators(void);
 
  void PrintComparators(const string aFileName);
 
 private:
 
  // GENERAL
  unsigned int verbosity;              
 
  // TRIGGERS
  ReadTriggers *triggers[COINC2_NT];   
  double toffset[COINC2_NT];           
  unsigned int nactive[COINC2_NT];     
  double snrmin[COINC2_NT];            
  double freqmin[COINC2_NT];           
  double freqmax[COINC2_NT];           
 
  // COINC
  double coinc_dt;                     
  vector <unsigned int> CoC[COINC2_NT];
  bool *CoTag;                         
  Segments *CoSeg;                     
 
  // COMPARATORS
  TH1D *hc_snrfrac[2];        
  TH1D *hc_freqfrac[2];       
  TGraph *gc_snrtime[4];      
  TGraph *gc_freqtime[4];     
  TGraph *gc_snrsnr;          
  TGraph *gc_freqfreq;        
 
  void SelectClusters(const unsigned int aSampleIndex);
  
  ClassDef(Coinc2,0)  
};
 
#endif