From f15a4108ba4e62f8f02c68c55e925edba39726c3 Mon Sep 17 00:00:00 2001
From: "Matthew F. McEneaney" <mfmceneaney@gmail.com>
Date: Wed, 4 Feb 2026 15:25:10 -0500
Subject: [PATCH 1/2] fix: Add to vectors instead of accessing elements.

---
 saga/include/hbanalysis.h | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/saga/include/hbanalysis.h b/saga/include/hbanalysis.h
index a3ebee0..210396b 100644
--- a/saga/include/hbanalysis.h
+++ b/saga/include/hbanalysis.h
@@ -1225,8 +1225,8 @@ void getKinBinnedHB(
             epsilon = eps_bg_pdf;
             epsilon_err = eps_bg_pdf_err;
             for (int idx=0; idx<nparams; idx++) {
-                ys_sb[idx] = asymfit_result_sb[k2++];
-                eys_sb[idx] = asymfit_result_sb[k2++];
+                ys_sb.push_back(asymfit_result_sb[k2++]);
+                eys_sb.push_back(asymfit_result_sb[k2++]);
                 ys[idx]  = (ys[idx] - epsilon * ys_sb[idx]) / (1.0 - epsilon);
                 eys[idx] = TMath::Sqrt(eys[idx]*eys[idx] + epsilon * epsilon * eys_sb[idx]*eys_sb[idx]) / (1.0 - epsilon);
             }

From b3f4c4bcc254fdc3aec8279550542a5fceaff140 Mon Sep 17 00:00:00 2001
From: "Matthew F. McEneaney" <mfmceneaney@gmail.com>
Date: Wed, 4 Feb 2026 17:06:05 -0500
Subject: [PATCH 2/2] fix: Updated debug messages.

---
 saga/include/hbanalysis.h | 3 +--
 1 file changed, 1 insertion(+), 2 deletions(-)

diff --git a/saga/include/hbanalysis.h b/saga/include/hbanalysis.h
index 210396b..7900b1b 100644
--- a/saga/include/hbanalysis.h
+++ b/saga/include/hbanalysis.h
@@ -997,7 +997,7 @@ void getKinBinnedHB(
 
         // Create signal region dataset for sideband subtraction
         if (use_sb_subtraction) {
-            LOG_DEBUG(Form("[%s]: Creating sideband dataset...", method_name.c_str()));
+            LOG_DEBUG(Form("[%s]: Creating signal dataset...", method_name.c_str()));
             data::createDataset(
                 binframe_sg,
                 ws_sg, //NOTE: Use separate sideband workspace for dataset, variable, and pdf name uniqueness.
@@ -1221,7 +1221,6 @@ void getKinBinnedHB(
         if (use_sb_subtraction) {
             LOG_DEBUG(Form("[%s]: Applying sideband subtraction...", method_name.c_str()));
             int k2 = 1 + 2 * (binvars.size() + depolvars.size() + rawasymvars.size());
-            LOG_DEBUG(Form("[%s]: k2 = %d AND asymfit_result_sb.size() = %d", method_name.c_str(), k2, asymfit_result_sb.size()));
             epsilon = eps_bg_pdf;
             epsilon_err = eps_bg_pdf_err;
             for (int idx=0; idx<nparams; idx++) {