xAODUncalibMeasSurfAcc.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
#include "ActsCalibrators/xAODUncalibMeasSurfAcc.h"
#include "ActsCalibrators/xAODUncalibMeasCalibrator.h"
#include "ActsGeometry/DetectorElementToActsGeometryIdMap.h"
#include "ActsGeometry/SurfaceOfMeasurementUtil.h"
#include "xAODAuxiliaryMeasurement/AuxiliaryMeasurement.h"
#include "xAODMuonPrepData/UtilFunctions.h"
 
namespace ActsTrk::detail{
    xAODUncalibMeasSurfAcc::xAODUncalibMeasSurfAcc(const ActsTrk::ITrackingGeometrySvc* trackGeoSvc):
        m_actsTrackingGeometry{trackGeoSvc->trackingGeometry().get()},
        m_detectorElementToGeometryIdMap{trackGeoSvc->surfaceIdMap()}{}
           
    const Acts::Surface* xAODUncalibMeasSurfAcc::operator()(const Acts::SourceLink& sourceLink) const {
        return get(xAODUncalibMeasCalibrator::unpack(sourceLink));
    }
    const Acts::Surface* xAODUncalibMeasSurfAcc::get(const xAOD::UncalibratedMeasurement* meas) const {
        switch (meas->type()) {
            using enum xAOD::UncalibMeasType;
            case PixelClusterType:
            case StripClusterType:
            case HGTDClusterType:{
                assert(m_detectorElementToGeometryIdMap);
                assert(m_actsTrackingGeometry);
                const auto geoKey = makeDetectorElementKey(meas->type(), meas->identifierHash());
                const auto geoid_iter = m_detectorElementToGeometryIdMap->find(geoKey);
                if (geoid_iter == m_detectorElementToGeometryIdMap->end()) {
                    return nullptr;
                }
                if (const Acts::Surface* surface = DetectorElementToActsGeometryIdMap::getSurface(*geoid_iter)) {
                    return surface;
                }
                return m_actsTrackingGeometry->findSurface( DetectorElementToActsGeometryIdMap::getValue(*geoid_iter));
            }
            case MdtDriftCircleType:
            case RpcStripType:
            case TgcStripType:
            case sTgcStripType:
            case MMClusterType:{
                return &xAOD::muonSurface(meas);
                break;
            }case Other: {
                const auto* pMeas = static_cast<const xAOD::AuxiliaryMeasurement*>(meas);
                return pMeas->surface().get();
            }
 
            default:
                break;
              
        }
#ifndef NDEBUG
        throw std::domain_error("xAODUncalibMeasSurfAcc() - Cannot decode surface type");
#endif        
        return nullptr;
    }
}