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/Geometry/GeometryContext.hpp"
#include "Acts/Geometry/GeometryIdentifier.hpp"
#include "Acts/MagneticField/MagneticFieldContext.hpp"
#include "Acts/Material/MaterialInteraction.hpp"
#include "Acts/Material/SurfaceMaterialMapper.hpp"
#include "Acts/Material/VolumeMaterialMapper.hpp"
#include "Acts/Utilities/Logger.hpp"
#include "ActsExamples/Framework/DataHandle.hpp"
#include "ActsExamples/Framework/IAlgorithm.hpp"
#include "ActsExamples/Framework/ProcessCode.hpp"
#include "ActsExamples/MaterialMapping/IMaterialWriter.hpp"

#include <climits>
#include <cstddef>
#include <cstdint>
#include <functional>
#include <map>
#include <memory>
#include <mutex>
#include <string>
#include <unordered_map>
#include <utility>
#include <vector>

namespace ActsExamples {
class IMaterialWriter;
struct AlgorithmContext;
}  // namespace ActsExamples

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 IAlgorithm {
 public:
  struct Config {
    // 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;

    std::string collection = "material_tracks";

    std::string mappingMaterialCollection = "mapped_material_tracks";

    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;
  };

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

  ~MaterialMapping() override;

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

  std::vector<std::pair<double, int>> scoringParameters(uint64_t surfaceID);

  const Config& config() const { return m_cfg; }

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

  ReadDataHandle<std::unordered_map<size_t, Acts::RecordedMaterialTrack>>
      m_inputMaterialTracks{this, "InputMaterialTracks"};
  WriteDataHandle<std::unordered_map<size_t, Acts::RecordedMaterialTrack>>
      m_outputMaterialTracks{this, "OutputMaterialTracks"};
};

}  // namespace ActsExamples