TrkMeasurementCalibrator.icc

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
#include "ActsGeometry/ATLASSourceLink.h"
#include "TrkSurfaces/Surface.h"
#include "TrkSurfaces/SurfaceBounds.h"
 
#include "InDetMeasurementUtilities/ClusterConversionUtilities.h"
#include "TrkRIO_OnTrack/RIO_OnTrack.h"
#include "InDetRIO_OnTrack/SCT_ClusterOnTrack.h"
#include "InDetRIO_OnTrack/PixelClusterOnTrack.h"
 
#include "HGTD_RIO_OnTrack/HGTD_ClusterOnTrack.h"
 
#include <format>
namespace ActsTrk::detail {
 
  template <typename trajectory_t>
  void TrkMeasurementCalibrator::calibrate(const Acts::GeometryContext& /*gctx*/,
                                           const Acts::CalibrationContext& /*cctx*/,
                                           const Acts::SourceLink& sl,
                                           TrackState_t<trajectory_t> trackState) const {
    const SourceLink_t sourceLink = unpack(sl);
 
    auto detEl = dynamic_cast<const InDetDD::SiDetectorElement *>(sourceLink->associatedSurface().associatedDetectorElement()); 
 
    if(auto pixelCluster = dynamic_cast<const InDet::PixelClusterOnTrack *>(sourceLink); pixelCluster != nullptr) {
      auto [locPos, locCov] = TrackingUtilities::convertPix_LocalPosCov(*pixelCluster->prepRawData());
      setState<2, trajectory_t>(ProjectorType::e2DimNoTime, locPos, locCov, sl, trackState);
    } else if(auto stripCluster = dynamic_cast<const InDet::SCT_ClusterOnTrack *>(sourceLink); stripCluster != nullptr) {
      auto [locPos, locCov] = TrackingUtilities::convertSCT_LocalPosCov(*stripCluster->prepRawData());
      if( detEl->isBarrel() ) {
        setState<1, trajectory_t>(ProjectorType::e1DimNoTime, locPos, locCov,
                                  sl, trackState);
      } else {
        setState<1, trajectory_t>(ProjectorType::e1DimRotNoTime, locPos, locCov,
                                  sl, trackState);
      }
    } else if(auto hgtdCluster = dynamic_cast<const HGTD_ClusterOnTrack*>(sourceLink); hgtdCluster != nullptr ) {
      auto [locPos, locCov] = TrackingUtilities::convertHGTD_LocalPosCov(*hgtdCluster->prepRawData());
      setState<3, trajectory_t>(ProjectorType::e2DimWithTime, locPos, locCov, sl, trackState);
     } else {
      std::stringstream linkDump{};
      sourceLink->dump(linkDump); 
      throw std::domain_error(
        std::format(
          "TrkMeasurementCalibrator::calibrate() - Invalid dimension {:d} for measurement {:} ",
          sourceLink->localParameters().dimension(), linkDump.str()
        )
      );
    }
  }
} // namespace ActsTrk::detail