ConeSurface.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file ConeSurface.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/Geometry/GeometryContext.hpp"
#include "Acts/Geometry/Polyhedron.hpp"
#include "Acts/Surfaces/ConeBounds.hpp"
#include "Acts/Surfaces/Surface.hpp"
#include "Acts/Utilities/Definitions.hpp"
#include "Acts/Utilities/ParameterDefinitions.hpp"
#include "Acts/Utilities/detail/RealQuadraticEquation.hpp"

namespace Acts {


class ConeSurface : public Surface {
  friend Surface;

 protected:
  ConeSurface(const Transform3D& transform, double alpha,
              bool symmetric = false);

  ConeSurface(const Transform3D& transform, double alpha, double zmin,
              double zmax, double halfPhi = M_PI);

  ConeSurface(const Transform3D& transform,
              const std::shared_ptr<const ConeBounds>& cbounds);

  ConeSurface(const ConeSurface& other);

  ConeSurface(const GeometryContext& gctx, const ConeSurface& other,
              const Transform3D& shift);

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

  ConeSurface& operator=(const ConeSurface& other);

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

  SurfaceType type() const override;

  RotationMatrix3D referenceFrame(const GeometryContext& gctx,
                                  const Vector3D& position,
                                  const Vector3D& momentum) const final;

  Vector3D normal(const GeometryContext& gctx,
                  const Vector2D& lposition) const final;

  Vector3D normal(const GeometryContext& gctx,
                  const Vector3D& position) const final;

  using Surface::normal;

  // Return method for the rotational symmetry axis
  // @return This returns the local z axis
  virtual Vector3D rotSymmetryAxis(const GeometryContext& gctx) const;

  const ConeBounds& bounds() const final;

  Vector3D localToGlobal(const GeometryContext& gctx, const Vector2D& lposition,
                         const Vector3D& momentum) const final;

  Result<Vector2D> globalToLocal(const GeometryContext& gctx,
                                 const Vector3D& position,
                                 const Vector3D& momentum) const final;

  SurfaceIntersection intersect(const GeometryContext& gctx,
                                const Vector3D& position,
                                const Vector3D& direction,
                                const BoundaryCheck& bcheck) const final;

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

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

  std::string name() const override;

  LocalCartesianToBoundLocalMatrix localCartesianToBoundLocalDerivative(
      const GeometryContext& gctx, const Vector3D& position) const final;

 protected:
  std::shared_ptr<const ConeBounds> m_bounds;  

 private:
  detail::RealQuadraticEquation intersectionSolver(
      const GeometryContext& gctx, const Vector3D& position,
      const Vector3D& direction) const;
};

#include "Acts/Surfaces/detail/ConeSurface.ipp"

}  // namespace Acts