CovarianceEngine.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file CovarianceEngine.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2019-2020 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

// Workaround for building on clang+libstdc++
#include "Acts/Utilities/detail/ReferenceWrapperAnyCompat.hpp"

#include "Acts/Definitions/Algebra.hpp"
#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/EventData/TrackParameters.hpp"
#include "Acts/EventData/detail/CorrectedTransformationFreeToBound.hpp"
#include "Acts/Geometry/GeometryContext.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Result.hpp"

#include <cmath>
#include <functional>
#include <tuple>
#include <variant>

namespace Acts {

namespace detail {

Result<std::tuple<BoundTrackParameters, BoundMatrix, double>> boundState(
    const GeometryContext& geoContext, BoundSquareMatrix& covarianceMatrix,
    BoundMatrix& jacobian, FreeMatrix& transportJacobian,
    FreeVector& derivatives, BoundToFreeMatrix& jacToGlobal,
    FreeVector& parameters, const ParticleHypothesis& particleHypothesis,
    bool covTransport, double accumulatedPath, const Surface& surface,
    const FreeToBoundCorrection& freeToBoundCorrection =
        FreeToBoundCorrection(false));

std::tuple<CurvilinearTrackParameters, BoundMatrix, double> curvilinearState(
    BoundSquareMatrix& covarianceMatrix, BoundMatrix& jacobian,
    FreeMatrix& transportJacobian, FreeVector& derivatives,
    BoundToFreeMatrix& jacToGlobal, const FreeVector& parameters,
    const ParticleHypothesis& particleHypothesis, bool covTransport,
    double accumulatedPath);

void transportCovarianceToBound(
    const GeometryContext& geoContext, BoundSquareMatrix& boundCovariance,
    BoundMatrix& fullTransportJacobian, FreeMatrix& freeTransportJacobian,
    FreeVector& freeToPathDerivatives, BoundToFreeMatrix& boundToFreeJacobian,
    FreeVector& freeParameters, const Surface& surface,
    const FreeToBoundCorrection& freeToBoundCorrection =
        FreeToBoundCorrection(false));

void transportCovarianceToCurvilinear(BoundSquareMatrix& boundCovariance,
                                      BoundMatrix& fullTransportJacobian,
                                      FreeMatrix& freeTransportJacobian,
                                      FreeVector& freeToPathDerivatives,
                                      BoundToFreeMatrix& boundToFreeJacobian,
                                      const Vector3& direction);

}  // namespace detail
}  // namespace Acts