/Users/jmonk/Physics/ForIA/ForIA/Analysis.hh

Go to the documentation of this file.

#ifndef FORIA_ANALYSIS_HH
#define FORIA_ANALYSIS_HH

#include "ForIA/Units.hh"

#include "boost/smart_ptr.hpp"

#include <string>
#include <vector>

namespace ForIA {
  
  class Event;
  class IHistogrammer;
  using std::string;
  using std::vector;
  

  class Analysis {
    
  public:
    Analysis(IHistogrammer *histogrammer, const string &name);
    
    virtual ~Analysis(){};
    
    const string &name()const;
    
    virtual bool initialise()=0;
    void executeClear(const Event &evt);
    
    virtual bool finalise()=0;
    
    bool storesEvent()const;
    
  protected:
    
    virtual bool execute(const Event &evt)=0;
    
    void storeEvent(bool doStore=true);
    
    bool bookHistogram1D(const string &path, const string &name, const string &title, 
                         int nBins, double xMin, double xMax, bool divideBinWidths=true);
    
    bool bookHistogram1D(const string &name, const string &title, 
                         int nBins, double xMin, double xMax, bool divideBinWidths=true);
    
    bool bookHistogram1D(const string &path, const string &name, const string &title, 
                         const vector<double> &bins, bool divideBinWidths=true);

    bool bookHistogram1D(const string &name, const string &title, 
                         const vector<double> &bins, bool divideBinWidths=true);
    
    
    void bookProfile1D(const string &path, const string &name, const string &title,
                       int nBins, double xMin, double xMax);
    
    void bookProfile1D(const string &name, const string &title,
                       int nBins, double xMin, double xMax);
    
    bool bookHistogram2D(const string &path, const string &name, const string &title, 
                         int nBinsX, double xMin, double xMax,
                         int nBinsY, double yMin, double yMax,
                         bool divideBinWidths=true);
    
    bool bookHistogram2D(const string &name, const string &title, 
                         int nBinsX, double xMin, double xMax,
                         int nBinsY, double yMin, double yMax, 
                         bool divideBinWidths=true);
    
    bool bookHistogram2D(const string &path, const string &name, const string &title,
                         const vector<double> &xbins, const vector<double> &ybins,
                         bool divideBinWidths=true);
    
    bool bookHistogram2D(const string &name, const string &title,
                         const vector<double> &xbins, const vector<double> &ybins,
                         bool divideBinWidths=true);
    
    
    void bookHistogram3D(const string &path, const string &name, const string &title,
                         int nBinsX, double xMin, double xMax,
                         int nBinsY, double yMin, double yMax,
                         int nBinsZ, double zMin, double zMax,
                         bool divideBinWidths=true);
    
    void bookHistogram3D(const string &name, const string &title,
                         int nBinsX, double xMin, double xMax,
                         int nBinsY, double yMin, double yMax,
                         int nBinsZ, double zMin, double zMax,
                         bool divideBinWidths=true);
    
    
    bool fillHistogram1D(const string &path, const string &name, double value, double weight);
    
    bool fillHistogram1D(const string &name, double value, double weight);
    
    void fillProfile1D(const string &path, const string &name, double xVal, double yVal, double weight);
    
    void fillProfile1D(const string &name, double xVal, double yVal, double weight);
    
    bool fillHistogram2D(const string &path, const string &name, double xVal, double yVal, double weight);

    bool fillHistogram2D(const string &name, double xVal, double yVal, double weight);
    
    void fillHistogram3D(const string &path, const string &name, double xVal, double yVal, double zVal, double weight);
    
    void fillHistogram3D(const string &name, double xVal, double yVal, double zVal, double weight);
    
    void normaliseHistogram1D(const string &path, const string &name, double intgl);

    void normaliseHistogram1D(const string &name, double intgl);
    
    void scaleHistogram1D(const string &path, const string &name, double scale);
    
    void scaleHistogram1D(const string &name, double scale);
    
    IHistogrammer *histograms();
    
    string analysisPath(const std::string path){
      return "/" + m_name + path;
    }
    
  private:
    
    Analysis(){}
    
    IHistogrammer *m_histogrammer;
    string m_name;
    bool m_storeEvent;
    
  };
  
  typedef boost::shared_ptr<Analysis> AnalysisPtr;
  typedef boost::shared_ptr<Analysis const> AnalysisConstPtr;
  
}

#endif