---

H1DCnv.cpp

Go to the documentation of this file.

#define ROOTHISTCNV_H1DCNV_CPP


// Include files
#include <iostream>
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/CnvFactory.h"
#include "GaudiKernel/GenericAddress.h"
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/IHistogram1D.h"
#include "GaudiKernel/IAxis.h"

#include "H1DCnv.h"
#include "RootDirFcn.h"

#include "TH1.h"
#include "TFile.h"
#include <list>

//------------------------------------------------------------------------------
//
// Implementation of class :  RootHistCnv::H1DCnv
//
// Author :                   Charles Leggett
//
//------------------------------------------------------------------------------


// Instantiation of a static factory class used by clients to create
// instances of this service
static CnvFactory<RootHistCnv::H1DCnv> s_factory;
const ICnvFactory& RootHistH1DCnvFactory = s_factory;


namespace RootHistCnv {

  StatusCode H1DCnv::createRep(DataObject* pObject, IOpaqueAddress*& refpAddress)  {

    IHistogram1D* h = dynamic_cast<IHistogram1D*>(pObject);
    refpAddress = 0;

    MsgStream log(msgSvc(), "RootHistCnv::H1DCnv");

    if ( 0 != h ) {

      std::string temp = h->title();
      //      log << MSG::DEBUG << "temp: " << temp << endreq;

      int sep = temp.find('|');
      if ( 0 >= sep ) {
        // No histogram ID, no histogram name
        log << MSG::ERROR << "Cannot separate histogram title and identifier" << endreq;
        log << MSG::ERROR << "while converting the histogram: " << temp << endreq;
        return StatusCode::FAILURE;
      }

      std::string histoID( temp, 0, sep );
      std::string title( temp, sep+1, temp.length() );

      int    nBinX    = h->xAxis()->bins();        // Number of bins
      float  xLow     = h->xAxis()->lowerEdge();   // Histogram lower edge
      float  xHigh    = h->xAxis()->upperEdge();   // Histogram upper edge
      float  x_uflow  = xLow  - 1.0;      // Underflow coordinate
      float  x_oflow  = xHigh + 1.0;      // Overflow coordinate

      int    entries  = 0;                // Number of entries in a bin
      float  height   = 0.;               // Height if a bin
      float  centre   = 0.;               // Centre of a bin
      float  weight   = 0.;               // height / entries
      int    i        = 0;                // Bin index
      int    j        = 0;                // Fill index
//        float *errs     = new float[nBinX]; // Array of bin errors
      Stat_t *errs    = new Stat_t[nBinX];


      std::string full = pObject->fullpath();
      std::string dir = pObject->location();
      RootTrimLeadingDir(dir,"/stat");
      RootMkdir(dir);
      RootCd(dir);

      log << MSG::INFO << "creating TH1F id: " << histoID 
          << "  title: " << title << " in " << dir << endreq;

      //      cdToPath(full);
      TH1F *rhist = new TH1F(histoID.c_str(),title.c_str(),nBinX,xLow,xHigh);

      for ( i = 0; i < nBinX; i++ )  {
        errs[i] = h->binError(i);
        height = h->binHeight(i);
        if( 0 != height ) {
          centre  = h->xAxis()->binCentre(i);
          entries = h->binEntries(i);
          weight  = height / entries;
          for( j = 0; j < entries; j++ ) {
            rhist->Fill(centre,weight);
          }
        }
      }
//        log << MSG::DEBUG << "filled hist" << endreq;

      rhist->SetError( errs );

//        log << MSG::DEBUG << "set errors" << endreq;

      delete [] errs;


      // Underflow bins
      entries          = h->binEntries( IHistogram::UNDERFLOW_BIN );
      height           = h->binHeight( IHistogram::UNDERFLOW_BIN );
      weight  = height / entries;
      for( j = 0; j < entries; j++ ) {
        rhist->Fill(x_uflow,weight);
      }
//        log << MSG::DEBUG << "filled underflow" << endreq;
      // Overflow bins
      entries          = h->binEntries( IHistogram::OVERFLOW_BIN );
      height           = h->binHeight( IHistogram::OVERFLOW_BIN );
      weight  = height / entries;
      for( j = 0; j < entries; j++ ) {
        rhist->Fill(x_oflow,weight);
      }
//        log << MSG::DEBUG << "filled overflow" << endreq;

      // Set the correct MEAN and RMS vales in the HBOOK histogram
//        float mean          = (float) h->mean();
//        float rms           = (float) h->rms();
//        float sumOfHeights  = (float) h->sumBinHeights();
//        float equivEntries  = (float) h->equivalentBinEntries();
      //      SET_MEAN_RMS( histoID, mean, rms, sumOfHeights, equivEntries );

      write(histoID, h, refpAddress);
//        log << MSG::DEBUG << "called write()" << endreq;
      return StatusCode::SUCCESS;
    }
    log << MSG::ERROR << "no histograms found" << endreq;
    return StatusCode::FAILURE;
  }

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

  StatusCode H1DCnv::write(std::string /*histoID*/, IHistogram1D* h, IOpaqueAddress*& refpAddress)   {

    MsgStream log(msgSvc(), "RootHistCnv::H1DCnv");

    DataObject* pObject = dynamic_cast<DataObject*>(h);
//      IDataDirectory* par = pObject->directory()->parent();
    std::string full(pObject->location());
//      std::string lpath(pObject->localPath());
//      std::string fpath(pObject->fullpath()); 

//      DataObject* pObjPar = pObject->parent();
//      std::string pPath = pObjPar->fullpath();

//      std::list<std::string> dirlst;
//      std::list<std::string>::const_iterator diritr;

//      while (pObjPar != NULL) {
//        std::string dir = pObjPar->localPath();
//        std::string pth = pObjPar->location();

//        dirlst.push_front(dir);

//        pObjPar = pObjPar->parent();
//      }
    
//      for (diritr=dirlst.begin(); diritr!=dirlst.end(); ++diritr) {
//        log << MSG::DEBUG << "pth: " << *diritr << endreq;
//      }
    

//      int icycle = 0;
    full.insert(0,"/");

//      // I know this is a big hack...but there is no other way to find out
//      // If the RZ directory exists or not!
//      if ( par != 0 && par->address() == 0 )    {
//        HMDIR(full,"S");
//      }
//      else  {
//        HCDIR(full);
//      }

//      log << MSG::INFO << "In R:H1DCnv::write" << endreq;
//      log << MSG::DEBUG << "location(): " << full << endreq;
//      log << MSG::DEBUG << "localPath(): " << lpath << endreq;
//      log << MSG::DEBUG << "fullpath(): " << fpath << endreq;
//      log << MSG::DEBUG << "parentPath: " << pPath << endreq;

//      log << MSG::DEBUG;
//      if( log.isActive() ) {
//        //      HPRINT(histoID);
//      }

    GenericAddress* addr = new GenericAddress( repSvcType() , objType() );
    addr->setContainerName(full);
    addr->setObjectName(h->title());
    refpAddress = addr;
    return StatusCode::SUCCESS;
  }

}    // namespace RootHistCnv

---