BGHistogrammer Class Reference

Inheritance diagram for BGHistogrammer:

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

Public Attributes
list	hist = [TH1F('h' + str(i), self.simHits[i], 400, -20, 20) for i in range(self.n)]
	list of histograms

Static Public Attributes
list	simHits
	list of SimHits to be histogramed
	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()

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()

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()

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()

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

list hist = [TH1F('h' + str(i), self.simHits[i], 400, -20, 20) for i in range(self.n)]

list of histograms

Definition at line 52 of file validate.py.

n

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