Prototype of new alignment configuration

Overview

This MR contains a prototype for configuring alignment jobs using a Moore-like configuration. Its main purpose is documenting the work done during past hackathons and to serve as basis for the next steps. It includes track and vertex selections, adding constraints and defining the elements to be aligned. The fundamental change is making the Alignment project depend on Moore instead of Phys.

Some documentation how to use the new configuration is provided in Alignment/Humboldt/README.md

Next steps and open points are summarised in #17 and #20

Configurations

The configurations added are

Alignment/Humboldt/options/AlignVPHalvesModules.py #example Velo alignment job
Alignment/Humboldt/options/AlignVPHalvesModules_PrKalman.py #example Velo job using PrKalman (see #16)
Alignment/Humboldt/options/AlignFTStationsLayers.py #example SciFi job
Alignment/Humboldt/options/AlignVPHalvesModules_Analyzer.py #run analyser part of example alignment job
Alignment/Humboldt/options/AlignVPHalvesModules_Iterator.py #run iterator part of example alignment job

Helper functions to set up new configuration are in Alignment/Humboldt/python/Humboldt/utils.py.

Central track and particle selections are started to be defined in

Alignment/Humboldt/python/Humboldt/TrackSelections.py
Alignment/Humboldt/python/Humboldt/ParticleSelections.py

New functionality (WIP)

• add AlignIterator taking parts from OnlineAlignment to perform iterator part offline

Tests

• add Alignment/Humboldt/tests/qmtest/humboldt.qms/align-vp-halves-modules.qmt to test new configuration
• add Alignment/Humboldt/tests/qmtest/humboldt.qms/analyzer-vp-halves-modules.qmt to test Analyzer part
• add Alignment/Humboldt/tests/qmtest/humboldt.qms/iterator-vp-halves-modules.qmt to test Iterator part
• add Alignment/Humboldt/tests/qmtest/humboldt.qms/align-vp-halves-modules-splititer.qmt to test gaudisplititer.py (WIP)

Uses Rec!2495 (merged)

Alignment/Humboldt/python/Humboldt/utils.py

+###############################################################################
+# (c) Copyright 2021 CERN for the benefit of the LHCb Collaboration           #
+#                                                                             #
+# This software is distributed under the terms of the GNU General Public      #
+# Licence version 3 (GPL Version 3), copied verbatim in the file "COPYING".   #
+#                                                                             #
+# In applying this licence, CERN does not waive the privileges and immunities #
+# granted to it by virtue of its status as an Intergovernmental Organization  #
+# or submit itself to any jurisdiction.                                       #
+###############################################################################
+
+from PyConf.control_flow import CompositeNode, NodeLogic
+from PyConf.tonic import configurable
+from PyConf.Algorithms import AlignAlgorithm
+from PyConf.Tools import AlignUpdateTool, EigenDiagSolvTool, AlignChisqConstraintTool
+from PyConf.application import configure
+
+
+@configurable
+def createAlignAlgorithm(pvs,
+                         tracks,
+                         odin,
+                         particles,
+                         elements,
+                         updatetool,
+                         name='HumboldtAlignAlgorithm',
+                         updateFinalize=False,
+                         outputlevel=3,
+                         nIterations=1,
+                         usecorrelations=True,
+                         chi2Outlier=10000,
+                         histoPrint=False,
+                         outputDataFile='',
+                         onlineMode=False,
+                         runList=[],
+                         useLocalFrame=True,
+                         xmlWriters=[],
+                         inputDataFiles=[]):
+    alignment = AlignAlgorithm(
+        name=name,
+        VertexLocation=pvs,
+        ODINLocation=odin,
+        TracksLocation=tracks,
+        ParticleLocation=particles,
+        UpdateTool=updatetool,
+        Elements=list(elements),
+        UpdateInFinalize=updateFinalize,
+        OutputLevel=outputlevel,
+        NumberOfIterations=nIterations,
+        UseCorrelations=usecorrelations,
+        Chi2Outlier=chi2Outlier,
+        HistoPrint=histoPrint,
+        OutputDataFile=outputDataFile,
+        OnlineMode=onlineMode,
+        RunList=runList,
+        UseLocalFrame=useLocalFrame,
+        XmlWriters=xmlWriters,
+        InputDataFiles=inputDataFiles)
+    return alignment
+
+
+@configurable
+def createAlignUpdateTool(name='updateTool',
+                          surveyconstraintTool=AlignChisqConstraintTool(),
+                          matrixSolverTool=EigenDiagSolvTool(public=True),
+                          constraints=[],
+                          minNumberOfHits=10,
+                          usePreconditioning=True,
+                          logFile="alignlog.txt",
+                          useWeightedAverage=False):
+    updatetool = AlignUpdateTool(
+        name=name,
+        MatrixSolverTool=matrixSolverTool,
+        SurveyConstraintTool=surveyconstraintTool,
+        Constraints=constraints,
+        MinNumberOfHits=minNumberOfHits,
+        UsePreconditioning=usePreconditioning,
+        LogFile=logFile,
+        UseWeightedAverage=useWeightedAverage,
+        public=True)
+    return updatetool
+
+
+def runAlignment(options,
+                 chisqConstraintTool,
+                 lagrangeConstrains,
+                 alignmentTracks,
+                 alignmentPVs,
+                 particles,
+                 odin,
+                 elementsToAlign,
+                 monitorList=[]):
+    updatetool = createAlignUpdateTool(
+        surveyconstraintTool=chisqConstraintTool,
+        constraints=lagrangeConstrains)
+    alignment = createAlignAlgorithm(
+        pvs=alignmentPVs,
+        tracks=alignmentTracks,
+        odin=odin,
+        particles=particles,
+        elements=elementsToAlign,
+        updatetool=updatetool)
+    top_node = CompositeNode(
+        "alignment",
+        monitorList + [alignment],
+        combine_logic=NodeLogic.NONLAZY_OR,
+        force_order=True)
+    configure(options, top_node, public_tools=[])