Noise transforms for LArCAF

This brings back in new configuration the transforms computing noise bursts in LArCAF jobs. The code is not part of standard reco, is used only in specialised LAr calibration jobs.

LArCalorimeter/LArCafJobs/python/LArNoiseConfig.py

+# Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
+
+from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
+from AthenaConfiguration.ComponentFactory import CompFactory 
+
+def LArNoiseCfg(flags):
+
+    result=ComponentAccumulator()
+
+    from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
+    result.merge(PoolReadCfg(flags))
+    
+    from LArGeoAlgsNV.LArGMConfig import LArGMCfg
+    result.merge(LArGMCfg(flags))
+    from TileGeoModel.TileGMConfig import TileGMCfg
+    result.merge(TileGMCfg(flags))
+
+    #Setup cabling
+    from LArCabling.LArCablingConfig import LArOnOffIdMappingCfg
+    result.merge(LArOnOffIdMappingCfg(flags))
+    # setup bad chan and missing febs
+    from LArBadChannelTool.LArBadChannelConfig import LArBadChannelCfg,LArBadFebCfg
+    LArBadChannelCfg(flags)
+    LArBadFebCfg(flags)
+
+    from LumiBlockComps.BunchCrossingCondAlgConfig import BunchCrossingCondAlgCfg
+    result.merge(BunchCrossingCondAlgCfg(flags))
+    
+    from LArCellRec.LArCollisionTimeConfig import LArCollisionTimeCfg
+    result.merge(LArCollisionTimeCfg(flags))
+
+    from AthenaMonitoring.TriggerInterface import getTrigDecisionTool
+    result.merge(getTrigDecisionTool(flags))
+    result.getPublicTool("TrigDecisionTool").TrigConfigSvc="TrigConf::TrigConfigSvc/TrigConfigSvc"
+
+    
+    noiseAlg=CompFactory.LArNoiseBursts("LArNoiseBursts")
+    noiseAlg.SigmaCut = flags.LArNoise.SigmaCut
+    noiseAlg.NumberOfBunchesInFront = flags.LArNoise.NumberOfBunchesInFront
+    noiseAlg.KeepOnlyCellID = flags.LArNoise.KeepOnlyCellID
+    result.addEventAlgo(noiseAlg)
+
+    if (flags.LArNoise.outNtupLAr!=""):
+        result.addService(CompFactory.THistSvc(Output=["TTREE DATAFILE='"+flags.LArNoise.outNtupLAr+"' OPT='RECREATE'",]))
+
+    if (flags.LArNoise.HECNoiseNtup!=""):
+        hecAlg=CompFactory.LArHECNoise("LArHECNoise")
+        hecAlg.MinDigitADC = flags.LArNoise.MinDigitADC
+        hecAlg.MaxDeltaT = flags.LArNoise.MaxDeltaT
+        from IOVDbSvc.IOVDbSvcConfig import addFolders
+        result.merge(addFolders(flags, "/LAR/ElecCalibFlat/Pedestal", "LAR_ONL", className="CondAttrListCollection", db="CONDBR2"))
+        LArPedestalCondAlg =  CompFactory.getComp("LArFlatConditionsAlg<LArPedestalFlat>")
+        result.addCondAlgo(LArPedestalCondAlg(ReadKey="/LAR/ElecCalibFlat/Pedestal", WriteKey="LArPedestal"))
+
+        result.addEventAlgo(hecAlg)
+        if result.getService("THistSvc") is not None:
+           result.getService("THistSvc").Output += ["HEC DATAFILE='"+flags.LArNoise.HECNoiseNtup+"' OPT='RECREATE'",]
+        else:   
+           result.addService(CompFactory.THistSvc(Output=["HEC DATAFILE='"+flags.LArNoise.HECNoiseNtup+"' OPT='RECREATE'",]))
+
+    return result
+
+def LArNoiseFromRawCfg(flags):
+
+    result=ComponentAccumulator()
+    result.debugMode = "trackCA trackEventAlgo"
+
+    from LArByteStream.LArRawDataReadingConfig import LArRawDataReadingCfg
+    result.merge(LArRawDataReadingCfg(flags))
+
+    #Setup cabling
+    from LArCabling.LArCablingConfig import LArOnOffIdMappingCfg
+    result.merge(LArOnOffIdMappingCfg(flags))
+    # setup bad chan and missing febs
+    from LArBadChannelTool.LArBadChannelConfig import LArBadChannelCfg,LArBadFebCfg
+    LArBadChannelCfg(flags)
+    LArBadFebCfg(flags)
+
+    from CaloRec.CaloRecoConfig import CaloRecoCfg
+    result.merge(CaloRecoCfg(flags))
+
+    from LumiBlockComps.BunchCrossingCondAlgConfig import BunchCrossingCondAlgCfg
+    result.merge(BunchCrossingCondAlgCfg(flags))
+    
+    from CaloTools.CaloNoiseCondAlgConfig import CaloNoiseCondAlgCfg
+    result.merge(CaloNoiseCondAlgCfg(flags,"totalNoise"))
+
+    from LArROD.LArFebErrorSummaryMakerConfig import LArFebErrorSummaryMakerCfg
+    result.merge(LArFebErrorSummaryMakerCfg(flags))
+    result.getEventAlgo("LArFebErrorSummaryMaker").CheckAllFEB=False
+
+    if (flags.LArNoise.outNtupLAr != "" or flags.LArNoise.HECNoiseNtup!=""):
+       from LArCellRec.LArTimeVetoAlgConfig import LArTimeVetoAlgCfg
+       result.merge(LArTimeVetoAlgCfg(flags))
+
+       if (flags.LArNoise.outHistLAr == ""):
+          from AthenaMonitoring.TriggerInterface import getTrigDecisionTool
+          result.merge(getTrigDecisionTool(flags))
+          result.getPublicTool("TrigDecisionTool").TrigConfigSvc="TrigConf::TrigConfigSvc/TrigConfigSvc"
+
+
+    if (flags.LArNoise.outNtupLAr != ""):
+       noiseAlg=CompFactory.LArNoiseBursts("LArNoiseBursts")
+       noiseAlg.SigmaCut = flags.LArNoise.SigmaCut
+       noiseAlg.NumberOfBunchesInFront = flags.LArNoise.NumberOfBunchesInFront
+       noiseAlg.KeepOnlyCellID = flags.LArNoise.KeepOnlyCellID
+       result.addEventAlgo(noiseAlg)
+
+       result.addService(CompFactory.THistSvc(Output=["TTREE DATAFILE='"+flags.LArNoise.outNtupLAr+"' OPT='RECREATE'",]))
+
+    if (flags.LArNoise.HECNoiseNtup!=""):
+       hecAlg=CompFactory.LArHECNoise("LArHECNoise")
+       hecAlg.MinDigitADC = flags.LArNoise.MinDigitADC
+       hecAlg.MaxDeltaT = flags.LArNoise.MaxDeltaT
+       from IOVDbSvc.IOVDbSvcConfig import addFolders
+       result.merge(addFolders(flags, "/LAR/ElecCalibFlat/Pedestal", "LAR_ONL", className="CondAttrListCollection", db="CONDBR2"))
+       LArPedestalCondAlg =  CompFactory.getComp("LArFlatConditionsAlg<LArPedestalFlat>")
+       result.addCondAlgo(LArPedestalCondAlg(ReadKey="/LAR/ElecCalibFlat/Pedestal", WriteKey="LArPedestal"))
+
+       result.addEventAlgo(hecAlg)
+       if result.getService("THistSvc") is not None:
+           result.getService("THistSvc").Output += ["HEC DATAFILE='"+flags.LArNoise.HECNoiseNtup+"' OPT='RECREATE'",]
+       else:   
+           result.addService(CompFactory.THistSvc(Output=["HEC DATAFILE='"+flags.LArNoise.HECNoiseNtup+"' OPT='RECREATE'",]))
+
+    if (flags.LArNoise.outHistLAr != ""):
+       from LArCellRec.LArNoisyROSummaryConfig import LArNoisyROSummaryCfg
+       result.merge(LArNoisyROSummaryCfg(flags))
+       from LArMonitoring.LArNoisyROMonAlg import LArNoisyROMonConfig 
+       result.merge(LArNoisyROMonConfig(flags, inKey="LArNoisyROSummary"))
+
+    return result
+
+
+
+if __name__=="__main__":
+    
+    from AthenaConfiguration.AllConfigFlags import ConfigFlags
+    from LArNoiseFlags import addNoiseFlags
+    addNoiseFlags(ConfigFlags)
+    #ConfigFlags.Input.Files=['/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MetadataTests/data18/data18_13TeV.00363979.physics_Main.daq.ESD.0750._0001.pool.root']
+    ConfigFlags.Input.Files=['/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/RecExOnline/data16_13TeV.00302347.express_express.merge.RAW._lb0432._SFO-ALL._0001.1']
+
+    from AthenaConfiguration.MainServicesConfig import MainServicesCfg
+    
+    cfg=MainServicesCfg(ConfigFlags)
+    #cfg.merge(LArNoiseCfg(ConfigFlags))
+    cfg.merge(LArNoiseFromRawCfg(ConfigFlags))
+
+    cfg.run(10)