MaterialMapping.hpp

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

// This file is part of the Acts project.
//
// Copyright (C) 2017-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/Material/SurfaceMaterialMapper.hpp"
#include "Acts/Material/VolumeMaterialMapper.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/Framework/BareAlgorithm.hpp"
#include "ActsExamples/MaterialMapping/IMaterialWriter.hpp"

#include <climits>
#include <memory>
#include <mutex>

namespace Acts {

class TrackingGeometry;
class ISurfaceMaterial;
class IVolumeMaterial;

using SurfaceMaterialMap =
    std::map<GeometryIdentifier, std::shared_ptr<const ISurfaceMaterial>>;

using VolumeMaterialMap =
    std::map<GeometryIdentifier, std::shared_ptr<const IVolumeMaterial>>;

using DetectorMaterialMaps = std::pair<SurfaceMaterialMap, VolumeMaterialMap>;
}  // namespace Acts

namespace ActsExamples {

class MaterialMapping : public ActsExamples::BareAlgorithm {
 public:
  struct Config {
    Config(std::reference_wrapper<const Acts::GeometryContext> gctx,
           std::reference_wrapper<const Acts::MagneticFieldContext> mctx)
        : geoContext(gctx), magFieldContext(mctx) {}

    std::string collection = "material-tracks";

    std::string mappingMaterialCollection = "MappedMaterialTracks";

    std::shared_ptr<Acts::SurfaceMaterialMapper> materialSurfaceMapper =
        nullptr;

    std::shared_ptr<Acts::VolumeMaterialMapper> materialVolumeMapper = nullptr;

    std::vector<std::shared_ptr<IMaterialWriter>> materialWriters;

    std::shared_ptr<const Acts::TrackingGeometry> trackingGeometry = nullptr;

    // Geometry context for the state creation
    std::reference_wrapper<const Acts::GeometryContext> geoContext;

    // MagneticField  context for the state creation
    std::reference_wrapper<const Acts::MagneticFieldContext> magFieldContext;
  };

  MaterialMapping(const Config& cfg,
                  Acts::Logging::Level level = Acts::Logging::INFO);

  ~MaterialMapping();

  ActsExamples::ProcessCode execute(
      const AlgorithmContext& context) const final override;

 private:
  Config m_cfg;  
  Acts::SurfaceMaterialMapper::State
      m_mappingState;  
  Acts::VolumeMaterialMapper::State
      m_mappingStateVol;  
};

}  // namespace ActsExamples