Learn R Programming

Luminescence (version 0.4.1)

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

Description

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.

Usage

calc_SourceDoseRate(measurement.date, calib.date, calib.dose.rate, 
    calib.error, source.type = "Sr-90", dose.rate.unit = "Gy/s")

Arguments

measurement.date
character (required): date of measurement in "YYYY-MM-DD"
calib.date
character (required): date of source calibration in "YYYY-MM-DD"
calib.dose.rate
numeric (required): dose rate at date of calibration in Gy/s or Gy/min
calib.error
numeric (required): error of dose rate at date of calibration Gy/s or Gy/min
source.type
character (with default): specify irrdiation source (Sr-90 or Co-60 or Am-214), see details for further information
dose.rate.unit
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 Se

Value

  • 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 $ parameters (list) .. $ source.type .. $ halflife .. $ dose.rate.unit

Function version

0.1 (2014-12-18 15:30:27)

Details

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 rllll{ # 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 [3] Co-60 5.274 y NNDC, Brookhaven National Laboratory }

References

NNDC, Brookhaven National Laboratory (http://www.nndc.bnl.gov/)

See Also

Second2Gray

Examples

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

##show results
get_RLum.Results(dose.rate)

##(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, dose.rate)

Run the code above in your browser using DataLab