python.physvalPostProcessingTools.ResolutionComputation Class Reference

Collaboration diagram for python.physvalPostProcessingTools.ResolutionComputation:

Public Member Functions
def	__init__ (self, histogram_name)
def	__call__ (self, root_file)

Static Public Member Functions
def	from_yaml (fragment)

Public Attributes
	histogram_name

Private Member Functions
def	_compute_s (self, hist, cum, prob)

Detailed Description

Definition at line 161 of file physvalPostProcessingTools.py.

Constructor & Destructor Documentation

__init__()

def python.physvalPostProcessingTools.ResolutionComputation.__init__	(		self,
			histogram_name
	)

Definition at line 163 of file physvalPostProcessingTools.py.

    def __init__(self, histogram_name):
        self.histogram_name = histogram_name
 

Member Function Documentation

__call__()

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

Definition at line 166 of file physvalPostProcessingTools.py.

    def __call__(self, root_file):
        hist = find_histo(root_file, self.histogram_name)
        if not hist:
            if crash_on_error:
                raise RuntimeError(f"Histogram '{self.histogram_name}' not found.")
            else:
                logging.error(f"Histogram '{self.histogram_name}' not found.")
                return
 
        # Step 1: Compute normalized cumulative sum
        cum = [hist.GetBinContent(i) for i in range(0, hist.GetNbinsX() + 2)]
        for i in range(1, len(cum)):
            cum[i] += cum[i - 1]
        total = cum[-1] or 1.0
        cum = [x / total for x in cum]
 
        # Step 2: Compute smallest intervals
        s68 = self._compute_s(hist, cum, 0.683)
        s90 = self._compute_s(hist, cum, 0.90)
 
        logging.info(f"Resolution widths for {self.histogram_name} -> s68: {s68:.4f}, s90: {s90:.4f}")
        logging.info("Todo: Take a 2D histo, take the quantity say xaxis, make its small intervals %s", "and then compute the Resolution widths and then save as a 1D histo with quantity on xaxis and Resolution widths on yaxis")
 

_compute_s()

def python.physvalPostProcessingTools.ResolutionComputation._compute_s (   self,   hist,   cum,   prob  )

private

Definition at line 189 of file physvalPostProcessingTools.py.

    def _compute_s(self, hist, cum, prob):
        # Compute the smallest interval that contains the given probability.
        min_width = float("inf")
        best = (0, 0)
        for i in range(len(cum)):
            target = cum[i] + prob
            if target > 1:
                break
            for j in range(i + 1, len(cum)):
                if cum[j] >= target:
                    width = hist.GetBinCenter(j) - hist.GetBinCenter(i)
                    if width < min_width:
                        min_width = width
                        best = (i, j)
                    break
        return 0.5 * (hist.GetBinCenter(best[1]) - hist.GetBinCenter(best[0]))
 

from_yaml()

def python.physvalPostProcessingTools.ResolutionComputation.from_yaml (   fragment )

static

Definition at line 207 of file physvalPostProcessingTools.py.

    def from_yaml(fragment):
        return ResolutionComputation(
            histogram_name=fragment["histogram"]
        )
 

Member Data Documentation

histogram_name

python.physvalPostProcessingTools.ResolutionComputation.histogram_name

Definition at line 164 of file physvalPostProcessingTools.py.

---

The documentation for this class was generated from the following file:

• physvalPostProcessingTools.py