Data.cpp

Go to the documentation of this file.

// $Id: Data.cpp,v 1.11 2006/06/07 15:34:35 fmachefe Exp $

// include files
#include "TH1F.h"
#include "TH2F.h"
#include "Object.h"
#include "Data.h"

//-----------------------------------------------------------------------------
// Implementation file for class : Processus
//
// 2004-07-23 : Frederic Machefert
//-----------------------------------------------------------------------------

//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
Data::Data(  ){
}

//=============================================================================
// Destructor
//=============================================================================
Data::~Data() {
}


void Data::addDataStream (std::string n,  std::string t){
  m_names .push_back(n);
  m_titles.push_back(t);
  m_values.push_back(new std::vector<double>);
}

void Data::addHisto1d (TH1D *h){
  m_h1d.push_back( h );
}

void Data::addHisto2d (TH2D *h){
  m_h2d.push_back( h );
}

void Data::purge (){
  m_names .clear();
  m_titles.clear();
  std::vector< std::vector<double>* >::iterator iter;
  for (iter=m_values.begin(); iter!=m_values.end(); ++iter){
    delete (*iter);
  }
  m_values.clear();
  m_h1d.clear();
  m_h2d.clear();
  m_hist1d.clear();
  debug("Data have been purged.","Data.purge");
}

StatusCode Data::clear( int i ){
  if (vectorPtr(i)!=0){
    vectorPtr(i)->clear();
    return StatusCode::SUCCESS;
  }
  else {
    return StatusCode::FAILURE;
  }
}

StatusCode Data::dataFill(int i, double val){
  if (vectorPtr(i)!=0){
    vectorPtr(i)->push_back(val);
    return StatusCode::SUCCESS;
  }
  else {
    return StatusCode::FAILURE;
  }
}

void Data::buildHistos(){
  std::vector<TH1D*>::iterator iter1d;
  for (iter1d=m_h1d.begin(); iter1d!=m_h1d.end(); ++iter1d){
    m_hist1d.push_back(Histo1D(*iter1d));
  }
  std::vector<TH2D*>::iterator iter2d;
  for (iter2d=m_h2d.begin(); iter2d!=m_h2d.end(); ++iter2d){
    m_hist2d.push_back(Histo2D(*iter2d));
  }
}

std::string Data::name(unsigned int row){
  if (row>m_names.size()){
    warning("Try to access a data stream value not defined ("+itos(row)+")","Data::name");
    return 0;
  }
  return m_names[row];
}

std::string Data::title(unsigned int row){
  if (row>m_titles.size()){
    warning("Try to access a data stream value not defined ("+itos(row)+")","Data::title");
    return 0;
  }
  return m_titles[row];
}

std::vector<double>* Data::vectorPtr(unsigned int row){
  if (row>m_values.size()){
    warning("Try to access a data stream value not defined ("+itos(row)+")",
        "Data::vecData");
    return 0;
  }
  return m_values[row];
}

std::vector<double> Data::vector(unsigned int row){
  std::vector<double>* ptr=vectorPtr(row);
  if ((ptr)!=0){
    return *ptr;
  }
  std::vector<double> a;
  return (a);
}

int Data::rowFromName(std::string name){
  std::vector<std::string>::const_iterator iter;
  for (iter=m_names.begin(); iter!=m_names.end(); ++iter){
    if ((*iter).compare(name)==0){
      return ((int)(iter-m_names.begin()));
    }
  }
  return -1;
}

std::vector<double> Data::vector(std::string name){
  int row=rowFromName(name);
  if (row>-1){
    std::vector<double>* ptr=vectorPtr(row);
    if ((ptr)!=0){
      return *ptr;
    }
  }
  std::vector<double> a;
  return (a);
}

void Data::print(){
  char line[100];
  sprintf(line,"Number of data streams : %i", m_names.size());
  info(line,"Data::print");
  std::vector<std::string>::const_iterator iter;
  for (iter=m_names.begin(); iter!=m_names.end(); ++iter){
    unsigned int index = (unsigned int)(iter-m_names.begin());
    sprintf(line," (%3i) %10s - %40s (%8i)", index,
            (*iter).c_str(),m_titles[index].c_str(),
            vectorPtr(index)->size());
    info(line,"Data::print");
  }
  sprintf(line,"Number of 1D histograms : %i", m_h1d.size());
  info(line,"Data::print");
  std::vector<TH1D*>::const_iterator h1diter;
  for (h1diter=m_h1d.begin(); h1diter!=m_h1d.end(); ++h1diter){
    unsigned int index = (unsigned int)(h1diter-m_h1d.begin());
    sprintf(line," (%3i) %10s - %50s", index,
            (*h1diter)->GetName(),(*h1diter)->GetTitle());
    info(line,"Data::print");
  }
  sprintf(line,"Number of 2D histograms : %i", m_h2d.size());
  info(line,"Data::print");
  std::vector<TH2D*>::const_iterator h2diter;
  for (h2diter=m_h2d.begin(); h2diter!=m_h2d.end(); ++h2diter){
    unsigned int index = (unsigned int)(h2diter-m_h2d.begin());
    sprintf(line," (%3i) %10s - %50s", index,
            (*h2diter)->GetName(),(*h2diter)->GetTitle());
    info(line,"Data::print");
  }
}