PerigeeSurface.cpp

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

// 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/.

#include "Acts/Surfaces/PerigeeSurface.hpp"

#include <iomanip>
#include <iostream>
#include <utility>

Acts::PerigeeSurface::PerigeeSurface(const Vector3D& gp)
    : LineSurface(Transform3D(Translation3D(gp.x(), gp.y(), gp.z())), nullptr) {
}

Acts::PerigeeSurface::PerigeeSurface(const Transform3D& transform)
    : GeometryObject(), LineSurface(transform) {}

Acts::PerigeeSurface::PerigeeSurface(const PerigeeSurface& other)
    : GeometryObject(), LineSurface(other) {}

Acts::PerigeeSurface::PerigeeSurface(const GeometryContext& gctx,
                                     const PerigeeSurface& other,
                                     const Transform3D& shift)
    : GeometryObject(), LineSurface(gctx, other, shift) {}

Acts::PerigeeSurface& Acts::PerigeeSurface::operator=(
    const PerigeeSurface& other) {
  if (this != &other) {
    LineSurface::operator=(other);
  }
  return *this;
}

Acts::Surface::SurfaceType Acts::PerigeeSurface::type() const {
  return Surface::Perigee;
}

std::string Acts::PerigeeSurface::name() const {
  return "Acts::PerigeeSurface";
}

std::ostream& Acts::PerigeeSurface::toStream(const GeometryContext& gctx,
                                             std::ostream& sl) const {
  sl << std::setiosflags(std::ios::fixed);
  sl << std::setprecision(7);
  sl << "Acts::PerigeeSurface:" << std::endl;
  const Vector3D& sfCenter = center(gctx);
  sl << "     Center position  (x, y, z) = (" << sfCenter.x() << ", "
     << sfCenter.y() << ", " << sfCenter.z() << ")";
  sl << std::setprecision(-1);
  return sl;
}

Acts::Polyhedron Acts::PerigeeSurface::polyhedronRepresentation(
    const GeometryContext& gctx, size_t /*lseg*/) const {
  // Prepare vertices and faces
  std::vector<Vector3D> vertices;
  std::vector<Polyhedron::FaceType> faces;
  std::vector<Polyhedron::FaceType> triangularMesh;

  const Transform3D& ctransform = transform(gctx);
  Vector3D left(0, 0, -100.);
  Vector3D right(0, 0, 100.);

  // The central wire/straw
  vertices.push_back(ctransform * left);
  vertices.push_back(ctransform * right);
  faces.push_back({0, 1});
  vertices.push_back(ctransform * Vector3D(0., 0., 0.));
  triangularMesh.push_back({0, 2, 1});

  return Polyhedron(vertices, faces, triangularMesh);
}