PPMSimBSMonitorAlgorithm.py

Go to the documentation of this file.

#
#  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
#
import sys 
 
 
def PPMSimBSMonitoringConfig(flags):
    '''Function to configure LVL1 PPM Sim algorithm in the monitoring system.'''
 
 
    # get the component factory - used for getting the algorithms
    from AthenaConfiguration.ComponentFactory import CompFactory
    from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
    result = ComponentAccumulator()
 
    # make the athena monitoring helper
    from AthenaMonitoring import AthMonitorCfgHelper
    helper = AthMonitorCfgHelper(flags,'PPMSimBSMonitoringCfg')
 
    # get any algorithms
    PPMSimBSMonAlg = helper.addAlgorithm(CompFactory.PPMSimBSMonitorAlgorithm,'PPMSimBSMonAlg')
    
    # import tools
    from DetDescrCnvSvc.DetDescrCnvSvcConfig import DetDescrCnvSvcCfg # Needed for L1TTIDTools (ATR-21865)
    helper.result().merge(DetDescrCnvSvcCfg(flags))
 
    from TrigT1CaloCondSvc.L1CaloCondConfig import L1CaloCondAlgCfg 
    helper.result().merge(L1CaloCondAlgCfg(flags,Physics=True, Calib1=False, Calib2=False)) 
 
    from TrigConfxAOD.TrigConfxAODConfig import getxAODConfigSvc
    helper.result().getPrimaryAndMerge(getxAODConfigSvc(flags))
 
    # add any steering
    groupName = 'PPMSimBSMonitor' # the monitoring group name is also used for the package name
    PPMSimBSMonAlg.PackageName = groupName
    
    # Histogram paths
    mainDir = 'L1Calo'
    trigPath = 'PPM/'
 
    # Trigger tower plots: eta-phi granularity
    etabins = [-4.9,-4.475,-4.050,-3.625,-3.2,-3.1,-2.9,
                     -2.7,-2.5,-2.4,-2.3,-2.2,-2.1,-2.0,-1.9,
                     -1.8,-1.7,-1.6,-1.5,-1.4,-1.3,-1.2,-1.1,
                     -1.0,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,
                     -0.2,-0.1,0.0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,
                     0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,
                     1.8,1.9,2.0,2.1,2.2,2.3,2.4,2.5,2.7,2.9,
                     3.1,3.2,3.625,4.050,4.475,4.9]
 
 
 
 
    phibins = 64
    phimin = 0
 
    # Common location 
    trigPath = "/PPM/Errors/Data_Simulation"
    
    
    
   
    group_names = {"peakf": "Peakfinder", "satBcid": "SaturatedBCID"}
    mismatch_map= {"NonZeroMatches": "Data/Simulation Non-Zero Matches", "ZeroMatches":  "Data/Simulation Zero Matches", "DataNoSim":"Data but no Simulation", "SimNoData":"Simulation but no Data"}
 
 
    # Eta-phi maps
    histPath = trigPath+'/PPMBcidSim'
    for names in group_names:
        for elem in mismatch_map:
            group_bcidsim = helper.addGroup(PPMSimBSMonAlg,  'group_Mismatch_{0}_{1}'.format(names, elem),  mainDir)
            group_bcidsim.defineHistogram('etaTT_2D,phiTT_2D;ppm_2d_etaPhi_tt_'+names+'_'+elem, title='PPM '+group_names[names]+' ' + mismatch_map[elem],  type='TH2D', path=histPath, xbins=etabins, ybins=phibins, ymin=phimin, ymax=phibins, opt='kAlwaysCreate')
    
 
 
    
    
    histPath_CP = trigPath+'/PPMLUTSim'
    histPath_JEP = trigPath+'/PPMLUTSim'
    lut_map =  { "SimEqData": "Data/Simulation Non-zero Matches", "SimNeData": "Data/Simulation Non-zero Mismatches", "SimNoData":  "Simulation but no Data", "DataNoSim": "Data but no Simulation" }
    group_names_lut = {"lutCp": ["CP", histPath_CP], "lutJep": ["JEP", histPath_JEP] }
    
    # EM layer
    for names in group_names_lut:
        for elem in lut_map:
            group_lut_em = helper.addGroup(PPMSimBSMonAlg,  'groupLUT{0}_EM_{1}'.format(group_names_lut[names][0],elem),  mainDir)
            group_lut_em.defineHistogram('etaTT_2D,phiTT_2D;ppm_em_2d_etaPhi_tt_'+names+'_'+elem, title='PPM LUT-'+group_names_lut[names][0]+' EM ' + lut_map[elem],  type='TH2D', path=group_names_lut[names][1], xbins=etabins, ybins=phibins, ymin=phimin, ymax=phibins, opt='kAlwaysCreate')
    
    # HAD layer 
    for names in group_names_lut:
        for elem in lut_map:
            group_lut_had = helper.addGroup(PPMSimBSMonAlg,  'groupLUT{0}_HAD_{1}'.format(group_names_lut[names][0],elem),  mainDir)
            group_lut_had.defineHistogram('etaTT_2D,phiTT_2D;ppm_had_2d_etaPhi_tt_'+names+'_'+elem, title='PPM LUT-'+group_names_lut[names][0]+' HAD ' + lut_map[elem],  type='TH2D', path=group_names_lut[names][1], xbins=etabins, ybins=phibins, ymin=phimin, ymax=phibins, opt='kAlwaysCreate')
 
    
  
    # Mismatch errors 
    nCrates          = 8
    nModulesPerCrate = 16
 
    histPath = trigPath+'/MismatchEventNumbers'
 
    y = 2 * (nModulesPerCrate + 16 * (nCrates % 2) )
 
    yErrorLabels = []
    
 
    crate_map= [ "cr0cr1", "cr2cr3", "cr4cr5", "cr6cr7"]
    
    for elem in crate_map:
        yErrorLabels =  BinErrors(y, elem)
        group_Error= helper.addGroup(PPMSimBSMonAlg, 'group_Error_{}'.format(elem), mainDir)     
        group_Error.defineHistogram('eventMonitor,y_2D;ppm_2d_LUT_MismatchEvents_'+ elem, title='PPM LUT Mismatch Event Numbers;Events with Error/Mismatch;Crate/Module', type='TH2I', path=histPath,xbins=10,xmin=0, xmax=10, ybins=y, ymin=0, ymax=y,ylabels=yErrorLabels,opt='kAlwaysCreate')
        
    # Finish up
 
    acc = helper.result()
    result.merge(acc)
    return result
 
def BinErrors(y, crs):
 
    cr0 = crs.split("cr")[1]
    cr1 = crs.split("cr")[2]
 
 
    yErrorLabels = []
    for i in range(0, int(y/2)):
        yErrorLabels.append(str(cr0)+'/'+str(i))
    for i in range(0, int(y/2)):
        yErrorLabels.append(str(cr1)+'/'+str(i))
        
    return yErrorLabels
        
 
 
 
 
if __name__=='__main__':
   
  
    from AthenaConfiguration.AllConfigFlags import initConfigFlags
 
    flags = initConfigFlags()
    flags.Exec.MaxEvents = 10
    flags.IOVDb.GlobalTag = 'CONDBR2-BLKPA-2022-02'
    flags.Input.Files =  ["/eos/atlas/atlastier0/rucio/data22_13p6TeV/express_express/00423433/data22_13p6TeV.00423433.express_express.recon.ESD.x653/data22_13p6TeV.00423433.express_express.recon.ESD.x653._lb0015._SFO-ALL._0001.2"]
 
    flags.Output.HISTFileName = 'ExampleMonitorOutput_LVL1.root'
 
    
    flags.lock()
 
    from AthenaConfiguration.MainServicesConfig import MainServicesCfg
    from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
    acc = MainServicesCfg(flags)
    acc.merge(PoolReadCfg(flags))
    
 
    PPMSimBSMonitorCfg = PPMSimBSMonitoringConfig(flags)
    acc.merge(PPMSimBSMonitorCfg)
 
 
    acc.printConfig(withDetails=True, summariseProps = True)
 
    sys.exit(acc.run().isFailure())