dune::HDF5Utils Namespace Reference

Classes
struct	HDFFileInfo
struct	HeaderInfo

Typedefs
typedef std::vector< raw::RawDigit >	RawDigits
typedef std::vector< raw::RDTimeStamp >	RDTimeStamps
typedef std::unique_ptr< HDFFileInfo >	HDFFileInfoPtr
typedef std::vector< Fragment >	Fragments

Functions
HDFFileInfoPtr	openFile (const std::string &fileName)
void	closeFile (HDFFileInfoPtr hdfFileInfoPtr)
std::deque< std::string >	getTopLevelGroupNames (HDFFileInfoPtr &hdfFileInfoPtr)
std::deque< std::string >	getMidLevelGroupNames (hid_t grp)
uint64_t	formatTrigTimeStamp (uint64_t trigTimeStamp)
bool	attrExists (hid_t object, const std::string &attrname)
hid_t	getGroupFromPath (hid_t fd, const std::string &path)
void	getHeaderInfo (hid_t the_group, const std::string &det_type, HeaderInfo &info)

Typedef Documentation

typedef std::vector<Fragment> dune::HDF5Utils::Fragments

Definition at line 57 of file HDF5Utils.h.

typedef std::unique_ptr<HDFFileInfo> dune::HDF5Utils::HDFFileInfoPtr

Definition at line 46 of file HDF5Utils.h.

typedef std::vector<raw::RawDigit> dune::HDF5Utils::RawDigits

Definition at line 23 of file HDF5Utils.h.

typedef std::vector<raw::RDTimeStamp> dune::HDF5Utils::RDTimeStamps

Definition at line 24 of file HDF5Utils.h.

Function Documentation

bool dune::HDF5Utils::attrExists	(	hid_t	object,
		const std::string &	attrname
	)

Definition at line 71 of file HDF5Utils.cc.

                                                           {
  // Save old error handler 
  H5E_auto_t old_func;
  void *old_client_data;
  H5Eget_auto(H5E_DEFAULT,&old_func, &old_client_data);
  
  // Turn off error handling */
  H5Eset_auto(H5E_DEFAULT,NULL, NULL);
  
  // Probe. On failure, retval is supposed to be negative
  
  hid_t retval = H5Aopen_name(object, attrname.data());
  
  // Restore previous error handler 
  H5Eset_auto(H5E_DEFAULT,old_func, old_client_data);
  
  bool result = (retval >= 0);
  return result;
}

void dune::HDF5Utils::closeFile

(

HDFFileInfoPtr

hdfFileInfoPtr

)

Definition at line 31 of file HDF5Utils.cc.

                                              {
  H5Fclose(hdfFileInfoPtr->filePtr);
  hdfFileInfoPtr->filePtr = 0;
}

uint64_t dune::HDF5Utils::formatTrigTimeStamp

(

uint64_t

trigTimeStamp

)

Definition at line 58 of file HDF5Utils.cc.

  {
    int clock_period = 20; //20 ns with 50MHz DAQClock //FIXME don't hardcode
    uint64_t trigtime = trigTimeStamp * clock_period; //time in ns 
    uint64_t million = 1e9;
    uint64_t trigtime_sec =  trigtime/million;
    //std::cout << "trigtime_sec: " << trigtime_sec << std::endl;
    uint64_t trigtime_nanosec = trigtime % million;
    uint64_t toReturn = (trigtime_sec << 32) | trigtime_nanosec;
    return toReturn;
  }

hid_t dune::HDF5Utils::getGroupFromPath	(	hid_t	fd,
		const std::string &	path
	)

Definition at line 91 of file HDF5Utils.cc.

                                                        {
  hid_t grp = H5Gopen(fd, path.data(), H5P_DEFAULT);
  return grp;
}

void dune::HDF5Utils::getHeaderInfo	(	hid_t	the_group,
		const std::string &	det_type,
		HeaderInfo &	info
	)

Definition at line 96 of file HDF5Utils.cc.

                                {
  hid_t datasetid = H5Dopen(the_group, det_type.data(), H5P_DEFAULT);
  hsize_t ds_size = H5Dget_storage_size(datasetid);
  //std::cout << "      Data Set Size (bytes): " << ds_size << std::endl;
  // todo -- check for zero size
  if (ds_size < 64) {
    //std::cout << "TriggerRecordHeader datset too small" << std::endl; 
  }
  
  size_t narray = ds_size / sizeof(char);
  size_t rdr = ds_size % sizeof(char);
  if (rdr > 0 || narray == 0) narray++;
  char *ds_data = new char[narray];
  H5Dread(datasetid, H5T_STD_I8LE, H5S_ALL, H5S_ALL,
                       H5P_DEFAULT, ds_data);
  //std::cout << std::hex << "      Retrieved data: ecode: " << ecode <<
  //"  first byte: " << firstbyte << " last byte: " <<
  //lastbyte << std::dec << std::endl;
  H5Dclose(datasetid);
  
  //int magic_word = 0;
  memcpy(&info.magicWord, &ds_data[0],4);
  //std::cout << "   Magic word: 0x" << std::hex << info.magicWord << std::dec <<
  //std::endl;
  
  //int version = 0;
  memcpy(&info.version, &ds_data[4],4);
  //std::cout << "   Version: " << std::dec << info.version << std::dec <<
  //std::endl;
  
  //uint64_t trignum=0;
  memcpy(&info.trigNum, &ds_data[8],8);
  //std::cout << "   Trig Num: " << std::dec << info.trigNum << std::dec <<
  //std::endl;
  
  //uint64_t trig_timestamp=0;
  memcpy(&info.trigTimestamp, &ds_data[16],8);
  //std::cout << "   Trig Timestamp: " << std::dec << info.trigTimestamp <<
  //std::dec << std::endl;
  
  //uint64_t nreq=0;
  memcpy(&info.nReq, &ds_data[24],8);
  //std::cout << "   No. of requested components:   " << std::dec << info.nReq <<
  //std::dec << std::endl;
  
  //int runno=0;
  memcpy(&info.runNum, &ds_data[32], 4);
  //std::cout << "   Run Number: " << std::dec << info.runNum << std::endl;
  //run_id = info.runNum;
  
  //int errbits=0;
  memcpy(&info.errBits, &ds_data[36], 4);
  //std::cout << "   Error bits: " << std::dec << info.errBits << std::endl;
  
  //short triggertype=0;
  memcpy(&info.triggerType, &ds_data[40], 2);
  //std::cout << "   Trigger type: " << std::dec << info.triggerType << std::endl;
  
  //delete[] ds_data;  // free up memory
} 

std::deque< std::string > dune::HDF5Utils::getMidLevelGroupNames

(

hid_t

grp

)

Definition at line 43 of file HDF5Utils.cc.

                                                     {
  std::deque<std::string> theList;
  hsize_t nobj = 0;
  H5Gget_num_objs(grp, &nobj);
  for (hsize_t idx = 0; idx < nobj; ++idx) {
    hsize_t len = H5Gget_objname_by_idx(grp, idx, NULL, 0 );
    char *memb_name = new char(len+1);
    H5Gget_objname_by_idx(grp, idx, memb_name, len+1 );
    theList.emplace_back(memb_name);
    delete memb_name;
  }
  return theList;
}

std::deque< std::string > dune::HDF5Utils::getTopLevelGroupNames

(

HDFFileInfoPtr &

hdfFileInfoPtr

)

Definition at line 36 of file HDF5Utils.cc.

                                                                          {
  hid_t grp = H5Gopen(hdfFileInfoPtr->filePtr, "/", H5P_DEFAULT);
  std::deque<std::string> theList = getMidLevelGroupNames(grp);
  H5Gclose(grp);
  return theList;
}

HDFFileInfoPtr dune::HDF5Utils::openFile

(

const std::string &

fileName

)

Definition at line 12 of file HDF5Utils.cc.

                                                   {
  HDFFileInfoPtr hdfFileInfoPtr(new HDFFileInfo());
  hdfFileInfoPtr->filePtr = H5Fopen(fileName.data(), H5F_ACC_RDONLY, H5P_DEFAULT);
  hdfFileInfoPtr->bytesWritten = 0;
  hdfFileInfoPtr->fileName = fileName;
  hdfFileInfoPtr->runNumber = 0;   // don't know it yet.  Fill it when we get fragments.
  hdfFileInfoPtr->dataFormatVersion = 0;

  hid_t grp = H5Gopen(hdfFileInfoPtr->filePtr,"/", H5P_DEFAULT);
  if (attrExists(grp, "data_format_version"))
    {
      hid_t ga = H5Aopen_name(grp, "data_format_version");
      H5Aread(ga, H5Aget_type(ga), &hdfFileInfoPtr->dataFormatVersion);
      H5Aclose(ga);
    }
  H5Gclose(grp);
  return hdfFileInfoPtr;
}