/afs/in2p3.fr/throng/virgo/virgoDev/SoftFlorent/GWOLLUM/v2r1/src/Coinc/Coinc2.h

Go to the documentation of this file.

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

#include "GwollumPlot.h"
#include "ReadTriggers.h"
#include "MakeTriggers.h"
#include "TriggerSelect.h"

using namespace std;

class Coinc2 : public GwollumPlot {

 public:

  Coinc2(const int aVerbose=0);

  Coinc2(ReadTriggers *aTriggerSample0, ReadTriggers *aTriggerSample1, const int aVerbose=0);

  virtual ~Coinc2(void);
  bool SetTriggerSamples(ReadTriggers *aTriggerSample0, ReadTriggers *aTriggerSample1, const double aTimeOffset0=0.0, const double aTimeOffset1=0.0, Segments *aValidSegments=NULL);

  int MakeCoinc(void);

  int MakeCoincWithSmartShift(void);

  void PrintCoinc(void);

  bool WriteCoinc(const string aOutDir, const string aOutFileName);

  inline void SetCoincDt(const double aDt){ coinc_dt=aDt; return; };

  inline void SetCoincDf(const double aDf){ coinc_df=aDf; return; };

  inline void SetCoincAf(const double aAf){ coinc_af=aAf; return; };

  inline int GetNCoinc(void){ return (int)CoC[0].size(); };

  inline int GetClusterIndex(const int aSampleIndex, const int aCoincIndex){ 
    if(aCoincIndex < 0 || aCoincIndex >= GetNCoinc()) return -1;
    if(aSampleIndex < 0 || aSampleIndex >= 2) return -1;
    return CoC[aSampleIndex][aCoincIndex];
  };

  inline void SetCoincTag(const int aCoincIndex, const bool aNewTag){ 
    if(aCoincIndex>=0&&aCoincIndex<GetNCoinc()) CoTag[aCoincIndex]=aNewTag;
  };

  inline void SetCoincTag(const bool aNewTag){ 
    for(int c=0; c<(int)CoC[0].size(); c++) CoTag[c]=aNewTag;
  };

  inline bool GetCoincTag(const int aCoincIndex){ 
    if(aCoincIndex>=0&&aCoincIndex<GetNCoinc()) return CoTag[aCoincIndex];
    return false;
  };

  bool MakeMonitors(const bool aWeighted = false);

  bool WriteMonitors(const string aOutDir, const string aOutFileName);

  bool MakeComparators(const string aTriggerId0, const string aTriggerId1);
  bool PrintComparator(const string aComparatorName);

 protected:

  // GENERAL
  int verbosity;              

  // INPUT
  ReadTriggers *triggers[2];  
  double toffset[2];          

  // COINC
  vector <int> CoC[2];        
  bool *CoTag;                
  Segments *CoSeg;            
  bool smartshift;            
  TTree *CoTree;              

  // PARAMETERS
  double coinc_dt;            
  double coinc_df;            
  double coinc_af;            

  // MONITORS
  TH2D *h_n[2];               
  TH2D *h_ncum[2];            
  TH2D *h_ct[2];              
  TH2D *h_ctcum[2];           
  TH2D *h_n_coinc[2];         
  TH2D *h_ncum_coinc[2];      
  TH2D *h_n_nocoinc[2];       
  TH2D *h_ncum_nocoinc[2];    

  // COMPARATORS
  TH2D *hc_asnr;              
  TH2D *hc_snrsnr;            
  TH2D *hc_af;                
  TH2D *hc_ff;                


 private:
  
  void SelectClusters(const int aSampleIndex);  

  inline double GetSmartShiftStart(const int aSampleIndex, const int aClIndex){
    return triggers[aSampleIndex]->GetClusterTimeEnd(aClIndex)+toffset[aSampleIndex]+triggers[aSampleIndex]->GetClusterizeDt()/2.0;
  };

  inline double GetSmartShiftEnd(const int aSampleIndex, const int aClIndex){
    double te=triggers[aSampleIndex]->GetClusterTimeEnd(aClIndex)+toffset[aSampleIndex]+triggers[aSampleIndex]->GetClusterizeDt()/2.0+triggers[aSampleIndex]->GetClusterDuration(aClIndex);
    if(aClIndex==triggers[aSampleIndex]->GetNClusters()-1) return te;// last cluster
    double temax=triggers[aSampleIndex]->GetClusterTimeStart(aClIndex+1)+toffset[aSampleIndex]-triggers[aSampleIndex]->GetClusterizeDt()/2.0;
    if(te<temax) return te;
    return temax;
  };

  ClassDef(Coinc2,0)  
};

#endif

---