R/read.bundletree.R

In moritzlindner/PatchR: An environment for working with ephys (PatchClamp) recordings

#' read.bundletree
#'
#' function inherited from ephys2 (/github/tdanker/ephys2/). Undocumented.
#' @noRd
read.bundletree<-function (myfile, bundlename = ".pul", con = NA,encoding=getOption("encoding"))
{
  name<-""
  rm(name)
  finally_close_con = is.na(con)
  if (is.na(con)) {
    con <- file(myfile, "rb", encoding=encoding)
  }
  seek(con, 0)
  signature <- readChar(con, 8)
  if (signature == "DAT2") {
    seek(con, where = 48)
    nitems <- readBin(con, "int", size = 1)
    seek(con, where = 64)
    bundleitems <- do.call(rbind, (lapply(0:nitems, function(item) {
      start <- readBin(con, "int", size = 4)
      end <- readBin(con, "int", size = 4)
      name <- readChar(con, 8)
      data.frame(start = start, end = end, name = name,
                 stringsAsFactors = F)
    })))
    start <- subset(bundleitems, name == bundlename)$start
    seek(con, where = start)
  }
  else {
    start <- 0
  }
  seek(con, where = start)
  magic <- readChar(con, nchars = 4)
  stopifnot(magic == "eerT")
  nLevels <- readBin(con, "int", size = 4)
  lvl_sizes <- lapply(1:nLevels, function(i) {
    readBin(con, "int", size = 4)
  })
  tree <- pm_load_nodes(con, nLevels, lvl_sizes, 1)
  if (finally_close_con)
    close(con)
  tree
}