Skip to content

Commit 465229f

Browse files
authored
Fix saving fit results for other levels and per_ptjet true (alisw#1006)
1 parent f4c720c commit 465229f

File tree

1 file changed

+6
-5
lines changed

1 file changed

+6
-5
lines changed

machine_learning_hep/analysis/analyzer_jets.py

Lines changed: 6 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -114,7 +114,6 @@ def __init__(self, datap, case, typean, period):
114114
self.h_effnew_pthf = {"pr": None, "np": None}
115115
self.hfeeddown_det = {"mc": {}, "data": {}}
116116
self.h_reflcorr = create_hist("h_reflcorr", ";#it{p}_{T}^{HF} (GeV/#it{c})", self.bins_candpt)
117-
self.h_fit_results = {}
118117
self.h_fit_results = {
119118
level: {
120119
param: create_hist(
@@ -127,7 +126,7 @@ def __init__(self, datap, case, typean, period):
127126
("#it{#mu}", "#it{#sigma}", "significance", "#it{#chi}^{2}"),
128127
)
129128
}
130-
for level in ("data", "mc")
129+
for level in ("data", "mc", "mcsig", "mcrefl")
131130
}
132131
self.n_events = {}
133132
self.n_colls_read = {}
@@ -394,7 +393,8 @@ def _roofit_mass(self, level, hist, ipt, pdfnames, param_names, fitcfg, roows=No
394393
c = TCanvas()
395394

396395
chi2 = frame.chiSquare()
397-
self.h_fit_results[level]["chi2"].SetBinContent(ipt + 1, chi2)
396+
if "ptjet" not in filename:
397+
self.h_fit_results[level]["chi2"].SetBinContent(ipt + 1, chi2)
398398
if chi2 > 5.0:
399399
self.logger.error(
400400
"Roofit fit is too bad: %s, ipt: %d, pthf: %g-%g, Chi2 = %g",
@@ -416,8 +416,9 @@ def _roofit_mass(self, level, hist, ipt, pdfnames, param_names, fitcfg, roows=No
416416
ws, res, pdfnames, param_names, mean_sgn, sigma_sgn
417417
)
418418
add_text_info_perf(textInfoLeft, sig, sig_err, bkg, bkg_err, s_over_b, s_over_b_err, signif, signif_err)
419-
self.h_fit_results[level]["significance"].SetBinContent(ipt + 1, signif)
420-
self.h_fit_results[level]["significance"].SetBinError(ipt + 1, signif_err)
419+
if "ptjet" not in filename:
420+
self.h_fit_results[level]["significance"].SetBinContent(ipt + 1, signif)
421+
self.h_fit_results[level]["significance"].SetBinError(ipt + 1, signif_err)
421422

422423
frame.Draw()
423424
textInfoRight.Draw()

0 commit comments

Comments
 (0)