In This Package:

PositionerTool.cc

Go to the documentation of this file.

#include "PositionerTool.h"

#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/MsgStream.h"

#include "DetDesc/LVolume.h"
#include "DetDesc/PVolume.h"
#include "DetDesc/DetectorElement.h"
#include "GaudiKernel/Point3DTypes.h"
#include "GaudiKernel/Transform3DTypes.h"

#include <string>

PositionerTool::PositionerTool(const std::string& type,
                               const std::string& name,
                               const IInterface* parent)
  : GaudiTool(type, name, parent)
  , m_detSvc(0)
{
  declareProperty("PhysicalVolume",m_physVolName="",
                  "Name of the physical volume to add to geometry");
  declareProperty("LogicalVolume",m_logVolName="",
                  "Name of the logical volume to associate with the physical");
  declareProperty("MotherVolume",m_motherLVolName="",
                  "Name of the mother logical volume where the physical volume is placed");
  declareProperty("Position",m_position=std::vector<double>(3,0.),
                  "Position of the new volume in the parent coordinates");
  double identity[] = {1., 0., 0., 
                       0., 1., 0.,
                       0., 0., 1.};
  declareProperty("Rotation",m_rotation=std::vector<double>(identity,
                                                            identity+9),
                  "Rotation of the new volume in the parent coordinates");
  declareProperty("Element",m_detElemPath="",
                  "Name of the detector element to place in the TDS");
  declareProperty("Support",m_support="",
                  "Path+Name of the supporting detector element in the TDS");
}

PositionerTool::~PositionerTool()
{
}



StatusCode PositionerTool::initialize()
{
    this->GaudiTool::initialize();

    debug() << "PositionerTool::initialize()" << endreq;
    StatusCode sc = service("DetectorDataSvc",m_detSvc,true);
    if (sc.isFailure()) return sc;

    return sc;
}

StatusCode PositionerTool::reinitialize()
{
    return StatusCode::SUCCESS;
}

StatusCode PositionerTool::finalize()
{
    m_detSvc->release();
    return this->GaudiTool::finalize();
}

StatusCode PositionerTool::queryInterface(const InterfaceID& riid, 
                                          void** ppvInterface)
{
  StatusCode sc = StatusCode::FAILURE;
  if (ppvInterface) {
    *ppvInterface = 0;
    
    if (IPositionerTool::interfaceID().versionMatch(riid)) {
      *ppvInterface = static_cast<IPositionerTool*>(this);
      sc = StatusCode::SUCCESS;
      addRef();
    }
    else sc = GaudiTool::queryInterface( riid, ppvInterface );    
  }
  return sc;
}



StatusCode PositionerTool::placeVolume()
{

  // Find parent logical volume for physical volume
  SmartDataPtr<LVolume> parentLogVol(m_detSvc, m_motherLVolName);
  if(!parentLogVol){
    error() << "PositionerTool::placeVolume(): "
            << "could not find parent volume of name: " << m_motherLVolName
            << endreq;
    return StatusCode::FAILURE;
  }

  // Check if physical volume already exists in parent logical volume
  ILVolume::PVolumes::const_iterator pvIter = parentLogVol->pvBegin();
  for(;pvIter != parentLogVol->pvEnd(); pvIter++){
    if((*pvIter)->name() == m_physVolName){
      warning() << "PositionerTool::placeVolume(): "
                << "physical volume " << m_physVolName
                << " already exists in mother volume " << m_motherLVolName
                << endreq;
      return StatusCode::SUCCESS;
    }
  }

  // Create the physical volume
  Gaudi::Rotation3D rotation(m_rotation.begin(), m_rotation.end());
  Gaudi::XYZVector position(m_position[0],m_position[1],m_position[2]);
  IPVolume* pv = 
    parentLogVol->createPVolume(m_physVolName, m_logVolName, 
                                Gaudi::Transform3D(rotation, 
                                                   Gaudi::XYZVector(position)));
  
  if( 0 == pv ) {
    error() << "PositionerTool::placeVolume(): "
            << "could not create volume: " << m_physVolName
            << endreq;
    return StatusCode::FAILURE;
  }

  StatusCode sc;
  
  if(m_detElemPath != ""){
    // Add the corresponding detector element
    DetectorElement* detElem = new DetectorElement;
    sc = m_detSvc->registerObject( m_support, m_detElemPath, detElem);
    if( sc.isFailure() ){
      error() << "Failed to register object with support / path: " 
              << m_support << " / " << m_detElemPath << endreq;
    }
    detElem->createGeometryInfo( m_logVolName, m_support, m_physVolName);
  }

  return sc;
}

---