Inheritance diagram for calibrateGlobalT0Offline:

Collaboration diagram for calibrateGlobalT0Offline:

Public Member Functions
def	initialize (self)

def	event (self)

def	getHistogramToFit (self)

def	terminate (self)

Public Attributes
	file
	Output root file.

	bunchTimeSep
	bunch separation time

	h1
	histogram of current offset

	h2
	histogram of current offset vs event number

	h1a
	histogram of current offset, wrap-around into [-1/2, 1/2] of bunch sep. More...

	h2a
	histogram of current offset vs event number, wrap-around into [-1/2, 1/2]

	h1b
	histogram of current offset, wrap-around into [0, 1] of bunch sep. More...

	h2b
	histogram of current offset vs event number, wrap-around into [0, 1]

	expNo
	experiment number

	run
	run number formatted as string with leading zeros

	t0
	common T0 used for the calibration of the input file

	t0Err
	common T0 uncertainty

Detailed Description

** Description **
Module to calculate the commont time offset of the TOP detector.
It selects two-track events, saves the bunch finder time offset in an histrogram and
fits it to extract the global time offset with respect to the RF clock.
It should be run in an individual job for each run, but in future it may be further
authomatize to run on multiple runs in the same job.

**Output**
Root file containing:
* The histogram of the bunch finder time offset
* The fit function
* The fit parameters ina tree format (usefuly when merging several files)

**Contributors**
* Marko Staric
* Umberto Tamponi

Definition at line 76 of file cdst_calibrateCommonT0.py.

Member Function Documentation

◆ event()

def event ( self )

Fills the histogram of the time offset

Definition at line 166 of file cdst_calibrateCommonT0.py.

     def event(self):
         """
         Fills the histogram of the time offset
         """
         recBunch = Belle2.PyStoreObj('TOPRecBunch')
         evtMetaData = Belle2.PyStoreObj('EventMetaData')
         if not recBunch.isValid():
             return
         if recBunch.isReconstructed() and recBunch.getNumTracks() == 2:
             offset = recBunch.getCurrentOffset()
             evtNum = evtMetaData.getEvent()
             self.h1.Fill(offset)
             self.h2.Fill(evtNum, offset)
             # wrap-around into [-1/2, 1/2] of bunch cycle
             a = offset - round(offset / self.bunchTimeSep, 0) * self.bunchTimeSep
             self.h1a.Fill(a)
             self.h2a.Fill(evtNum, a)
             # wrap-around into [0, 1] of bunch cycle
             b = offset - round(offset / self.bunchTimeSep - 0.5, 0) * self.bunchTimeSep
             self.h1b.Fill(b)
             self.h2b.Fill(evtNum, b)
  

◆ getHistogramToFit()

def getHistogramToFit ( self )

Selects a histogram with the peak closest to histogram center

Definition at line 188 of file cdst_calibrateCommonT0.py.

◆ initialize()

def initialize ( self )

Creates the histogram used for the commoT0 calculation, and
takes the necessary objects from the DataStore

Definition at line 96 of file cdst_calibrateCommonT0.py.

◆ terminate()

def terminate ( self )

Performs the fit of the bunch finder offset distribution,
using a line + gaussian function.
The results of the fit are saved in a tree

Definition at line 199 of file cdst_calibrateCommonT0.py.

Member Data Documentation

◆ h1a

h1a

histogram of current offset, wrap-around into [-1/2, 1/2] of bunch sep.

time

Definition at line 124 of file cdst_calibrateCommonT0.py.

◆ h1b

h1b

histogram of current offset, wrap-around into [0, 1] of bunch sep.