MaterialGridHelper.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file MaterialGridHelper.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/Material/AccumulatedVolumeMaterial.hpp"
#include "Acts/Material/MaterialSlab.hpp"
#include "Acts/Utilities/BinUtility.hpp"
#include "Acts/Utilities/Definitions.hpp"
#include "Acts/Utilities/detail/Axis.hpp"
#include "Acts/Utilities/detail/Grid.hpp"

#include <stdexcept>

namespace Acts {

using RecordedMaterialVolumePoint =
    std::vector<std::pair<Acts::MaterialSlab, std::vector<Acts::Vector3D>>>;

using EAxis = Acts::detail::EquidistantAxis;
using Grid2D =
    Acts::detail::Grid<Acts::AccumulatedVolumeMaterial, EAxis, EAxis>;
using Grid3D =
    Acts::detail::Grid<Acts::AccumulatedVolumeMaterial, EAxis, EAxis, EAxis>;
using MaterialGrid2D =
    Acts::detail::Grid<Acts::Material::ParametersVector, EAxis, EAxis>;
using MaterialGrid3D =
    Acts::detail::Grid<Acts::Material::ParametersVector, EAxis, EAxis, EAxis>;

Grid2D createGrid(std::array<double, 3> gridAxis1,
                  std::array<double, 3> gridAxis2);

Grid3D createGrid(std::array<double, 3> gridAxis1,
                  std::array<double, 3> gridAxis2,
                  std::array<double, 3> gridAxis3);

std::function<double(Acts::Vector3D)> globalToLocalFromBin(
    Acts::BinningValue& type);

Grid2D createGrid2D(
    const BinUtility& bins,
    std::function<Acts::Vector2D(Acts::Vector3D)>& transfoGlobalToLocal);

Grid3D createGrid3D(
    const BinUtility& bins,
    std::function<Acts::Vector3D(Acts::Vector3D)>& transfoGlobalToLocal);

MaterialGrid2D mapMaterialPoints(
    Grid2D& grid, const Acts::RecordedMaterialVolumePoint& mPoints,
    std::function<Acts::Vector2D(Acts::Vector3D)>& transfoGlobalToLocal);

MaterialGrid3D mapMaterialPoints(
    Grid3D& grid, const Acts::RecordedMaterialVolumePoint& mPoints,
    std::function<Acts::Vector3D(Acts::Vector3D)>& transfoGlobalToLocal);

}  // namespace Acts