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

Description Usage Arguments Details Value Author(s) References Examples

Description

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

Usage

1
2
3
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 curbach@gmx.de

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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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 Jan. 13, 2021, 8 a.m.