Luminescence (version 0.3.1)

Analyse_SAR.OSLdata: Analyse SAR CW-OSL measurements.

Description

The function analyses SAR CW-OSL curve data and provides a summary of the measured data for every position. The output of the function is optimized for SAR OSL measurements on quartz.

Usage

Analyse_SAR.OSLdata(input.data, signal.integral, background.integral, 
    position, run, set, info.measurement = "unkown measurement", 
    log = "", output.plot = FALSE, output.plot.single = FALSE, 
    cex.global = 1)

Arguments

input.data

Risoe.BINfileData-class (required): input data from a Risoe BIN file, produced by the function readBIN2R.

signal.integral

vector (required): channels used for the signal integral, e.g. signal.integral=c(1:2)

background.integral

vector (required): channels used for the background integral, e.g. background.integral=c(85:100)

position

vector (optional): reader positions that want to be analysed (e.g. position=c(1:48). Empty positions are automatically omitted. If no value is given all positions are analysed by default.

run

vector (optional): range of runs used for the analysis. If no value is given the range of the runs in the sequence is deduced from the Risoe.BINfileData object.

set

vector (optional): range of sets used for the analysis. If no value is given the range of the sets in the sequence is deduced from the Risoe.BINfileData object.

info.measurement

character (with default): option to provide information about the measurement on the plot output (e.g. name of the BIN or BINX file).

log

character (with default): a character string which contains "x" if the x axis is to be logarithmic, "y" if the y axis is to be logarithmic and "xy" or "yx" if both axes are to be logarithmic. See

output.plot

logical (with default): plot output (TRUE/FALSE)

output.plot.single

logical (with default): single plot output (TRUE/FALSE) to allow for plotting the results in single plot windows. Requires output.plot = TRUE.

cex.global

numeric (with default): global scaling factor.

Value

• A plot (optional) and list is returned containing the following elements:
• LnLxTnTx
  data.frame of all calculated Lx/Tx values including signal, background counts and the dose points.
• RejectionCriteria
  data.frame with values that might by used as rejection criteria. NA is produced if no R0 dose point exists.
• SARParameters
  data.frame of additional measurement parameters obtained from the BIN file, e.g. preheat or read temperature (not valid for all types of measurements).

Details

The function works only for standard SAR protocol measurements introduced by Murray and Wintle (2000) with CW-OSL curves. For the calculation of the Lx/Tx value the function calc_OSLLxTxRatio is used. Provided rejection criteria recyling ratio: calculated for every repeated regeneration dose point. recuperation: recuperation rate calculated by comparing the Lx/Tx values of the zero regeneration point with the Ln/Tn value (the Lx/Tx ratio of the natural signal). For methodological background see Aitken and Smith (1988) IRSL/BOSL: the integrated counts (signal.integral) of a IRSL curve are compared with the integrated counts of the first regenerated dose point. It is assumed that IRSL curves got the same dose as the first regenerated dose point. Note: This is not the IR depletation ratio described by Duller (2003).

References

Aitken, M.J. & Smith, B.W., 1988. Optical dating: recuperation after bleaching. Quaternary Science Reviews, 7, pp. 387-393. Duller, G., 2003. Distinguishing quartz and feldspar in single grain luminescence measurements. Radiation Measurements, 37 (2), pp. 161-165. Murray, A.S. & Wintle, A.G., 2000. Luminescence dating of quartz using an improved single-aliquot regenerative-dose protocol. Radiation Measurements, 32, pp. 57-73.

See Also

calc_OSLLxTxRatio, Risoe.BINfileData-class, readBIN2R and for further analysis plot_GrowthCurve

Examples

##load data
data(ExampleData.BINfileData, envir = environment())

##analyse data
output <- Analyse_SAR.OSLdata(input.data = CWOSL.SAR.Data, 
                              signal.integral = c(1:5), 
                              background.integral = c(900:1000),
                              position = c(1:1),
                              output.plot = TRUE)

##combine results relevant for further analysis
output.SAR <- data.frame(Dose = output$LnLxTnTx[[1]]$Dose, 
                         LxTx = output$LnLxTnTx[[1]]$LxTx, 
                         LxTx.Error = output$LnLxTnTx[[1]]$LxTx.Error)
output.SAR