testMuonGen1.py

Go to the documentation of this file.

#!/usr/bin/env python
# Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
"""
  Run ACTS geometry construction for Muon
"""
from AthenaConfiguration.AllConfigFlags import initConfigFlags
flags = initConfigFlags()
 
from AthenaConfiguration.TestDefaults import defaultGeometryTags, defaultConditionsTags
flags.GeoModel.AtlasVersion = defaultGeometryTags.RUN4
flags.IOVDb.GlobalTag = defaultConditionsTags.RUN4_MC
flags.GeoModel.Align.Dynamic = False
flags.Input.Files = ['/cvmfs/atlas-nightlies.cern.ch/repo/data/data-art/SimCoreTests/valid1.410000.PowhegPythiaEvtGen_P2012_ttbar_hdamp172p5_nonallhad.evgen.EVNT.e4993.EVNT.08166201._000012.pool.root.1']
 
flags.Detector.GeometryITkPixel = False
flags.Detector.GeometryITkStrip = False
flags.Detector.GeometryBpipe = False
flags.Detector.GeometryCalo = False
flags.Detector.GeometryMuon = True
 
flags.Concurrency.NumThreads = 64
flags.Concurrency.NumConcurrentEvents = 64
 
flags.Exec.MaxEvents = 100
 
from MuonGeoModelTestR4.testGeoModel import geoModelFileDefault, configureDefaultTagsCfg
flags.GeoModel.SQLiteDBFullPath = geoModelFileDefault(useR4Layout = True)
 
flags.GeoModel.SQLiteDB = True
configureDefaultTagsCfg(flags)
flags.lock()
flags.dump()
 
from AthenaConfiguration.MainServicesConfig import MainServicesCfg
acc = MainServicesCfg( flags )
 
from ActsConfig.ActsGeometryConfig import ActsExtrapolationAlgCfg, ActsTrackingGeometrySvcCfg
 
from AthenaCommon.Constants import INFO
tgSvc = ActsTrackingGeometrySvcCfg(flags,
                                   OutputLevel=INFO,
                                   RunConsistencyChecks=False,
                                   #  ConsistencyCheckOutput="trk_geo_check.csv", # enable debug output writing
                                   ObjDebugOutput=True)
 
acc.merge(tgSvc)
 
alg = ActsExtrapolationAlgCfg(flags,
                              OutputLevel=INFO,
                              NParticlesPerEvent = int(100),
                              WritePropStep = True,
                              EtaRange = [-5, 5],
                              PtRange = [20, 100])
 
acc.merge(alg)
acc.printConfig()
sc = acc.run()
 
if sc.isFailure():
    import sys
    sys.exit(1)