CylinderSurface.hpp

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

#include <cmath>

namespace Acts {


class CylinderSurface : public Surface {
  friend Surface;

 protected:
  CylinderSurface(std::shared_ptr<const CylinderBounds> cbounds,
                  const DetectorElementBase& detelement);

  CylinderSurface(const Transform3D& transform, double radius, double halfz,
                  double halfphi = M_PI, double avphi = 0.);

  CylinderSurface(const Transform3D& transform,
                  const std::shared_ptr<const CylinderBounds>& cbounds);

  CylinderSurface(const CylinderSurface& other);

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

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

  CylinderSurface& operator=(const CylinderSurface& other);

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

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

  SurfaceType type() const override;

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

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

  using Surface::normal;

  virtual Vector3D rotSymmetryAxis(const GeometryContext& gctx) const;

  const CylinderBounds& 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;

  std::string name() const override;

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

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

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

 private:
  detail::RealQuadraticEquation intersectionSolver(
      const Transform3D& transform, const Vector3D& position,
      const Vector3D& direction) const;
};

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

}  // namespace Acts