StandaloneGeometrySetup.cxx

Go to the documentation of this file.

#include "larcorealg/Geometry/StandaloneGeometrySetup.h"

// LArSoft libraries
#include "larcorealg/Geometry/ChannelMapAlg.h"
#include "larcorealg/Geometry/GeometryCore.h"

// CET libraries
#include "fhiclcpp/ParameterSet.h"

// ROOT libraries
#include "TGeoManager.h"

// C/C++ standard libraries
#include <string>
#include <utility> // std::forward()
#include <memory> // std::make_unique(), std::make_shared()


//------------------------------------------------------------------------------
std::unique_ptr<geo::GeometryCore>
lar::standalone::SetupGeometryWithChannelMapping
  (
    fhicl::ParameterSet const& pset,
    std::shared_ptr<geo::ChannelMapAlg> channelMap
  )
{
  auto const bForceReload = true;
  
  //
  // create the geometry object
  //
  auto geom = std::make_unique<geo::GeometryCore>(pset);
  
  //
  // extract of relevant configuration parameters
  //
  std::string relPath              = pset.get<std::string>("RelativePath",     ""   );
  const bool disableWiresInG4      = pset.get<bool>       ("DisableWiresInG4", false);
  const std::string GDMLFileName   = pset.get<std::string>("GDML"                   );
//  const std::string ROOTFileName   = pset.get<std::string>("ROOT"                   );
  
  // add a final directory separator ("/") to relPath if not already there
  if (!relPath.empty() && (relPath.back() != '/')) relPath += '/';
  
  // We are going to find files now.
  // cet::search_path constructor decides if the constructor argument is a path
  // or an environment variable (in this case, the latter)
  cet::search_path sp("FW_SEARCH_PATH");
  
  //
  // "GDML" file (for GEANT4)
  //
  // this is our hint for the path; start with the relative path:
  std::string GDMLFilePathHint = relPath + GDMLFileName;
  
  // special if geometry with no wires is used for GEANT4 simulation
  if(disableWiresInG4) {
    GDMLFilePathHint.insert(
      std::min(GDMLFilePathHint.rfind(".gdml"), GDMLFilePathHint.length()),
      "_nowires"
      );
  } // if disable wires
  
  std::string GDMLFilePath;
  if( !sp.find_file(GDMLFilePathHint, GDMLFilePath) ) {
    throw cet::exception("StaticLoadGeometry")
      << "Can't find geometry file '" << GDMLFilePathHint
      << "' (for GEANT4)!\n";
  }
  
  //
  // "ROOT" file (for geometry)
  //
  // this is our hint for the path; start with the relative path:
  std::string ROOTFilePathHint = relPath + GDMLFileName;
  
  std::string ROOTFilePath;
  if( !sp.find_file(ROOTFilePathHint, ROOTFilePath) ) {
    throw cet::exception("StaticLoadGeometry")
      << "Can't find geometry file '" << ROOTFilePathHint
      << "' (for geometry)!\n";
  }
  
  //
  // initialize the geometry with the files we have found
  //
  geom->LoadGeometryFile(GDMLFilePath, ROOTFilePath, bForceReload);
  
  //
  // create and apply channel mapping
  //
  
  geom->ApplyChannelMap(channelMap);
  
  return geom;
} // lar::standalone::SetupGeometryWithChannelMapping()


//------------------------------------------------------------------------------