InDetGlobalPrimaryVertexMonAlgCfg.py

Go to the documentation of this file.

#
#  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
#
 
"""
@file InDetGlobalPrimaryVertexMonAlgCfg.py
@author Leonid Serkin and Per Johansson
@date April 2020
@brief Configuration for Run 3 based on InDetGlobalPrimaryVertexMonTool.cxx
"""
 
from AthenaConfiguration.ComponentFactory import CompFactory
 
def HistoInDetGlobalPrimaryVertexMonAlgCfg(helper, alg):
 
    # this creates a "pvGroup" called "alg" which will put its histograms into the subdirectory "PrimaryVertex"
    pvGroup = helper.addGroup(alg, 'PrimaryVertex')
    pathpv = '/InDetGlobal/PrimaryVertex'
 
    # begin histogram definitions
    varName = 'm_PvX;pvX' #done
    title = 'Primary vertex: x;x (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=500, xmin=-1.0, xmax=1.0)
 
    varName = 'm_PvY;pvY' #done
    title = 'Primary vertex: y;y (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=500, xmin=-1.0, xmax=1.0)
 
    varName = 'm_PvZ;pvZ' #done
    title = 'Primary vertex: z;z (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=100, xmin=-200., xmax=200.)
 
    varName = 'm_PvN;pvN' #done
    title = 'Total number of vertices (primary and pile up);Total number of vertices;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=60, xmin=-0.5, xmax=59.5)
 
    varName = 'm_nPriVtx;pvNPriVtx' #done
    title = 'Number of primary vertices;Number of primary vertices;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=3, xmin=-0.5, xmax=2.5)
 
    varName = 'm_nPileupVtx;pvNPileupVtx' #done
    title   = 'Number of pileup vertices;Number of pile up vertices;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=50, xmin=-0.5, xmax=49.5)
    
    varName = 'm_PvErrX;pvErrX'  #done
    title   = 'Primary vertex: #sigma_{x}; #sigma_{x} (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=200, xmin=0., xmax=0.03)
  
    varName = 'm_PvErrY;pvErrY'  #done
    title   = 'Primary vertex: #sigma_{y}; #sigma_{y} (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=200, xmin=0., xmax=0.03)
 
    varName = 'm_PvErrZ;pvErrZ'   #done
    title   = 'Primary vertex: #sigma_{z}; #sigma_{z} (mm);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=100, xmin=0., xmax=0.1)
 
    varName = 'm_PvChiSqDoF;pvChiSqDof'   #done
    title   = 'Primary vertex: #Chi^{2}/DoF of vertex fit;#Chi^{2}/DoF;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=100, xmin=0., xmax=5.)
 
    varName = 'm_PvN;pvNTracks'  #done
    title = 'Number of tracks in primary vertex;Number of tracks;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=300, xmin=0., xmax=300.)
 
    varName = 'm_PvTrackPt;pvTrackPt'  #done
    title   = 'Primary vertex: original track p_{t};p_{t} (GeV);Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=100, xmin=0., xmax=20.)
 
    varName = 'm_PvTrackEta;pvTrackEta'  #done
    title   = 'Primary vertex: original track #eta; #eta;Events'
    pvGroup.defineHistogram(varName, type='TH1F', path=pathpv, title=title, xbins=100, xmin=-3, xmax=3)
 
# end histograms
 
 
def InDetGlobalPrimaryVertexMonAlgCfg(helper, acc,
                                      flags, name="InDetGlobalPrimaryVertexMonAlg",
                                      **kwargs):
    kwargs.setdefault("doEnhancedMonitoring", True)
 
    from AthenaMonitoring.FilledBunchFilterToolConfig import FilledBunchFilterToolCfg
    from AthenaMonitoring.AtlasReadyFilterConfig import AtlasReadyFilterCfg
 
    monAlg = helper.addAlgorithm(
        CompFactory.InDetGlobalPrimaryVertexMonAlg, name,
        addFilterTools = [FilledBunchFilterToolCfg(flags), AtlasReadyFilterCfg(flags)],
        **kwargs)
 
    HistoInDetGlobalPrimaryVertexMonAlgCfg(helper, monAlg)
    return