OxEventResults.h

Go to the documentation of this file.

 
#ifndef __OxEventResults__
#define __OxEventResults__
 
#include <TGraphErrors.h>
#include "OxEventUtils.h"
#include "Oomicron.h"
 
using namespace std;
 
class OxEventResults{
 
public:
  
  OxEventResults(const string aName, OxEventParam* aParam, Omicron *aOmicron);
 
  virtual ~OxEventResults(void);
 
  void Reset(void);
 
  void SaveEvent(const unsigned int aStartTime, const unsigned int aEndTime);
 
  void Process(const unsigned int aEventsNorm);
 
  double RateToRank(const double aRate);
  
  inline void SetLoudestN(const unsigned int aN=0){
    g_loudest->Set(aN);
    for(int p=0; p<g_loudest->GetN(); p++) g_loudest->SetPoint(p, -1.0, -1.0);
  };
 
  inline unsigned int GetEventN(void){
    return (unsigned int)h1_rk->GetSumOfWeights();
  };
 
  inline double GetEventRate(const double aRank){
    if(segments->GetLiveTime()>0.0) return h1_rk_cum->Interpolate(aRank)/segments->GetLiveTime();
    return 0.0;
  };
 
  inline unsigned int GetLoudestN(void){ return (unsigned int)g_loudest->GetN(); };
 
  inline Long64_t GetLoudest(double &aRank, const unsigned int aLoudestIndex){
    if(aLoudestIndex>=GetLoudestN()) return -1;
    aRank = g_loudest->GetX()[aLoudestIndex];
    return (Long64_t)g_loudest->GetY()[aLoudestIndex];
  };
 
  inline Segments* GetSegments(void){ return segments; };
 
  inline unsigned int GetCumulativeLiveTime(void){ return livetime; };
 
  inline TH1D* GetRankPlot(void){ return h1_rk; };
 
  inline TH1D* GetCumulativeRankPlot(void){ return h1_rk_cum; };
 
  inline TH1D* GetFrequencyPlot(void){ return h1_frequency; };
 
  inline TH1D* GetXiqPlot(void){ return h1_Xiq; };
 
  inline TH1D* GetXiqMaxPlot(void){ return h1_Xiq_max; };
 
  inline TH2D* GetRankFrequencyPlot(void){ return h2_rk_frequency; };
 
  inline TH2D* GetSliceVeto2dPlot(void){ return h2_sliceveto; };
 
  inline TH1D* GetSliceVeto1dPlot(const unsigned int aQindex){ return h1_sliceveto[aQindex]; };
 
  inline TH1D* GetTimeVariancePlot(const unsigned int aQindex){ return h1_tvar[aQindex]; };
 
  TGraphErrors* GetIfarPlot(OxEventResults *aBackground);
 
private:
 
  // ANALYSIS CONTAINERS
  OxEventParam *param;           
  Segments *segments;            
  unsigned int livetime;         
  TH1D *h1_rk;                   
  TH1D *h1_rk_cum;               
  TH1D *h1_frequency;            
  TH1D *h1_Xiq;                  
  TH1D *h1_Xiq_max;              
  TH1D **h1_sliceveto;           
  TH1D **h1_tvar;                
  TH2D *h2_rk_frequency;         
  TH2D *h2_sliceveto;            
  TGraphErrors* g_loudest;       
 
  // OMICRON PARAMETERS
  unsigned int omicron_nq;       
  unsigned int* omicron_nf;      
  
  ClassDef(OxEventResults,0)  
};
 
 
#endif