extractMaxGeometryElements.h

Go to the documentation of this file.

#ifndef LARCOREALG_GEOMETRY_DETAILS_EXTRACTMAXGEOMETRYELEMENTS_H
#define LARCOREALG_GEOMETRY_DETAILS_EXTRACTMAXGEOMETRYELEMENTS_H

// LArSoft libraries
#include "larcorealg/Geometry/CryostatGeo.h"
#include "larcorealg/Geometry/PlaneGeo.h"
#include "larcorealg/Geometry/TPCGeo.h"
#include "larcorealg/Geometry/WireReadoutGeom.h"

// C/C++ standard libraries
#include <array>
#include <cstddef> // std::size_t
#include <vector>

namespace geo::details {

  template <std::size_t Levels = 4U>
  static std::array<unsigned int, Levels> extractMaxGeometryElements(
    std::vector<CryostatGeo> const& Cryostats,
    WireReadoutGeom const& wireReadoutGeom);

} // namespace geo::details

// -----------------------------------------------------------------------------
// --- template implementation
// -----------------------------------------------------------------------------
template <std::size_t Levels /* = 4U */>
std::array<unsigned int, Levels> geo::details::extractMaxGeometryElements(
  std::vector<CryostatGeo> const& Cryostats,
  WireReadoutGeom const& wireReadoutGeom)
{
  static_assert(Levels > 0U);
  static_assert(Levels <= 4U);

  std::array<unsigned int, Levels> maxElements;
  maxElements.fill(0U);

  auto setMax = [&maxElements](std::size_t index, unsigned int value) {
    if (maxElements[index] < value) maxElements[index] = value;
  };

  setMax(0U, Cryostats.size());
  if constexpr (Levels > 1U) {
    for (CryostatGeo const& cryo : Cryostats) {
      setMax(1U, cryo.NTPC());
      if constexpr (Levels > 2U) {
        for (TPCGeo const& TPC : cryo.IterateTPCs()) {
          setMax(2U, wireReadoutGeom.Nplanes(TPC.ID()));
          if constexpr (Levels > 3U) {
            for (PlaneGeo const& plane : wireReadoutGeom.Iterate<PlaneGeo>(TPC.ID())) {
              setMax(3U, plane.Nwires());
            } // for planes
          }   // if do wires
        }     // for TPCs
      }       // if do planes
    }         // for cryostats
  }           // if do TPCs

  return maxElements;
} // geo::details::extractMaxGeometryElements()

// -----------------------------------------------------------------------------

#endif // LARCOREALG_GEOMETRY_DETAILS_EXTRACTMAXGEOMETRYELEMENTS_H