CorrectedTransformationFreeToBound.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2021-2022 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/Definitions/Algebra.hpp"
#include "Acts/Definitions/Direction.hpp"
#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/EventData/detail/TransformationFreeToBound.hpp"
#include "Acts/Geometry/GeometryContext.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Logger.hpp"

#include <optional>
#include <tuple>

namespace Acts {

struct FreeToBoundCorrection {
  bool apply = false;

  ActsScalar alpha = 0.1;
  ActsScalar beta = 2;

  ActsScalar cosIncidentAngleMinCutoff = 1e-5;
  ActsScalar cosIncidentAngleMaxCutoff = 0.99500417;

  FreeToBoundCorrection() = default;

  FreeToBoundCorrection(bool apply_, ActsScalar alpha_, ActsScalar beta_);

  explicit FreeToBoundCorrection(bool apply_);

  operator bool() const;
};

namespace detail {

struct CorrectedFreeToBoundTransformer {
  CorrectedFreeToBoundTransformer(ActsScalar alpha, ActsScalar beta,
                                  ActsScalar cosIncidentAngleMinCutoff,
                                  ActsScalar cosIncidentAngleMaxCutoff);

  CorrectedFreeToBoundTransformer(
      const FreeToBoundCorrection& freeToBoundCorrection);

  CorrectedFreeToBoundTransformer() = default;
  CorrectedFreeToBoundTransformer(const CorrectedFreeToBoundTransformer&) =
      default;
  CorrectedFreeToBoundTransformer(CorrectedFreeToBoundTransformer&&) = default;
  CorrectedFreeToBoundTransformer& operator=(
      const CorrectedFreeToBoundTransformer&) = default;
  CorrectedFreeToBoundTransformer& operator=(
      CorrectedFreeToBoundTransformer&&) = default;

  std::optional<std::tuple<BoundVector, BoundSquareMatrix>> operator()(
      const FreeVector& freeParams, const FreeSquareMatrix& freeCovariance,
      const Surface& surface, const GeometryContext& geoContext,
      Direction navDir = Direction::Forward,
      const Logger& logger = getDummyLogger()) const;

 private:
  ActsScalar m_alpha = 0.1;
  ActsScalar m_beta = 2;

  ActsScalar m_cosIncidentAngleMinCutoff = 1e-5;

  ActsScalar m_cosIncidentAngleMaxCutoff = 0.99500417;
};

}  // namespace detail
}  // namespace Acts