WaveformAnalyzer Class Reference

Inheritance diagram for WaveformAnalyzer:

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

Public Attributes
	feProps
	output ntuple
	nWaveForms
	number of waveforms processed
	nWaveFormsOutOfOrder
	number of waveforms out of order
	f
	object of output file
	plotCounter
	counts how many plots were created.

Detailed Description

Analyzes waveform and FE data

Definition at line 120 of file TOPFE_qualityPlots.py.

Member Function Documentation

event()

def event

(

self

)

Event processor: get data and print to screen

Definition at line 149 of file TOPFE_qualityPlots.py.

    def event(self):
        '''
        Event processor: get data and print to screen
        '''
        evtMetaData = Belle2.PyStoreObj('EventMetaData')
        event = evtMetaData.obj().getEvent()
        run = evtMetaData.obj().getRun()
        waveforms = Belle2.PyStoreArray('TOPRawWaveforms')
        for waveform in waveforms:
            chan = waveform.getChannel()
            self.nWaveForms += 1
            # declare some counters for our ntuple
            nInFirstWindow, nNarrowPeaks, nHeight150, nOscillations = 0, 0, 0, 0
 
            # waveform.areWindowsInOrder is a bit too strict at the moment
            wins = np.array(waveform.getStorageWindows())
            if not np.all(wins[:-1] <= wins[1:]):  # change false for pics
                self.nWaveFormsOutOfOrder += 1
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "windowOrder")
                    self.plotCounter += 1
 
            wf = np.array(waveform.getWaveform())
            if False and args.plotWaveforms:
                wf_display(waveform, run, event, "general")
                self.plotCounter += 1
 
            # If TOPWaveformFeatureExtractor has been run then there are extra
            # TOPRawDigits related to this waveform
            rawDigits = waveform.getRelationsWith("TOPRawDigits")
            nTOPRawDigits = len(rawDigits)
 
            if False and args.plotWaveforms and nTOPRawDigits > 2:
                # There are at least 3 Top raw digits
                wf_display(waveform, run, event, "manyPeaks")
                self.plotCounter += 1
            if False and args.plotWaveforms and nTOPRawDigits > 3:
                # There are at least 4 TOPRawDigits
                wf_display(waveform, run, event, "tooManyPeaks")
                self.plotCounter += 1
            raw = rawDigits[0]
            fePeakHt = raw.getValuePeak()
            fePeakTDC = raw.getSamplePeak()
            fePeakWd = raw.getSampleFall() - raw.getSampleRise()
            fePeakIntegral = raw.getIntegral()
            # we are sorting the peaks. The FE should always point to the highest one
 
            # the channel we get from the waveform encodes the asic, carrier, and the actual channel
            if 0 < fePeakTDC < 65 and chan % 8 == 0:
                # To see if cal pulse is in first window
                nInFirstWindow += 1
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "calPuls_firstWin")
                    self.plotCounter += 1
 
            # there are strange bumps in the ADC distribtion of the cal channels
            if (140 < fePeakHt < 220) and chan % 8 == 0:
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "strangeADCBump_1")
                    self.plotCounter += 1
            if (300 < fePeakHt < 410) and chan % 8 == 0:
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "strangeADCBump_2")
                    self.plotCounter += 1
 
            fePeak1Ht = -1
            fePeak1TDC = -1
            fePeak1Wd = -1
            # fePeak1Integral = -1
            if nTOPRawDigits > 1:
                # if there is a second TOPRawDigit then get it for our ntuple
                fePeak1Ht = rawDigits[1].getValuePeak()
                fePeak1TDC = rawDigits[1].getSamplePeak()
                fePeak1Wd = rawDigits[1].getSampleFall() - rawDigits[1].getSampleRise()
                # fePeak1Integral = rawDigits[1].getIntegral()
 
            if nTOPRawDigits > 1 and fePeak1TDC < 64:
                # counting the times that the second digit is found in the first window
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "secondPeakInFirstWindow")
                    self.plotCounter += 1
 
            if fePeakWd < 5 or nTOPRawDigits > 1 and fePeak1Wd < 5:
                # counting thin peaks
                nNarrowPeaks += 1
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "thinpeak")
                    self.plotCounter += 1
 
            if 145 < fePeak1Ht < 155 and chan % 8 == 0:
                # Jan needs to write a comment about why we're doing this test
                nHeight150 += 1
                if False and args.plotWaveforms:
                    wf_display(waveform, run, event, "peak1Ht_is_150")
                    self.plotCounter += 1
 
            if nTOPRawDigits > 5 and (fePeakHt - fePeak1Ht) < 5:
                if np.mean(wf) < 10:
                    # want to count the portion of times we see this oscillation behavior
                    # guard the call to numpy.mean by the previous if statements
                    nOscillations += 1
                    if False and args.plotWaveforms:
                        wf_display(waveform, run, event, "oscillations")
                        self.plotCounter += 1
 
            self.feProps.Fill(
                event,
                run,
                fePeak1Ht,
                fePeak1TDC,
                fePeak1Wd,
                fePeakHt,
                fePeakTDC,
                fePeakWd,
                fePeakIntegral,
                nTOPRawDigits,
                chan,
                nNarrowPeaks,
                nInFirstWindow,
                nHeight150,
                nOscillations,
            )
 
            # only plot the first 10 figures.
            if args.plotWaveforms and self.plotCounter >= 10:
                evtMetaData.obj().setEndOfData()
 

initialize()

def initialize

(

self

)

 Initialize the Module: open the canvas. 

Definition at line 125 of file TOPFE_qualityPlots.py.

    def initialize(self):
        ''' Initialize the Module: open the canvas. '''
 
        
        self.feProps = TNtuple(
            "feProps",
            "feProps",
            "event:run"
            ":fePeak1Ht:fePeak1TDC:fePeak1Wd"
            ":fePeakHt:fePeakTDC:fePeakWd:fePeakIntegral"
            ":nTOPRawDigits:ch:nNarrowPeaks:nInFirstWindow:nHeight150:nOscillations")
 
        
        self.nWaveForms = 0
 
        
        self.nWaveFormsOutOfOrder = 0
 
        
        self.f = TFile.Open(args.output, "RECREATE")
 
        
        self.plotCounter = 0
 

terminate()

def terminate

(

self

)

End of run

Definition at line 276 of file TOPFE_qualityPlots.py.

    def terminate(self):
        '''End of run'''
        self.f.WriteTObject(self.feProps)
        self.f.Close()
        print("# Waveforms:", self.nWaveForms)
        print("out of order:", self.nWaveFormsOutOfOrder)
 
 
# Set the log level to show only error and fatal messages
# b2.set_log_level(b2.LogLevel.ERROR)
 
# Database

Member Data Documentation

f

f

object of output file

Definition at line 144 of file TOPFE_qualityPlots.py.

feProps

feProps

output ntuple

Definition at line 129 of file TOPFE_qualityPlots.py.

nWaveForms

nWaveForms

number of waveforms processed

Definition at line 138 of file TOPFE_qualityPlots.py.

nWaveFormsOutOfOrder

nWaveFormsOutOfOrder

number of waveforms out of order

Definition at line 141 of file TOPFE_qualityPlots.py.

plotCounter

plotCounter

counts how many plots were created.

Stop after 10

Definition at line 147 of file TOPFE_qualityPlots.py.

---

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

• top/examples/TOPFE_qualityPlots.py