RootMaterialDecorator.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file RootMaterialDecorator.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 "ActsExamples/Framework/ProcessCode.hpp"
#include <Acts/Definitions/Algebra.hpp>
#include <Acts/Geometry/GeometryIdentifier.hpp>
#include <Acts/Geometry/TrackingVolume.hpp>
#include <Acts/Material/IMaterialDecorator.hpp>
#include <Acts/Material/ISurfaceMaterial.hpp>
#include <Acts/Material/IVolumeMaterial.hpp>
#include <Acts/Surfaces/Surface.hpp>
#include <Acts/Utilities/Logger.hpp>

#include <map>
#include <memory>
#include <mutex>
#include <string>
#include <utility>

class TFile;

namespace Acts {
class ISurfaceMaterial;
class IVolumeMaterial;

using SurfaceMaterialMap =
    std::map<GeometryIdentifier, std::shared_ptr<const ISurfaceMaterial>>;
using VolumeMaterialMap =
    std::map<GeometryIdentifier, std::shared_ptr<const IVolumeMaterial>>;
}  // namespace Acts

namespace ActsExamples {

class RootMaterialDecorator : public Acts::IMaterialDecorator {
 public:
  class Config {
   public:
    std::string folderSurfaceNameBase = "SurfaceMaterial";
    std::string folderVolumeNameBase = "VolumeMaterial";
    std::string voltag = "_vol";
    std::string boutag = "_bou";
    std::string laytag = "_lay";
    std::string apptag = "_app";
    std::string sentag = "_sen";
    std::string ntag = "n";
    std::string vtag = "v";
    std::string otag = "o";
    std::string mintag = "min";
    std::string maxtag = "max";
    std::string ttag = "t";
    std::string x0tag = "x0";
    std::string l0tag = "l0";
    std::string atag = "A";
    std::string ztag = "Z";
    std::string rhotag = "rho";
    std::string fileName = "material-maps.root";
  };

  RootMaterialDecorator(const Config& config, Acts::Logging::Level level);

  ~RootMaterialDecorator() override;

  void decorate(Acts::Surface& surface) const final {
    // Null out the material for this surface
    if (m_clearSurfaceMaterial) {
      surface.assignSurfaceMaterial(nullptr);
    }
    // Try to find the surface in the map
    auto sMaterial = m_surfaceMaterialMap.find(surface.geometryId());
    if (sMaterial != m_surfaceMaterialMap.end()) {
      surface.assignSurfaceMaterial(sMaterial->second);
    }
  }

  void decorate(Acts::TrackingVolume& volume) const final {
    // Null out the material for this volume
    if (m_clearSurfaceMaterial) {
      volume.assignVolumeMaterial(nullptr);
    }
    // Try to find the surface in the map
    auto vMaterial = m_volumeMaterialMap.find(volume.geometryId());
    if (vMaterial != m_volumeMaterialMap.end()) {
      volume.assignVolumeMaterial(vMaterial->second);
    }
  }

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

 private:
  Config m_cfg;

  std::unique_ptr<const Acts::Logger> m_logger{nullptr};

  TFile* m_inputFile{nullptr};

  Acts::SurfaceMaterialMap m_surfaceMaterialMap;

  Acts::VolumeMaterialMap m_volumeMaterialMap;

  bool m_clearSurfaceMaterial{true};

  const Acts::Logger& logger() const { return *m_logger; }
};

}  // namespace ActsExamples