RpcGeomInfo.cc

Go to the documentation of this file.

#include "RpcGeomInfo.h"

#include "DetDesc/IDetectorElement.h"
#include "DetDesc/IGeometryInfo.h"

#include "GaudiKernel/Point3DTypes.h"
#include "GaudiKernel/Vector3DTypes.h"

#include "CLHEP/Vector/ThreeVector.h"

using namespace CLHEP;

RpcGeomInfo::RpcGeomInfo(unsigned int rpcid,
                         IDetectorElement* me,
                         IDetectorElement* parent)
    : m_id(rpcid)
    , m_me(me)
    , m_parent(parent)
    , m_gp(0), m_lp(0), m_gd(0), m_ld(0)
{
}

RpcGeomInfo::~RpcGeomInfo()
{
    if (m_gp) delete m_gp; m_gp = 0;
    if (m_lp) delete m_lp; m_lp = 0;
    if (m_gd) delete m_gd; m_gd = 0;
    if (m_ld) delete m_ld; m_ld = 0;
}

unsigned int RpcGeomInfo::rpcid() const
{
    return m_id;
}

const Hep3Vector& RpcGeomInfo::globalPosition() const
{
    if (m_gp) return *m_gp;

    Gaudi::XYZPoint zero(0,0,0);
    Gaudi::XYZPoint gp = m_me->geometry()->toGlobal(zero);

    m_gp = new Hep3Vector(gp.x(),gp.y(),gp.z());
    return *m_gp;
}

const Hep3Vector& RpcGeomInfo::localPosition() const
{
    if (m_lp) return *m_lp;

    Gaudi::XYZPoint zero(0,0,0);
    Gaudi::XYZPoint gp = m_me->geometry()->toGlobal(zero);
    Gaudi::XYZPoint lp = m_parent->geometry()->toLocal(gp);
    m_lp = new Hep3Vector(lp.x(),lp.y(),lp.z());
    return *m_lp;
}

const Hep3Vector& RpcGeomInfo::globalDirection() const
{
    if (m_gd) return *m_gd;

    Gaudi::XYZVector norm(0,0,1);
    Gaudi::XYZVector gd = m_me->geometry()->toGlobal(norm);
    m_gd = new Hep3Vector(gd.x(),gd.y(),gd.z());
    return *m_gd;
}

const Hep3Vector& RpcGeomInfo::localDirection() const
{
    if (m_ld) return *m_ld;

    Gaudi::XYZVector norm(0,0,1);
    Gaudi::XYZVector gd = m_me->geometry()->toGlobal(norm);
    Gaudi::XYZVector ld = m_parent->geometry()->toLocal(gd);
    m_ld = new Hep3Vector(ld.x(),ld.y(),ld.z());
    return *m_ld;
}

const IDetectorElement& RpcGeomInfo::detectorElement() const
{
    return *m_me;
}

const IDetectorElement& RpcGeomInfo::parentDetector() const
{
    return *m_parent;
}