Histogrammer Class Reference

Inheritance 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
	nTOPev
	numnber of events
	minTOPTime
	min TOP hit time
	maxTOPTime
	max TOP hit time
	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 171 of file validateOverlay.py.

    def event(self):
        ''' Event processor: fill histograms '''
 
        for h in self.hist:
            digits = Belle2.PyStoreArray(h.GetName())
            if h.GetName() == 'CDCHits':
                continue
            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')
        nCDC_in = 0
        nCDC_out = 0
        for cdcHit in cdcHits:  # count wires of inner/outer layers
            if cdcHit.getISuperLayer() == 0 and cdcHit.getADCCount() > 15:
                nCDC_in += 1
            if cdcHit.getISuperLayer() != 0 and cdcHit.getADCCount() > 18:
                nCDC_out += 1
        self.hist2[4].Fill(nCDC_in/self.nCDCin / self.CDCITIn)
        self.hist2[5].Fill(nCDC_out/self.nCDCout / self.CDCITout)
        self.hist2[6].Fill((nCDC_in+nCDC_out) / (self.nCDCin * self.CDCITIn + self.nCDCout * self.CDCITout))
        self.hist[3].Fill(nCDC_in+nCDC_out)
 
        # TOP
        self.nTOPev += 1
        for TOPDigit in Belle2.PyStoreArray('TOPDigits'):
            if TOPDigit.getHitQuality() != Belle2.TOPDigit.c_Good:
                continue
            self.hist2[7].Fill(TOPDigit.getTime())
            if self.minTOPTime:
                self.minTOPTime = min(self.minTOPTime, TOPDigit.getTime())
                self.maxTOPTime = max(self.maxTOPTime, TOPDigit.getTime())
            else:
                self.minTOPTime = TOPDigit.getTime()
                self.maxTOPTime = TOPDigit.getTime()
 

initialize()

def initialize

(

self

)

 Initialize the Module: book histograms and set descriptions and checks

Definition at line 39 of file validateOverlay.py.

    def initialize(self):
        ''' Initialize the Module: book histograms and set descriptions and checks'''
 
        
        self.hist = []
        self.hist.append(TH1F('PXDDigits', 'PXDDigits (no data reduction)',
                              200, 0, 100000))
        self.hist.append(TH1F('SVDShaperDigits', 'SVDShaperDigits', 100, 0, 10000))
        self.hist.append(TH1F('SVDShaperDigitsZS5', 'ZS5 SVDShaperDigits', 100, 0, 8000))
        self.hist.append(TH1F('CDCHits', 'CDCHits above threshold', 100, 0, 6000))
        self.hist.append(TH1F('TOPDigits', 'TOPDigits', 100, 0, 2000))
        self.hist.append(TH1F('ARICHDigits', 'ARICHDigits', 100, 0, 300))
        self.hist.append(TH1F('ECLDigits', 'ECLDigits, m_Amp > 500 (roughly 25 MeV)',
                              100, 0, 400))
        self.hist.append(TH1F('KLMDigits', 'KLMDigits', 150, 0, 150))
 
        
        self.hist2 = []
        self.hist2.append(TH1F('PXDL1Occupancy', 'PXD L1 Occupancy', 100, 0, 4))
        self.hist2.append(TH1F('PXDL2Occupancy', 'PXD L2 Occupancy', 100, 0, 4))
        self.hist2.append(TH1F('PXDOccupancy', 'PXD Occupancy', 100, 0, 4))
        self.hist2.append(TH1F('SVDOccupancy', 'SVD L3 Occupancy', 100, 0, 15))
        self.hist2.append(TH1F('CDCInnerHitRate', 'CDC Hit Rate per wire for the Inner layers', 100, 0, 1000))
        self.hist2.append(TH1F('CDCOuterHitRate', 'CDC Hit Rate per wire for the Outer layers', 100, 0, 1000))
        self.hist2.append(TH1F('CDCHitRate', 'CDC Hit Rate per wire', 100, 0, 1000))
        self.hist2.append(TH1F("TOPHitTime", "TOP Hit Time Distribution; time [ns]", 300, -150, 150))
        
        self.nTOPev = 0
        
        self.minTOPTime = None
        
        self.maxTOPTime = None
 
        
        self.nPXDL1 = 250*768*2 * 8
        
        self.nPXDL2 = 250*768*2 * 12
        if len(sys.argv) == 3 and int(sys.argv[2]) == 1003:
            self.nPXDL2 = 250*768*2 * 2
        
        self.nPXD = self.nPXDL1 + self.nPXDL2
 
        
        self.nSVDL3 = 768 * 7 * 2 * 2
 
        
        self.nCDCin = 160 * 8
        
        self.CDCITIn = 408.7 * 1e-6
        
        self.nCDCout = 160*6 + 192*6 + 224*6 + 256*6 + 288*6 + 320*6 + 352*6 + 384 * 6
        
        self.CDCITout = 754.6 * 1e-6
 
        
        self.allh = self.hist+self.hist2
 
        for h in self.allh:
            if h.GetName() == 'PXDOccupancy':
                h.SetXTitle('PXD occupancy [%]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'PXD L1+L2 Occupancy per event - no data reduction (with BG overlay only and no event generator)')
            elif h.GetName() == 'PXDL1Occupancy':
                h.SetXTitle('PXD L1 occupancy [%]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'PXD L1 Occupancy per event - no data reduction (with BG overlay only and no event generator)')
            elif h.GetName() == 'PXDL2Occupancy':
                h.SetXTitle('PXD L2 occupancy [%]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'PXD L2 Occupancy per event - no data reduction (with BG overlay only and no event generator)')
            elif h.GetName() == 'SVDOccupancy':
                h.SetXTitle('SVD L3 ZS5 occupancy - average u,v [%]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'SVD ZS5 L3 Occupancy per event average u,v (with BG overlay only and no event generator)')
            elif h.GetName() == 'CDCHitRate':
                h.SetXTitle('CDC Hit Rate [kHz/wire]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'CDC Hit Rate, averaged inner/outer wires (with BG overlay only and no event generator)')
            elif h.GetName() == 'CDCInnerHitRate':
                h.SetXTitle('CDC Hit Rate for Inner Layers [kHz/wire]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'CDC Hit Rate per wire, for the Inner layers (with BG overlay only and no event generator)')
            elif h.GetName() == 'CDCOuterHitRate':
                h.SetXTitle('CDC Hit Rate for Outer Layers [kHz/wire]')
                option = 'shifter'
                descr = TNamed(
                    'Description',
                    'CDC Hit Rate per wire, for the Outer layers (with BG overlay only and no event generator)')
            else:
                h.SetXTitle('number of digits in event')
                option = 'shifter'
                descr = TNamed('Description', 'Number of background ' + h.GetName() +
                               ' per event (with BG overlay only and no event generator)')
 
            if h.GetName() == 'PXDDigits':
                option = 'expert'
            elif h.GetName() == 'SVDShaperDigits':
                option = 'expert'
            elif h.GetName() == 'SVDShaperDigitsZS5':
                option = 'expert'
            elif h.GetName() == 'CDCHits':
                option = 'expert'
            elif h.GetName() == 'TOPHitTime':
                option = 'expert'
 
            h.SetYTitle('entries per bin')
            h.GetListOfFunctions().Add(descr)
            check = TNamed('Check', 'Distribution must agree with its reference')
            h.GetListOfFunctions().Add(check)
            contact = TNamed('Contact', 'marko.staric@ijs.si')
            h.GetListOfFunctions().Add(contact)
            options = TNamed('MetaOptions', option)
            h.GetListOfFunctions().Add(options)
 
            ntplvar = "PXDL1_occupancy:PXDL2_occupancy:PXD_occupancy:SVDL3ZS5_occupancy"
            ntplvar += ":CDCInner_hitRate:CDCOuter_hitRate:CDC_hitRate:TOP_hitRate"
 
            
            self.ntpl = TNtuple("ntpl", "Average Background Levels", (ntplvar))
 

terminate()

def terminate

(

self

)

 Write histograms to file.

Definition at line 235 of file validateOverlay.py.

    def terminate(self):
        """ Write histograms to file."""
 
        TOPbg = self.hist2[7].Integral(101, 150)/self.nTOPev/50e-3/16/32
 
        # ntpl content "PXDL1:PXDL2:PXD:SVDL3ZS5:CDCInner:CDCOuter:CDC:TOP");
        values = [self.hist2[0].GetMean(), self.hist2[1].GetMean(), self.hist2[2].GetMean(),
                  self.hist2[3].GetMean(), self.hist2[4].GetMean(), self.hist2[5].GetMean(),
                  self.hist2[6].GetMean(), float(TOPbg)]
 
        self.ntpl.Fill(array.array("f", values),)
        tfile = TFile('overlayPlots.root', 'recreate')
        self.ntpl.Write()
        for h in self.allh:
            h.Write()
        tfile.Close()
 
 

Member Data Documentation

allh

allh

merged list of all histograms

Definition at line 94 of file validateOverlay.py.

CDCITIn

CDCITIn

integration time of the inner CDC layers [ms]

Definition at line 87 of file validateOverlay.py.

CDCITout

CDCITout

integration time of the outer CDC layers [ms]

Definition at line 91 of file validateOverlay.py.

hist

hist

list of digits histograms

Definition at line 43 of file validateOverlay.py.

hist2

hist2

list of occupancy-like histograms

Definition at line 56 of file validateOverlay.py.

maxTOPTime

maxTOPTime

max TOP hit time

Definition at line 70 of file validateOverlay.py.

minTOPTime

minTOPTime

min TOP hit time

Definition at line 68 of file validateOverlay.py.

nCDCin

nCDCin

number of inner CDC wires

Definition at line 85 of file validateOverlay.py.

nCDCout

nCDCout

number of outer CDC wires

Definition at line 89 of file validateOverlay.py.

nPXD

nPXD

number of PXD pixels (L1+L2)

Definition at line 79 of file validateOverlay.py.

nPXDL1

nPXDL1

number of PXD pixels (L1)

Definition at line 73 of file validateOverlay.py.

nPXDL2

nPXDL2

number of PXD pixels (L2)

Definition at line 75 of file validateOverlay.py.

nSVDL3

nSVDL3

number of L3 strips (u+v)

Definition at line 82 of file validateOverlay.py.

nTOPev

nTOPev

numnber of events

Definition at line 66 of file validateOverlay.py.

ntpl

ntpl

ntuple to store background levels

Definition at line 169 of file validateOverlay.py.

---

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

• background/validation/validateOverlay.py