fit.effectivemass: Fits a constant to effective mass data
In hadron: Analysis Framework for Monte Carlo Simulation Data in Physics

fit.effectivemass

R 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 `NA`s 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

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.