Luminescence (version 0.3.4)

CW2pLM: Transform a CW-OSL curve into a pLM-OSL curve

Description

Transforms a conventionally measured continuous-wave (CW) curve into a pseudo linearly modulated (pLM) curve using the equations given in Bulur (2000).

Usage

CW2pLM(values)

Arguments

values

RLum.Data.Curve or data.frame (required): RLum.Data.Curve data object. Alternatively, a data.frame of the measured curve data o

Value

• The function returns the same data type as the input data type with the transformed curve values.
• data.frame
  generic R data structure
• RLum.Data.Curve
  package RLum object

Function version

0.4 (2014-04-13 14:29:18)

Details

According to Bulur (2000) the curve data are transformed by introducing two new parameters P (stimulation period) and u (transformed time): $$P=2*max(t)$$ $$u=\sqrt{(2*t*P)}$$ The new count values are then calculated by $$ctsNEW = cts(u/P)$$ and the returned data.frame is produced by: data.frame(u,ctsNEW)

References

Bulur, E., 2000. A simple transformation for converting CW-OSL curves to LM-OSL curves. Radiation Measurements, 32, 141-145. Further Reading Bulur, E., 1996. An Alternative Technique For Optically Stimulated Luminescence (OSL) Experiment. Radiation Measurements, 26, 701-709.

See Also

CW2pHMi, CW2pLMi, CW2pPMi, fit_LMCurve, lm, RLum.Data.Curve The output of the function can be further used for LM-OSL fitting: CW2pLMi, CW2pHMi, CW2pPMi, fit_LMCurve, RLum.Data.Curve, plot_RLum

Examples

##read curve from CWOSL.SAR.Data transform curve and plot values
data(ExampleData.BINfileData, envir = environment())

##read id for the 1st OSL curve
id.OSL <- CWOSL.SAR.Data@METADATA[CWOSL.SAR.Data@METADATA[,"LTYPE"] == "OSL","ID"]

##produce x and y (time and count data for the data set)
x<-seq(CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"],
       CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"],
       by = CWOSL.SAR.Data@METADATA[id.OSL[1],"HIGH"]/CWOSL.SAR.Data@METADATA[id.OSL[1],"NPOINTS"])
y <- unlist(CWOSL.SAR.Data@DATA[id.OSL[1]])
values <- data.frame(x,y)

##transform values
values.transformed <- CW2pLM(values)

##plot
plot(values.transformed)