Skip to content

Commit 11d4ece

Browse files
DelloStrittoLuigi Dello Stritto
andauthored
fix processer and hadron analyzer (alisw#972)
* fix processer and hadron analyzer * fix * fix * fix --------- Co-authored-by: Luigi Dello Stritto <[email protected]>
1 parent 0e7a4e3 commit 11d4ece

File tree

2 files changed

+13
-6
lines changed

2 files changed

+13
-6
lines changed

machine_learning_hep/analysis/analyzerdhadrons.py

Lines changed: 9 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -93,6 +93,7 @@ def __init__(self, datap, case, typean, period):
9393
self.include_reflection = datap["analysis"][self.typean].get(
9494
"include_reflection", False)
9595

96+
self.p_nevents = datap["analysis"][self.typean]["nevents"]
9697
self.p_sigmamb = datap["analysis"]["sigmamb"]
9798
self.p_br = datap["ml"]["opt"]["BR"]
9899

@@ -494,9 +495,14 @@ def makenormyields(self): # pylint: disable=import-outside-toplevel, too-many-b
494495
filemass = TFile.Open(self.n_filemass)
495496
hevents = filemass.Get("all_events")
496497
hselevents = filemass.Get("sel_events")
497-
norm, selnorm = self.calculate_norm(self.logger, hevents, hselevents)
498-
histonorm.SetBinContent(1, selnorm)
499-
self.logger.warning("Number of events %d", norm)
498+
499+
if self.p_nevents is not None:
500+
selnorm = self.p_nevents
501+
else:
502+
norm, selnorm = self.calculate_norm(self.logger, hevents, hselevents)
503+
histonorm.SetBinContent(1, selnorm)
504+
self.logger.warning("Number of events %d", norm)
505+
500506
self.logger.warning("Number of events after event selection %d", selnorm)
501507

502508
if self.p_dobkgfromsideband:

machine_learning_hep/processer.py

Lines changed: 4 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -250,7 +250,7 @@ def __init__(self, case, datap, run_param, mcordata, p_maxfiles, # pylint: disab
250250
for ipt in range(self.p_nptbins)]
251251
self.f_evt_count_ml = os.path.join(self.d_pkl_ml, self.n_evt_count_ml)
252252

253-
self.lpt_gensk_sl = [self.n_gen_sl.replace(".p", "_%s%d_%d.p" %
253+
self.lpt_gensk_sl = [self.n_gen_sl.replace(".p", "_%s%d_%d.p" %
254254
(self.v_var_binning, self.lpt_anbinmin[i], self.lpt_anbinmax[i]))
255255
for i in range(self.p_nptbins)] if self.n_gen_sl else None
256256

@@ -282,8 +282,9 @@ def __init__(self, case, datap, run_param, mcordata, p_maxfiles, # pylint: disab
282282
self.lpt_gensk[ipt]) for ipt in range(self.p_nptbins)]
283283
self.lpt_gendecmerged = [os.path.join(self.d_pkl_decmerged, self.lpt_gensk[ipt])
284284
for ipt in range(self.p_nptbins)]
285-
self.mptfiles_gensk_sl = [createlist(self.d_pklsk, self.l_path,
286-
self.lpt_gensk_sl[ipt]) for ipt in range(self.p_nptbins)] if self.lpt_gensk_sl else None
285+
self.mptfiles_gensk_sl = [createlist(self.d_pklsk, self.l_path, self.lpt_gensk_sl[ipt])
286+
for ipt in range(self.p_nptbins)] \
287+
if self.lpt_gensk_sl else None
287288

288289
# self.triggerbit = datap["analysis"][self.typean]["triggerbit"]
289290
self.runlistrigger = runlisttrigger

0 commit comments

Comments
 (0)