eta_c with wide priors

We started with the following definition of parameters:

E[0] = 1.640 +- 1.640
log(E[n]-E[n-1]) = -0.868 +- 4.000
log(Z[n,local]) = 0.300 +- 4.000
log(Z[n,1S]) = 1.700 +- 4.000
log(Z[n,2S]) = 0.500 +- 4.000
Nstates = 3
Tmin,Tmax = 1,32
All SVD cuts disabled
Preconditioning of data covariance matrix enabled (rescale to have unit diagonal prior to inversion)
Blocksize = 1
Datasets are:
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_d_d_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_1S_1S_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_2S_2S_0.119_0.119_p000.bz2
Functional form of fit is Sum_{elevel} Z^2 * exp(-E[elevel]*Tmid) * 2*cosh(E[elevel]*(Tmid-t)) where
- Tmid is the middle-timeslice, for our case, Tmid=32.
- Z[n] = exp(alpha[n])
- elevel ranges from 0 to Nlevels
- E[0] is the ground state energy
- E[n] (n>0) is defined via (E[n]-E[n-1]) = exp(beta[n])
- alpha[n] and beta[n] are our parameters

Fit results for local-local correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 110188.724823
X²/dof 26.929 / 32 26.9/32
E0 1.64243e+00 +- 2.42319e-04 1.6424310 +- 0.0002423
log(dE1) -6.92602e-01 +- 3.36279e-02 -0.6925723 +- 0.0336258
log(dE2) -4.33833e-01 +- 3.88573e-02 -0.4338020 +- 0.0388634
log(Z0) 3.06220e-01 +- 2.44496e-03 0.3062213 +- 0.0024449
log(Z1) 5.90928e-01 +- 3.95491e-02 0.5909641 +- 0.0395471
log(Z2) 7.83968e-01 +- 2.27519e-02 0.7839473 +- 0.0227529
Bootstrap fit results for local-local correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof 26.929 / 32 26.9/32
E0 1.642431 0.000292 0.000209 1.6424 +0.0003 -0.0002
E1 2.142703 0.015135 0.020308 2.1427 +0.0151 -0.0203
E2 2.790724 0.038943 0.041389 2.7908 +0.0388 -0.0419
log(Z0) 0.306220 0.001900 0.003334 0.3062 +0.0019 -0.0033
log(Z1) 0.590928 0.035580 0.047408 0.5910 +0.0355 -0.0475
log(Z2) 0.783968 0.026264 0.021616 0.7839 +0.0263 -0.0216
Fit results for 1S-1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 121008.558143
X²/dof 30.594 / 32 30.6/32 8.47368966673 / 16 Paul's got the wrong number of dof
E0 1.64229e+00 +- 2.35762e-04 1.6422894 +- 0.0002358 1.64228 +- 0.00022 Good agreement from all three
log(dE1) -2.05145e-01 +- 7.17224e-01 -0.1973006 +- 0.6519374 -0.174906 +- 0.61
log(dE2) -1.01862e+00 +- 3.38070e+00 -0.9337536 +- 3.4554010 0.152791 +- 3.7
log(Z0) 1.69205e+00 +- 2.08989e-03 1.6920431 +- 0.0020910 1.69302 +- 0.0017
log(Z1) 6.03994e-01 +- 1.63954e+00 0.6326816 +- 1.3834059 0.648704 +- 0.82
log(Z2) 3.29008e-01 +- 2.89145e+00 0.2787477 +- 2.8781121 0.553633 +- 0.98
Fit results for 2S-2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 8506.473993
X²/dof 28.696 / 32 28.7/32
E0 1.64206e+00 +- 6.62302e-04 1.6420620 +- 0.0006623
log(dE1) -9.03575e-01 +- 5.32491e-02 -0.9035746 +- 0.0532489
log(dE2) -7.41408e-02 +- 2.87793e-01 -0.0741373 +- 0.2877932
log(Z0) 4.66482e-01 +- 7.38988e-03 0.4664818 +- 0.0073899
log(Z1) 1.17649e+00 +- 5.22354e-02 1.1764946 +- 0.0522352
log(Z2) 1.09187e+00 +- 5.05955e-02 1.0918740 +- 0.0505959
Fit results for local+1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 280480.736939
X²/dof 62.431 / 64
E0 1.64240e+00 +- 2.25414e-04
log(dE1) -7.37958e-01 +- 2.60752e-02
log(dE2) -4.68157e-01 +- 2.65329e-02
log(Z0,l) 3.04420e-01 +- 2.25991e-03
log(Z0,1S) 1.69265e+00 +- 1.94607e-03
log(Z1,l) 5.39848e-01 +- 3.06599e-02
log(Z1,1S) -1.22190e-01 +- 2.73177e-01
log(Z2,l) 8.11370e-01 +- 1.49740e-02
log(Z2,1S) 7.87498e-01 +- 7.03348e-02
Fit results for local+1S+2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 705307.249906
X²/dof 130.305 / 96
E0 1.64251e+00 +- 2.21393e-04
log(dE1) -8.60766e-01 +- 2.06709e-02
log(dE2) -5.15699e-01 +- 1.55669e-02
log(Z0,l) 2.98840e-01 +- 2.15222e-03
log(Z0,1S) 1.68889e+00 +- 1.87951e-03
log(Z0,2S) 4.68844e-01 +- 3.22214e-03
log(Z1,l) 4.08479e-01 +- 2.37557e-02
log(Z1,1S) 2.26353e-01 +- 1.12677e-01
log(Z1,2S) 1.17910e+00 +- 1.91798e-02
log(Z2,l) 8.66326e-01 +- 7.88017e-03
log(Z2,1S) 6.15040e-01 +- 7.78124e-02
log(Z2,2S) 9.47982e-01 +- 4.30701e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	110188.724823
X²/dof	26.929 / 32	26.9/32
E0	1.64243e+00 +- 2.42319e-04	1.6424310 +- 0.0002423
log(dE1)	-6.92602e-01 +- 3.36279e-02	-0.6925723 +- 0.0336258
log(dE2)	-4.33833e-01 +- 3.88573e-02	-0.4338020 +- 0.0388634
log(Z0)	3.06220e-01 +- 2.44496e-03	0.3062213 +- 0.0024449
log(Z1)	5.90928e-01 +- 3.95491e-02	0.5909641 +- 0.0395471
log(Z2)	7.83968e-01 +- 2.27519e-02	0.7839473 +- 0.0227529

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²
X²/dof	26.929 / 32	26.9/32
E0	1.642431 0.000292 0.000209	1.6424 +0.0003 -0.0002
E1	2.142703 0.015135 0.020308	2.1427 +0.0151 -0.0203
E2	2.790724 0.038943 0.041389	2.7908 +0.0388 -0.0419
log(Z0)	0.306220 0.001900 0.003334	0.3062 +0.0019 -0.0033
log(Z1)	0.590928 0.035580 0.047408	0.5910 +0.0355 -0.0475
log(Z2)	0.783968 0.026264 0.021616	0.7839 +0.0263 -0.0216

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	121008.558143
X²/dof	30.594 / 32	30.6/32	8.47368966673 / 16	Paul's got the wrong number of dof
E0	1.64229e+00 +- 2.35762e-04	1.6422894 +- 0.0002358	1.64228 +- 0.00022	Good agreement from all three
log(dE1)	-2.05145e-01 +- 7.17224e-01	-0.1973006 +- 0.6519374	-0.174906 +- 0.61
log(dE2)	-1.01862e+00 +- 3.38070e+00	-0.9337536 +- 3.4554010	0.152791 +- 3.7
log(Z0)	1.69205e+00 +- 2.08989e-03	1.6920431 +- 0.0020910	1.69302 +- 0.0017
log(Z1)	6.03994e-01 +- 1.63954e+00	0.6326816 +- 1.3834059	0.648704 +- 0.82
log(Z2)	3.29008e-01 +- 2.89145e+00	0.2787477 +- 2.8781121	0.553633 +- 0.98

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	8506.473993
X²/dof	28.696 / 32	28.7/32
E0	1.64206e+00 +- 6.62302e-04	1.6420620 +- 0.0006623
log(dE1)	-9.03575e-01 +- 5.32491e-02	-0.9035746 +- 0.0532489
log(dE2)	-7.41408e-02 +- 2.87793e-01	-0.0741373 +- 0.2877932
log(Z0)	4.66482e-01 +- 7.38988e-03	0.4664818 +- 0.0073899
log(Z1)	1.17649e+00 +- 5.22354e-02	1.1764946 +- 0.0522352
log(Z2)	1.09187e+00 +- 5.05955e-02	1.0918740 +- 0.0505959

eta_c with narrow priors

We started with the following definition of parameters:

E[0] = 1.65 +- .05
log(E[n]-E[n-1]) = -9.16291e-01 +- 4.81212e-01
log(Z[n,local]) = 1 +- 2 for n=0, 0 +- 2 for n!=0
log(Z[n,1S]) = 1 +- 2 for n=0, 0 +- 2 for n!=0
log(Z[n,2S]) = 1 +- 2 for n=1, 0 +- 2 for n!=1
Nstates = 3
Tmin,Tmax = 1,32
All SVD cuts disabled
Preconditioning of data covariance matrix enabled (rescale to have unit diagonal prior to inversion)
Blocksize = 1
Datasets are:
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_d_d_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_1S_1S_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/pi_2S_2S_0.119_0.119_p000.bz2
Functional form of fit is Sum_{elevel} Z^2 * exp(-E[elevel]*Tmid) * 2*cosh(E[elevel]*(Tmid-t)) where
- Tmid is the middle-timeslice, for our case, Tmid=32.
- Z[n] = exp(alpha[n])
- elevel ranges from 0 to Nlevels
- E[0] is the ground state energy
- E[n] (n>0) is defined via (E[n]-E[n-1]) = exp(beta[n])
- alpha[n] and beta[n] are our parameters

Fit results for local-local correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 748340.813688
X²/dof 28.466 / 32
E0 1.64243e+00 +- 2.42362e-04
log(dE1) -6.95824e-01 +- 3.37434e-02
log(dE2) -4.37890e-01 +- 3.80152e-02
log(Z0) 3.06129e-01 +- 2.44777e-03
log(Z1) 5.86868e-01 +- 3.96370e-02
log(Z2) 7.86163e-01 +- 2.25763e-02
Fit results for 1S-1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 67798.050633
X²/dof 31.096 / 32
E0 1.64229e+00 +- 2.34800e-04
log(dE1) -7.30577e-01 +- 3.50252e-01
log(dE2) -8.07327e-01 +- 3.86047e-01
log(Z0) 1.69189e+00 +- 2.03337e-03
log(Z1) -7.53091e-01 +- 1.73075e+00
log(Z2) 7.97421e-01 +- 7.18224e-02
Fit results for 2S-2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 18490.688657
X²/dof 31.288 / 32
E0 1.64192e+00 +- 6.75241e-04
log(dE1) -9.36941e-01 +- 5.92885e-02
log(dE2) -2.73010e-01 +- 2.32551e-01
log(Z0) 4.64351e-01 +- 7.63830e-03
log(Z1) 1.13783e+00 +- 6.13099e-02
log(Z2) 1.08709e+00 +- 4.63178e-02
Fit results for local+1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 1481863.986694
X²/dof 63.777 / 64
E0 1.64240e+00 +- 2.25414e-04
log(dE1) -7.39510e-01 +- 2.61005e-02
log(dE2) -4.70080e-01 +- 2.62710e-02
log(Z0,l) 3.04398e-01 +- 2.25811e-03
log(Z0,1S) 1.69265e+00 +- 1.94259e-03
log(Z1,l) 5.37747e-01 +- 3.06855e-02
log(Z1,1S) -1.25101e-01 +- 2.72673e-01
log(Z2,l) 8.12302e-01 +- 1.49119e-02
log(Z2,1S) 7.87235e-01 +- 6.98235e-02
Fit results for local+1S+2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 2107200.233945
X²/dof 131.597 / 96 132/96
E0 1.64251e+00 +- 2.21383e-04 1.6425092 +- 0.0002214
log(dE1) -8.61037e-01 +- 2.06612e-02 -0.8610995 +- 0.0206630
log(dE2) -5.16194e-01 +- 1.55377e-02 -0.5162205 +- 0.0155348
log(Z0,l) 2.98858e-01 +- 2.15167e-03 0.2988578 +- 0.0021517
log(Z0,1S) 1.68891e+00 +- 1.87883e-03 1.6889089 +- 0.0018789
log(Z0,2S) 4.68848e-01 +- 3.22131e-03 0.4688481 +- 0.0032214
log(Z1,l) 4.07945e-01 +- 2.37473e-02 0.4078702 +- 0.0237487
log(Z1,1S) 2.25013e-01 +- 1.12793e-01 0.2248930 +- 0.1128157
log(Z1,2S) 1.17887e+00 +- 1.91715e-02 1.1787981 +- 0.0191731
log(Z2,l) 8.66457e-01 +- 7.86905e-03 0.8664808 +- 0.0078680
log(Z2,1S) 6.15265e-01 +- 7.76388e-02 0.6153493 +- 0.0776209
log(Z2,2S) 9.48182e-01 +- 4.30152e-02 0.9483335 +- 0.0429986
Bootstrap fit results for local+1S+2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof 130.305 / 96 132/96
E0 1.642509 0.000342 0.000155 1.6425 +0.0003 -0.0002
E1 2.065339 0.010821 0.011540 2.0652 +0.0108 -0.0115
E2 2.662416 0.020659 0.018248 2.6620 +0.0207 -0.0182
log(Z0,l) 0.298840 0.001987 0.003326 0.2989 +0.0020 -0.0033
log(Z0,1S) 1.688890 0.001653 0.002767 1.6889 +0.0016 -0.0028
log(Z0,2S) 0.468843 0.003345 0.004490 0.4688 +0.0034 -0.0045
log(Z1,l) 0.408456 0.033272 0.031342 0.4079 +0.0331 -0.0308
log(Z1,1S) 0.226329 0.135369 0.138897 0.2249 +0.1360 -0.1408
log(Z1,2S) 1.179087 0.022289 0.026506 1.1788 +0.0222 -0.0263
log(Z2,l) 0.866333 0.010156 0.010428 0.8665 +0.0098 -0.0103
log(Z2,1S) 0.615048 0.086083 0.113912 0.6153 +0.0863 -0.1143
log(Z2,2S) 0.948009 0.051848 0.055409 0.9483 +0.0514 -0.0560

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	2107200.233945
X²/dof	131.597 / 96	132/96
E0	1.64251e+00 +- 2.21383e-04	1.6425092 +- 0.0002214
log(dE1)	-8.61037e-01 +- 2.06612e-02	-0.8610995 +- 0.0206630
log(dE2)	-5.16194e-01 +- 1.55377e-02	-0.5162205 +- 0.0155348
log(Z0,l)	2.98858e-01 +- 2.15167e-03	0.2988578 +- 0.0021517
log(Z0,1S)	1.68891e+00 +- 1.87883e-03	1.6889089 +- 0.0018789
log(Z0,2S)	4.68848e-01 +- 3.22131e-03	0.4688481 +- 0.0032214
log(Z1,l)	4.07945e-01 +- 2.37473e-02	0.4078702 +- 0.0237487
log(Z1,1S)	2.25013e-01 +- 1.12793e-01	0.2248930 +- 0.1128157
log(Z1,2S)	1.17887e+00 +- 1.91715e-02	1.1787981 +- 0.0191731
log(Z2,l)	8.66457e-01 +- 7.86905e-03	0.8664808 +- 0.0078680
log(Z2,1S)	6.15265e-01 +- 7.76388e-02	0.6153493 +- 0.0776209
log(Z2,2S)	9.48182e-01 +- 4.30152e-02	0.9483335 +- 0.0429986

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²
X²/dof	130.305 / 96	132/96
E0	1.642509 0.000342 0.000155	1.6425 +0.0003 -0.0002
E1	2.065339 0.010821 0.011540	2.0652 +0.0108 -0.0115
E2	2.662416 0.020659 0.018248	2.6620 +0.0207 -0.0182
log(Z0,l)	0.298840 0.001987 0.003326	0.2989 +0.0020 -0.0033
log(Z0,1S)	1.688890 0.001653 0.002767	1.6889 +0.0016 -0.0028
log(Z0,2S)	0.468843 0.003345 0.004490	0.4688 +0.0034 -0.0045
log(Z1,l)	0.408456 0.033272 0.031342	0.4079 +0.0331 -0.0308
log(Z1,1S)	0.226329 0.135369 0.138897	0.2249 +0.1360 -0.1408
log(Z1,2S)	1.179087 0.022289 0.026506	1.1788 +0.0222 -0.0263
log(Z2,l)	0.866333 0.010156 0.010428	0.8665 +0.0098 -0.0103
log(Z2,1S)	0.615048 0.086083 0.113912	0.6153 +0.0863 -0.1143
log(Z2,2S)	0.948009 0.051848 0.055409	0.9483 +0.0514 -0.0560

h_c

We started with the following definition of parameters:

E[0] = 2.0 +- 0.6
log(E[n]-E[n-1]) = -9.16291e-01 +- 4.81212e-01
log(Z[n=0]) = 1 +- 2
log(Z[n!=0]) = 0 +- 2
Nstates = 3
Tmin,Tmax[local] = 2,16
Tmin,Tmax[1S] = 2,15
Tmin,Tmax[2S] = 2,10
All SVD cuts disabled
Preconditioning of data covariance matrix enabled (rescale to have unit diagonal prior to inversion)
Blocksize = 1
Datasets are:
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_d_d_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_1S_1S_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_2S_2S_0.119_0.119_p000.bz2
Functional form of fit is Sum_{elevel} Z^2 * exp(-E[elevel]*Tmid) * 2*cosh(E[elevel]*(Tmid-t)) where
- Tmid is the middle-timeslice, for our case, Tmid=32.
- Z[n] = exp(alpha[n])
- elevel ranges from 0 to Nlevels
- E[0] is the ground state energy
- E[n] (n>0) is defined via (E[n]-E[n-1]) = exp(beta[n])
- alpha[n] and beta[n] are our parameters

Fit results for local-local correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 3792418.458006
X²/dof 14.404 / 15
E0 2.00799e+00 +- 1.78115e-02 2.0079519 +- 0.0178319
log(dE1) -9.84754e-01 +- 4.53713e-01 -0.9855137 +- 0.4537165
log(dE2) -7.56732e-01 +- 3.04517e-01 -0.7561813 +- 0.3042313
log(Z0) -4.67671e-01 +- 1.00567e-01 -0.4678958 +- 0.1007287
log(Z1) -2.48589e-01 +- 4.80517e-01 -0.2489305 +- 0.4798682
log(Z2) 2.82869e-01 +- 1.35049e-01 0.2830501 +- 0.1346708
Bootstrap fit results for local-local correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof
E0 2.007989 0.016420 0.208186 2.0080 +0.0164 -0.2048
E1 2.381520 0.029891 0.283130 2.3812 +0.0300 -0.2791
E2 2.850717 0.060360 0.095620 2.8507 +0.0603 -0.0955
log(Z0) -0.467671 0.079717 1.690148 -0.4679 +0.0784 -1.6145
log(Z1) -0.248589 0.081088 0.346861 -0.2489 +0.0836 -0.3402
log(Z2) 0.282869 0.086369 0.000114 0.2831 +0.0845 -0.0006 Errorbars very asymmetric, negative one doesn't agree
Fit results for 1S-1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 148.548588
X²/dof 10.308 / 14 10.3/14
E0 1.97789e+00 +- 7.27617e-03 1.9778980 +- 0.0072704
log(dE1) -9.57117e-01 +- 4.80528e-01 -0.9569252 +- 0.4805310
log(dE2) -9.14724e-01 +- 4.81075e-01 -0.9146911 +- 0.4810755
log(Z0) 9.79912e-01 +- 2.44640e-02 0.9799476 +- 0.0244354
log(Z1) -7.64942e-01 +- 1.39159e+00 -0.7661252 +- 1.3927066
log(Z2) -9.63830e-01 +- 1.79032e+00 -0.9649399 +- 1.7907030
Bootstrap fit results for 1S-1S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof
E0 1.977886 0.004394 0.007071 1.9779 +0.0049 -0.0070
E1 2.361885 0.016620 0.037581 2.3620 +0.0201 -0.0376 bootstrap errors disagree in first digit
E2 2.762512 0.020956 0.041088 2.7626 +0.0243 -0.0407 bootstrap errors disagree at 10% level
log(Z0) 0.979912 0.009472 0.033256 0.9799 +0.0102 -0.0339
log(Z1) -0.764942 0.635831 0.384350 -0.7661 +0.6283 -0.4123 bootstrap errors disagree at 5% level
log(Z2) -0.963830 0.184103 0.240689 -0.9649 +0.1587 -0.3506 bootstrap errors disagree in first digit
Fit results for 2S-2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 31515.761941
X²/dof 4.517 / 9
E0 2.02943e+00 +- 1.09314e-02 2.0293920 +- 0.0109838
log(dE1) -7.30246e-01 +- 4.45969e-01 -0.7303296 +- 0.4460627
log(dE2) -7.68428e-01 +- 4.55936e-01 -0.7681543 +- 0.4560984
log(Z0) 6.41525e-01 +- 3.75962e-02 0.6413890 +- 0.0378600
log(Z1) -7.60076e-01 +- 1.88499e+00 -0.7512787 +- 1.8811065
log(Z2) 5.68546e-01 +- 2.44629e-01 0.5677745 +- 0.2463036
Bootstrap fit results for 2S-2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof
E0 2.029427 0.006929 0.008641 2.0294 +0.0069 -0.0084
E1 2.511218 0.064767 0.068641 2.5111 +0.0644 -0.0724 errors disagree at 5% level
E2 2.974959 0.113414 0.124898 2.9750 +0.1157 -0.1284
log(Z0) 0.641525 0.014640 0.033044 0.6414 +0.0152 -0.0330
log(Z1) -0.760076 0.397152 0.138882 -0.7513 +0.3950 -0.1514 errors disagree at 10% level
log(Z2) 0.568546 0.133709 0.171539 0.5678 +0.1343 -0.1737

Parameter	Damian's result	Jim's result
Initial X²	3792418.458006
X²/dof	14.404 / 15
E0	2.00799e+00 +- 1.78115e-02	2.0079519 +- 0.0178319
log(dE1)	-9.84754e-01 +- 4.53713e-01	-0.9855137 +- 0.4537165
log(dE2)	-7.56732e-01 +- 3.04517e-01	-0.7561813 +- 0.3042313
log(Z0)	-4.67671e-01 +- 1.00567e-01	-0.4678958 +- 0.1007287
log(Z1)	-2.48589e-01 +- 4.80517e-01	-0.2489305 +- 0.4798682
log(Z2)	2.82869e-01 +- 1.35049e-01	0.2830501 +- 0.1346708

Parameter	Damian's result	Jim's result	Comments
Initial X²
X²/dof
E0	2.007989 0.016420 0.208186	2.0080 +0.0164 -0.2048
E1	2.381520 0.029891 0.283130	2.3812 +0.0300 -0.2791
E2	2.850717 0.060360 0.095620	2.8507 +0.0603 -0.0955
log(Z0)	-0.467671 0.079717 1.690148	-0.4679 +0.0784 -1.6145
log(Z1)	-0.248589 0.081088 0.346861	-0.2489 +0.0836 -0.3402
log(Z2)	0.282869 0.086369 0.000114	0.2831 +0.0845 -0.0006	Errorbars very asymmetric, negative one doesn't agree

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	148.548588
X²/dof	10.308 / 14	10.3/14
E0	1.97789e+00 +- 7.27617e-03	1.9778980 +- 0.0072704
log(dE1)	-9.57117e-01 +- 4.80528e-01	-0.9569252 +- 0.4805310
log(dE2)	-9.14724e-01 +- 4.81075e-01	-0.9146911 +- 0.4810755
log(Z0)	9.79912e-01 +- 2.44640e-02	0.9799476 +- 0.0244354
log(Z1)	-7.64942e-01 +- 1.39159e+00	-0.7661252 +- 1.3927066
log(Z2)	-9.63830e-01 +- 1.79032e+00	-0.9649399 +- 1.7907030

Parameter	Damian's result	Jim's result	Comments
Initial X²
X²/dof
E0	1.977886 0.004394 0.007071	1.9779 +0.0049 -0.0070
E1	2.361885 0.016620 0.037581	2.3620 +0.0201 -0.0376	bootstrap errors disagree in first digit
E2	2.762512 0.020956 0.041088	2.7626 +0.0243 -0.0407	bootstrap errors disagree at 10% level
log(Z0)	0.979912 0.009472 0.033256	0.9799 +0.0102 -0.0339
log(Z1)	-0.764942 0.635831 0.384350	-0.7661 +0.6283 -0.4123	bootstrap errors disagree at 5% level
log(Z2)	-0.963830 0.184103 0.240689	-0.9649 +0.1587 -0.3506	bootstrap errors disagree in first digit

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	31515.761941
X²/dof	4.517 / 9
E0	2.02943e+00 +- 1.09314e-02	2.0293920 +- 0.0109838
log(dE1)	-7.30246e-01 +- 4.45969e-01	-0.7303296 +- 0.4460627
log(dE2)	-7.68428e-01 +- 4.55936e-01	-0.7681543 +- 0.4560984
log(Z0)	6.41525e-01 +- 3.75962e-02	0.6413890 +- 0.0378600
log(Z1)	-7.60076e-01 +- 1.88499e+00	-0.7512787 +- 1.8811065
log(Z2)	5.68546e-01 +- 2.44629e-01	0.5677745 +- 0.2463036

Parameter	Damian's result	Jim's result	Comments
Initial X²
X²/dof
E0	2.029427 0.006929 0.008641	2.0294 +0.0069 -0.0084
E1	2.511218 0.064767 0.068641	2.5111 +0.0644 -0.0724	errors disagree at 5% level
E2	2.974959 0.113414 0.124898	2.9750 +0.1157 -0.1284
log(Z0)	0.641525 0.014640 0.033044	0.6414 +0.0152 -0.0330
log(Z1)	-0.760076 0.397152 0.138882	-0.7513 +0.3950 -0.1514	errors disagree at 10% level
log(Z2)	0.568546 0.133709 0.171539	0.5678 +0.1343 -0.1737

h_c

We started with the following definition of parameters:

E[0] = 2.0 +- 0.6
log(E[n]-E[n-1]) = -9.16291e-01 +- 4.81212e-01
log(Z[n=0]) = 1 +- 2
log(Z[n!=0]) = 0 +- 2
Nstates = 3
Tmin,Tmax = 2,15
All SVD cuts disabled
Preconditioning of data covariance matrix enabled (rescale to have unit diagonal prior to inversion)
Blocksize = 1
Datasets are:
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_d_d_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_1S_1S_0.119_0.119_p000.bz2
- lqcd:/oldhome/simone/qcd/onia/l2064f21b676m010m050/Ds_k0.119_csw1.72-t0/fits/corrs/b1_2S_2S_0.119_0.119_p000.bz2
Functional form of fit is Sum_{elevel} Z^2 * exp(-E[elevel]*Tmid) * 2*cosh(E[elevel]*(Tmid-t)) where
- Tmid is the middle-timeslice, for our case, Tmid=32.
- Z[n] = exp(alpha[n])
- elevel ranges from 0 to Nlevels
- E[0] is the ground state energy
- E[n] (n>0) is defined via (E[n]-E[n-1]) = exp(beta[n])
- alpha[n] and beta[n] are our parameters

Fit results for local+1S+2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X² 6099781.788249
X²/dof 35.886 / 42 = 0.854 35.9/42
E0 1.97855e+00 +- 9.24241e-03 1.9785355 +- 0.0093057
log(dE1) -1.51010e+00 +- 3.19966e-01 -1.5097784 +- 0.3198652
log(dE2) -4.96693e-01 +- 7.51517e-02 -0.4966108 +- 0.0751879
log(Z0,l) -6.97943e-01 +- 8.42495e-02 -0.6980366 +- 0.0846268
log(Z0,1S) 9.77174e-01 +- 4.09081e-02 0.9770397 +- 0.0412248
log(Z0,2S) 3.79396e-01 +- 6.73019e-02 0.3792892 +- 0.0676448
log(Z1,l) -3.57252e-01 +- 1.46158e-01 -0.3569833 +- 0.1461035
log(Z1,1S) -6.97609e-01 +- 1.48453e+00 -0.6898565 +- 1.4697194 1% disagreement, but huge errors
log(Z1,2S) 3.33134e-01 +- 1.12546e-01 0.3335169 +- 0.1125486
log(Z2,l) 3.53934e-01 +- 2.48176e-02 0.3539141 +- 0.0248267
log(Z2,1S) -1.05817e+00 +- 1.80858e+00 -1.0730568 +- 1.8201065
log(Z2,2S) 3.23981e-01 +- 1.88384e-01 0.3235137 +- 0.1885131
Bootstrap fit results for local+1S+2S correlator:

Parameter Damian's result Jim's result Paul's result Comments
Initial X²
X²/dof 35.886 / 42 = 0.854 35.9/42
E0 1.978553 0.004804 0.014121 1.9785 +0.0041 -0.0135
E1 2.199442 0.044283 0.062380 2.1995 +0.0398 -0.0517 Damian's errors 10% larger than Jim's
E2 2.807982 0.079793 0.045456 2.8081 +0.0658 -0.0383 Damian's errors 10% larger than Jim's
log(Z0,l) -0.697943 0.048281 0.111731 -0.6980 +0.0339 -0.1058
log(Z0,1S) 0.977174 0.013388 0.061790 0.9770 +0.0134 -0.0574
log(Z0,2S) 0.379396 0.034993 0.108507 0.3793 +0.0335 -0.1064
log(Z1,l) -0.357252 0.111480 0.059415 -0.3570 +0.1046 -0.0468
log(Z1,1S) -0.697609 0.770996 1.083922 -0.6899 +0.7531 -1.0880
log(Z1,2S) 0.333134 0.143688 0.043497 0.3335 +0.1412 -0.0406
log(Z2,l) 0.353934 0.025183 0.016970 0.3539 +0.0249 -0.0175
log(Z2,1S) -1.058168 0.409742 0.689854 -1.0731 +0.4040 -0.7245
log(Z2,2S) 0.323981 0.137795 0.195635 0.3235 +0.1175 -0.1853 Damian's errors 10% larger than Jim's

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	280480.736939
X²/dof	62.431 / 64
E0	1.64240e+00 +- 2.25414e-04
log(dE1)	-7.37958e-01 +- 2.60752e-02
log(dE2)	-4.68157e-01 +- 2.65329e-02
log(Z0,l)	3.04420e-01 +- 2.25991e-03
log(Z0,1S)	1.69265e+00 +- 1.94607e-03
log(Z1,l)	5.39848e-01 +- 3.06599e-02
log(Z1,1S)	-1.22190e-01 +- 2.73177e-01
log(Z2,l)	8.11370e-01 +- 1.49740e-02
log(Z2,1S)	7.87498e-01 +- 7.03348e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	705307.249906
X²/dof	130.305 / 96
E0	1.64251e+00 +- 2.21393e-04
log(dE1)	-8.60766e-01 +- 2.06709e-02
log(dE2)	-5.15699e-01 +- 1.55669e-02
log(Z0,l)	2.98840e-01 +- 2.15222e-03
log(Z0,1S)	1.68889e+00 +- 1.87951e-03
log(Z0,2S)	4.68844e-01 +- 3.22214e-03
log(Z1,l)	4.08479e-01 +- 2.37557e-02
log(Z1,1S)	2.26353e-01 +- 1.12677e-01
log(Z1,2S)	1.17910e+00 +- 1.91798e-02
log(Z2,l)	8.66326e-01 +- 7.88017e-03
log(Z2,1S)	6.15040e-01 +- 7.78124e-02
log(Z2,2S)	9.47982e-01 +- 4.30701e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	748340.813688
X²/dof	28.466 / 32
E0	1.64243e+00 +- 2.42362e-04
log(dE1)	-6.95824e-01 +- 3.37434e-02
log(dE2)	-4.37890e-01 +- 3.80152e-02
log(Z0)	3.06129e-01 +- 2.44777e-03
log(Z1)	5.86868e-01 +- 3.96370e-02
log(Z2)	7.86163e-01 +- 2.25763e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	67798.050633
X²/dof	31.096 / 32
E0	1.64229e+00 +- 2.34800e-04
log(dE1)	-7.30577e-01 +- 3.50252e-01
log(dE2)	-8.07327e-01 +- 3.86047e-01
log(Z0)	1.69189e+00 +- 2.03337e-03
log(Z1)	-7.53091e-01 +- 1.73075e+00
log(Z2)	7.97421e-01 +- 7.18224e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	18490.688657
X²/dof	31.288 / 32
E0	1.64192e+00 +- 6.75241e-04
log(dE1)	-9.36941e-01 +- 5.92885e-02
log(dE2)	-2.73010e-01 +- 2.32551e-01
log(Z0)	4.64351e-01 +- 7.63830e-03
log(Z1)	1.13783e+00 +- 6.13099e-02
log(Z2)	1.08709e+00 +- 4.63178e-02

Parameter	Damian's result	Jim's result	Paul's result	Comments
Initial X²	1481863.986694
X²/dof	63.777 / 64
E0	1.64240e+00 +- 2.25414e-04
log(dE1)	-7.39510e-01 +- 2.61005e-02
log(dE2)	-4.70080e-01 +- 2.62710e-02
log(Z0,l)	3.04398e-01 +- 2.25811e-03
log(Z0,1S)	1.69265e+00 +- 1.94259e-03
log(Z1,l)	5.37747e-01 +- 3.06855e-02
log(Z1,1S)	-1.25101e-01 +- 2.72673e-01
log(Z2,l)	8.12302e-01 +- 1.49119e-02
log(Z2,1S)	7.87235e-01 +- 6.98235e-02

Parameter	Damian's result	Jim's result	Comments
Initial X²	6099781.788249
X²/dof	35.886 / 42 = 0.854	35.9/42
E0	1.97855e+00 +- 9.24241e-03	1.9785355 +- 0.0093057
log(dE1)	-1.51010e+00 +- 3.19966e-01	-1.5097784 +- 0.3198652
log(dE2)	-4.96693e-01 +- 7.51517e-02	-0.4966108 +- 0.0751879
log(Z0,l)	-6.97943e-01 +- 8.42495e-02	-0.6980366 +- 0.0846268
log(Z0,1S)	9.77174e-01 +- 4.09081e-02	0.9770397 +- 0.0412248
log(Z0,2S)	3.79396e-01 +- 6.73019e-02	0.3792892 +- 0.0676448
log(Z1,l)	-3.57252e-01 +- 1.46158e-01	-0.3569833 +- 0.1461035
log(Z1,1S)	-6.97609e-01 +- 1.48453e+00	-0.6898565 +- 1.4697194	1% disagreement, but huge errors
log(Z1,2S)	3.33134e-01 +- 1.12546e-01	0.3335169 +- 0.1125486
log(Z2,l)	3.53934e-01 +- 2.48176e-02	0.3539141 +- 0.0248267
log(Z2,1S)	-1.05817e+00 +- 1.80858e+00	-1.0730568 +- 1.8201065
log(Z2,2S)	3.23981e-01 +- 1.88384e-01	0.3235137 +- 0.1885131

etac with wide priors

etac with narrow priors

hc

hc

eta_c with wide priors

eta_c with narrow priors

h_c

h_c