GaudiPI.h

Go to the documentation of this file.

#ifndef GAUDISVC_GAUDIPI_H
#define GAUDISVC_GAUDIPI_H

#include "AIDA_visibility_hack.h"

#include <vector>
#include <utility>
#include "AIDA/IHistogram1D.h"
#include "AIDA/IHistogram2D.h"
#include "AIDA/IHistogram3D.h"
#include "AIDA/IProfile1D.h"
#include "AIDA/IProfile2D.h"
#include "GaudiKernel/HistogramBase.h"

class DataObject;
class TH2D;

namespace Gaudi {
  typedef std::vector<double> Edges;

  template <class Q, class T> T* getRepresentation(const Q& hist)  {
    T* result = 0;
    const HistogramBase *p = dynamic_cast<const HistogramBase*>(&hist);
    if ( p )  {
      result = dynamic_cast<T*>(p->representation());
    }
    return result;
  }
  std::pair<DataObject*,AIDA::IHistogram1D*> createH1D(const AIDA::IHistogram1D& hist);
  std::pair<DataObject*,AIDA::IHistogram1D*> createH1D(const std::string& title, int nBins=10, double lowerEdge=0., double upperEdge=1.);
  std::pair<DataObject*,AIDA::IHistogram1D*> createH1D(const std::string& title, const Edges& e);

  std::pair<DataObject*,AIDA::IHistogram2D*> createH2D(const AIDA::IHistogram2D& hist);
  std::pair<DataObject*,AIDA::IHistogram2D*> createH2D(TH2D* rep);
  std::pair<DataObject*,AIDA::IHistogram2D*> createH2D(const std::string& title, int binsX, double iminX, double imaxX, int binsY, double iminY, double imaxY);
  std::pair<DataObject*,AIDA::IHistogram2D*> createH2D(const std::string& title, const Edges& eX, const Edges& eY);
  std::pair<DataObject*,AIDA::IHistogram1D*> slice1DX(const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);
  std::pair<DataObject*,AIDA::IProfile1D*>   profile1DX(const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);
  std::pair<DataObject*,AIDA::IHistogram1D*> project1DX(const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);
  std::pair<DataObject*,AIDA::IHistogram1D*> slice1DY  (const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);
  std::pair<DataObject*,AIDA::IProfile1D*>   profile1DY(const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);
  std::pair<DataObject*,AIDA::IHistogram1D*> project1DY(const std::string& name, const AIDA::IHistogram2D& h, int firstbin, int lastbin);

  std::pair<DataObject*,AIDA::IHistogram3D*> createH3D(const AIDA::IHistogram3D& hist);
  std::pair<DataObject*,AIDA::IHistogram3D*> createH3D(const std::string& title, int nBinsX, double xlow, double xup, int nBinsY, double ylow, double yup, int nBinsZ, double zlow, double zup);
  std::pair<DataObject*,AIDA::IHistogram3D*> createH3D(const std::string& title, const Edges  & eX,const Edges& eY, const Edges& eZ);

  std::pair<DataObject*,AIDA::IProfile1D*> createProf1D(const AIDA::IProfile1D& hist);
  std::pair<DataObject*,AIDA::IProfile1D*> createProf1D
  (const std::string& title, int nBins, double xlow, double xup, double ylow, double yup, const std::string& opt="" );
  std::pair<DataObject*,AIDA::IProfile1D*> createProf1D
  (const std::string& title, const Edges& e, double ylow, double yup, const std::string& opt="");

  std::pair<DataObject*,AIDA::IProfile2D*> createProf2D(const AIDA::IProfile2D& hist);
  std::pair<DataObject*,AIDA::IProfile2D*> createProf2D(const std::string& title, int binsX, double iminX, double imaxX, int binsY, double iminY, double imaxY, double lowerValue, double upperValue);
  std::pair<DataObject*,AIDA::IProfile2D*> createProf2D(const std::string& title, const Edges& eX, const Edges& eY, double lowerValue, double upperValue);
}
#endif // GAUDIPI_GAUDIPI_H