Luminescence (version 0.8.6)

calc_SourceDoseRate: Calculation of the source dose rate via the date of measurement


Calculating the dose rate of the irradiation source via the date of measurement based on: source calibration date, source dose rate, dose rate error. The function returns a data.frame that provides the input argument dose_rate for the function Second2Gray.


calc_SourceDoseRate( = Sys.Date(),,
  calib.dose.rate, calib.error, source.type = "Sr-90",
  dose.rate.unit = "Gy/s", predict = NULL)


character or Date (with default): Date of measurement in "YYYY-MM-DD". If no value is provided, the date will be set to today. The argument can be provided as vector.

character or Date (required): date of source calibration in "YYYY-MM-DD"


numeric (required): dose rate at date of calibration in Gy/s or Gy/min


numeric (required): error of dose rate at date of calibration Gy/s or Gy/min


character (with default): specify irrdiation source (Sr-90 or Co-60 or Am-214), see details for further information


character (with default): specify dose rate unit for input (Gy/min or Gy/s), the output is given in Gy/s as valid for the function Second2Gray


integer (with default): option allowing to predicit the dose rate of the source over time in days set by the provided value. Starting date is the value set with, e.g., calc_SourceDoseRate(..., predict = 100) calculates the source dose rate for the next 100 days.


Returns an S4 object of type '>RLum.Results. Slot data contains a list with the following structure:

$ dose.rate (data.frame)
.. $ dose.rate
.. $ dose.rate.error
.. $ date (corresponding measurement date)
$ parameters (list)
.. $ source.type
.. $ halflife
.. $ dose.rate.unit
$ call (the original function call)

The output should be accessed using the function get_RLum. A plot method of the output is provided via plot_RLum

Function version

0.3.1 (2018-01-25 16:46:06)

How to cite

Fuchs, M.C., Kreutzer, S. (2018). calc_SourceDoseRate(): Calculation of the source dose rate via the date of measurement. Function version 0.3.1. In: Kreutzer, S., Burow, C., Dietze, M., Fuchs, M.C., Schmidt, C., Fischer, M., Friedrich, J. (2018). Luminescence: Comprehensive Luminescence Dating Data Analysis. R package version 0.8.6.


Calculation of the source dose rate based on the time elapsed since the last calibration of the irradiation source. Decay parameters assume a Sr-90 beta source. $$dose.rate = D0 * exp(-log(2) / T.1/2 * t)$$ with: D0 <- calibration dose rate T.1/2 <- half-life of the source nuclide (here in days) t <- time since source calibration (in days) log(2) / T.1/2 equals the decay constant lambda

Information on the date of measurements may be taken from the data's original .BIN file (using e.g., BINfile <- readBIN2R() and the slot BINfile@METADATA$DATE)

Allowed source types and related values

# Source type T.1/2 Reference [1]
Sr-90 28.90 y NNDC, Brookhaven National Laboratory [2] Am-214
432.6 y NNDC, Brookhaven National Laboratory # Source type T.1/2


NNDC, Brookhaven National Laboratory

See Also

Second2Gray, get_RLum, plot_RLum



##(1) Simple function usage
##Basic calculation of the dose rate for a specific date
dose.rate <-  calc_SourceDoseRate( = "2012-01-27",
                         = "2014-12-19",
                                  calib.dose.rate = 0.0438,
                                  calib.error = 0.0019)

##show results

##(2) Usage in combination with another function (e.g., Second2Gray() )
## load example data
data(ExampleData.DeValues, envir = environment())

## use the calculated variable dose.rate as input argument
## to convert De(s) to De(Gy)
Second2Gray(ExampleData.DeValues$BT998, dose.rate)

##(3) source rate prediction and plotting
dose.rate <-  calc_SourceDoseRate( = "2012-01-27",
                         = "2014-12-19",
                                  calib.dose.rate = 0.0438,
                                  calib.error = 0.0019,
                                  predict = 1000)

##(4) export output to a LaTeX table (example using the package 'xtable')
# }

# }
# }