TileTBPulseMonitorAlgorithm Namespace Reference

Functions
	getPMT (partition, channel)
	TileTBPulseMonitoringConfig (flags, timeRange=[-100, 100], fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], useDemoCabling=2018, useFELIX=False, **kwargs)

Detailed Description

@file TileTBPulseMonitorAlgorithm.py
@brief Python configuration of TileTBPulseMonitorAlgorithm algorithm for the Run III

Function Documentation

getPMT()

TileTBPulseMonitorAlgorithm.getPMT	(		partition,
			channel )

Function to get PMT number: 0,1 

Definition at line 13 of file TileTBPulseMonitorAlgorithm.py.

def getPMT(partition, channel):
    ''' Function to get PMT number: 0,1 '''
 
    channelToPMT = [
        # LBA,LBC:
        0,  1,  0,  1,  0,  1,  0,  1,  0,  1,  0,  1,
        0,  1,  0,  1,  0,  1,  0,  1,  0,  1,  0,  1,
        0,  1,  0,  1,  0,  1, -1, -1,  0,  1,  0,  1,
        0,  1,  0,  1,  0,  1,  0, -1,  0,  1,  0,  1,
        # EBA,EBC:
        0,  0,  0,  1,  0,  1,  0,  1,  0,  1,  0,  1,
        0,  0,  0,  1,  0,  1, -1, -1,  0,  1,  0,  1,
        -1, -1, -1, -1, -1, -1,  0,  0,  1, -1, -1,  1,
        1,  1,  0,  0,  1,  0, -1, -1, -1, -1, -1, -1 ]
 
    # In gap scintillators (E3,E4,E1,E2) there is only one pmt per cell
    if partition in ['EBA', 'EBC'] and channel in [0, 1, 12, 13]:
        pmt = 0
    else:
        pmt = channelToPMT[channel+48] if partition in ['EBA', 'EBC'] else channelToPMT[channel]
 
        # Mirroring of odd/even numbers in negative side
        # (central symmetry of negative/positive drawers)
        if (pmt != -1 and partition in ['LBC', 'EBC']):
            pmt = 1 - pmt
 
    return pmt
 
 

TileTBPulseMonitoringConfig()

TileTBPulseMonitorAlgorithm.TileTBPulseMonitoringConfig	(		flags,
			timeRange = [-100, 100],
			fragIDs = [0x100, 0x101, 0x200, 0x201, 0x402],
			useDemoCabling = 2018,
			useFELIX = False,
		**	kwargs )

Function to configure TileTBPulseMonitorAlgorithm algorithm in the monitoring system.

Definition at line 42 of file TileTBPulseMonitorAlgorithm.py.

def TileTBPulseMonitoringConfig(flags, timeRange=[-100, 100], fragIDs=[0x100, 0x101, 0x200, 0x201, 0x402], useDemoCabling=2018, useFELIX=False, **kwargs):
 
    ''' Function to configure TileTBPulseMonitorAlgorithm algorithm in the monitoring system.'''
 
    suffix = "Flx" if useFELIX else ""
    topPath = 'TestBeam/' + ('Felix' if useFELIX else 'Legacy')
 
    from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
    result = ComponentAccumulator()
 
    from TileGeoModel.TileGMConfig import TileGMCfg
    result.merge(TileGMCfg(flags))
 
    from TileConditions.TileCablingSvcConfig import TileCablingSvcCfg
    result.merge(TileCablingSvcCfg(flags))
 
    from TileConditions.TileInfoLoaderConfig import TileInfoLoaderCfg
    result.merge(TileInfoLoaderCfg(flags))
 
    from AthenaMonitoring import AthMonitorCfgHelper
    helper = AthMonitorCfgHelper(flags, f'TileTBPulse{suffix}Monitoring')
 
    from AthenaConfiguration.ComponentFactory import CompFactory
    tileTBPulseMonAlg = helper.addAlgorithm(CompFactory.TileTBPulseMonitorAlgorithm, f'TileTBPulse{suffix}MonAlg')
 
    tileTBPulseMonAlg.TriggerChain = ''
 
    from TileCalibBlobObjs.Classes import TileCalibUtils as Tile
 
    modules = []
    if fragIDs:
        for fragID in fragIDs:
            ros = fragID >> 8
            drawer = fragID & 0x3F
            modules += [Tile.getDrawerString(ros, drawer)]
    else:
        for ros in range(1, Tile.MAX_ROS):
            for drawer in range(0, Tile.MAX_DRAWER):
                fragIDs += [(ros << 8) | drawer]
                modules += [Tile.getDrawerString(ros, drawer)]
 
    tileTBPulseMonAlg.TileFragIDs = fragIDs
 
    kwargs.setdefault('TileDigitsContainer', f'TileDigits{suffix}Cnt')
    kwargs.setdefault('TileRawChannelContainer', flags.Tile.RawChannelContainer.replace('TileRawChannel', f'TileRawChannel{suffix}'))
    for k, v in kwargs.items():
        setattr(tileTBPulseMonAlg, k, v)
 
    run = str(flags.Input.RunNumbers[0])
 
    # Configure histogram with TileTBPulseMonAlg algorithm execution time
    executeTimeGroup = helper.addGroup(tileTBPulseMonAlg, 'TileTBPulseMonExecuteTime', topPath)
    executeTimeGroup.defineHistogram('TIME_execute', path='PulseShape', type='TH1F',
                                     title='Time for execute TileTBPulseMonAlg algorithm;time [#mus]',
                                     xbins=100, xmin=0, xmax=10000)
 
    from TileMonitoring.TileMonitoringCfgHelper import getCellName
 
    def addPulseShapeHistogramsArray(helper, modules, algorithm, name, title, path, type='TH2D',
                                     xbins=100, xmin=-100, xmax=100, ybins=100, ymin=-0.2, ymax=1.5,
                                     run='', value='', aliasPrefix='', useDemoCabling=2018):
        ''' This function configures 2D (or 1D Profile) histograms with Tile pulse shape per module, channel, gain  '''
 
        pulseShapeArray = helper.addArray([modules], algorithm, name, topPath=path)
        for postfix, tool in pulseShapeArray.Tools.items():
            moduleName = postfix[1:]
            partition = moduleName[:3]
            module = int(moduleName[3:]) - 1
            for channel in range(0, Tile.MAX_CHAN):
                legacyChannel = getLegacyChannelForDemonstrator(useDemoCabling, partition, module, channel)
                pmt = getPMT(partition, legacyChannel)
                pmtName = f'Channel_{channel}' if pmt < 0 else {0 : 'PMT_Up', 1 : 'PMT_Down'}[pmt]
                cell = getCellName(partition, legacyChannel)
                cellName = cell.replace('B', 'BC') if (partition in ['LBA','LBC'] and cell and cell[0] == 'B' and cell != 'B9') else cell
                for gain in range(0, Tile.MAX_GAIN):
                    gainName = {0 : 'lo', 1 : 'hi'}[gain]
                    fullPath = f'{partition}/{moduleName}'
                    name = f'time_{channel}_{gain},amplitude_{channel}_{gain};{aliasPrefix}{cellName}_{moduleName}_{pmtName}_{gainName}'
                    fullTitle = f'Run {run} {moduleName} Channel {channel} Gain {gainName}: {title};time [ns];Normalized Units'
                    tool.defineHistogram(name, title=fullTitle, path=fullPath, type=type,
                                         xbins=xbins, xmin=xmin, xmax=xmax, ybins=ybins, ymin=ymin, ymax=ymax)
        return pulseShapeArray
 
    addPulseShapeHistogramsArray(helper, modules, tileTBPulseMonAlg, name='TilePulseShape', title='Pulse shape',
                                 path=f'{topPath}/PulseShape', xbins=abs(timeRange[1]), xmin=timeRange[0], xmax=timeRange[1],
                                 run=run, aliasPrefix='pulseShape_', useDemoCabling=useDemoCabling)
 
    addPulseShapeHistogramsArray(helper, modules, tileTBPulseMonAlg, name='TilePulseShapeProfile',
                                 title='Pulse shape profile', path=f'{topPath}/PulseShape', type='TProfile',
                                 xbins=abs(timeRange[1]), xmin=timeRange[0], xmax=timeRange[1],
                                 ybins=None, ymin=None, ymax=None, run=run, aliasPrefix='pulseShapeProfile_',
                                 useDemoCabling=useDemoCabling)
 
    accumalator = helper.result()
    result.merge(accumalator)
    return result