FacesHelper.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 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/Polyhedron.hpp"
#include "Acts/Utilities/Definitions.hpp"

#include <numeric>
#include <utility>
#include <vector>

namespace Acts {

namespace detail {

struct FacesHelper {
  using FaceVector = std::vector<Polyhedron::FaceType>;

  static std::pair<FaceVector, FaceVector> convexFaceMesh(
      const std::vector<Vector3D>& vertices, bool centerLast = false) {
    FaceVector faces;
    FaceVector triangularMesh;
    // Write the face
    unsigned int offset = centerLast ? 1 : 0;
    std::vector<size_t> face(vertices.size() - offset);
    std::iota(face.begin(), face.end(), 0);
    faces.push_back(face);
    unsigned int anker = centerLast ? vertices.size() - 1 : 0;
    for (unsigned int it = 2 - offset; it < vertices.size() - offset; ++it) {
      triangularMesh.push_back({anker, it - 1, it});
    }
    // Close for centered reference point
    if (centerLast) {
      triangularMesh.push_back({anker, vertices.size() - 2, 0});
    }
    return {faces, triangularMesh};
  }

  static std::pair<FaceVector, FaceVector> cylindricalFaceMesh(
      const std::vector<Vector3D>& vertices, bool fullTwoPi = true) {
    FaceVector faces;
    FaceVector triangularMesh;
    size_t nqfaces = 0.5 * vertices.size();
    size_t reduce = (not fullTwoPi) ? 1 : 0;
    for (size_t iface = 0; iface < nqfaces - reduce; ++iface) {
      size_t p2 = (iface + 1 == nqfaces) ? 0 : iface + 1;
      std::vector<size_t> face = {iface, p2, p2 + nqfaces, nqfaces + iface};
      faces.push_back(face);
      std::vector<size_t> triA = {iface, p2, p2 + nqfaces};
      triangularMesh.push_back(triA);
      std::vector<size_t> triB = {p2 + nqfaces, nqfaces + iface, iface};
      triangularMesh.push_back(triB);
    }
    return {faces, triangularMesh};
  }
};

}  // namespace detail

}  // namespace Acts