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 <Acts/Geometry/GeometryID.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/Definitions.hpp>
#include <Acts/Utilities/Logger.hpp>
#include <map>
#include <mutex>

#include "ACTFW/Framework/ProcessCode.hpp"

class TFile;

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

namespace FW {

class RootMaterialDecorator : public Acts::IMaterialDecorator {
 public:
  class Config {
   public:
    std::string folderNameBase = "Material";
    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";
    std::shared_ptr<const Acts::Logger> logger;
    // The name of the writer
    std::string name = "";

    Config(const std::string& lname = "MaterialReader",
           Acts::Logging::Level lvl = Acts::Logging::INFO)
        : logger(Acts::getDefaultLogger(lname, lvl)), name(lname) {}
  };

  RootMaterialDecorator(const Config& cfg);

  ~RootMaterialDecorator();

  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.geoID());
    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.geoID());
    if (vMaterial != m_volumeMaterialMap.end()) {
      volume.assignVolumeMaterial(vMaterial->second);
    }
  }

 private:
  Config m_cfg;

  TFile* m_inputFile{nullptr};

  Acts::SurfaceMaterialMap m_surfaceMaterialMap;

  Acts::VolumeMaterialMap m_volumeMaterialMap;

  bool m_clearSurfaceMaterial{true};

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

}  // namespace FW