dEmcGeaTowerEvalWrapper.C

Go to the documentation of this file.

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

ClassImp(dEmcGeaTowerEvalWrapper)

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

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

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

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

DEMCGEATOWEREVAL_ST*
dEmcGeaTowerEvalWrapper::TableData()
{
  return fTableData;
}

DEMCGEATOWEREVAL_ST&
dEmcGeaTowerEvalWrapper::operator[](const size_t& row)
{
  return fTableData[row];
}

const DEMCGEATOWEREVAL_ST&
dEmcGeaTowerEvalWrapper::operator[](const size_t& row) const
{
  return fTableData[row];
}

void
dEmcGeaTowerEvalWrapper::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) {
     DEMCGEATOWEREVAL_ST* newData = new DEMCGEATOWEREVAL_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
dEmcGeaTowerEvalWrapper::SetRowCount(const size_t& n)
{
  if (n > fTableHeader->maxlen) {
     fTableHeader->nok = fTableHeader->maxlen;
  }
  else if (n >= 0) {
     fTableHeader->nok = n;
  }
}

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

void
dEmcGeaTowerEvalWrapper::Streamer(TBuffer& b)
{
   // Stream an object of class dEmcGeaTowerEvalWrapper.
   // 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(DEMCGEATOWEREVAL_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 DEMCGEATOWEREVAL_ST[max_rows];
           fTableHeader->data_pointer = (long)fTableData;

     SetMaxRowCount(max_rows);
     SetType("dEmcGeaTowerEvalWrapper");

     for (long i=0; i<RowCount(); i++) {
        b >> fTableData[i].id;
        b >> fTableData[i].twrkey;
        b >> fTableData[i].input;
        b >> fTableData[i].arm;
        b >> fTableData[i].sector;
        b >> fTableData[i].indy;
        b >> fTableData[i].indz;
        b.ReadStaticArray(fTableData[i].trkno);
        b.ReadStaticArray(fTableData[i].pid);
        b.ReadStaticArray(fTableData[i].ptot);
        b.ReadStaticArray((float *)fTableData[i].vertex);
        b.ReadStaticArray(fTableData[i].ancestry);
        b.ReadStaticArray(fTableData[i].itparent);
        b.ReadStaticArray(fTableData[i].idparent);
        b.ReadStaticArray(fTableData[i].edep);
        b >> fTableData[i].mease;
        b >> fTableData[i].tof;
        b.ReadStaticArray(fTableData[i].toffirst);
        b >> fTableData[i].dist_nom;
        b >> fTableData[i].dist_act;
        b >> fTableData[i].sinthe_nom;
        b >> fTableData[i].sinthe_act;
     }
   }
   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].twrkey;
        b << fTableData[i].input;
        b << fTableData[i].arm;
        b << fTableData[i].sector;
        b << fTableData[i].indy;
        b << fTableData[i].indz;
        b.WriteArray(fTableData[i].trkno,3);
        b.WriteArray(fTableData[i].pid,3);
        b.WriteArray(fTableData[i].ptot,3);
        b.WriteArray((float *)fTableData[i].vertex,9);
        b.WriteArray(fTableData[i].ancestry,3);
        b.WriteArray(fTableData[i].itparent,3);
        b.WriteArray(fTableData[i].idparent,3);
        b.WriteArray(fTableData[i].edep,3);
        b << fTableData[i].mease;
        b << fTableData[i].tof;
        b.WriteArray(fTableData[i].toffirst,3);
        b << fTableData[i].dist_nom;
        b << fTableData[i].dist_act;
        b << fTableData[i].sinthe_nom;
        b << fTableData[i].sinthe_act;
     }
   }

}