add fill_hist funtion

hbt/config/hist_hooks.py

@@ -1,6 +1,7 @@
 import order as od
 import numpy as np
 from columnflow.util import maybe_import
+from columnflow.columnar_util import fill_hist
 
 hist = maybe_import("hist")
 
@@ -82,16 +83,11 @@ def ana_ABCD(task, hists):
                 else:
                     create_regions(rate, rate_int_dic, rate_int, cat_id, regions, mc_hist, data_hist, cat_name, cat_parent)
 
-            # initialise empty QCD histograms
-            hists[cfg.processes.n.qcd] = rate.copy()
-            hists[cfg.processes.n.qcd].reset()
+            # NOTE: the ABCD regions have the same entries for both parent categories (incl and 2j) why ?
 
-            # for cat_parent in regions.keys():
-            # rates[cat_parent] = {}
-            # per bin conversion factor
-            # rates[cat_parent]["scale_factor"] = rates[cat_parent]["C"].values() / rates[cat_parent]["D"].values()
-            # replace NaN with zero in case of division by zero
-            #scale_factor = np.nan_to_num(rates["C"].values() / rates["D"].values())
+            # initialise empty QCD histograms
+            qcd_hist = rate.copy()
+            qcd_hist = qcd_hist.reset()
 
             scale_factor = {}
             for cat_parent in regions.keys():
@@ -102,7 +98,9 @@ def ana_ABCD(task, hists):
             embed()
             quit()
 
-            # hists[cfg.processes.n.qcd].axes[0][0] = rates[cat_parent]["B"] * rates[cat_parent]["scale_factor"]
+            # add QCD histogram to hists
+            # TODO: BUG in fill_hist function (ValueError: got multi-dimensional hist but only one dimensional data)
+            # hists[cfg.processes.n.qcd] = fill_hist(qcd_hist, rate) * scale_factor["incl"]
 
             return hists