R/TC-scan_block.R

In bdcaf/rawTof: rawTof

#' @include TofClass.R
NULL


#' read a block with specified indices
#' @param indrange  list with two elements where the first is the start index and the second the number of scans
#' @return matrix of counts
#' @export
#' @rdname TofClass-class
setGeneric(name = "scan_block",
           def = function(.Object, indrange) {
             standardGeneric("scan_block")
                     })

#' read a block with specified indices
#' @import rhdf5
#' @return matrix of counts
#' @export
#' @rdname TofClass-class
setMethod("scan_block", "TofClass",
          function(.Object, indrange){
            dd <- .Object@.dims
            H5Sselect_hyperslab(.Object@.scanspace,
                                start = c(indrange[[1]], 1, 1, 1),
                                count = c(indrange[[2]],
                                          dd[[2]], dd[[3]], dd[[4]]) )
            h5spacemem <- H5Screate_simple(
                       c(indrange[[2]], dd[[2]], dd[[3]], dd[[4]]))
            block <- matrix(H5Dread(h5dataset = .Object@.scandata,
                                    h5spaceFile = .Object@.scanspace,
                                    h5spaceMem = h5spacemem),
                            nrow = indrange[[2]])
          H5Sclose(h5spacemem)
          block / .Object@single_ion_signal
          })