Public Member Functions
	__init__ (self, exp, run, histName, pdfName)

	makeGraph (self, x, y)

	makeText (self, x, y, s)

	initialize (self)

	terminate (self)

	beginRun (self)

	endRun (self)

	event (self, eventHits, tt_ctime, raw_time)

Public Attributes
	exp = exp
	internal copy of experiment number

	run = run
	internal copy of run number

	histName = histName
	internal copy of the pathname of the output histogram ROOT file

	pdfName = pdfName
	internal copy of the pathname of the output histogram PDF file

	histogramFile = ROOT.TFile.Open(self.histName, "RECREATE")
	Output ROOT TFile that will contain the histograms/scatterplots.

	hist_ttc_trigtime
	histogram of the tt_ctime relative to triggertime

	hist_rawKLMlane
	histogram of the hit's lane

	hist_rawKLMsizeMultihit
	histogram of number of hits, including multiple entries on one readout channel

	hist_rawKLMsize
	histogram of number of hits, at most one entry per readout channel

	hist_PerChannelMultiplicity
	scatterplot of multiplicity of entries in one readout channel vs lane/axis

	hist_RPCLaneAxisOccupancy
	scatterplot of number of mapped RPC hits by lane/axis vs sector, at most one entry per readout channel

	hist_ScintLaneAxisOccupancy
	scatterplot of number of mapped scint hits by lane/axis vs sector, at most one entry per readout channel

list	hist_ChannelOccupancy = [0, 0]
	scatterplots of channel occupancy (1 hit per readout channel) for each axis

list	hist_ChannelOccupancyAL
	histograms of channel occupancy (1 hit per readout channel), indexed by axis/lane

	hist_RPCTimeLowBitsBySector
	scatterplot of RPC TDC low-order bits vs sector (should be 0 since granularity is 4 ns)

	hist_RPCTime
	histogram of RPC TDC value relative to event's REVO9 trigger time in last word of event

	hist_RPCTime2
	histogram of RPC TDC value relative to event's ctime in event header

list	hist_RPCTimePerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
	histograms of RPC TDC value relative to event's trigger time for axis 0, indexed by lane

list	hist_RPCTimePerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
	histograms of RPC TDC value relative to event's trigger time for axis 1, indexed by lane

list	hist_RPCTime2PerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
	histograms of RPC TDC value relative to event's ctime for axis 0, indexed by lane

list	hist_RPCTime2PerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
	histograms of RPC TDC value relative to event's ctime for axis 1, indexed by lane

	hist_RPCTdcRange = ROOT.TH1F('RPCTdcRange', expRun + 'RPC TDC-range in event;TDCMax - TDCMin (ns)', 512, -0.5, 2047.5)
	histogram of RPC TDC range in event

	hist_RPCRevotimeRange
	histogram of RPC REVO9 range in event

	hist_revotimeRPCtdc
	scatterplot of RPC REVO9 range vs TDC value in event

	hist_revotimeRPCtdc2
	scatterplot of RPC REVO9 range vs TDC value corrected for DC-processing delay in event

	hist_jRPCtdc
	scatterplot of RPC calibrated time vs hit's index

	hist_jRPCtdc2
	scatterplot of RPC calibrated time vs hit's index, corrected for DC-processing delay

	hist_ScintTimeLowBitsBySector
	scatterplot of scint TDC low-order bits vs sector

	hist_ScintTime
	histogram of scint TDC value relative to event's trigger time

	hist_ScintCtime
	histogram of scint CTIME value relative to event's ctime

	hist_ScintCtime0
	histogram of scint CTIME value relative to event's trigger time

	hist_ScintCtimeRange
	histogram of scint CTIME range in event

Static Public Attributes
int	BKLM_ID = 0x07000000
	COPPER base identifier for BKLM readout.

int	EKLM_ID = 0x08000000
	COPPER base identifier for EKLM readout.

int	BKLM_STRIP_BIT = 0
	bit position for strip-1 [0..47]

int	BKLM_PLANE_BIT = 6
	bit position for plane-1 [0..1]; 0 is inner-plane

int	BKLM_LAYER_BIT = 7
	bit position for layer-1 [0..14]; 0 is innermost

int	BKLM_SECTOR_BIT = 11
	bit position for sector-1 [0..7]; 0 is on the +x axis and 2 is on the +y axis

int	BKLM_SECTION_BIT = 14
	bit position for section [0..1]; forward is 0

int	BKLM_MAXSTRIP_BIT = 15
	bit position for maxStrip-1 [0..47]

int	BKLM_STRIP_MASK = 0x3f
	bit mask for strip-1 [0..47]

tuple	BKLM_PLANE_MASK = (1 << BKLM_PLANE_BIT)
	bit mask for plane-1 [0..1]; 0 is inner-plane

tuple	BKLM_LAYER_MASK = (15 << BKLM_LAYER_BIT)
	bit mask for layer-1 [0..15]; 0 is innermost and 14 is outermost

tuple	BKLM_SECTOR_MASK = (7 << BKLM_SECTOR_BIT)
	bit mask for sector-1 [0..7]; 0 is on the +x axis and 2 is on the +y axis

tuple	BKLM_SECTION_MASK = (1 << BKLM_SECTION_BIT)
	bit mask for section [0..1]; forward is 0

tuple	BKLM_MAXSTRIP_MASK = (63 << BKLM_MAXSTRIP_BIT)
	bit mask for maxStrip-1 [0..47]

tuple	BKLM_MODULEID_MASK = (BKLM_SECTION_MASK \| BKLM_SECTOR_MASK \| BKLM_LAYER_MASK)
	bit mask for unique module identifier (end, sector, layer)

Detailed Description

Fill BKLM histograms of values from RawKLMs, KLMDigits, BKLMHit1ds, and BKLMHit2ds;
(optionally) draw event displays from these data-objects.

Definition at line 18 of file EventInspectorPocketDAQ.py.

Constructor & Destructor Documentation

◆ init()

__init__	(	self,
		exp,
		run,
		histName,
		pdfName )

Constructor

Arguments:
    exp (str): formatted experiment number
    run (str): formatter run number
    histName (str): path name of the output histogram ROOT file
    pdfName (str): path name of the output histogram PDF file
    eventPdfName (str): path name of the output event-display PDF file
    maxDisplays (int): max # of events displays to write
    minRPCHits (int): min # of RPC BKLMHit2ds in any sector for event display

Definition at line 53 of file EventInspectorPocketDAQ.py.

    def __init__(self, exp, run, histName, pdfName):
        """Constructor
 
        Arguments:
            exp (str): formatted experiment number
            run (str): formatter run number
            histName (str): path name of the output histogram ROOT file
            pdfName (str): path name of the output histogram PDF file
            eventPdfName (str): path name of the output event-display PDF file
            maxDisplays (int): max # of events displays to write
            minRPCHits (int): min # of RPC BKLMHit2ds in any sector for event display
        """
        super().__init__()
        
        self.exp = exp
        
        self.run = run
        
        self.histName = histName
        
        self.pdfName = pdfName
 

Member Function Documentation

◆ beginRun()

beginRun ( self )

Handle begin of run: print diagnostic message

Definition at line 341 of file EventInspectorPocketDAQ.py.

    def beginRun(self):
        """Handle begin of run: print diagnostic message"""
        print('beginRun', self.run)
 

◆ endRun()

endRun ( self )

Handle end of run: print diagnostic message

Definition at line 345 of file EventInspectorPocketDAQ.py.

    def endRun(self):
        """Handle end of run: print diagnostic message"""
        print('endRun', self.run)
 

◆ event()

event	(	self,
		eventHits,
		tt_ctime,
		raw_time )

Process one event: fill histograms

Definition at line 349 of file EventInspectorPocketDAQ.py.

    def event(self, eventHits, tt_ctime, raw_time):
        """Process one event: fill histograms"""
 
        sectorFB = 0
        n = len(eventHits)
        countAllMultihit = 2 * n + 1 if (n > 0) else 0
        countAll = 1
        channelMultiplicity = {}
        minCtime = 99999
        minRPCTdc = 99999
        maxRPCTdc = 0
        minRPCCtime = 99999
        maxRPCCtime = 0
        minScintCtime = 99999
        maxScintCtime = 0
        self.hist_ttc_trigtime.Fill((tt_ctime - raw_time) & 0x0fff)
        for j in range(0, n):
            items = eventHits[j]
            lane = items[0]
            channel = items[1]
            axis = items[2]
            ctime = items[3]  # this is revo9 time for RPCs, ctime for scintillators
            flag = 1 if (lane > 2) else 2
            isRPC = (flag == 1)
            isScint = (flag == 2)  # noqa (F841) kept for completeness
            laneAxisChannel = (((lane << 1) + axis) << 7) + channel
            if laneAxisChannel not in channelMultiplicity:
                countAll = countAll + 2
                channelMultiplicity[laneAxisChannel] = 0
            channelMultiplicity[laneAxisChannel] = channelMultiplicity[laneAxisChannel] + 1
            if ctime < minCtime:
                minCtime = ctime
            self.hist_rawKLMlane.Fill(lane)
        for j in range(0, n):
            items = eventHits[j]
            lane = items[0]
            channel = items[1]
            axis = items[2]
            ctime = items[3]  # this is revo9 time for RPCs, ctime for scintillators
            tdc = items[4]
            charge = items[5]  # noqa (F841) kept for completeness
            flag = 1 if (lane > 2) else 2
            isRPC = (flag == 1)
            isScint = (flag == 2)  # noqa (F841) kept for completeness
            laneAxisChannel = (((lane << 1) + axis) << 7) + channel
            laneAxis = axis if ((lane < 1) or (lane > 20)) else ((lane << 1) + axis)
            if laneAxisChannel in channelMultiplicity:
                if channelMultiplicity[laneAxisChannel] > 1:
                    self.hist_PerChannelMultiplicity.Fill(channelMultiplicity[laneAxisChannel], laneAxis)
                # DIVOT del channelMultiplicity[laneAxisChannel] # consider only first hit in the channel/axis/lane of this dc
                t = (tdc - raw_time) & 0x03ff  # in ns, range is 0..1023
                t2 = (tdc - tt_ctime) & 0x03ff  # in ns, range is 0..1023
                ct = ((ctime << 3) - tt_ctime) & 0x3ff
                if isRPC:
                    if tdc < minRPCTdc:
                        minRPCTdc = tdc
                    if tdc > maxRPCTdc:
                        maxRPCTdc = tdc
                    if ctime < minRPCCtime:
                        minRPCCtime = ctime
                    if ctime > maxRPCCtime:
                        maxRPCCtime = ctime
                    self.hist_RPCTimeLowBitsBySector.Fill(sectorFB, (tdc & 3))
                    self.hist_RPCLaneAxisOccupancy.Fill(sectorFB, laneAxis)
                    self.hist_RPCTime.Fill(t)
                    self.hist_RPCTime2.Fill(t2)
                    if axis == 0:
                        self.hist_RPCTimePerLayerA0[lane].Fill(t)
                        self.hist_RPCTime2PerLayerA0[lane].Fill(t2)
                    else:
                        self.hist_RPCTimePerLayerA1[lane].Fill(t)
                        self.hist_RPCTime2PerLayerA1[lane].Fill(t2)
                    if n > 0:
                        t0j = 0.75 * j
                        self.hist_revotimeRPCtdc.Fill(t, ctime - minCtime)
                        self.hist_revotimeRPCtdc2.Fill(t - t0j, ctime - minCtime)
                        self.hist_jRPCtdc.Fill(t, j)
                        self.hist_jRPCtdc2.Fill(t - t0j, j)
                else:
                    if ctime < minScintCtime:
                        minScintCtime = ctime
                    if ctime > maxScintCtime:
                        maxScintCtime = ctime
                    self.hist_ScintTimeLowBitsBySector.Fill(sectorFB, (tdc & 3))
                    self.hist_ScintLaneAxisOccupancy.Fill(sectorFB, laneAxis)
                    self.hist_ScintTime.Fill(t)
                    self.hist_ScintCtime.Fill(ct)
                    self.hist_ScintCtime0.Fill(((ctime << 3) - raw_time) & 0x3ff)
                self.hist_ChannelOccupancy[axis].Fill(lane, channel)
                self.hist_ChannelOccupancyAL[lane][axis].Fill(channel)
        if n > 1:
            if maxRPCTdc > 0:
                self.hist_RPCTdcRange.Fill(maxRPCTdc - minRPCTdc)
            if maxRPCCtime > 0:
                self.hist_RPCRevotimeRange.Fill((maxRPCCtime - minRPCCtime) << 3)
            if maxScintCtime > 0:
                self.hist_ScintCtimeRange.Fill((maxScintCtime - minScintCtime) << 3)
        self.hist_rawKLMsizeMultihit.Fill(countAllMultihit)
        self.hist_rawKLMsize.Fill(countAll)

◆ initialize()

initialize ( self )

Handle job initialization: create histograms

Definition at line 105 of file EventInspectorPocketDAQ.py.

    def initialize(self):
        """Handle job initialization: create histograms"""
 
        expRun = f'e{int(self.exp):02d}r{int(self.run)}: '
 
        
        self.histogramFile = ROOT.TFile.Open(self.histName, "RECREATE")
 
        # create the rawKLM histograms
 
        
        self.hist_ttc_trigtime = ROOT.TH1F('ttc_trigtime',
                                           expRun + 'tt_ctime relative to triggertime;' +
                                           'tt_ctime - triggertime (ns)',
                                           256, -0.5, 2047.5)
        
        self.hist_rawKLMlane = ROOT.TH1F('rawKLMlane',
                                         expRun + 'RawKLM lane;Lane (scint: 1..7, RPC: 8..20)',
                                         21, -0.5, 20.5)
        
        self.hist_rawKLMsizeMultihit = ROOT.TH1F('rawKLMsizeMultihit',
                                                 expRun + 'RawKLM word count (N/channel)',
                                                 200, -0.5, 199.5)
        
        self.hist_rawKLMsize = ROOT.TH1F('rawKLMsize',
                                         expRun + 'RawKLM word count (1/channel)',
                                         200, -0.5, 199.5)
        
        self.hist_PerChannelMultiplicity = ROOT.TH2F(
            'PerChannelMultiplicity',
            expRun + 'Per-channel multiplicity (N/channel > 1);' +
                     'Per-channel multiplicity;' +
                     '(Lane #) * 2 + (Axis #)',
            30, -0.5, 29.5, 42, -0.5, 41.5)
        
        self.hist_RPCLaneAxisOccupancy = ROOT.TH2F(
            'RPCLaneAxisOccupancy',
            expRun + 'Lane/axis occupancy of RPC channels (1/channel);' +
                     'Sector # (always 0);' +
                     '(Lane #) * 2 + (Axis #)',
            3, -1.5, 1.5, 42, -0.5, 41.5)
        
        self.hist_ScintLaneAxisOccupancy = ROOT.TH2F(
            'ScintLaneAxisOccupancy',
            expRun + 'Lane/axis occupancy of scint channels (1/channel);' +
                     'Sector # (always 0);' +
                     '(Lane #) * 2 + (Axis #)',
            3, -1.5, 1.5, 42, -0.5, 41.5)
        
        self.hist_ChannelOccupancy = [0, 0]
        self.hist_ChannelOccupancy[0] = ROOT.TH2F('ChannelOccupancy_a0',
                                                  expRun + 'Channel occupancy for axis 0;lane;channel',
                                                  42, -0.25, 20.75, 128, -0.25, 63.75)
        self.hist_ChannelOccupancy[1] = ROOT.TH2F('ChannelOccupancy_a1',
                                                  expRun + 'Channel occupancy for axis 1;lane;channel',
                                                  42, -0.25, 20.75, 128, -0.25, 63.75)
        
        self.hist_ChannelOccupancyAL = [
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0],
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0],
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0]]
        for lane in range(0, 21):
            nChannels = 64 if (lane > 2) else 128
            label = f'ChannelOccupancy_A0L{lane}'
            title = f'{expRun}Channel occupancy for axis 0 lane {lane};channel'
            self.hist_ChannelOccupancyAL[lane][0] = ROOT.TH1F(label, title, nChannels, -0.5, nChannels - 0.5)
            label = f'ChannelOccupancy_A1L{lane}'
            title = f'{expRun}Channel occupancy for axis 1 lane {lane};channel'
            self.hist_ChannelOccupancyAL[lane][1] = ROOT.TH1F(label, title, nChannels, -0.5, nChannels - 0.5)
        
        self.hist_RPCTimeLowBitsBySector = ROOT.TH2F('RPCTimeLowBitsBySector',
                                                     expRun + 'RPC TDC lowest-order bits;' +
                                                     'Sector # (always 0);' +
                                                     'TDC % 4 (ns) [should be 0]',
                                                     3, -1.5, 1.5, 8, -0.25, 3.75)
        
        self.hist_RPCTime = ROOT.TH1F('RPCTime',
                                      expRun + 'RPC tdc relative to event trigtime;tdc - triggerTime (ns)',
                                      256, -0.5, 1023.5)
        
        self.hist_RPCTime2 = ROOT.TH1F('RPCTime_Ctime',
                                       expRun + 'RPC tdc relative to event ctime;tdc - trigCtime (ns)',
                                       256, -0.5, 1023.5)
        
        self.hist_RPCTimePerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
        for lane in range(0, 21):
            label = f'RPCTimeA0L{lane:02d}'
            title = f'{expRun}RPC axis 0 lane {lane} time relative to trigtime;t - triggerTime (ns)'
            self.hist_RPCTimePerLayerA0[lane] = ROOT.TH1F(label, title, 256, -0.5, 1023.5)
        
        self.hist_RPCTimePerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
        for lane in range(0, 21):
            label = f'RPCTimeA1L{lane:02d}'
            title = f'{expRun}RPC axis 1 lane {lane} time relative to trigtime;t - triggerTime (ns)'
            self.hist_RPCTimePerLayerA1[lane] = ROOT.TH1F(label, title, 256, -0.5, 1023.5)
        
        self.hist_RPCTime2PerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
        for lane in range(0, 21):
            label = f'RPCTime2A0L{lane:02d}'
            title = f'{expRun}RPC axis 0 lane {lane} time relative to trigCtime;t - trigCtime (ns)'
            self.hist_RPCTime2PerLayerA0[lane] = ROOT.TH1F(label, title, 256, -0.5, 1023.5)
        
        self.hist_RPCTime2PerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
        for lane in range(0, 21):
            label = f'RPCTime2A1L{lane:02d}'
            title = f'{expRun}RPC axis 1 lane {lane} time relative to trigCtime;t - trigCtime (ns)'
            self.hist_RPCTime2PerLayerA1[lane] = ROOT.TH1F(label, title, 256, -0.5, 1023.5)
        
        self.hist_RPCTdcRange = ROOT.TH1F('RPCTdcRange', expRun + 'RPC TDC-range in event;TDCMax - TDCMin (ns)', 512, -0.5, 2047.5)
        
        self.hist_RPCRevotimeRange = ROOT.TH1F('RPCRevotimeRange',
                                               expRun + 'RPC revotime-range in event;revotimeMax - revotimeMin (ns)',
                                               128, -0.5, 8191.5)
        
        self.hist_revotimeRPCtdc = ROOT.TH2F('revotimeRPCtdc',
                                             expRun + 'RPC TDC vs revotime;tdc (ns);revotime - minRevotime',
                                             64, 767.5, 1023.5, 10, -0.5, 159.5)
        
        self.hist_revotimeRPCtdc2 = ROOT.TH2F(
            'revotimeRPCtdc2',
            expRun + 'RPC TDC vs revotime;' +
            'tdc - dt(index) (ns);' +
            'revotime - minRevotime',
            64, 767.5, 1023.5, 10, -0.5, 159.5)
        
        self.hist_jRPCtdc = ROOT.TH2F('jRPCtdc',
                                      expRun + 'RPC TDC vs hit index;tdc (ns);Hit index',
                                      64, 767.5, 1023.5, 60, -0.5, 59.5)
        
        self.hist_jRPCtdc2 = ROOT.TH2F('jRPCtdc2',
                                       expRun + 'RPC TDC vs hit index;tdc - dt(index) (ns);Hit index',
                                       64, 767.5, 1023.5, 60, -0.5, 59.5)
        
        self.hist_ScintTimeLowBitsBySector = ROOT.TH2F(
            'ScintTimeLowBitsBySector',
            expRun + 'Scint TDC lowest-order bits;' +
            'Sector # (always 0);' +
            'TDC % 4 (ns)',
            3, -1.5, 1.5, 8, -0.25, 3.75)
        
        self.hist_ScintTime = ROOT.TH1F('ScintTime',
                                        expRun + 'Scint tdc distribution;tdc - triggerTime (ns)',
                                        256, -0.5, 1023.5)
        
        self.hist_ScintCtime = ROOT.TH1F('ScintCtime',
                                         expRun + 'Scint ctime distribution;ctime - triggerCtime (ns)',
                                         32, -0.5, 1023.5)
        
        self.hist_ScintCtime0 = ROOT.TH1F('ScintCtime0',
                                          expRun + 'Scint ctime distribution;ctime - triggerTime (ns)',
                                          32, -0.5, 1023.5)
        
        self.hist_ScintCtimeRange = ROOT.TH1F(
            'ScintCtimeRange', expRun + 'Scint ctime-range in event;ctimeMax - ctimeMin (ns)', 128, -0.5, 8191.5)
 

◆ makeGraph()

makeGraph	(	self,
		x,
		y )

Create and return a ROOT TGraph

Arguments:
  x[] (real): x coordinates
  y[] (real): y coordinates

Definition at line 75 of file EventInspectorPocketDAQ.py.

    def makeGraph(self, x, y):
        """Create and return a ROOT TGraph
 
        Arguments:
          x[] (real): x coordinates
          y[] (real): y coordinates
        """
        graph = ROOT.TGraph()
        for i in range(0, len(x)):
            graph.SetPoint(i, x[i], y[i])
        graph.SetLineColor(2)
        graph.SetLineWidth(1)
        return graph
 

◆ makeText()

makeText	(	self,
		x,
		y,
		s )

Create and return a ROOT TLatex with the following properties:
size = 0.04, color = red, alignment = middle centre, angle = 90 degrees

Arguments:
  x (real): x coordinate
  y (real): y coordinate
  s (str):  character string

Definition at line 89 of file EventInspectorPocketDAQ.py.

    def makeText(self, x, y, s):
        """Create and return a ROOT TLatex with the following properties:
        size = 0.04, color = red, alignment = middle centre, angle = 90 degrees
 
        Arguments:
          x (real): x coordinate
          y (real): y coordinate
          s (str):  character string
        """
        text = ROOT.TLatex(x, y, s)
        text.SetTextSize(0.04)
        text.SetTextColor(2)
        text.SetTextAlign(22)
        text.SetTextAngle(90)
        return text
 

◆ terminate()

terminate ( self )

Handle job termination: draw histograms, close output files

Definition at line 260 of file EventInspectorPocketDAQ.py.

    def terminate(self):
        """Handle job termination: draw histograms, close output files"""
 
        canvas = ROOT.TCanvas("canvas", self.pdfName, 1600, 1600)
        title = f'{self.pdfName}['
        canvas.SaveAs(title)
        canvas.Clear()
        canvas.GetPad(0).SetGrid(1, 1)
        canvas.GetPad(0).Update()
        self.hist_rawKLMlane.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_rawKLMlane.GetName()}")
        self.hist_rawKLMsizeMultihit.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_rawKLMsizeMultihit.GetName()}")
        self.hist_rawKLMsize.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_rawKLMsize.GetName()}")
        self.hist_PerChannelMultiplicity.Draw("box")
        canvas.Print(self.pdfName, f"Title:{self.hist_PerChannelMultiplicity.GetName()}")
        canvas.Clear()
        canvas.Divide(2, 1)
        canvas.GetPad(0).SetGrid(1, 1)
        canvas.GetPad(1).SetGrid(1, 1)
        canvas.GetPad(2).SetGrid(1, 1)
        for sectorFB in range(0, 1):
            canvas.cd(1)
            self.hist_ChannelOccupancy[0].Draw("colz")
            canvas.cd(2)
            self.hist_ChannelOccupancy[1].Draw("colz")
            canvas.Print(self.pdfName, f"Title:{self.hist_ChannelOccupancy[0].GetName()}")
        for lane in range(0, 21):
            n0 = self.hist_ChannelOccupancyAL[lane][0].GetEntries()
            n1 = self.hist_ChannelOccupancyAL[lane][1].GetEntries()
            if n0 + n1 > 0:
                canvas.cd(1)
                self.hist_ChannelOccupancyAL[lane][0].Draw()
                canvas.cd(2)
                self.hist_ChannelOccupancyAL[lane][1].Draw()
                canvas.Print(self.pdfName, f"Title:{self.hist_ChannelOccupancyAL[lane][0].GetName()}")
        canvas.Clear()
        canvas.Divide(1, 1)
        self.hist_ttc_trigtime.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_ttc_trigtime.GetName()}")
        self.hist_RPCTimeLowBitsBySector.Draw("box")
        canvas.Print(self.pdfName, f"Title:{self.hist_RPCTimeLowBitsBySector.GetName()}")
        self.hist_RPCTime.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_RPCTime.GetName()}")
        for lane in range(0, 21):
            if self.hist_RPCTimePerLayerA0[lane].GetEntries() > 0:
                self.hist_RPCTimePerLayerA0[lane].Draw()
                canvas.Print(self.pdfName, f"Title:{self.hist_RPCTimePerLayerA0[lane].GetName()}")
        for lane in range(0, 21):
            if self.hist_RPCTimePerLayerA1[lane].GetEntries() > 0:
                self.hist_RPCTimePerLayerA1[lane].Draw()
                canvas.Print(self.pdfName, f"Title:{self.hist_RPCTimePerLayerA1[lane].GetName()}")
        self.hist_RPCTdcRange.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_RPCTdcRange.GetName()}")
        self.hist_RPCRevotimeRange.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_RPCRevotimeRange.GetName()}")
        self.hist_revotimeRPCtdc.Draw("colz")
        canvas.Print(self.pdfName, f"Title:{self.hist_revotimeRPCtdc.GetName()}")
        # self.hist_revotimeRPCtdc2.Draw("colz")
        # canvas.Print(self.pdfName, "Title:{0}".format(self.hist_revotimeRPCtdc2.GetName()))
        self.hist_jRPCtdc.Draw("colz")
        canvas.Print(self.pdfName, f"Title:{self.hist_jRPCtdc.GetName()}")
        # self.hist_jRPCtdc2.Draw("colz")
        # canvas.Print(self.pdfName, "Title:{0}".format(self.hist_jRPCtdc2.GetName()))
        self.hist_ScintTimeLowBitsBySector.Draw("box")
        canvas.Print(self.pdfName, f"Title:{self.hist_ScintTimeLowBitsBySector.GetName()}")
        self.hist_ScintTime.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_ScintTime.GetName()}")
        self.hist_ScintCtime.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_ScintCtime.GetName()}")
        self.hist_ScintCtime0.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_ScintCtime0.GetName()}")
        self.hist_ScintCtimeRange.Draw()
        canvas.Print(self.pdfName, f"Title:{self.hist_ScintCtimeRange.GetName()}")
        pdfNameLast = f'{self.pdfName}]'
        canvas.Print(pdfNameLast, f"Title:{self.hist_ScintCtimeRange.GetName()}")
        self.histogramFile.Write()
        self.histogramFile.Close()
        print('Goodbye')
 

Member Data Documentation

◆ BKLM_ID

int BKLM_ID = 0x07000000

static

COPPER base identifier for BKLM readout.

Definition at line 23 of file EventInspectorPocketDAQ.py.

◆ BKLM_LAYER_BIT

int BKLM_LAYER_BIT = 7

static

bit position for layer-1 [0..14]; 0 is innermost

Definition at line 31 of file EventInspectorPocketDAQ.py.

◆ BKLM_LAYER_MASK

tuple BKLM_LAYER_MASK = (15 << BKLM_LAYER_BIT)

static

bit mask for layer-1 [0..15]; 0 is innermost and 14 is outermost

Definition at line 43 of file EventInspectorPocketDAQ.py.

◆ BKLM_MAXSTRIP_BIT

int BKLM_MAXSTRIP_BIT = 15

static

bit position for maxStrip-1 [0..47]

Definition at line 37 of file EventInspectorPocketDAQ.py.

◆ BKLM_MAXSTRIP_MASK

tuple BKLM_MAXSTRIP_MASK = (63 << BKLM_MAXSTRIP_BIT)

static

bit mask for maxStrip-1 [0..47]

Definition at line 49 of file EventInspectorPocketDAQ.py.

◆ BKLM_MODULEID_MASK

tuple BKLM_MODULEID_MASK = (BKLM_SECTION_MASK | BKLM_SECTOR_MASK | BKLM_LAYER_MASK)

static

bit mask for unique module identifier (end, sector, layer)

Definition at line 51 of file EventInspectorPocketDAQ.py.

◆ BKLM_PLANE_BIT

int BKLM_PLANE_BIT = 6

static

bit position for plane-1 [0..1]; 0 is inner-plane

Definition at line 29 of file EventInspectorPocketDAQ.py.

◆ BKLM_PLANE_MASK

tuple BKLM_PLANE_MASK = (1 << BKLM_PLANE_BIT)

static

bit mask for plane-1 [0..1]; 0 is inner-plane

Definition at line 41 of file EventInspectorPocketDAQ.py.

◆ BKLM_SECTION_BIT

int BKLM_SECTION_BIT = 14

static

bit position for section [0..1]; forward is 0

Definition at line 35 of file EventInspectorPocketDAQ.py.

◆ BKLM_SECTION_MASK

tuple BKLM_SECTION_MASK = (1 << BKLM_SECTION_BIT)

static

bit mask for section [0..1]; forward is 0

Definition at line 47 of file EventInspectorPocketDAQ.py.

◆ BKLM_SECTOR_BIT

int BKLM_SECTOR_BIT = 11

static

bit position for sector-1 [0..7]; 0 is on the +x axis and 2 is on the +y axis

Definition at line 33 of file EventInspectorPocketDAQ.py.

◆ BKLM_SECTOR_MASK

tuple BKLM_SECTOR_MASK = (7 << BKLM_SECTOR_BIT)

static

bit mask for sector-1 [0..7]; 0 is on the +x axis and 2 is on the +y axis

Definition at line 45 of file EventInspectorPocketDAQ.py.

◆ BKLM_STRIP_BIT

int BKLM_STRIP_BIT = 0

static

bit position for strip-1 [0..47]

Definition at line 27 of file EventInspectorPocketDAQ.py.

◆ BKLM_STRIP_MASK

int BKLM_STRIP_MASK = 0x3f

static

bit mask for strip-1 [0..47]

Definition at line 39 of file EventInspectorPocketDAQ.py.

◆ EKLM_ID

int EKLM_ID = 0x08000000

static

COPPER base identifier for EKLM readout.

Definition at line 25 of file EventInspectorPocketDAQ.py.

◆ exp

exp = exp

internal copy of experiment number

Definition at line 67 of file EventInspectorPocketDAQ.py.

◆ hist_ChannelOccupancy

list hist_ChannelOccupancy = [0, 0]

scatterplots of channel occupancy (1 hit per readout channel) for each axis

Definition at line 154 of file EventInspectorPocketDAQ.py.

◆ hist_ChannelOccupancyAL

list hist_ChannelOccupancyAL

Initial value:

=  [
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0],
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0], [0, 0],
            [0, 0], [0, 0], [0, 0], [0, 0], [0, 0]]

histograms of channel occupancy (1 hit per readout channel), indexed by axis/lane

Definition at line 162 of file EventInspectorPocketDAQ.py.

◆ hist_jRPCtdc

hist_jRPCtdc

Initial value:

=  ROOT.TH2F('jRPCtdc',
                                      expRun + 'RPC TDC vs hit index;tdc (ns);Hit index',
                                      64, 767.5, 1023.5, 60, -0.5, 59.5)

scatterplot of RPC calibrated time vs hit's index

Definition at line 230 of file EventInspectorPocketDAQ.py.

◆ hist_jRPCtdc2

hist_jRPCtdc2

Initial value:

=  ROOT.TH2F('jRPCtdc2',
                                       expRun + 'RPC TDC vs hit index;tdc - dt(index) (ns);Hit index',
                                       64, 767.5, 1023.5, 60, -0.5, 59.5)

scatterplot of RPC calibrated time vs hit's index, corrected for DC-processing delay

Definition at line 234 of file EventInspectorPocketDAQ.py.

◆ hist_PerChannelMultiplicity

hist_PerChannelMultiplicity

Initial value:

=  ROOT.TH2F(
            'PerChannelMultiplicity',
            expRun + 'Per-channel multiplicity (N/channel > 1);' +
                     'Per-channel multiplicity;' +
                     '(Lane #) * 2 + (Axis #)',
            30, -0.5, 29.5, 42, -0.5, 41.5)

scatterplot of multiplicity of entries in one readout channel vs lane/axis

Definition at line 133 of file EventInspectorPocketDAQ.py.

◆ hist_rawKLMlane

hist_rawKLMlane

Initial value:

=  ROOT.TH1F('rawKLMlane',
                                         expRun + 'RawKLM lane;Lane (scint: 1..7, RPC: 8..20)',
                                         21, -0.5, 20.5)

histogram of the hit's lane

Definition at line 121 of file EventInspectorPocketDAQ.py.

◆ hist_rawKLMsize

hist_rawKLMsize

Initial value:

=  ROOT.TH1F('rawKLMsize',
                                         expRun + 'RawKLM word count (1/channel)',
                                         200, -0.5, 199.5)

histogram of number of hits, at most one entry per readout channel

Definition at line 129 of file EventInspectorPocketDAQ.py.

◆ hist_rawKLMsizeMultihit

hist_rawKLMsizeMultihit

Initial value:

=  ROOT.TH1F('rawKLMsizeMultihit',
                                                 expRun + 'RawKLM word count (N/channel)',
                                                 200, -0.5, 199.5)

histogram of number of hits, including multiple entries on one readout channel

Definition at line 125 of file EventInspectorPocketDAQ.py.

◆ hist_revotimeRPCtdc

hist_revotimeRPCtdc

Initial value:

=  ROOT.TH2F('revotimeRPCtdc',
                                             expRun + 'RPC TDC vs revotime;tdc (ns);revotime - minRevotime',
                                             64, 767.5, 1023.5, 10, -0.5, 159.5)

scatterplot of RPC REVO9 range vs TDC value in event

Definition at line 219 of file EventInspectorPocketDAQ.py.

◆ hist_revotimeRPCtdc2

hist_revotimeRPCtdc2

Initial value:

=  ROOT.TH2F(
            'revotimeRPCtdc2',
            expRun + 'RPC TDC vs revotime;' +
            'tdc - dt(index) (ns);' +
            'revotime - minRevotime',
            64, 767.5, 1023.5, 10, -0.5, 159.5)

scatterplot of RPC REVO9 range vs TDC value corrected for DC-processing delay in event

Definition at line 223 of file EventInspectorPocketDAQ.py.

◆ hist_RPCLaneAxisOccupancy

hist_RPCLaneAxisOccupancy

Initial value:

=  ROOT.TH2F(
            'RPCLaneAxisOccupancy',
            expRun + 'Lane/axis occupancy of RPC channels (1/channel);' +
                     'Sector # (always 0);' +
                     '(Lane #) * 2 + (Axis #)',
            3, -1.5, 1.5, 42, -0.5, 41.5)

scatterplot of number of mapped RPC hits by lane/axis vs sector, at most one entry per readout channel

Definition at line 140 of file EventInspectorPocketDAQ.py.

◆ hist_RPCRevotimeRange

hist_RPCRevotimeRange

Initial value:

=  ROOT.TH1F('RPCRevotimeRange',
                                               expRun + 'RPC revotime-range in event;revotimeMax - revotimeMin (ns)',
                                               128, -0.5, 8191.5)

histogram of RPC REVO9 range in event

Definition at line 215 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTdcRange

hist_RPCTdcRange = ROOT.TH1F('RPCTdcRange', expRun + 'RPC TDC-range in event;TDCMax - TDCMin (ns)', 512, -0.5, 2047.5)

histogram of RPC TDC range in event

Definition at line 213 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTime

hist_RPCTime

Initial value:

=  ROOT.TH1F('RPCTime',
                                      expRun + 'RPC tdc relative to event trigtime;tdc - triggerTime (ns)',
                                      256, -0.5, 1023.5)

histogram of RPC TDC value relative to event's REVO9 trigger time in last word of event

Definition at line 181 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTime2

hist_RPCTime2

Initial value:

=  ROOT.TH1F('RPCTime_Ctime',
                                       expRun + 'RPC tdc relative to event ctime;tdc - trigCtime (ns)',
                                       256, -0.5, 1023.5)

histogram of RPC TDC value relative to event's ctime in event header

Definition at line 185 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTime2PerLayerA0

list hist_RPCTime2PerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

histograms of RPC TDC value relative to event's ctime for axis 0, indexed by lane

Definition at line 201 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTime2PerLayerA1

list hist_RPCTime2PerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

histograms of RPC TDC value relative to event's ctime for axis 1, indexed by lane

Definition at line 207 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTimeLowBitsBySector

hist_RPCTimeLowBitsBySector

Initial value:

=  ROOT.TH2F('RPCTimeLowBitsBySector',
                                                     expRun + 'RPC TDC lowest-order bits;' +
                                                     'Sector # (always 0);' +
                                                     'TDC % 4 (ns) [should be 0]',
                                                     3, -1.5, 1.5, 8, -0.25, 3.75)

scatterplot of RPC TDC low-order bits vs sector (should be 0 since granularity is 4 ns)

Definition at line 175 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTimePerLayerA0

list hist_RPCTimePerLayerA0 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

histograms of RPC TDC value relative to event's trigger time for axis 0, indexed by lane

Definition at line 189 of file EventInspectorPocketDAQ.py.

◆ hist_RPCTimePerLayerA1

list hist_RPCTimePerLayerA1 = [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

histograms of RPC TDC value relative to event's trigger time for axis 1, indexed by lane

Definition at line 195 of file EventInspectorPocketDAQ.py.

◆ hist_ScintCtime

hist_ScintCtime

Initial value:

=  ROOT.TH1F('ScintCtime',
                                         expRun + 'Scint ctime distribution;ctime - triggerCtime (ns)',
                                         32, -0.5, 1023.5)

histogram of scint CTIME value relative to event's ctime

Definition at line 249 of file EventInspectorPocketDAQ.py.

◆ hist_ScintCtime0

hist_ScintCtime0

Initial value:

=  ROOT.TH1F('ScintCtime0',
                                          expRun + 'Scint ctime distribution;ctime - triggerTime (ns)',
                                          32, -0.5, 1023.5)

histogram of scint CTIME value relative to event's trigger time

Definition at line 253 of file EventInspectorPocketDAQ.py.

◆ hist_ScintCtimeRange

hist_ScintCtimeRange

Initial value:

= ROOT.TH1F(

'ScintCtimeRange', expRun + 'Scint ctime-range in event;ctimeMax - ctimeMin (ns)', 128, -0.5, 8191.5)

histogram of scint CTIME range in event

Definition at line 257 of file EventInspectorPocketDAQ.py.

◆ hist_ScintLaneAxisOccupancy

hist_ScintLaneAxisOccupancy

Initial value:

=  ROOT.TH2F(
            'ScintLaneAxisOccupancy',
            expRun + 'Lane/axis occupancy of scint channels (1/channel);' +
                     'Sector # (always 0);' +
                     '(Lane #) * 2 + (Axis #)',
            3, -1.5, 1.5, 42, -0.5, 41.5)

scatterplot of number of mapped scint hits by lane/axis vs sector, at most one entry per readout channel

Definition at line 147 of file EventInspectorPocketDAQ.py.

◆ hist_ScintTime

hist_ScintTime

Initial value:

=  ROOT.TH1F('ScintTime',
                                        expRun + 'Scint tdc distribution;tdc - triggerTime (ns)',
                                        256, -0.5, 1023.5)

histogram of scint TDC value relative to event's trigger time

Definition at line 245 of file EventInspectorPocketDAQ.py.

◆ hist_ScintTimeLowBitsBySector

hist_ScintTimeLowBitsBySector

Initial value:

=  ROOT.TH2F(
            'ScintTimeLowBitsBySector',
            expRun + 'Scint TDC lowest-order bits;' +
            'Sector # (always 0);' +
            'TDC % 4 (ns)',
            3, -1.5, 1.5, 8, -0.25, 3.75)

scatterplot of scint TDC low-order bits vs sector

Definition at line 238 of file EventInspectorPocketDAQ.py.

◆ hist_ttc_trigtime

hist_ttc_trigtime

Initial value:

=  ROOT.TH1F('ttc_trigtime',
                                           expRun + 'tt_ctime relative to triggertime;' +
                                           'tt_ctime - triggertime (ns)',
                                           256, -0.5, 2047.5)

histogram of the tt_ctime relative to triggertime

Definition at line 116 of file EventInspectorPocketDAQ.py.

◆ histName

histName = histName

internal copy of the pathname of the output histogram ROOT file

Definition at line 71 of file EventInspectorPocketDAQ.py.

◆ histogramFile

histogramFile = ROOT.TFile.Open(self.histName, "RECREATE")

Output ROOT TFile that will contain the histograms/scatterplots.

Definition at line 111 of file EventInspectorPocketDAQ.py.

◆ pdfName

pdfName = pdfName

internal copy of the pathname of the output histogram PDF file

Definition at line 73 of file EventInspectorPocketDAQ.py.

◆ run

run = run

internal copy of run number

Definition at line 69 of file EventInspectorPocketDAQ.py.

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

klm/bklm/scripts/EventInspectorPocketDAQ.py

Public Member Functions

Public Attributes

Static Public Attributes

Detailed Description

Constructor & Destructor Documentation

◆ __init__()

Member Function Documentation

◆ beginRun()

◆ endRun()

◆ event()

◆ initialize()

◆ makeGraph()

◆ makeText()

◆ terminate()

Member Data Documentation

◆ BKLM_ID

◆ BKLM_LAYER_BIT

◆ BKLM_LAYER_MASK

◆ BKLM_MAXSTRIP_BIT

◆ BKLM_MAXSTRIP_MASK

◆ BKLM_MODULEID_MASK

◆ BKLM_PLANE_BIT

◆ BKLM_PLANE_MASK

◆ BKLM_SECTION_BIT

◆ BKLM_SECTION_MASK

◆ BKLM_SECTOR_BIT

◆ BKLM_SECTOR_MASK

◆ BKLM_STRIP_BIT

◆ BKLM_STRIP_MASK

◆ EKLM_ID

◆ exp

◆ hist_ChannelOccupancy

◆ hist_ChannelOccupancyAL

◆ hist_jRPCtdc

◆ hist_jRPCtdc2

◆ hist_PerChannelMultiplicity

◆ hist_rawKLMlane

◆ hist_rawKLMsize

◆ hist_rawKLMsizeMultihit

◆ hist_revotimeRPCtdc

◆ hist_revotimeRPCtdc2

◆ hist_RPCLaneAxisOccupancy

◆ hist_RPCRevotimeRange

◆ hist_RPCTdcRange

◆ hist_RPCTime

◆ hist_RPCTime2

◆ hist_RPCTime2PerLayerA0

◆ hist_RPCTime2PerLayerA1

◆ hist_RPCTimeLowBitsBySector

◆ hist_RPCTimePerLayerA0

◆ hist_RPCTimePerLayerA1

◆ hist_ScintCtime

◆ hist_ScintCtime0

◆ hist_ScintCtimeRange

◆ hist_ScintLaneAxisOccupancy

◆ hist_ScintTime

◆ hist_ScintTimeLowBitsBySector

◆ hist_ttc_trigtime

◆ histName

◆ histogramFile

◆ pdfName

◆ run

◆ init()