fit.effectivemass: Fits a constant to effective mass data

fit.effectivemassR Documentation

Fits a constant to effective mass data

Description

Performs a correlated fit of a constant to data generated with bootstrap.effectivemass.

Usage

fit.effectivemass(cf, t1, t2, useCov = FALSE, replace.na = TRUE,
  boot.fit = TRUE, autoproceed = FALSE, every)

Arguments

cf

An object of class effectivemass generated by a call to bootstrap.effectivemass.

t1, t2

The fit range. If several correlators are fitted, this is automatically replicated accordingly. The fit range is adjusted such that NAs are removed from the fit. They must fulfill t1<t2. For symmetric correlators, they must both run from 0 to Time/2-1, otherwise from 0 to Time-1.

useCov

Use the correlated chisquare. This works only for not too noisy data.

replace.na

The functions inverted to determine the effective mass values might, due to fluctuations, return NA. If replace.na=TRUE, these are reaplaced in the bootstrap samples by randomly chosen values from the distribution that are not NA. Otherwise the fits in which the NA values occur will fail.

boot.fit

If set to FALSE, the effective mass fit is not bootstrapped, even though bootstrap samples are still used to estimate the variance-covariance matrix for the correlated fit. This is a useful time-saver if error information is not strictly necessary. Of course, this affects the return values related to the bootstrap, which are set to NA.

autoproceed

When the inversion of the variance-covariance matrix fails, the default behaviour is to abort the fit. Setting this to TRUE means that the fit is instead continued with a diagonal inverse of the variance-covariance matrix.

every

Fit only a part of the data points. Indices that are not multiples of every are skipped. If no value is provided, all points are taken into account.

Details

A correlated chisquare minimisation is performed on the original data as well as on all bootstrap samples generated by bootstrap.effectivemass. The inverse covariance matrix is generated as described in hep-lat/9412087 in case of too little data to relibably estimate it.

Value

An object with class effectivemassfit is returned. It contains all the data of the input object effMass with the following additional member objects:

opt.res: the object returned by the optim on the original data.

massfit.tsboot: the bootstrap values of the mass and the chisquare function.

ii: the index array of data used in the fit.

invCovMatrix: the inverse covariance matrix.

dof: the degrees of freedom of the fit.

t1,t2: the fit range.

Author(s)

Carsten Urbach, curbach@gmx.de

References

C.Michael, A.McKerrell, Phys.Rev. D51 (1995) 3745-3750, hep-lat/9412087

See Also

bootstrap.effectivemass, bootstrap.gevp, gevp2cf, invertCovMatrix

Examples


data(samplecf)
samplecf <- bootstrap.cf(cf=samplecf, boot.R=99, boot.l=2, seed=1442556)
effmass <- fit.effectivemass(bootstrap.effectivemass(cf=samplecf), t1=15, t2=23)
summary(effmass)
plot(effmass, ylim=c(0.14,0.15))


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