TRRegionXMLHandler.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/
 
// class header
#include "TRRegionXMLHandler.h"
 
// package includes
#include "TRTRadiatorParameters.h"
#include "TRTTransitionRadiation.h"
 
// Athena includes
#include "IdDictDetDescr/IdDictManager.h"
#include "IdDict/IdDictMgr.h"
 
// Gaudi includes
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/Bootstrap.h"
 
// Geant4 includes
#include "G4LogicalVolumeStore.hh" // For logical volume setup
 
// For XML parsigin
#include "boost/property_tree/xml_parser.hpp"
#include "boost/property_tree/ptree.hpp"
 
// STL includes
#include <iostream>
#include <fstream>
#include <string>
#include <stdexcept>
 
TRRegionXMLHandler::TRRegionXMLHandler(TRTTransitionRadiation *tr) :
  AthMessaging("TRRegionXMLHandler"),
  m_theProcess(tr)
{}
 
void TRRegionXMLHandler::Process(const std::string& name)
{
  ISvcLocator * svcLocator = Gaudi::svcLocator(); // from Bootstrap
 
  ATH_MSG_DEBUG("This is TRRegionXMLHandler. Handler called");
 
  SmartIF<StoreGateSvc> smartSG{svcLocator->service("StoreGateSvc")};
  if (smartSG) { m_storeGate = smartSG.get(); }
  if( !m_storeGate ) {
    ATH_MSG_ERROR("Unable to locate StoreGate! Stopping!");
    throw std::runtime_error("Unable to locate StoreGate!");
  }
  SmartIF<StoreGateSvc> detStore{svcLocator->service( "DetectorStore")};
  if( !detStore ) {
    ATH_MSG_ERROR("Unable to locate DetectorStore! Leaving!");
    throw std::runtime_error("Unable to locate DetectorStore!");
  }
 
  const IdDictManager * idDictMgr{};
  if (StatusCode::SUCCESS == detStore->retrieve(idDictMgr, "IdDict")) {
    if (idDictMgr) {
      std::string tag = idDictMgr->manager()->tag();
      m_initialLayoutIdDict =
        (tag == "initial_layout" || tag == "destaged_layout");
    }
  } else {
    ATH_MSG_FATAL("Could not retrieve geometry layout. TR process is not to be trusted in the following");
    throw std::runtime_error("Could not retrieve geometry layout!");
  }
 
  // Crack open the XML file
  std::filebuf fb;
  if (!fb.open(name,std::ios::in)){
    ATH_MSG_FATAL("Could not open file " << name << " bombing out");
    throw std::runtime_error("Could not open file!");
  }
  std::istream is(&fb);
  boost::property_tree::ptree pt;
  read_xml(is, pt);
 
  for( boost::property_tree::ptree::value_type const& v : pt.get_child("FADS") ) {
    if( v.first == "TRRegionParameters" ) {
 
      std::string volName=v.second.get<std::string>("<xmlattr>.RadiatorName");
      double foilThickness=v.second.get<double>("<xmlattr>.RadiatorFoilThickness");
      double gasThickness=v.second.get<double>("<xmlattr>.RadiatorGasThickness");
      int regionFlag=v.second.get<int>("<xmlattr>.RadiatorBARRELorENDCAP");
      std::string detectorPart=v.second.get<std::string>("<xmlattr>.DetectorPart");
    
      G4LogicalVolumeStore *g4lvs = G4LogicalVolumeStore::GetInstance();
      unsigned int numberOfVolumes = 0;
      for (const auto log_vol : *g4lvs){
        if (volName == static_cast<const std::string&>(log_vol->GetName())) {
          TRTRadiatorParameters rad( log_vol,
                                     foilThickness,gasThickness,
                                     (BEflag)regionFlag );
          m_theProcess->AddRadiatorParameters(rad);
          ++numberOfVolumes;
        }
      }
    
      if( numberOfVolumes == 0 ) {
        // In case of destaged geometry (no C wheels) any missing volume
        // except the radiators corresponding to the C wheels will cause
        // the code to abort
        if ( m_initialLayoutIdDict != 0 ) {
          if  ( ( volName!="TRT::MainRadiatorC" ) &&
                ( volName!="TRT::ThinRadiatorC" )    ) {
            ATH_MSG_FATAL(" Volume-name " << volName
                          <<" not found! Geometry layout "
                          << m_initialLayoutIdDict);
            std::abort();
          }
        }
      } // Didn't get any volumes
    } // Loop over region parameters objects
  } // Loop over FADS things in the XML file
 
}