Histogrammer Class Reference

Inheritance diagram for Histogrammer:

Collaboration diagram for Histogrammer:

Public Member Functions
def	initialize (self)
def	event (self)
def	terminate (self)

Public Attributes
	hist
	list of digits histograms
	hist2
	list of occupancy-like histograms
	nPXDL1
	number of PXD pixels (L1)
	nPXDL2
	number of PXD pixels (L2)
	nPXD
	number of PXD pixels (L1+L2)
	nSVDL3
	number of L3 strips (u+v)
	nCDCin
	number of inner CDC wires
	CDCITIn
	integration time of the inner CDC layers [ms]
	nCDCout
	number of outer CDC wires
	CDCITout
	integration time of the outer CDC layers [ms]
	allh
	merged list of all histograms
	ntpl
	ntuple to store background levels

Detailed Description

Make validation histograms for BG overlay.
This validation assumes that all channels are active
(no masked strip, dead pixels, ...)

Definition at line 31 of file validateOverlay.py.

Member Function Documentation

event()

def event

(

self

)

 Event processor: fill histograms 

Definition at line 143 of file validateOverlay.py.

    def event(self):
        ''' Event processor: fill histograms '''
 
        for h in self.hist:
            digits = Belle2.PyStoreArray(h.GetName())
            if h.GetName() == 'ECLDigits':
                n = 0
                for digit in digits:
                    if digit.getAmp() > 500:  # roughly 25 MeV
                        n += 1
                h.Fill(n)
            else:
                h.Fill(digits.getEntries())
 
        # PXD
        pxdDigits = Belle2.PyStoreArray('PXDDigits')
        self.hist2[2].Fill(pxdDigits.getEntries()/self.nPXD * 100)
        nL1 = 0
        nL2 = 0
        for pxdDigit in pxdDigits:
            if pxdDigit.getSensorID().getLayerNumber() == 1:  # check L1/L2
                nL1 += 1
            else:
                nL2 += 1
        self.hist2[0].Fill(nL1/self.nPXDL1 * 100)
        self.hist2[1].Fill(nL2/self.nPXDL2 * 100)
 
        # SVD
        svdDigits = Belle2.PyStoreArray('SVDShaperDigitsZS5')
        nL3 = 0
        for svdDigit in svdDigits:
            if svdDigit.getSensorID().getLayerNumber() == 3:  # only check SVD L3
                nL3 += 1
        self.hist2[3].Fill(nL3/self.nSVDL3 * 100)
 
        # CDC
        cdcHits = Belle2.PyStoreArray('CDCHits')
        self.hist2[6].Fill(cdcHits.getEntries() / (self.nCDCin * self.CDCITIn + self.nCDCout * self.CDCITout))
        nCDC_in = 0
        nCDC_out = 0
        for cdcHit in cdcHits:
            if cdcHit.getICLayer() < 8:  # count wires of inner/outer layers
                nCDC_in += 1
            else:
                nCDC_out += 1
        self.hist2[4].Fill(nCDC_in/self.nCDCin / self.CDCITIn)
        self.hist2[5].Fill(nCDC_out/self.nCDCout / self.CDCITout)
 

initialize()

def initialize

(

self

)

 Initialize the Module: book histograms and set descriptions and checks

Definition at line 39 of file validateOverlay.py.

terminate()

def terminate

(

self

)

 Write histograms to file.

Definition at line 191 of file validateOverlay.py.

---

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

• background/validation/validateOverlay.py