| # Gaudi Framework, version v21r7p1 |

Home |  | Generated: 15 Feb 2010

# P2D.cpp

Go to the documentation of this file.

    00001 #include "GaudiKernel/DataObject.h"
    00002 #include "GaudiKernel/ObjectFactory.h"
    00003 #include "GaudiPI.h"
    00004 #include "Generic2D.h"
    00005 #include "TProfile2D.h"
    00006 #include "TH2D.h"
    00007
    00008 namespace Gaudi {
    00009
    00016   class GAUDI_API Profile2D : public DataObject, public Generic2D<AIDA::IProfile2D,TProfile2D>   {
    00017   public:
    00019     Profile2D() {
    00020       m_classType = "IProfile2D";
    00021       m_rep = new TProfile2D();
    00022       m_rep->SetErrorOption("s");
    00023       m_rep->SetDirectory(0);
    00024       m_sumEntries = 0;
    00025     }
    00027     Profile2D(TProfile2D* rep);
    00029     virtual ~Profile2D()  {}
    00031     bool fill(double x,double y,double z,double weight)  {
    00032       m_rep->Fill(x,y,z,weight);
    00033       return true;
    00034     }
    00036     virtual const CLID& clID() const { return Gaudi::Profile2D::classID(); }
    00037     static const CLID& classID()     { return CLID_ProfileH2; }
    00038   };
    00039 }
    00040 typedef Gaudi::Profile2D P2D;
    00041 DECLARE_DATAOBJECT_FACTORY(P2D)
    00042
    00043 namespace Gaudi {
    00044   template <>
    00045   void* Generic2D<AIDA::IProfile2D,TProfile2D>::cast(const std::string& className) const  {
    00046     if (className == "AIDA::IProfile2D")
    00047       return const_cast<AIDA::IProfile2D*>((AIDA::IProfile2D*)this);
    00048     else if (className == "AIDA::IProfile")
    00049       return const_cast<AIDA::IProfile*>((AIDA::IProfile*)this);
    00050     else if (className == "AIDA::IBaseHistogram")
    00051       return const_cast<AIDA::IBaseHistogram*>((AIDA::IBaseHistogram*)this);
    00052     return 0;
    00053   }
    00054
    00055   template <>
    00056   int  Generic2D<AIDA::IProfile2D,TProfile2D>::binEntries(int idX, int idY) const  {
    00057     int rBin = m_rep->GetBin(rIndexX(idX),rIndexY(idY));
    00058     return int(m_rep->GetBinEntries(rBin)+0.5);
    00059   }
    00060
    00061   template <>
    00062   void Generic2D<AIDA::IProfile2D,TProfile2D>::adoptRepresentation(TObject* rep) {
    00063     TProfile2D* imp = dynamic_cast<TProfile2D*>(rep);
    00064     if ( imp )  {
    00065       if ( m_rep ) delete m_rep;
    00066       m_rep = imp;
    00067       m_xAxis.initialize(m_rep->GetXaxis(),true);
    00068       m_yAxis.initialize(m_rep->GetYaxis(),true);
    00069       setTitle(m_rep->GetTitle());
    00070       return;
    00071     }
    00072     throw std::runtime_error("Cannot adopt native histogram representation.");
    00073   }
    00074 }
    00075
    00076 std::pair<DataObject*,AIDA::IProfile2D*> Gaudi::createProf2D(const std::string& title, const Edges& eX, const Edges& eY, double /* zlow */ , double /* zup */) {
    00077   // Not implemented in ROOT! Can only use TProfile2D with no z-limits
    00078   Profile2D* p = new Profile2D(new TProfile2D(title.c_str(),title.c_str(),eX.size()-1,&eX.front(), eY.size()-1,&eY.front()/*,zlow,zup */));
    00079   return std::pair<DataObject*,AIDA::IProfile2D*>(p,p);
    00080 }
    00081
    00082 std::pair<DataObject*,AIDA::IProfile2D*>
    00083 Gaudi::createProf2D(const std::string& title,int binsX,double xlow,double xup,int binsY,double ylow,double yup,double zlow,double zup) {
    00084   Profile2D* p = new Profile2D(new TProfile2D(title.c_str(),title.c_str(),binsX,xlow,xup,binsY,ylow,yup,zlow,zup));
    00085   return std::pair<DataObject*,AIDA::IProfile2D*>(p,p);
    00086 }
    00087
    00088 std::pair<DataObject*,AIDA::IProfile2D*> Gaudi::createProf2D(const AIDA::IProfile2D& hist)  {
    00089   TProfile2D *h = getRepresentation<AIDA::IProfile2D,TProfile2D>(hist);
    00090   Profile2D *n = h ? new Profile2D(new TProfile2D(*h)) : 0;
    00091   return std::pair<DataObject*,AIDA::IProfile2D*>(n,n);
    00092 }
    00093
    00094 Gaudi::Profile2D::Profile2D(TProfile2D* rep)    {
    00095   m_rep = 0;
    00096   m_classType = "IProfile2D";
    00097   rep->SetDirectory(0);
    00098   adoptRepresentation(rep);
    00099   m_sumEntries = 0;
    00100 }