gaussConvol: Gauss convolution

Description Usage Arguments Value See Also Examples

View source: R/supl.R

Description

Function to calculate the convolution of a given function with the Gaussian function

Usage

1
2
		   
gaussConvol(SQ, Q, Qdamp=0.0457, err=1e-6)

Arguments

SQ

numeric vector containing function to be convoluted.

Q

numeric vector containing grid points corresponding to SQ.

Qdamp

numeric indicating standard deviation parameter for the Gaussian function.

err

numeric, relative accuracy requested.

Value

numeric vector of function values.

See Also

calcTotalScatt

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## simulate a particle 
Cu <- createAtom("Cu")
part <- simPart(atoms=list(Cu), atomsShell=list(Cu), r=20, 
    rcore=14, latticep=4.08, latticepShell=3.89)

## calculate total scattering function
gQ <- calcTotalScatt(part, type="neutron",  minQ=0.771, maxQ=18,
    dQ=0.01, scatterLength=c(4.87, 7.97), sigma=c(.01, .01))
plot(gQ$Q, gQ$gQ, type="l")
						  
## simulate instrumental resolution effect
Q <- gQ$Q
gQ <- gaussConvol(SQ=gQ$gQ, Q=gQ$Q, Qdamp=0.061, err=1e-5)
lines(Q, gQ, col=2)		   

nanop documentation built on May 30, 2017, 7:22 a.m.

Search within the nanop package
Search all R packages, documentation and source code