DQTLumiMonAlg.py

Go to the documentation of this file.

#
#  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
#
 
def DQTLumiMonAlgConfig(flags):
    from AthenaMonitoring import AthMonitorCfgHelper
    from AthenaConfiguration.ComponentFactory import CompFactory
    DQTLumiMonAlg = CompFactory.DQTLumiMonAlg
 
    helper = AthMonitorCfgHelper(flags, 'DQTLumiMonAlgCfg')
    # Three instances of the algorithm. One using any trigger, another using only muon
    # triggers, and the final using only electron triggers.
    DQTLumiMonAlgConfigByTriggerChain(helper, DQTLumiMonAlg, flags)
    if flags.DQ.useTrigger:
        DQTLumiMonAlgConfigByTriggerChain(helper, DQTLumiMonAlg, flags,
                                        'CATEGORY_monitoring_muonIso','EF_muX')
        DQTLumiMonAlgConfigByTriggerChain(helper, DQTLumiMonAlg, flags,
                                        'CATEGORY_primary_single_ele','EF_eX')
    return helper.result()
 
 
def DQTLumiMonAlgConfigByTriggerChain(helper, algConfObj, flags, triggerChain='', triggerPath=''):
    monAlg = helper.addAlgorithm(algConfObj, 'DQTLumiMonAlg'+triggerPath)
 
    if triggerChain:
        monAlg.TriggerChain = triggerChain
        if not triggerPath:
            triggerPath = triggerChain
    else:
        triggerPath = 'AnyTrigger'
 
    groupPath = 'GLOBAL/Luminosity/'+triggerPath
    group = helper.addGroup(monAlg, 'default', groupPath, 'run')
    pixelgroup = helper.addGroup(monAlg, 'pixel', groupPath, 'run')
    nLooseT = 'Number of Loose Vertices'
    nTightT = 'Number of Tight Vertices'
    nPixClusterT = 'Number of Pixel Clusters'
 
    # Vertex parameters. A vertex
    monAlg.TightTrackWeight = 0.01
    monAlg.TightNTracks = 4
 
    lbdict = { 'xmin': min(flags.Input.LumiBlockNumbers)-0.5 if flags.Input.LumiBlockNumbers else 0.5,
               'xmax': max(flags.Input.LumiBlockNumbers)+0.5 if flags.Input.LumiBlockNumbers else 1.5 }
    lbdict['xbins'] = int(lbdict['xmax']-lbdict['xmin'])
 
    # Raw plots of lumi variables
    group.defineHistogram('avgIntPerXing', type='TH1F', xbins=100, xmin=0, xmax=100,
                          title='Average number of interactions per bunch Xing')
    group.defineHistogram('intPerXing', type='TH1F', xbins=100, xmin=0, xmax=100,
                          title='Actual number of interactions per bunch Xing')
 
    # Plots of luminosity variables as a function of lb number
    group.defineHistogram('LB,avgLumi', type='TProfile',
                          title='Average Lumi vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,avgIntPerXing', type='TProfile',
                          title='Number of interactions per event;LB;<#mu>_{LB}',
                          opt='kAddBinsDynamically', merge='merge', **lbdict)
    group.defineHistogram('LB,lumiPerBCID', type='TProfile',
                          title='Instantaneous Luminosity vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,intPerXing', type='TProfile',
                          title='Instaneous interactions vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,duration', type='TProfile',
                          title='Lumi Block time in sec vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,avgLiveFrac', type='TProfile',
                          title='Average live fraction lumi over all BCIDs vs LB',
                          opt='kAddBinsDynamically', merge='merge', **lbdict)
    group.defineHistogram('LB,liveFracPerBCID', type='TProfile',
                          title='Current BCID lumi vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,lumiWeight', type='TProfile',
                          title='Current BCID lumi vs LB', opt='kAddBinsDynamically',
                          merge='merge', **lbdict)
 
    # Vertex counts vs. lb number (loose and tight)
    group.defineHistogram('LB,nLooseVtx', type='TProfile',
                          title=nLooseT+' per event;LB;<NlooseVtx/event>_{LB}',
                          opt='kAddBinsDynamically', merge='merge', **lbdict)
    group.defineHistogram('LB,nTightVtx', type='TProfile',
                          title=nTightT+' per event;LB;<NtightVtx/event>_{LB}',
                          opt='kAddBinsDynamically', merge='merge', **lbdict)
 
    # Vertex counts per avg mu vs. lb number. (Same as above, just weighted.)
    group.defineHistogram('LB,nLooseVtxPerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                          title=nLooseT+' per event per Mu;LB;<NlooseVtx/event/#mu>_{LB}',
                          merge='merge', **lbdict)
    group.defineHistogram('LB,nTightVtxPerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                          title=nTightT+' per event per Mu;LB;<NtightVtx/event/#mu>_{LB}',
                          merge='merge', **lbdict)
 
    # Vertex counts vs. avg mu.
    group.defineHistogram('avgIntPerXing,nLooseVtx', type='TProfile', xbins=250, xmin=0, xmax=100,
                          title=nLooseT+' per event;#mu;NlooseVtx/event')
    group.defineHistogram('avgIntPerXing,nTightVtx', type='TProfile', xbins=250, xmin=0, xmax=100,
                          title=nTightT+' per event;#mu;NtightVtx/event')
 
    # Pixel clusters vs lb number. All, and then eparately for each endcap and barrel layer.
    pixelgroup.defineHistogram('LB,nClustersAll', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', all;LB;<NclustersAll/event>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersECA', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', endcap A;LB;<NlustersECA/event>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersECC', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', endcap C;LB;<NclustersECC/event>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB0', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', barrel layer 0;LB;<NclustersB0/event>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB1', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', barrel layer 1;LB;<NclustersB1/event>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB2', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+', barrel layer 2;LB;<NclustersB2/event>_{LB}',
                               merge='merge', **lbdict)
 
    # Pixel clusters per avg mu vs. lb number. (Same as above, just weighted.)
    pixelgroup.defineHistogram('LB,nClustersAllPerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, all;LB;<NclustersAll/event/#mu>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersECAPerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, endcap A;LB;<NlustersECA/event/#mu>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersECCPerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, endcap C;LB;<NclustersECC/event/#mu>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB0PerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, barrel layer 0;LB;<NclustersB0/event/#mu>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB1PerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, barrel layer 1;LB;<NclustersB1/event/#mu>_{LB}',
                               merge='merge', **lbdict)
    pixelgroup.defineHistogram('LB,nClustersB2PerAvgMu', type='TProfile', opt='kAddBinsDynamically',
                               title=nPixClusterT+' per Mu, barrel layer 2;LB;<NclustersB2/event/#mu>_{LB}',
                               merge='merge', **lbdict)
 
    # Pixel clusters vs. avg mu. All, then separately for each endcap and barrel layer.
    pixelgroup.defineHistogram('avgIntPerXing,nClustersAll', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', all;#mu;NclustersAll/event')
    pixelgroup.defineHistogram('avgIntPerXing,nClustersECA', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', endcap A;#mu;NclustersECA/event')
    pixelgroup.defineHistogram('avgIntPerXing,nClustersECC', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', endcap C;#mu;NclustersECC/event')
    pixelgroup.defineHistogram('avgIntPerXing,nClustersB0', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', barrel layer 0;#mu;NclustersB0/event')
    pixelgroup.defineHistogram('avgIntPerXing,nClustersB1', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', barrel layer 1;#mu;NclustersB1/event')
    pixelgroup.defineHistogram('avgIntPerXing,nClustersB2', type='TProfile',
                               xbins=250, xmin=0, xmax=100,
                               title=nPixClusterT+', barrel layer 2;#mu;NclustersB2/event')