BGHistogrammer Class Reference

Inheritance diagram for BGHistogrammer:

Public Member Functions
def	initialize (self)
def	event (self)
def	setErrors (self, hist)
def	terminate (self)

Public Attributes
	hist
	list of histograms

Static Public Attributes
list	simHits
	list of SimHits to be histogramed
len	n = len(simHits)
	length of simHits list

Detailed Description

Make validation histograms for BG mixer.

Definition at line 29 of file validate.py.

Member Function Documentation

event()

def event

(

self

)

 Event processor: fill histograms 

Definition at line 68 of file validate.py.

    def event(self):
        ''' Event processor: fill histograms '''
 
        for i in range(self.n):
            hits = Belle2.PyStoreArray(self.simHits[i])
            for hit in hits:
                time = hit.getGlobalTime() / 1000
                self.hist[i].Fill(time)
 

initialize()

def initialize

(

self

)

 Initialize the Module: set histogram axis titles, description and check

Definition at line 48 of file validate.py.

    def initialize(self):
        ''' Initialize the Module: set histogram axis titles, description and check'''
 
        
        self.hist = [TH1F('h' + str(i), self.simHits[i], 400, -20, 20) for i in range(self.n)]
 
        for i in range(self.n):
            self.hist[i].GetXaxis().SetTitle('time [#mus]')
            self.hist[i].GetYaxis().SetTitle('entries/bin')
            descr = TNamed('Description', 'Time distribution of ' + self.simHits[i] +
                           ' for mixed background')
            self.hist[i].GetListOfFunctions().Add(descr)
            check = TNamed('Check', 'Distribution must be flat in the time window. ' +
                           'Bin statistics is not Poissonian.')
            self.hist[i].GetListOfFunctions().Add(check)
            contact = TNamed('Contact', 'marko.staric@ijs.si')
            self.hist[i].GetListOfFunctions().Add(contact)
            options = TNamed('MetaOptions', 'shifter')
            self.hist[i].GetListOfFunctions().Add(options)
 

setErrors()

def setErrors	(		self,
			hist
	)

 Set bin errors - they are not Poissonian

Definition at line 77 of file validate.py.

    def setErrors(self, hist):
        ''' Set bin errors - they are not Poissonian'''
        hmax = hist.GetMaximum()
        temp = TH1F('temp', 'temporary', 100, 0, hmax * 1.1)
        for i in range(hist.GetNbinsX()):
            h = hist.GetBinContent(i + 1)
            if h > hmax * 0.2:
                temp.Fill(h)
        mean = temp.GetMean()
        rms = temp.GetStdDev()
        if mean > 0:
            for i in range(hist.GetNbinsX()):
                h = hist.GetBinContent(i + 1)
                e = 0
                if h > 0:
                    e = rms * math.sqrt(h / mean)
                hist.SetBinError(i + 1, e)
 

terminate()

def terminate

(

self

)

 Write histograms to the file.

Definition at line 95 of file validate.py.

    def terminate(self):
        """ Write histograms to the file."""
 
        tfile = TFile('timeDistributions.root', 'recreate')
        for i in range(self.n):
            self.setErrors(self.hist[i])
            self.hist[i].Write()
        tfile.Close()
 
 
# Create path

Member Data Documentation

hist

hist

list of histograms

Definition at line 52 of file validate.py.

n

len n = len(simHits)

static

length of simHits list

Definition at line 46 of file validate.py.

simHits

list simHits

static

Initial value:

=  [
        'PXDSimHits',
        'SVDSimHits',
        'CDCSimHits',
        'TOPSimHits',
        'ARICHSimHits',
        'ECLHits',
        'KLMSimHits',
    ]

list of SimHits to be histogramed

Definition at line 36 of file validate.py.

---

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

• background/validation/validate.py