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(std::shared_ptr<const Transform3D> htrans, double alpha,
              bool symmetric = false);

  ConeSurface(std::shared_ptr<const Transform3D> htrans, double alpha,
              double zmin, double zmax, double halfPhi = M_PI);

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

  ConeSurface(const ConeSurface& other);

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

 public:
  ~ConeSurface() override = default;

  ConeSurface() = delete;

  ConeSurface& operator=(const ConeSurface& other);

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

  SurfaceType type() const override;

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

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

  const 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 const Vector3D rotSymmetryAxis(const GeometryContext& gctx) const;

  const ConeBounds& bounds() const final;

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

  bool globalToLocal(const GeometryContext& gctx, const Vector3D& position,
                     const Vector3D& momentum, Vector2D& lposition) const final;

  Intersection intersectionEstimate(
      const GeometryContext& gctx, const Vector3D& position,
      const Vector3D& direction,
      const BoundaryCheck& bcheck = false) 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;

 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