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/AlignTrTools/src/TrackDoubleHitPruner.cpp

-/*****************************************************************************\
-* (c) Copyright 2000-2018 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.                                       *
-\*****************************************************************************/
-// Gaudi
-#include "GaudiKernel/ToolHandle.h"
-
-#include "GaudiAlg/GaudiAlgorithm.h"
-
-// track interfaces
-#include "TrackInterfaces/ITrackKalmanFilter.h"
-#include "TrackKernel/TrackFunctors.h"
-
-// event
-#include "Event/FitNode.h"
-#include "Event/Track.h"
-#include "Event/TrackFitResult.h"
-
-#include <string>
-
-class TrackDoubleHitPruner : public GaudiAlgorithm {
-
-public:
-  // Constructors and destructor
-  TrackDoubleHitPruner( const std::string& name, ISvcLocator* pSvcLocator );
-  virtual ~TrackDoubleHitPruner();
-
-  StatusCode execute() override;
-  StatusCode initialize() override;
-  StatusCode finalize() override;
-
-private:
-  std::string                    m_inputLocation;
-  ToolHandle<ITrackKalmanFilter> m_kalmanFilter;
-};
-
-DECLARE_COMPONENT( TrackDoubleHitPruner )
-
-TrackDoubleHitPruner::TrackDoubleHitPruner( const std::string& name, ISvcLocator* pSvcLocator )
-    : GaudiAlgorithm( name, pSvcLocator ), m_kalmanFilter( "TrackKalmanFilter", this ) {
-  // constructor
-  declareProperty( "TrackLocation", m_inputLocation = LHCb::TrackLocation::Default );
-}
-
-TrackDoubleHitPruner::~TrackDoubleHitPruner() {
-  // destructor
-}
-
-StatusCode TrackDoubleHitPruner::initialize() {
-  StatusCode sc = GaudiAlgorithm::initialize();
-  sc            = m_kalmanFilter.retrieve();
-  return sc;
-}
-
-StatusCode TrackDoubleHitPruner::finalize() {
-  m_kalmanFilter.release().ignore();
-  return GaudiAlgorithm::finalize();
-}
-
-bool sameOTMonoLayer( const LHCb::OTChannelID& rhs, const LHCb::OTChannelID& lhs ) {
-  return rhs.sequentialUniqueLayer() == lhs.sequentialUniqueLayer() && rhs.straw() / 32 == lhs.straw() / 32;
-}
-
-bool sameSTModule( const LHCb::STChannelID& rhs, const LHCb::STChannelID& lhs ) {
-  return rhs.station() == lhs.station() && rhs.layer() == lhs.layer() && rhs.detRegion() == lhs.detRegion() &&
-         rhs.sector() == lhs.sector();
-}
-
-bool sameVeloSensor( const LHCb::VeloChannelID& rhs, const LHCb::VeloChannelID& lhs ) {
-  return rhs.sensor() == lhs.sensor();
-}
-
-StatusCode TrackDoubleHitPruner::execute() {
-  // get them as a range .... hen const-cast:-)
-  LHCb::Track::Range tracks = get<LHCb::Track::Range>( m_inputLocation );
-  for ( const LHCb::Track* track : tracks ) {
-    // we assume that the nodes are sorted. we now check if there are
-    // neighbouring nodes in the same detector layer. this is rather
-    // detector specific.
-    if ( track->fitResult() ) {
-      const LHCb::FitNode* prevnode( 0 );
-      bool                 trackWasModified( false );
-      for ( const LHCb::FitNode* node : nodes( *track ) ) {
-        if ( node->type() == LHCb::FitNode::Type::HitOnTrack ) {
-          if ( prevnode ) {
-            LHCb::LHCbID idA = prevnode->measurement().lhcbID();
-            LHCb::LHCbID idB = node->measurement().lhcbID();
-            if ( idA.detectorType() == idB.detectorType() ) {
-              bool samelayer( false );
-              if ( idA.isVelo() )
-                samelayer = sameVeloSensor( idA.veloID(), idB.veloID() );
-              else if ( idA.isTT() || idA.isIT() )
-                samelayer = sameSTModule( idA.stID(), idB.stID() );
-              // if( idA.isOT() )
-              // samelayer = sameOTMonoLayer( idA.otID(), idB.otID()) ;
-
-              if ( samelayer ) {
-                if ( std::abs( node->residual() ) > std::abs( prevnode->residual() ) )
-                  const_cast<LHCb::FitNode*>( node )->setType( LHCb::FitNode::Type::Outlier );
-                else
-                  const_cast<LHCb::FitNode*>( prevnode )->setType( LHCb::FitNode::Type::Outlier );
-                trackWasModified = true;
-              }
-            }
-          }
-          prevnode = node;
-        } else {
-          prevnode = 0;
-        }
-      }
-      counter( "Fraction of modified tracks" ) += trackWasModified;
-      if ( trackWasModified ) {
-        LHCb::Track* nctrack = const_cast<LHCb::Track*>( track );
-        StatusCode   sc      = m_kalmanFilter->fit( *nctrack );
-        if ( sc.isFailure() ) nctrack->setFitStatus( LHCb::Track::FitStatus::FitFailed );
-      }
-    }
-  }
-  return StatusCode::SUCCESS;
-}