tests/testthat/test-mutator.R

In gamma: Dose Rate Estimation from in-Situ Gamma-Ray Spectrometry Measurements

# GammaSpectrum ================================================================
test_that("GammaSpectrum", {
  file <- system.file("extdata/LaBr.CNF", package = "gamma")
  spectrum <- read(file)

  expect_length(spectrum, 1024)
  # hash
  expect_equal(get_hash(spectrum), "e0a2c67173cf7f407db7148ae0058bbf")
  # names
  expect_equal(get_names(spectrum), "LaBr")
  set_names(spectrum) <- "X"
  expect_equal(get_names(spectrum), "X")
  expect_error(set_names(spectrum) <- c("A", "B"),
               "must be a character vector of length one, not 2")
  # live and real times
  expect_equal(get_livetime(spectrum), 3385.54)
  expect_equal(get_realtime(spectrum), 3403.67)
  # channels
  expect_length(get_channels(spectrum), 1024)
  expect_equal(range_channels(spectrum), c(1, 1024))
  # counts
  expect_length(get_counts(spectrum), 1024)
  # rates
  expect_length(get_rates(spectrum), 1024)
  # energy
  expect_length(get_energy(spectrum), 1024)
  expect_equal(range_energy(spectrum), c(-7.004032, 3124.914528))
  ## set energy calibration method
  lines <- data.frame(
    channel = c(115, 496, 870),
    energy = c(238, 1461, 2615)
  )
  set_energy_calibration(spectrum) <- lines
  expect_s4_class(spectrum, class = "GammaSpectrum")
  expect_true(has_calibration(spectrum))

  spectra <- as(list(spectrum, spectrum), "GammaSpectra")
  set_energy_calibration(spectra) <- lines
  expect_s4_class(spectra, class = "GammaSpectra")
  expect_true(all(has_calibration(spectra)))

  ## check the get method
  expect_s3_class(get_energy_calibration(spectrum), "lm")
  expect_s3_class(get_energy_calibration(spectra)[[2]], "lm")

})
# GammaSpectra =================================================================
test_that("GammaSpectra", {
  dir <- system.file("extdata/BDX_LaBr_1/calibration", package = "gamma")
  spectra <- read(dir)

  expect_length(spectra, 7)
  expect_length(names(spectra), 7)

  # subset
  expect_length(spectra[], 7)
  expect_length(spectra[NULL], 7)
  expect_length(spectra[1], 1)
  expect_length(spectra[-7], 6)
  expect_length(spectra[1:3], 3)
  expect_length(spectra[c(1, 3)], 2)
  expect_length(spectra["BRIQUE"], 1)
  expect_length(spectra[c("BRIQUE", "C347")], 2)
  expect_type(spectra[1:3, "energy"], "list")

  expect_s4_class(spectra[[1]], "GammaSpectrum")
  expect_s4_class(spectra[["BRIQUE"]], "GammaSpectrum")

  # lengths
  expect_true(all(lengths(spectra) == 1024))
  # hash
  expect_true(all(nchar(get_hash(spectra)) == 32))
  # names
  expect_true(all(nchar(get_names(spectra)) > 0))
  set_names(spectra) <- LETTERS[1:7]
  expect_true(all(get_names(spectra) == LETTERS[1:7]))
  expect_true(all(get_names(spectra) == names(spectra)))
  # live and real times
  expect_length(get_livetime(spectra), 7)
  expect_length(get_realtime(spectra), 7)
  # channels
  expect_true(all(lengths(get_channels(spectra)) == 1024))
  expect_equal(dim(range_channels(spectra)), c(7, 2))
  # counts
  expect_true(all(lengths(get_counts(spectra)) == 1024))
  # rate
  expect_true(all(lengths(get_rates(spectra)) == 1024))
  # energy
  expect_true(all(lengths(get_energy(spectra)) == 1024))
  expect_equal(dim(range_energy(spectra)), c(7, 2))

  expect_equal(dim(summarise(spectra)), c(7, 7))
})