We report our progress in data analysis on two-point correlation functions of the $B$ meson using sequential Bayesian method.
The data set of measurement is obtained using the Oktay-Kronfeld (OK) action for the bottom quarks (valence quarks) and the HISQ action for the light quarks on the MILC HISQ lattices.
We find that the old initial guess for the $\chi^2$ minimizer in the fitting code is poor enough to slow down the analysis somewhat.
In order to find a better initial guess, we adopt the Newton method.
We find that the Newton method provides a natural test to check whether the $\chi^2$ minimizer finds a local minimum or the global minimum, and it also reduces the number of iterations dramatically.