SingleCurvilinearTrackParameters.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2016-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

#include "Acts/EventData/SingleBoundTrackParameters.hpp"
#include "Acts/Surfaces/PlaneSurface.hpp"

namespace Acts {

template <typename charge_t>
class SingleCurvilinearTrackParameters
    : public SingleBoundTrackParameters<charge_t> {
  using Base = SingleBoundTrackParameters<charge_t>;

 public:
  using Scalar = BoundScalar;
  using ParametersVector = BoundVector;
  using CovarianceMatrix = BoundSymMatrix;

  SingleCurvilinearTrackParameters(
      const Vector4D& pos4, const Vector3D& dir, Scalar p, Scalar q,
      std::optional<CovarianceMatrix> cov = std::nullopt)
      : Base(Surface::makeShared<PlaneSurface>(pos4.segment<3>(ePos0), dir),
             detail::transformFreeToCurvilinearParameters(
                 pos4[eTime], dir, (q != Scalar(0)) ? (q / p) : (1 / p)),
             q, std::move(cov)) {
    assert((0 <= p) and "Absolute momentum must be positive");
  }

  template <typename T = charge_t,
            std::enable_if_t<std::is_default_constructible_v<T>, int> = 0>
  SingleCurvilinearTrackParameters(
      const Vector4D& pos4, const Vector3D& dir, Scalar qOverP,
      std::optional<CovarianceMatrix> cov = std::nullopt)
      : Base(Surface::makeShared<PlaneSurface>(pos4.segment<3>(ePos0), dir),
             detail::transformFreeToCurvilinearParameters(pos4[eTime], dir,
                                                          qOverP),
             std::move(cov)) {}

  SingleCurvilinearTrackParameters(
      const Vector4D& pos4, Scalar phi, Scalar theta, Scalar p, Scalar q,
      std::optional<CovarianceMatrix> cov = std::nullopt)
      : Base(Surface::makeShared<PlaneSurface>(
                 pos4.segment<3>(ePos0),
                 makeDirectionUnitFromPhiTheta(phi, theta)),
             detail::transformFreeToCurvilinearParameters(
                 pos4[eTime], phi, theta, (q != Scalar(0)) ? (q / p) : (1 / p)),
             q, std::move(cov)) {
    assert((0 <= p) and "Absolute momentum must be positive");
  }

  template <typename T = charge_t,
            std::enable_if_t<std::is_default_constructible_v<T>, int> = 0>
  SingleCurvilinearTrackParameters(
      const Vector4D& pos4, Scalar phi, Scalar theta, Scalar qOverP,
      std::optional<CovarianceMatrix> cov = std::nullopt)
      : Base(Surface::makeShared<PlaneSurface>(
                 pos4.segment<3>(ePos0),
                 makeDirectionUnitFromPhiTheta(phi, theta)),
             detail::transformFreeToCurvilinearParameters(pos4[eTime], phi,
                                                          theta, qOverP),
             std::move(cov)) {}

  // this class does not have a custom default constructor and thus should not
  // provide any custom default cstors, dstor, or assignment. see ISOCPP C.20.
};

}  // namespace Acts