SpacePointFormation2/detail/StripSpacePointCalibrationImpl.hpp

0001 // This file is part of the ACTS project.
0002 //
0003 // Copyright (C) 2016 CERN for the benefit of the ACTS project
0004 //
0005 // This Source Code Form is subject to the terms of the Mozilla Public
0006 // License, v. 2.0. If a copy of the MPL was not distributed with this
0007 // file, You can obtain one at https://mozilla.org/MPL/2.0/.
0008
0009 #pragma once
0010
0011 #include "Acts/EventData/StripSpacePointCalibrationDetails.hpp"
0012 #include "Acts/Utilities/detail/StdArrayLinalg.hpp"
0013
0014 #include <cmath>
0015
0016 namespace Acts::detail {
0017
0018 // Intentionally not using Eigen here as there is a noticeable performance
0019 // (10-20%) penalty for small vector sizes.
0020
0021 inline OuterStripSpacePointCalibrationDetailsDerived
0022 deriveOuterStripSpacePointCalibrationDetails(const std::array<float, 3>& ihv,
0023                                              const std::array<float, 3>& ohv,
0024                                              const std::array<float, 3>& iosv,
0025                                              const std::array<float, 3>& oc) {
0026   OuterStripSpacePointCalibrationDetailsDerived result{};
0027   result.innerCrossOuterHalfVector = stdArrayCross(ihv, ohv);
0028   result.innerToOuterSeparationCrossOuterHalfVector = stdArrayCross(iosv, ohv);
0029   result.innerToOuterSeparationCrossInnerHalfVector = stdArrayCross(iosv, ihv);
0030   result.outerCenter = oc;
0031   result.outerHalfVector = ohv;
0032   return result;
0033 }
0034
0035 inline OuterStripSpacePointCalibrationDetailsDerived
0036 deriveOuterStripSpacePointCalibrationDetails(
0037     const OuterStripSpacePointCalibrationDetails& sp) {
0038   return deriveOuterStripSpacePointCalibrationDetails(
0039       sp.innerHalfVector, sp.outerHalfVector, sp.innerToOuterSeparation,
0040       sp.outerCenter);
0041 }
0042
0043 inline bool calibrateOuterStripSpacePoint(
0044     const std::array<float, 3>& direction,
0045     const std::array<float, 3>& ihvCrossOhv,
0046     const std::array<float, 3>& iosvCrossOhv,
0047     const std::array<float, 3>& iosvCrossIhv, const std::array<float, 3>& oc,
0048     const std::array<float, 3>& ohv, std::array<float, 3>& calibrated,
0049     const float tolerance) {
0050   // scale = innerStripHalfVector dot (outerStripHalfVector cross direction)
0051   const float scale = stdArrayDot(direction, ihvCrossOhv);
0052
0053   // sInner = innerToOuterSeparationVector dot (outerStripHalfVector cross
0054   // direction) Check if direction is inside the inner detector element
0055   const float sInner = stdArrayDot(direction, iosvCrossOhv);
0056   if (std::abs(sInner) > std::abs(scale) * tolerance) {
0057     return false;
0058   }
0059
0060   // sOuter = innerToOuterSeparationVector dot (innerStripHalfVector cross
0061   // direction) Check if direction is inside the outer detector element
0062   const float sOuter = stdArrayDot(direction, iosvCrossIhv);
0063   if (std::abs(sOuter) > std::abs(scale) * tolerance) {
0064     return false;
0065   }
0066
0067   // Corrected position using the outer strip center and direction
0068   const float sOuterNorm = sOuter / scale;
0069   calibrated = stdArrayAddScaled(oc, ohv, sOuterNorm);
0070   return true;
0071 }
0072
0073 inline bool calibrateOuterStripSpacePoint(
0074     const std::array<float, 3>& direction,
0075     const OuterStripSpacePointCalibrationDetailsDerived& sp,
0076     std::array<float, 3>& calibrated, const float tolerance) {
0077   return detail::calibrateOuterStripSpacePoint(
0078       direction, sp.innerCrossOuterHalfVector,
0079       sp.innerToOuterSeparationCrossOuterHalfVector,
0080       sp.innerToOuterSeparationCrossInnerHalfVector, sp.outerCenter,
0081       sp.outerHalfVector, calibrated, tolerance);
0082 }
0083
0084 }  // namespace Acts::detail