Computes the pcac mass from the PP and the AP (PA) correlators and estimates the errors using the gamma method
pcac(psfilename, apfilename, pafilename, from = 3, to = 3, fit = F,
skip = 0, plotit = F, S = 1.5)filename of the file from which to read the PP correlator. It is supposed to be in GWC code format. mandatory.
filename of the file from which to read the AP correlator. It is supposed to be in GWC code format. Either PA or AP correlator (or both) must be given. If both are given, both are used by averaging.
filename of the file from which to read the PA correlator. It is supposed to be in GWC code format. Either PA or AP correlator (or both) must be given. If both are given, both are used by averaging.
the effective mass is computed starting with t=from
the effective mass is computed ending with t=to
logical. if TRUE a fit is performed to all t-values to determine the pcac mass.
no of measurements to skip at the beginning of the file
logical. if TRUE a plot is drawn.
passed to uwerr, see documentation of uwerr.
returns a data.frame with the results. The object is also of class
massfit which can be plotted using the generic function plot.
the symmetric difference operator is used.