MaterialMapUtils.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2019 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/Definitions/Algebra.hpp"
#include "Acts/Definitions/Units.hpp"
#include "Acts/Material/InterpolatedMaterialMap.hpp"
#include "Acts/Material/Material.hpp"
#include "Acts/Utilities/detail/AxisFwd.hpp"
#include "Acts/Utilities/detail/Grid.hpp"

#include <array>
#include <cstddef>
#include <functional>
#include <vector>

// Convenience functions to ease creation of and Acts::InterpolatedMaterialMap
// and to avoid code duplication. Currently implemented for the two most common
// formats: rz and xyz.

namespace Acts {

class Material;

MaterialMapper<detail::Grid<Material::ParametersVector, detail::EquidistantAxis,
                            detail::EquidistantAxis>>
materialMapperRZ(const std::function<size_t(std::array<size_t, 2> binsRZ,
                                            std::array<size_t, 2> nBinsRZ)>&
                     materialVectorToGridMapper,
                 std::vector<double> rPos, std::vector<double> zPos,
                 const std::vector<Acts::Material>& material,
                 double lengthUnit = UnitConstants::mm);

MaterialMapper<detail::Grid<Material::ParametersVector, detail::EquidistantAxis,
                            detail::EquidistantAxis, detail::EquidistantAxis>>
materialMapperXYZ(const std::function<size_t(std::array<size_t, 3> binsXYZ,
                                             std::array<size_t, 3> nBinsXYZ)>&
                      materialVectorToGridMapper,
                  std::vector<double> xPos, std::vector<double> yPos,
                  std::vector<double> zPos,
                  const std::vector<Material>& material,
                  double lengthUnit = UnitConstants::mm);

}  // namespace Acts