CbmRecoSetupUtils.cxx

Go to the documentation of this file.

/* Copyright (C) 2025 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Sergei Zharko [committer] */

 
#include "CbmRecoSetupUtils.h"
 
#include <Logger.h>
 
#include <TGeoArb8.h>
#include <TGeoManager.h>
#include <TGeoNode.h>
#include <TString.h>
 
#include <iomanip>
#include <limits>
 
#include <fmt/format.h>
using cbm::RecoSetupUtils;
 
// ---------------------------------------------------------------------------------------------------------------------
//
cbm::algo::GeoVolume RecoSetupUtils::ReadVolume(const TString& path)
{
  if (!gGeoManager->cd(path)) {
    std::stringstream msg;
    msg << "Node with path " << path << " is not found in the geo manager";
    throw std::runtime_error(msg.str());
  }
 
  // Get a shape of the bounding box
  const auto* bBox = dynamic_cast<const TGeoBBox*>(gGeoManager->GetCurrentVolume()->GetShape());
  if (!bBox) {
    throw std::runtime_error(fmt::format("The shape of the volume under path {} is not a TGeoBBox", path.Data()));
  }
 
  std::array<std::array<double, 3>, 8> locBoxVertices = {
    {// Vertex coordinates in the local frame of the box
     {bBox->GetOrigin()[0] + bBox->GetDX(), bBox->GetOrigin()[1] + bBox->GetDY(), bBox->GetOrigin()[2] + bBox->GetDZ()},
     {bBox->GetOrigin()[0] + bBox->GetDX(), bBox->GetOrigin()[1] + bBox->GetDY(), bBox->GetOrigin()[2] - bBox->GetDZ()},
     {bBox->GetOrigin()[0] + bBox->GetDX(), bBox->GetOrigin()[1] - bBox->GetDY(), bBox->GetOrigin()[2] + bBox->GetDZ()},
     {bBox->GetOrigin()[0] + bBox->GetDX(), bBox->GetOrigin()[1] - bBox->GetDY(), bBox->GetOrigin()[2] - bBox->GetDZ()},
     {bBox->GetOrigin()[0] - bBox->GetDX(), bBox->GetOrigin()[1] + bBox->GetDY(), bBox->GetOrigin()[2] + bBox->GetDZ()},
     {bBox->GetOrigin()[0] - bBox->GetDX(), bBox->GetOrigin()[1] + bBox->GetDY(), bBox->GetOrigin()[2] - bBox->GetDZ()},
     {bBox->GetOrigin()[0] - bBox->GetDX(), bBox->GetOrigin()[1] - bBox->GetDY(), bBox->GetOrigin()[2] + bBox->GetDZ()},
     {bBox->GetOrigin()[0] - bBox->GetDX(), bBox->GetOrigin()[1] - bBox->GetDY(),
      bBox->GetOrigin()[2] - bBox->GetDZ()}}};
 
  std::array<std::pair<double, double>, 3> bounds = {
    {// Bounds of the algo::GeoVolume
     {+std::numeric_limits<double>::max(), -std::numeric_limits<double>::max()},
     {+std::numeric_limits<double>::max(), -std::numeric_limits<double>::max()},
     {+std::numeric_limits<double>::max(), -std::numeric_limits<double>::max()}}};  // min/max for a dimension
 
  const auto* matrix{gGeoManager->GetCurrentMatrix()};
  for (int iVertex = 0; iVertex < 8; ++iVertex) {
    std::array<double, 3> glbVertex;
    matrix->LocalToMaster(locBoxVertices[iVertex].data(), glbVertex.data());
    for (int iDim = 0; iDim < 3; ++iDim) {
      bounds[iDim].first  = std::min(bounds[iDim].first, glbVertex[iDim]);
      bounds[iDim].second = std::max(bounds[iDim].second, glbVertex[iDim]);
    }
  }
 
  gGeoManager->cd();
  return algo::GeoVolume(bounds[0], bounds[1], bounds[2]);
}
 
 
// ---------------------------------------------------------------------------------------------------------------------
//
cbm::algo::GeoVolume RecoSetupUtils::ReadVolume(const std::string& path) { return ReadVolume(TString(path)); }