track_parametrization.hpp

// This file is part of the ACTS project.
//
// Copyright (C) 2016 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 https://mozilla.org/MPL/2.0/.
 
#pragma once
 
// Project include(s)
#include "detray/definitions/algebra.hpp"
#include "detray/utils/concepts.hpp"
 
namespace detray {

enum bound_indices : unsigned int {
  // Local position on the reference surface.
  // This is intentionally named different from the position components in
  // the other data vectors, to clarify that this is defined on a surface
  // while the others are defined in free space.
  e_bound_loc0 = 0u,
  e_bound_loc1 = 1u,
  // Direction angles
  e_bound_phi = 2u,
  e_bound_theta = 3u,
  // Global inverse-momentum-like parameter, i.e. q/p or 1/p
  // The naming is inconsistent for the case of neutral track parameters where
  // the value is interpreted as 1/p not as q/p. This is intentional to avoid
  // having multiple aliases for the same element and for lack of an
  // acceptable
  // common name.
  e_bound_qoverp = 4u,
  e_bound_time = 5u,
  // Last uninitialized value contains the total number of components
  e_bound_size,
};

enum free_indices : unsigned int {
  // Spatial position
  // The spatial position components must be stored as one continuous block.
  e_free_pos0 = 0u,
  e_free_pos1 = e_free_pos0 + 1u,
  e_free_pos2 = e_free_pos0 + 2u,
  // Time
  e_free_time = 3u,
  // (Unit) direction
  // The direction components must be stored as one continuous block.
  e_free_dir0 = 4u,
  e_free_dir1 = e_free_dir0 + 1u,
  e_free_dir2 = e_free_dir0 + 2u,
  // Global inverse-momentum-like parameter, i.e. q/p or 1/p
  // See BoundIndices for further information
  e_free_qoverp = 7u,
  // Last uninitialized value contains the total number of components
  e_free_size,
};

template <concepts::algebra algebra_t>
using free_vector = dmatrix<algebra_t, e_free_size, 1>;

template <concepts::algebra algebra_t>
using free_matrix = dmatrix<algebra_t, e_free_size, e_free_size>;

template <concepts::algebra algebra_t>
using bound_vector = dmatrix<algebra_t, e_bound_size, 1>;

template <concepts::algebra algebra_t>
using bound_matrix = dmatrix<algebra_t, e_bound_size, e_bound_size>;

template <concepts::algebra algebra_t>
using bound_to_free_matrix = dmatrix<algebra_t, e_free_size, e_bound_size>;

template <concepts::algebra algebra_t>
using bound_to_free_jacobian_submatrix = dmatrix<algebra_t, 3, 2>;

template <concepts::algebra algebra_t>
using free_to_bound_matrix = dmatrix<algebra_t, e_bound_size, e_free_size>;

template <concepts::algebra algebra_t>
using free_to_bound_jacobian_submatrix = dmatrix<algebra_t, 2, 3>;

template <concepts::algebra algebra_t>
using free_to_path_matrix = dmatrix<algebra_t, 1, e_free_size>;

template <concepts::algebra algebra_t>
using path_to_free_matrix = dmatrix<algebra_t, e_free_size, 1>;
 
}  // namespace detray