extract.obs: Extract One or More Gamma Combinations from am CMI Correlator

In hadron: Analysis Framework for Monte Carlo Simulation Data in Physics

Extract One or More Gamma Combinations from am CMI Correlator

Description

Extracts one or more gamma matrix combinations (observables) from a correlator stored in cmi format

Usage

extract.obs(cmicor, vec.obs = c(1), ind.vec = c(1, 2, 3, 4, 5), sym.vec,
  sign.vec, verbose = FALSE, symmetrise = TRUE)

Arguments

cmicor	an correlator object in cmi format
vec.obs	vector containing the numbers of observables to be extracted.
ind.vec	Index vector indexing the column numbers in cmicor to be used. The first must be the observable index, the second the smearing type index, the third the time, the fourth C(+t) and the fifth C(-t). Index vector indexing the column numbers in cmiloop to be used. The first must be the observable index, the second the smearing type index, the third the time, the fourth ReTL, the fifth ImTL, the sixth ReTF and the seventh ImTF.
sym.vec	a vector of bools of length equal to the number of observables indicating whether C(t) is symmetric in t, i.e. whether C(+t) and C(-t) should be added or subtracted. If not given C(+t) and C(-t) will be assumed to be symmetric.
sign.vec	a sign vector of length equal to the number of observables indicating whether the corresponding correlation function should be multiplied by +-1.
verbose	Increases verbosity of the function.
symmetrise	if set to TRUE, the correlation function will be averaged for t and Time-t, with the sign depending on the value of sym. Note that currently the correlator with t-values larger than Time/2 will be discarded.

Details

C(t) and C(-t) are averaged as indicated by sym.vec.

Value

returns a list containing

cf	for extract.obs: array containing the correlation function with dimension number of files times (nrObsnrStypes(Time/2+1)). C(t) and C(-t) are averaged according to sym.vec. for extract.loop: ReTL
icf	for extract.loop only: ImTL
scf	for extract.loop only: ReTF
sicf	for extract.loop only: ImTF
Time	The time extent of the correlation functions.
nrStypes	The number of smearing combinations.
nrObs	The number of observables.
nrSamples	for extrac.loop only: the number of samples found in the files.

Author(s)

Carsten Urbach, [e-mail]

See Also

readcmicor, readcmidatafiles,

Examples

files <- paste0(system.file(package="hadron"), "/extdata/outprcvn.dddd.00.0000")
X <- readcmifiles(files, skip=0,
                  colClasses=c("integer", "integer","integer","numeric","numeric"))
Y <- extract.obs(X)
Y

hadron documentation built on Sept. 9, 2022, 5:06 p.m.