cdhnew: finite size corrections a la Colangelo, Duerr, Haefeli, but...

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

finite size corrections a la Colangelo, Duerr, Haefeli, but re-expanded as series in the quark mass

Description

finite size corrections a la Colangelo, Duerr, Haefeli, but re-expanded as series in the quark mass

Usage

cdhnew(parm = rep(0, times = 6), rev = -1, aLamb1 = 0.055,
  aLamb2 = 0.58, aLamb3, aLamb4, ampiV, afpiV, aF0, a2B0mu, L,
  printit = FALSE, use.cimpl = TRUE)

Arguments

parm	m parameters
rev	rev=-1 corrects from L to L = infinity, rev=+1 the other way around
aLamb1	The four low energy Lambda1-Lambda4constants in lattice units.
aLamb2	see aLamb1.
aLamb3	see aLamb1.
aLamb4	see aLamb1.
ampiV	pseudo scalar mass values to be corrected
afpiV	pseudo scalar decay constant values to be corrected
aF0	af0 in lattice units
a2B0mu	2 B0 mu in lattice units, where mu is the quark mass and B0 a low energy constant
L	the lattice spatial extent
printit	if set to TRUE the corrections are printed
use.cimpl	use the four times faster direct c Implementation of the correction routine

Details

see reference for details. We use the simplyfied formulae for the S quantities, see eq. (59) in first reference.

Value

a list with the corrected values for mpi and fpi

Author(s)

Carsten Urbach [e-mail]

References

Gilberto Colangelo, Stephan Durr, Christoph Haefeli, Nucl.Phys.B721:136-174,2005. hep-lat/0503014

and

R. Frezzotti, V. Lubicz, S. Simula, arXiv:0812.4042 hep-lat

Examples

mu <- c(0.004, 0.006, 0.008, 0.010, 0.004)
L <- c(24, 24, 24, 24, 32)
mps <- c(0.14448, 0.17261, 0.19858, 0.22276, 0.14320)
fps <- c(0.06577, 0.07169, 0.07623, 0.07924, 0.06730)
aLamb1 <- 0.05
aLamb2 <- 0.5
aLamb3 <- 0.38
aLamb4 <- 0.66
aF0    <- 0.051
a2B    <- 5.64
cdhres <- cdhnew(rev=+1, aLamb1=aLamb1, aLamb2=aLamb2, aLamb3=aLamb3,
                 aLamb4=aLamb4, ampiV=mps, afpiV=fps, aF0=aF0,
                 a2B0mu=a2B*mu, L=L, printit=TRUE)
cdhres$mpiFV
cdhres$fpiFV

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