MuonPhaseII/MuonDetDescr/MuonGeoModelTestR4/python/testGeoModel.py

Go to the documentation of this file.

# Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
from AthenaConfiguration.ComponentFactory import CompFactory
 
 
 
class MuonPhaseIITestDefaults:
    
    EVGEN_PG = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/EVGEN_ParticleGun_FourMuon_Pt10to500.root"]
    
    HITS_PG_R3 = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/R3SimHits.pool.root"]
    
    HITS_PG_R4 = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/R4SimHits.pool.root"]
    
    HITS_PG_R3_MSOnly = []
    
    HITS_PG_R4_MSOnly = []
    
    DATA_BS = [
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-11._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-12._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-13._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-14._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-15._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-16._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-17._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-18._0001.data",
        "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonRecRTT/data25_13p6TeV/data25_13p6TeV.00508073.physics_Main.daq.RAW._lb0277._SFO-19._0001.data"]
    
    RDO_R3 = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/myRDO.R3.pool.root"]
    
    RDO_R4 = ["/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/myRDO.R4.pool.root"]
    
 
    
    GEODB_R3 = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/ATLAS-R3S-2021-03-02-00.db"
    
    GEODB_R3MSOnly = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/ATLAS-R3S-2021-03-02-00_MSOnly.db"
    
    GEODB_R4 = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/ATLAS-P2-RUN4-01-00-00.db"
    
    GEODB_R4MSOnly = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/ATLAS-P2-RUN4-01-00-00_MSOnly.db"
    
    GEODB_ITk_R3MS = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/ATLAS-P2-RUN4-01-00-00_R3MS.db"
    
    GEODB_TOROID = "/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/MuonGeomRTT/GeoDB/MUON_TOROID.db"
 
def SetupArgParser():
    from argparse import ArgumentParser
 
    parser = ArgumentParser()
    parser.add_argument("--threads", type=int, help="number of threads", default=1)
    parser.add_argument("--inputFile", "-i", default= MuonPhaseIITestDefaults.EVGEN_PG, 
                        help="Input file to run on ", nargs="+")
    parser.add_argument("--geoModelFile", default = MuonPhaseIITestDefaults.GEODB_R3, help="GeoModel SqLite file containing the muon geometry.")
    parser.add_argument("--defaultGeoFile", help="Use the  predefined GeoModel files on cvmfs", choices=["NONE", "RUN3", "RUN4", 
                                                                                                         "RUN3MSOnly", "RUN4MSOnly", 
                                                                                                         "ITkR3MS" ], default="NONE")
    parser.add_argument("--chambers", default=["all"], nargs="+", help="Chambers to check. If string is all, all chambers will be checked")
    parser.add_argument("--excludedChambers", default=[], nargs="+", help="Chambers to exclude. If string contains 'none', all chambers will be checked. Note: adding a chamber to --excludedChambers will overwrite it being in --chambers.")
    parser.add_argument("--outRootFile", default="NewGeoModelDump.root", help="Output ROOT file to dump the geomerty")
    parser.add_argument("--nEvents", help="Number of events to run", type = int ,default = 1)
    parser.add_argument("--skipEvents", help="Number of events to skip", type = int, default = 0)
    parser.add_argument("--noMdt", help="Disable the Mdts from the geometry", action='store_true', default = False)
    parser.add_argument("--noRpc", help="Disable the Rpcs from the geometry", action='store_true', default = False)
    parser.add_argument("--noTgc", help="Disable the Tgcs from the geometry", action='store_true', default = False)
    parser.add_argument("--noMM", help="Disable the MMs from the geometry", action='store_true', default = False)
    parser.add_argument("--noSTGC", help="Disable the sTgcs from the geometry", action='store_true', default = False)
    parser.add_argument("--eventPrintoutLevel", type=int, help="Interval of event heartbeat printouts from the loop manager", default = 1)
    parser.add_argument("--localMdtMezzJSON", default="", help="")
    parser.add_argument("--localMdtCablingJSON", default="", help="")
    return parser
 
def setupServicesCfg(flags):
    from AthenaConfiguration.MainServicesConfig import MainServicesCfg
    result = MainServicesCfg(flags)
    
    from AthenaConfiguration.Enums import Format
    if flags.Input.Format is Format.POOL:
        from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
        result.merge(PoolReadCfg(flags))
    elif flags.Input.Format == Format.BS:
        from ByteStreamCnvSvc.ByteStreamConfig import ByteStreamReadCfg
        result.merge(ByteStreamReadCfg(flags)) 
 
    from PerfMonComps.PerfMonCompsConfig import PerfMonMTSvcCfg
    result.merge(PerfMonMTSvcCfg(flags))
    from MuonConfig.MuonGeometryConfig import MuonIdHelperSvcCfg
    result.merge(MuonIdHelperSvcCfg(flags))
    return result
 
def GeoModelMdtTestCfg(flags, name = "GeoModelMdtTest", localMezzanineJSON="", localCablingJSON="",**kwargs):
    result = ComponentAccumulator()
    from MuonConfig.MuonCablingConfig import MDTCablingConfigCfg
    result.merge(MDTCablingConfigCfg(flags,MezzanineJSON=localMezzanineJSON, CablingJSON=localCablingJSON))
    the_alg = CompFactory.MuonGMR4.GeoModelMdtTest(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def GeoModelRpcTestCfg(flags, name = "GeoModelRpcTest", **kwargs):
    result = ComponentAccumulator()
    the_alg = CompFactory.MuonGMR4.GeoModelRpcTest(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def GeoModelTgcTestCfg(flags, name = "GeoModelTgcTest", **kwargs):
    result = ComponentAccumulator()
    the_alg = CompFactory.MuonGMR4.GeoModelTgcTest(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def GeoModelsTgcTestCfg(flags, name = "GeoModelsTgcTest", **kwargs):
    result = ComponentAccumulator()
    the_alg = CompFactory.MuonGMR4.GeoModelsTgcTest(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def GeoModelMmTestCfg(flags, name = "GeoModelMmTest", **kwargs):
    result = ComponentAccumulator()
    the_alg = CompFactory.MuonGMR4.GeoModelMmTest(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def NswGeoPlottingAlgCfg(flags, name="NswGeoPlotting", **kwargs):
    result = ComponentAccumulator()
    kwargs.setdefault("TestActsSurface", True)
    kwargs.setdefault("plotTgc", flags.Detector.GeometryTGC)
    kwargs.setdefault("plotStgc", flags.Detector.GeometrysTGC)
    kwargs.setdefault("plotMm", flags.Detector.GeometryMM)
    
    the_alg = CompFactory.MuonGMR4.NswGeoPlottingAlg(name, **kwargs)
    result.addEventAlgo(the_alg, primary = True)
    return result
 
def configureDefaultTagsCfg(flags):    
    from AthenaCommon.Logging import logging
    log = logging.getLogger('GeometryConfiguration')
 
    if not flags.GeoModel.SQLiteDB:
        raise ValueError("Default tag configuration only works for SQLite")
    
    from AthenaConfiguration.TestDefaults import defaultGeometryTags
    flags.GeoModel.AtlasVersion = defaultGeometryTags.RUN2    
    from AthenaConfiguration.Enums import LHCPeriod
    if flags.GeoModel.Run == LHCPeriod.Run3:   
        flags.GeoModel.AtlasVersion = defaultGeometryTags.RUN3
    elif flags.GeoModel.Run == LHCPeriod.Run4:
          flags.GeoModel.AtlasVersion = defaultGeometryTags.RUN4
    else:
        raise ValueError(f"Invalid run period {flags.GeoModel.Run}")
    from MuonConfig.MuonConfigUtils import configureCondTag
    configureCondTag(flags)
 
    log.info(f"Setup {flags.GeoModel.AtlasVersion} geometry loading {flags.GeoModel.SQLiteDBFullPath}")
    log.info(f"Use conditions tag {flags.IOVDb.GlobalTag}")
    
 
def setupGeoR4TestCfg(args,  flags = None):
    
    if flags is None:
        from AthenaConfiguration.AllConfigFlags import initConfigFlags
        flags = initConfigFlags()
    flags.Concurrency.NumThreads = args.threads
    flags.Concurrency.NumConcurrentEvents = args.threads
    flags.Exec.MaxEvents = args.nEvents
    flags.Exec.SkipEvents = args.skipEvents
    from os import path, system, listdir
    flags.Input.Files = []
    
    for fileArg in args.inputFile:
        if path.isdir(fileArg):
            flags.Input.Files += [ "{dir}/{file}".format(dir=fileArg, file=y) for y in listdir(fileArg) ]
        else:
            if fileArg[fileArg.rfind(".")+1 :]not in ["txt", "conf"]:
                 flags.Input.Files+=[fileArg]
            else:
                with open(fileArg) as inStream:
                   flags.Input.Files+=[ line.strip() for line in inStream if line[0]!='#'] 
 
    flags.Exec.FPE= 500
    flags.Exec.EventPrintoutInterval = 500
    
    if args.defaultGeoFile == "RUN3":
        flags.GeoModel.SQLiteDBFullPath =  MuonPhaseIITestDefaults.GEODB_R3
    elif args.defaultGeoFile == "RUN4":
        flags.GeoModel.SQLiteDBFullPath =  MuonPhaseIITestDefaults.GEODB_R4
    elif args.defaultGeoFile == "RUN3MSOnly":
        flags.GeoModel.SQLiteDBFullPath = MuonPhaseIITestDefaults.GEODB_R3MSOnly
    elif args.defaultGeoFile == "RUN4MSOnly":
        flags.GeoModel.SQLiteDBFullPath = MuonPhaseIITestDefaults.GEODB_R4MSOnly
    elif args.defaultGeoFile == "ITkR3MS":
        flags.GeoModel.SQLiteDBFullPath = MuonPhaseIITestDefaults.GEODB_ITk_R3MS
    elif args.geoModelFile.startswith("root://"):
        if not path.exists("Geometry/{geoTag}.db".format(geoTag=args.geoTag)):
            print ("Copy geometry file from EOS {source}".format(source = args.geoModelFile))
            system("mkdir Geometry/")
            system("xrdcp {source} Geometry/{geoTag}.db".format(source = args.geoModelFile,
                                                                geoTag=args.geoTag))
                                
        args.geoModelFile = "Geometry/{geoTag}.db".format(geoTag=args.geoTag)
    else:
        flags.GeoModel.SQLiteDBFullPath = args.geoModelFile
 
    flags.GeoModel.SQLiteDB = True
    configureDefaultTagsCfg(flags)
    
    flags.Detector.GeometryBpipe = False
    
    flags.Detector.GeometryBCM = False
    flags.Detector.GeometryPixel = False
    flags.Detector.GeometrySCT = False
    flags.Detector.GeometryTRT = False
    
    flags.Detector.GeometryPLR = False
    flags.Detector.GeometryBCMPrime = False
    flags.Detector.GeometryITkPixel = False
    flags.Detector.GeometryITkStrip = False
    
    flags.Detector.GeometryHGTD = False
    
    flags.Detector.GeometryLAr = False
    flags.Detector.GeometryTile = False
    flags.Detector.GeometryMBTS = False
    flags.Detector.GeometryCalo = False
    
    flags.Detector.GeometryCSC = False
    if args.noSTGC:
        flags.Detector.GeometrysTGC = False
    if args.noMM:
        flags.Detector.GeometryMM = False
    if args.noTgc:
        flags.Detector.GeometryTGC = False
    if args.noRpc:
        flags.Detector.GeometryRPC = False    
    if args.noMdt:
        flags.Detector.GeometryMDT = False
 
    flags.Acts.TrackingGeometry.UseBlueprint = True
 
    flags.Scheduler.CheckDependencies = True
    flags.Scheduler.ShowDataDeps = True
    flags.Scheduler.ShowDataFlow = True
    flags.Scheduler.ShowControlFlow = True
    flags.Scheduler.EnableVerboseViews = True
    flags.Scheduler.AutoLoadUnmetDependencies = True
    #flags.PerfMon.doFullMonMT = True
    flags.lock()
    flags.dump(evaluate = True)
    cfg = setupServicesCfg(flags)
 
    from MuonConfig.MuonGeometryConfig import MuonGeoModelCfg
    cfg.merge(MuonGeoModelCfg(flags))
    
    if not flags.Muon.usePhaseIIGeoSetup:
        print ("WARNING: New Muon plugin is not part of the Geometry file {geoDBFile}".format(geoDBFile=args.geoModelFile))
    else:
        from ActsAlignmentAlgs.AlignmentAlgsConfig import ActsGeometryContextAlgCfg
        cfg.merge(ActsGeometryContextAlgCfg(flags))
 
    cfg.getService("MessageSvc").verboseLimit = 10000000
    cfg.getService("MessageSvc").debugLimit = 10000000
 
    return flags, cfg
 
def executeTest(cfg):
    cfg.printConfig(withDetails=True, summariseProps=True)
    if not cfg.run().isSuccess(): exit(1)
 
if __name__=="__main__":
    args = SetupArgParser().parse_args()
    flags, cfg = setupGeoR4TestCfg(args)  
    from MuonConfig.MuonConfigUtils import setupHistSvcCfg
    cfg.merge(setupHistSvcCfg(flags, outFile = args.outRootFile, outStream="GEOMODELTESTER"))
    chambToTest =  args.chambers if len([x for x in args.chambers if x =="all"]) ==0 else []
    chambToExclude = args.excludedChambers
    
    
    if flags.Muon.usePhaseIIGeoSetup:
        cfg.getCondAlgo("MuonDetectorCondAlg").checkGeo = True
        cfg.getCondAlgo("MuonDetectorCondAlg").dumpGeo = True
        from TrackingGeometryCondAlg.AtlasTrackingGeometryCondAlgConfig import TrackingGeometryCondAlgCfg
        cfg.merge(TrackingGeometryCondAlgCfg(flags))
    
    
    cfg.getService("MessageSvc").setVerbose = []
 
    if flags.Detector.GeometryMDT:
        if not flags.Muon.usePhaseIIGeoSetup:
            from MuonGeoModelTest.testGeoModel import GeoModelMdtTestCfg as LegacyTestCfg
            cfg.merge(LegacyTestCfg(flags,
                                         TestStations = [ch for ch in chambToTest if ch[0] == "B" or ch[0] == "E"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "B" or ch[0] == "E"]))
        else:
            cfg.merge(GeoModelMdtTestCfg(flags, 
                                         TestStations = [ch for ch in chambToTest if ch[0] == "B" or ch[0] == "E"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "B" or ch[0] == "E"],
                                         localMezzanineJSON=args.localMdtMezzJSON,
                                         localCablingJSON=args.localMdtCablingJSON,
                                         ReadoutSideXML="ReadoutSides.xml",
                                         ExtraInputs=[( 'MuonGM::MuonDetectorManager' , 'ConditionStore+MuonDetectorManager' )]))
 
    if flags.Detector.GeometryRPC: 
        if not flags.Muon.usePhaseIIGeoSetup:
            from MuonGeoModelTest.testGeoModel import GeoModelRpcTestCfg as LegacyTestCfg
            cfg.merge(LegacyTestCfg(flags,
                                         TestStations = [ch for ch in chambToTest if ch[0] == "B"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "B"]))
        else:
            cfg.merge(GeoModelRpcTestCfg(flags, 
                                         TestStations = [ch for ch in chambToTest if ch[0] == "B"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "B"],
                                         ExtraInputs=[( 'MuonGM::MuonDetectorManager' , 'ConditionStore+MuonDetectorManager' )]))
 
    if flags.Detector.GeometryTGC:
        if not flags.Muon.usePhaseIIGeoSetup:
            from MuonGeoModelTest.testGeoModel import GeoModelTgcTestCfg as LegacyTestCfg
            cfg.merge(LegacyTestCfg(flags,
                                         TestStations = [ch for ch in chambToTest if ch[0] == "T"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "T"]))
        else:
            cfg.merge(GeoModelTgcTestCfg(flags, 
                                         TestStations = [ch for ch in chambToTest if ch[0] == "T"],
                                         ExcludeStations = [ch for ch in chambToExclude if ch[0] == "T"],
                                         ExtraInputs=[( 'MuonGM::MuonDetectorManager' , 'ConditionStore+MuonDetectorManager' )]))
 
    if flags.Detector.GeometryMM:
        if not flags.Muon.usePhaseIIGeoSetup:
            from MuonGeoModelTest.testGeoModel import GeoModelMmTestCfg as LegacyTestCfg
            cfg.merge(LegacyTestCfg(flags,
                                        TestStations = [ch for ch in chambToTest if ch[0] == "M"],
                                        ExcludeStations = [ch for ch in chambToExclude if ch[0] == "M"])) 
        else:
            cfg.merge(GeoModelMmTestCfg(flags, 
                                        TestStations = [ch for ch in chambToTest if ch[0] == "M"],
                                        ExcludeStations = [ch for ch in chambToExclude if ch[0] == "M"],
                                        ExtraInputs=[( 'MuonGM::MuonDetectorManager' , 'ConditionStore+MuonDetectorManager' )]))
    
    if flags.Detector.GeometrysTGC:
        if not flags.Muon.usePhaseIIGeoSetup:
            from MuonGeoModelTest.testGeoModel import GeoModelsTgcTestCfg as LegacyTestCfg
            cfg.merge(LegacyTestCfg(flags,
                                        TestStations = [ch for ch in chambToTest if ch[0] == "S"],
                                        ExcludeStations = [ch for ch in chambToExclude if ch[0] == "S"]))  
        else:
            cfg.merge(GeoModelsTgcTestCfg(flags, 
                                          TestStations = [ch for ch in chambToTest if ch[0] == "S"],
                                          ExcludeStations = [ch for ch in chambToExclude if ch[0] == "S"],
                                          ExtraInputs=[( 'MuonGM::MuonDetectorManager' , 'ConditionStore+MuonDetectorManager' )]))
    
    executeTest(cfg)