DiscSurface.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file DiscSurface.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/Definitions/Algebra.hpp"
#include "Acts/Definitions/Tolerance.hpp"
#include "Acts/Definitions/TrackParametrization.hpp"
#include "Acts/Geometry/GeometryContext.hpp"
#include "Acts/Geometry/Polyhedron.hpp"
#include "Acts/Surfaces/BoundaryCheck.hpp"
#include "Acts/Surfaces/DiscBounds.hpp"
#include "Acts/Surfaces/InfiniteBounds.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/BinningType.hpp"
#include "Acts/Utilities/Result.hpp"

#include <cmath>
#include <cstddef>
#include <memory>
#include <string>

namespace Acts {

class DetectorElementBase;
class DiscBounds;
class SurfaceBounds;

class DiscSurface : public Surface {
#ifndef DOXYGEN
  friend Surface;
#endif

 protected:
  DiscSurface(const Transform3& transform, double rmin, double rmax,
              double hphisec = M_PI);

  DiscSurface(const Transform3& transform, double minhalfx, double maxhalfx,
              double minR, double maxR, double avephi = 0., double stereo = 0.);

  DiscSurface(const Transform3& transform,
              std::shared_ptr<const DiscBounds> dbounds = nullptr);

  DiscSurface(std::shared_ptr<const DiscBounds> dbounds,
              const DetectorElementBase& detelement);

  DiscSurface(const DiscSurface& other);

  DiscSurface(const GeometryContext& gctx, const DiscSurface& other,
              const Transform3& shift);

 public:
  ~DiscSurface() override = default;
  DiscSurface() = delete;

  DiscSurface& operator=(const DiscSurface& other);

  SurfaceType type() const override;

  Vector3 normal(const GeometryContext& gctx,
                 const Vector2& lposition) const final;

  using Surface::normal;

  Vector3 binningPosition(const GeometryContext& gctx,
                          BinningValue bValue) const final;

  const SurfaceBounds& bounds() const final;

  Vector3 localToGlobal(const GeometryContext& gctx, const Vector2& lposition,
                        const Vector3& direction) const final;

  Result<Vector2> globalToLocal(
      const GeometryContext& gctx, const Vector3& position,
      const Vector3& direction,
      double tolerance = s_onSurfaceTolerance) const final;

  Vector2 localPolarToCartesian(const Vector2& lpolar) const;

  Vector2 localCartesianToPolar(const Vector2& lcart) const;

  Vector2 localPolarToLocalCartesian(const Vector2& locpol) const;

  Vector3 localCartesianToGlobal(const GeometryContext& gctx,
                                 const Vector2& lposition) const;

  Vector2 globalToLocalCartesian(const GeometryContext& gctx,
                                 const Vector3& position,
                                 double tol = 0.) const;

  BoundToFreeMatrix boundToFreeJacobian(
      const GeometryContext& gctx, const BoundVector& boundParams) const final;

  FreeToBoundMatrix freeToBoundJacobian(
      const GeometryContext& gctx, const FreeVector& parameters) const final;

  double pathCorrection(const GeometryContext& gctx, const Vector3& position,
                        const Vector3& direction) const final;

  SurfaceMultiIntersection intersect(
      const GeometryContext& gctx, const Vector3& position,
      const Vector3& direction, const BoundaryCheck& bcheck = false,
      ActsScalar tolerance = s_onSurfaceTolerance) const final;

  double binningPositionValue(const GeometryContext& gctx,
                              BinningValue bValue) const final;

  std::string name() const override;

  Polyhedron polyhedronRepresentation(const GeometryContext& gctx,
                                      size_t lseg) const override;

  ActsMatrix<2, 3> localCartesianToBoundLocalDerivative(
      const GeometryContext& gctx, const Vector3& position) const final;

 protected:
  std::shared_ptr<const DiscBounds> m_bounds;  
};
}  // end of namespace Acts