---

H1DVarCnv.cpp

Go to the documentation of this file.

#define ROOTHISTCNV_H1DVARCNV_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 "H1DVarCnv.h"
#include "RootDirFcn.h"

#include "TH1.h"

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


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


namespace RootHistCnv {

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

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

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

    if ( 0 != h ) {

      std::string temp = h->title();

      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 numStr( temp, 0, sep );
      int histoID = atoi(numStr.c_str()); // Histogram ID

      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
      Stat_t *errs     = new Stat_t[nBinX]; // Array of bin errors

      float *edge = new float[nBinX+1];
      for ( i = 0; i < nBinX; i++ )  {
        edge[i] = h->xAxis()->binLowerEdge(i);
      }
      edge[nBinX] = h->xAxis()->upperEdge();

      // Book the HBOOK 1D histogram with variable binning
      std::string title( temp, sep+1, temp.length() );

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

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

      TH1F *rhist = new TH1F(numStr.c_str(), title.c_str(), nBinX, edge);

      delete [] edge;

      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);
          }
        }
      }
      rhist->SetError( errs );
      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);
      }
      // 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);
      }

//        // 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);
      return StatusCode::SUCCESS;
    }
    return StatusCode::FAILURE;
  }

  StatusCode H1DVarCnv::write(int /*histoID*/, IHistogram1D* h, 
                              IOpaqueAddress*& refpAddress)   {
    MsgStream log(msgSvc(), "RootHistCnv::H1DVarCnv");

    DataObject* pObject = dynamic_cast<DataObject*>(h);
    //    IDataDirectory* par = pObject->directory()->parent();
    std::string full(pObject->location());
    //    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::DEBUG;
    if( log.isActive() ) {
//        HPRINT(histoID);
    }
//      HROUT(histoID, &icycle, "", 0);
//      HDELET(histoID);
    GenericAddress* addr = new GenericAddress( repSvcType() , objType() );
    addr->setContainerName(full);
    addr->setObjectName(h->title());
    refpAddress = addr;
    return StatusCode::SUCCESS;
  }

}    // namespace RootHistCnv

---