In HISKP-LQCD/hadron: Analysis Framework for Monte Carlo Simulation Data in Physics

knitr::opts_chunk$set(echo = TRUE)

devtools::load_all()

offset_val <- 50
offset_err <- 0.1
offset_matrix <- matrix(
  rnorm(nrow(samplecf$cf) * ncol(samplecf$cf), offset_val, offset_err),
  nrow = nrow(samplecf$cf),
  ncol = ncol(samplecf$cf))

samplecf_offset <- samplecf
samplecf_offset$cf <- samplecf$cf + offset_matrix

samplecf_boot <- bootstrap.cf(samplecf)
samplecf_offset_boot <- bootstrap.cf(samplecf_offset)

plot(samplecf_boot, col = 'blue')
plot(samplecf_offset_boot, col = 'red', rep = TRUE)

effmass <- bootstrap.effectivemass(samplecf_boot)
effmass_offset <- bootstrap.effectivemass(samplecf_offset_boot)

plot(effmass_offset, col = 'red')
plot(effmass, col = 'blue', rep = TRUE)

fit <- new_matrixfit(samplecf_boot, 8, 20)
plot(fit, log = 'y')
residual_plot(fit)

fit_offset <- new_matrixfit(samplecf_offset_boot, 8, 20)
plot(fit_offset, log = 'y')
residual_plot(fit_offset)

fit_offset2 <- new_matrixfit(samplecf_offset_boot, 8, 20, model = 'single_constant')
plot(fit_offset2, log = 'y')
residual_plot(fit_offset2)

HISKP-LQCD/hadron documentation built on Sept. 17, 2022, 10:03 a.m.