OnTrackCalibrator.icc

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#include "ActsCalibrators/xAODUncalibMeasCalibrator.h"
 
 
namespace ActsTrk::detail {
  
  template <typename traj_t>
  OnTrackCalibrator<traj_t> OnTrackCalibrator<traj_t>::NoCalibration(const ActsTrk::ITrackingGeometryTool* trackGeoTool)
  {
    ToolHandle<IOnTrackCalibratorTool<traj_t>> null(nullptr);
    return OnTrackCalibrator<traj_t>(trackGeoTool, null, null, null);
  }
  
  template <typename traj_t>
  OnTrackCalibrator<traj_t>::OnTrackCalibrator(const ActsTrk::ITrackingGeometryTool* trackGeoTool,
                 const ToolHandle<IOnTrackCalibratorTool<traj_t>>& pixelTool,
                 const ToolHandle<IOnTrackCalibratorTool<traj_t>>& stripTool,
                 const ToolHandle<IOnTrackCalibratorTool<traj_t>>& hgtdTool)
    : m_surfAcc{trackGeoTool} {
    if (pixelTool.isEnabled()) {
      pixelTool->connect(*this);
    } else {
      pixelCalibrator.template connect<&OnTrackCalibrator<traj_t>::passthrough<2, xAOD::PixelCluster>>(this);
    }
    
    if (stripTool.isEnabled()) {
      stripTool->connect(*this);
    } else {
      stripCalibrator.template connect<&OnTrackCalibrator<traj_t>::passthrough<1, xAOD::StripCluster>>(this);
    }
    if (hgtdTool.isEnabled()) {
      hgtdTool->connect(*this);
    } else {
  // cppcheck-suppress missingReturn; false positive
        hgtdCalibrator.template connect<&OnTrackCalibrator<traj_t>::passthrough<3, xAOD::HGTDCluster>>(this);
    }
  }
  
  template <typename traj_t>
  void OnTrackCalibrator<traj_t>::calibrate(const Acts::GeometryContext& geoctx,
                                            const Acts::CalibrationContext& cctx,
                                            const Acts::SourceLink& link,
                                            TrackStateProxy state) const {
    
    const xAOD::UncalibratedMeasurement& measurement{*xAODUncalibMeasCalibrator::unpack(link)};
    // @TODO in principle only need to lookup surface for strips
    switch (measurement.type()) {
        using enum xAOD::UncalibMeasType;
        using enum ProjectorType;
        case PixelClusterType: {
          assert(pixelCalibrator.connected());
          auto [pos, cov] = pixelCalibrator(geoctx, cctx,
                                            static_cast<const xAOD::PixelCluster&>(measurement), state);
          setState<2, traj_t>(e2DimNoTime, pos, cov, link, state);
          break;
        } case StripClusterType: {
          assert(stripCalibrator.connected());
          auto [pos, cov] = stripCalibrator(geoctx, cctx,
                                            static_cast<const xAOD::StripCluster&>(measurement), state);
          const Acts::Surface *surface = m_surfAcc.get(&measurement);
          if (!surface) {
            throw std::runtime_error("OnTrackCalibrator<traj_t>::calibrate() -- No surface for measurement.");
          }
          setState<1, traj_t>(surface->bounds().type() == Acts::SurfaceBounds::eAnnulus ? e1DimRotNoTime : e1DimNoTime,
                              pos, cov, link, state);
          break;
        } case HGTDClusterType: {
          assert(hgtdCalibrator.connected());
          auto [pos, cov] = hgtdCalibrator(geoctx, cctx,
                                           static_cast<const xAOD::HGTDCluster&>(measurement), state);
          setState<3, traj_t>(e2DimWithTime, pos, cov, link, state);
          break;
        } default:
            throw std::domain_error("OnTrackCalibrator can only handle Pixel or Strip measurements");
      }
  }
  
  template <typename traj_t>
  template <std::size_t Dim, typename Cluster>
  std::pair<xAOD::MeasVector<Dim>, xAOD::MeasMatrix<Dim>>
  OnTrackCalibrator<traj_t>::passthrough(const Acts::GeometryContext& /*gctx*/,
                                         const Acts::CalibrationContext& /*cctx*/,
                                         const Cluster& cluster,
                                         const TrackStateProxy& /*state*/) const {
    return std::make_pair(cluster.template localPosition<Dim>(),
                          cluster.template localCovariance<Dim>());
  }
} // namespace ActsTrk