Collaboration diagram for python.physvalPostProcessingTools.EfficiencyComputation:

Public Member Functions
def	__init__ (self, numerator, denominator, output)

def	__call__ (self, root_file)

Static Public Member Functions
def	from_yaml (fragment)

Public Attributes
	numerator

	denominator

	output

Detailed Description

Definition at line 109 of file physvalPostProcessingTools.py.

Constructor & Destructor Documentation

◆ init()

def python.physvalPostProcessingTools.EfficiencyComputation.__init__	(	self,
		numerator,
		denominator,
		output
	)

Definition at line 110 of file physvalPostProcessingTools.py.

     def __init__(self, numerator, denominator, output):
         self.numerator = numerator
         self.denominator = denominator
         self.output = output
  

Member Function Documentation

◆ call()

def python.physvalPostProcessingTools.EfficiencyComputation.__call__	(	self,
		root_file
	)

Definition at line 115 of file physvalPostProcessingTools.py.

     def __call__(self, root_file):
         num = find_histo(root_file, self.numerator)
         den = find_histo(root_file, self.denominator)
         if not num or not den:
             if crash_on_error:
                 raise RuntimeError(f"Missing histogram: {self.numerator} or {self.denominator}")
             else:
                 logging.error(f"Missing histogram: {self.numerator} or {self.denominator}")
                 return
  
         if not ROOT.TEfficiency.CheckConsistency(num, den):
             logging.warning("Inconsistency detected between numerator and denominator histograms.")
             logging.info("Attempting to fix...")
             for i in range(1, num.GetNbinsX() + 1):
                 num_val = num.GetBinContent(i)
                 den_val = den.GetBinContent(i)
                 if den_val == 0:
                     den.SetBinContent(i, 1e-6)
                     den_val = 1e-6  # avoid zero denominator by assigning a tiny epsilon
                     logging.info(f"Setting denominator to a smaller value for bin {i} as it is zero. Any other options?")
                 if num_val > den_val:
                     logging.info(f"For Bin {i}: Num ({num_val}) > Den ({den_val})! Adjusting by setting Den to Num. Any other options?")
                     den.SetBinContent(i, num_val)
             if not ROOT.TEfficiency.CheckConsistency(num, den):
                 logging.error("Unable to fix histogram inconsistencies. Aborting efficiency calculation.")
                 return
  
         eff = ROOT.TEfficiency(num, den)
         path_parts = self.output.strip("/").split("/")
         *dir_path, obj_name = path_parts
         directory_path = "/".join(dir_path)
         eff.SetName(obj_name)
         eff.SetTitle(f"{self.numerator}/{self.denominator};{num.GetXaxis().GetTitle()};Efficiency")
         ensure_directory_exists(root_file, directory_path)
         root_file.cd(directory_path)
         eff.Write()
         logging.info(f"Saved efficiency '{obj_name}' in '{directory_path}'.")
  

◆ from_yaml()

def python.physvalPostProcessingTools.EfficiencyComputation.from_yaml ( fragment )

static

Definition at line 154 of file physvalPostProcessingTools.py.

     def from_yaml(fragment):
         return EfficiencyComputation(
             numerator=fragment["numerator"],
             denominator=fragment["denominator"],
             output=fragment["output"]
         )
  