dEmcClusterWrapper.C

Go to the documentation of this file.

//INCLUDECHECKER: Removed this line: #include <cstdlib>
#include "dEmcClusterWrapper.h"

ClassImp(dEmcClusterWrapper)

dEmcClusterWrapper::dEmcClusterWrapper(const char* name, const size_t& max_rows)
  : PHTable(name,max_rows)
{
  size_t rowSize = sizeof(DEMCCLUSTER_ST);
  if (max_rows > 0) { // Avoid allocating a space of zero size!
     fTableData = new DEMCCLUSTER_ST[max_rows];
     SetMaxRowCount(max_rows);
  }
  else {
     fTableData = new DEMCCLUSTER_ST[1];
     SetMaxRowCount(1);
  }

  SetRowSize(rowSize);
  SetType("dEmcCluster");
  fTableHeader->dsl_pointer  = (long)this;
  fTableHeader->data_pointer = (long)fTableData;
}

dEmcClusterWrapper::~dEmcClusterWrapper()
{
  delete [] fTableData;
}

void*
dEmcClusterWrapper::RawTableData()
{
  return static_cast<void*>(fTableData);
}

DEMCCLUSTER_ST*
dEmcClusterWrapper::TableData()
{
  return fTableData;
}

DEMCCLUSTER_ST&
dEmcClusterWrapper::operator[](const size_t& row)
{
  return fTableData[row];
}

const DEMCCLUSTER_ST&
dEmcClusterWrapper::operator[](const size_t& row) const
{
  return fTableData[row];
}

void
dEmcClusterWrapper::SetMaxRowCount(const size_t& max_rows)
{
  // Avoid reallocing a space of zero size!
  if (max_rows <= 0) {
     return;
  }

  // Ensure that the current row count is not out of range.
  if (fTableHeader->nok > max_rows) {
     fTableHeader->nok = max_rows;
  }

  // If table needs to grow, allocate a new area for it.
  if (max_rows > fTableHeader->maxlen) {
     DEMCCLUSTER_ST* newData = new DEMCCLUSTER_ST[max_rows];
     if (fTableData) {
        for (long i = 0; i < fTableHeader->nok; i++) {
           newData[i] = fTableData[i];
        }
        delete [] fTableData;
     }
     fTableData = newData;
     fTableHeader->data_pointer = (long)fTableData;
  }

  fTableHeader->maxlen = max_rows;
}
void
dEmcClusterWrapper::SetRowCount(const size_t& n)
{
  if (n > fTableHeader->maxlen) {
     fTableHeader->nok = fTableHeader->maxlen;
  }
  else if (n >= 0) {
     fTableHeader->nok = n;
  }
}

void
dEmcClusterWrapper::SetRowSize(const size_t& row_size)
{
  if (row_size > 0) {
     fTableHeader->rbytes = row_size;
  }
}

void
dEmcClusterWrapper::Streamer(TBuffer& b)
{
   // Stream an object of class dEmcClusterWrapper.
   // What should be done on output if the table is empty?

   if (b.IsReading()) {
     Version_t v = b.ReadVersion();
     PHTable::Streamer(b);         // Read the table header.
     if (RowSize() != sizeof(DEMCCLUSTER_ST)) {
       // Sanity check failed.  Need some error message here.
       return;
     }

           // Reallocate the table explicitly here; the size of the data array
           // may be inconsistent with the max. row count variable in the header
           // (since the ROOT I/O default-constructs the former, and reads
           // the header for the latter).
           size_t max_rows = MaxRowCount();
           if (max_rows <= 0) { // Avoid allocating a space of zero size!
              max_rows = 1;
           }

           delete [] fTableData;
           fTableData = new DEMCCLUSTER_ST[max_rows];
           fTableHeader->data_pointer = (long)fTableData;

     SetMaxRowCount(max_rows);
     SetType("dEmcClusterWrapper");

     for (long i=0; i<RowCount(); i++) {
        b >> fTableData[i].id;
        b >> fTableData[i].runno;
        b >> fTableData[i].evno;
        b >> fTableData[i].clusno;
        b >> fTableData[i].method;
        b >> fTableData[i].type;
        b >> fTableData[i].arm;
        b >> fTableData[i].sector;
        b.ReadStaticArray(fTableData[i].xyz);
        b.ReadStaticArray(fTableData[i].dxyz);
        b >> fTableData[i].e;
        b >> fTableData[i].ecorr;
        b >> fTableData[i].de;
        b >> fTableData[i].tof;
        b >> fTableData[i].ecent;
        b >> fTableData[i].tofcorr;
        b >> fTableData[i].dtof;
        b >> fTableData[i].qual;
        b >> fTableData[i].pid;
        b >> fTableData[i].prob_photon;
        b >> fTableData[i].prob_neuhad;
        b >> fTableData[i].chi2;
        b >> fTableData[i].nsh;
        b >> fTableData[i].chi2_sh;
        b >> fTableData[i].prob_photon_sh;
        b.ReadStaticArray(fTableData[i].e_sh);
        b.ReadStaticArray(fTableData[i].ecorr_sh);
        b.ReadStaticArray(fTableData[i].de_sh);
        b.ReadStaticArray((float *)fTableData[i].xyz_sh);
        b.ReadStaticArray((float *)fTableData[i].dxyz_sh);
     }
   }
   else {
     b.WriteVersion(IsA());
     PHTable::Streamer(b);         // Write the table header.
     for (long i=0; i<RowCount(); i++) {
        b << fTableData[i].id;
        b << fTableData[i].runno;
        b << fTableData[i].evno;
        b << fTableData[i].clusno;
        b << fTableData[i].method;
        b << fTableData[i].type;
        b << fTableData[i].arm;
        b << fTableData[i].sector;
        b.WriteArray(fTableData[i].xyz,3);
        b.WriteArray(fTableData[i].dxyz,3);
        b << fTableData[i].e;
        b << fTableData[i].ecorr;
        b << fTableData[i].de;
        b << fTableData[i].tof;
        b << fTableData[i].ecent;
        b << fTableData[i].tofcorr;
        b << fTableData[i].dtof;
        b << fTableData[i].qual;
        b << fTableData[i].pid;
        b << fTableData[i].prob_photon;
        b << fTableData[i].prob_neuhad;
        b << fTableData[i].chi2;
        b << fTableData[i].nsh;
        b << fTableData[i].chi2_sh;
        b << fTableData[i].prob_photon_sh;
        b.WriteArray(fTableData[i].e_sh,2);
        b.WriteArray(fTableData[i].ecorr_sh,2);
        b.WriteArray(fTableData[i].de_sh,2);
        b.WriteArray((float *)fTableData[i].xyz_sh,6);
        b.WriteArray((float *)fTableData[i].dxyz_sh,6);
     }
   }

}