Geant4Mapping.cpp

// $Id: Geant4Mapping.cpp 588 2013-06-03 11:41:35Z markus.frank $
//====================================================================
//  AIDA Detector description implementation for LCD
//--------------------------------------------------------------------
//
//  Author     : M.Frank
//
//====================================================================

#include "DD4hep/Printout.h"
#include "DDG4/Geant4Mapping.h"
#include "G4PVPlacement.hh"
#include <stdexcept>

using namespace DD4hep::Simulation;
using namespace DD4hep::Geometry;
using namespace std;

/// Initializing Constructor
Geant4Mapping::Geant4Mapping(LCDD& lcdd) : m_lcdd(lcdd), m_dataPtr(0)
{
}

/// Standard destructor
Geant4Mapping::~Geant4Mapping()  {
  if ( m_dataPtr ) delete m_dataPtr;
  m_dataPtr = 0;
}

/// Possibility to define a singleton instance
Geant4Mapping& Geant4Mapping::instance()  {
  static Geant4Mapping inst(LCDD::getInstance());
  return inst;
}

/// When resolving pointers, we must check for the validity of the data block
void Geant4Mapping::checkValidity() const   {
  if ( m_dataPtr ) return;
  throw runtime_error("Geant4Mapping: Attempt to access an invalid data block!");
}

/// Create new data block. Delete old data block if present.
Geant4GeometryInfo& Geant4Mapping::init()  {
  Geant4GeometryInfo* p = detach();
  if ( p ) delete p;
  attach(new Geant4GeometryInfo());
  return data();
}

/// Release data and pass over the ownership
Geant4GeometryInfo* Geant4Mapping::detach()   {
  Geant4GeometryInfo* p = m_dataPtr;
  m_dataPtr = 0;
  return p;
}

/// Set a new data block
void Geant4Mapping::attach(Geant4GeometryInfo* data)   {
  m_dataPtr = data;
}

/// Access the volume manager
Geant4VolumeManager Geant4Mapping::volumeManager() const   {
  if ( m_dataPtr )  {
    if ( m_dataPtr->g4Paths.empty() )  {
      return Geant4VolumeManager(m_lcdd,m_dataPtr);
    }
    return Geant4VolumeManager(Geometry::Handle<Geant4GeometryInfo>(m_dataPtr));
  }
  throw runtime_error(format("Geant4Mapping","Cannot create volume manager without Geant4 geometry info [Invalid-Info]"));
}

/// Accessor to resolve geometry placements
PlacedVolume Geant4Mapping::placement(const G4VPhysicalVolume* node)  const   {
  checkValidity();
  const PlacementMap& m = m_dataPtr->g4Placements;
  for(PlacementMap::const_iterator i = m.begin(); i != m.end(); ++i)
    if ( (*i).second == node ) return PlacedVolume((*i).first);
  return PlacedVolume(0);
}