CreateMisalignmentITk.py

Go to the documentation of this file.

#!/usr/bin/env python
 
#Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
 
 
from AthenaConfiguration.AllConfigFlags import initConfigFlags
from AthenaConfiguration.MainServicesConfig import MainServicesCfg    
from InDetAlignGenAlgs.InDetAlignAlgsConfig import CreateITkMisalignAlgCfg
from AthenaConfiguration.DetectorConfigFlags import setupDetectorFlags
from InDetAlignGenTools.InDetAlignGenToolsConfig import ITkAlignDBTool
import sys
from AthenaConfiguration.ComponentFactory import CompFactory
 
 
 
def getFlags(**kwargs):
    flags=initConfigFlags()
    
    flags.Input.isMC             = True
 
    flags.Input.Files = []
    
    from AthenaConfiguration.TestDefaults import defaultGeometryTags,defaultConditionsTags
    flags.GeoModel.AtlasVersion = defaultGeometryTags.RUN4
    flags.GeoModel.Align.Dynamic = False
    
    flags.ITk.Geometry.isAlignable = True
    flags.ITk.Align.alignITkPixel = True
    flags.ITk.Align.alignITkStrip = True
 
 
 
    flags.ITk.Geometry.AllLocal = False
    detectors = [
    "ITkPixel",
    "ITkStrip",
    "Bpipe"
    ]
    setupDetectorFlags(flags, custom_list=detectors, toggle_geometry=True)
    flags.TrackingGeometry.MaterialSource = "Input"
 
 
    #Define the output database file name and add it to the flags
    if 'MisalignMode' not in kwargs.keys():
        MisalignMode = 11 # Radial
    else:
        MisalignMode=int(kwargs.get('MisalignMode',11))
    if MisalignMode == 2:
        translation_scale_str = kwargs.pop('Translation_Scale', '0.1,0.1,0.1')
        rotation_scale_str = kwargs.pop('Rotation_Scale', '0.1,0.1,0.1')
        Translation_Scale = [float(x) for x in translation_scale_str.split(',')]
        Rotation_Scale = [float(x) for x in rotation_scale_str.split(',')]
    if MisalignMode == 7:
        local_translation_str = kwargs.pop('Local_Translation', '0.0,0.0,0.0')
        local_rotation_str = kwargs.pop('Local_Rotation', '0.0,0.0,0.0')
        Local_Translation = [float(x) for x in local_translation_str.split(',')]
        Local_Rotation = [float(x) for x in local_rotation_str.split(',')] 
    databaseFilename     = 'MisalignmentSet%s.db' % (MisalignMode)
    flags.IOVDb.DBConnection="sqlite://;schema=%s;dbname=OFLCOND" % (databaseFilename) 
    flags.IOVDb.GlobalTag = defaultConditionsTags.RUN4_MC
 
 
    # This should run serially for the moment.
    flags.Concurrency.NumThreads = 1
    flags.Concurrency.NumConcurrentEvents = 1
    flags.lock()
 
    return flags
 
def CreateMis(flags,name="CreateITkMisalignAlg",**kwargs):
    #
    #--------------------------------------------------------------
    # Geometry section
    #--------------------------------------------------------------
    misalignmentOnTopOfExistingSet = False
    readDBPoolFile = False
 
    createFreshDB = not(readDBPoolFile or misalignmentOnTopOfExistingSet)
 
    MisalignMode = kwargs.pop('MisalignMode',11)
 
    shiftInMicrons = 100
 
    if MisalignMode in [11, 12, 31]:
        shiftInMicrons = 500
 
    outFiles = 'MisalignmentSet%s' % (MisalignMode)
    misalignModeMap = {0:'no Misalignment',
                    1: 'misalignment by 6 parameters',
                    2: 'random misalignment',
                    3: 'IBL-stave temperature dependent bowing',
                    7: 'misalignment according to module indices',
                    11: 'R deltaR (radial expansion)', 12: 'Phi deltaR (ellipse)',13: 'Z deltaR (funnel)',
                    21: 'R deltaPhi (curl)', 22: 'Phi deltaPhi (clamshell) ',23:'Z deltaPhi (twist)',
                    31: 'R deltaZ (telescope)',32:'Phi deltaZ (skew)',33:'Z deltaZ (z-expansion)'}
    
     
    local_translation_str = kwargs.pop('Local_Translation', '0.0,0.0,0.0')
    local_rotation_str = kwargs.pop('Local_Rotation', '0.0,0.0,0.0')
    Local_Translation = [float(x) for x in local_translation_str.split(',')]
    Local_Rotation = [float(x) for x in local_rotation_str.split(',')] 
    translation_scale_str = kwargs.pop('Translation_Scale', '0.1,0.1,0.1')
    rotation_scale_str = kwargs.pop('Rotation_Scale', '0.1,0.1,0.1')
    Translation_Scale = [float(x) for x in translation_scale_str.split(',')]
    Rotation_Scale = [float(x) for x in rotation_scale_str.split(',')]
    acc=MainServicesCfg(flags)
    print ("\n CreateMisalignAlg: Creation of misalignment mode %s: %s \n" % (int(MisalignMode),misalignModeMap.get(int(MisalignMode),'unknown')))
    kwargs.setdefault("ASCIIFilenameBase",outFiles)
    kwargs.setdefault("SQLiteTag",'MisalignmentMode_'+str(misalignModeMap.get(int(MisalignMode),'unknown')))
    kwargs.setdefault("MisalignMode",int(MisalignMode))
    kwargs.setdefault("Translation_Scale",Translation_Scale)
    kwargs.setdefault("Rotation_Scale",Rotation_Scale)
    kwargs.setdefault("Local_Translation",Local_Translation)
    kwargs.setdefault("Local_Rotation",Local_Rotation)
    kwargs.setdefault("Index", str(kwargs.pop('Index', '')))
    kwargs.setdefault("MaxShift",shiftInMicrons)
    kwargs.setdefault("CreateFreshDB",createFreshDB)
    #Create and configure the AlignDB tool
    
    outputAlignFolder="/Indet/AlignITk" #The folder name from which initial misalignments should be read - currently may clash if this is the same is the input folder!
 
    writeDBPoolFile=True   #Activate or deactivate writing to outFiles + '.pool.root'
    kargsTool={}
    kargsTool.setdefault("SCTTwoSide",True)
    kargsTool.setdefault("DBRoot",outputAlignFolder)
    kargsTool.setdefault("forceUserDBConfig",True)
    if writeDBPoolFile:
        print("Writing DB Pool File")
        InDetCondStream=CompFactory.AthenaOutputStreamTool("CondStream_write",OutputFile = outFiles+'.pool.root')
        InDetCondStream.PoolContainerPrefix="<type>"
        InDetCondStream.TopLevelContainerName=""
        InDetCondStream.SubLevelBranchName="<key>"
        kargsTool.setdefault("CondStream",InDetCondStream)
    dbTool = acc.popToolsAndMerge(ITkAlignDBTool(flags,**kargsTool))
 
    kwargs.setdefault("IDAlignDBTool",dbTool)
 
    cfg=CreateITkMisalignAlgCfg(flags,name=name,**kwargs)
 
    acc.merge(cfg)
    if writeDBPoolFile:
        print("To be writen DB pool File")
 
    return acc
 
if __name__ == "__main__":
    if(len(sys.argv[1:])):
        kwargs=dict(arg.split('=') for arg in sys.argv[1:])
        print(kwargs)
    else:
        kwargs={}
        print("no args")
    #Get the basic configuration flags
    flags=getFlags(**kwargs)
    flags.dump()
    #Add the tools and the algorithm to the accumulator
    acc=CreateMis(flags,**kwargs)
    acc.printConfig()
 
    #run
    sc=acc.run(10)
    if sc.isFailure():
        print("Failed to run the Misalignment Algorithm")
        sys.exit(-1)