Finds Peak Position from Gaussian — xrd.find.Peak • rigakuXRD

Fits a Gaussian form to a peak

xrd.find.Peak(
  TwoTheta,
  Intensity,
  peakPos = NA,
  Try.Sigma = c(0.1, 0.4, 0.2, 0.15),
  thetaDelta = 5,
  verbose = FALSE
)

Arguments

TwoTheta: angle
Intensity: intensity signal
peakPos: approximate angle of peak position
Try.Sigma: vector with peak widths used to start fitting
thetaDelta: width of angle to search peak
verbose: logical, if TRUE provides extra information

Value

fit peak position

Examples

filename = xrd.getSampleFiles(fileExt='asc')[1]
d = xrd.read.ASC(filename)
peak.pos = xrd.find.Peak(d$theta, d$I, 38.2)
plot(d$theta,d$I,log='y',col='red')
abline(v=peak.pos,col='blue')