LArPfoObjects.cc

Go to the documentation of this file.

 
#include "larpandoracontent/LArObjects/LArPfoObjects.h"
 
#include "Objects/CaloHit.h"
 
using namespace pandora;
 
namespace lar_content
{
 
LArTrackState::LArTrackState(const CartesianVector &position, const CartesianVector &direction, const CaloHit *const pCaloHit) :
    TrackState(position, direction),
    m_pCaloHit(pCaloHit)
{
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
LArTrackState::LArTrackState(const CartesianVector &position, const CartesianVector &direction) :
    TrackState(position, direction),
    m_pCaloHit(nullptr)
{
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArTrackState::GetDirection() const
{
    return this->GetMomentum();
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CaloHit *LArTrackState::GetCaloHit() const
{
    if (m_pCaloHit)
        return m_pCaloHit;
 
    throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
//------------------------------------------------------------------------------------------------------------------------------------------
 
LArShowerPCA::LArShowerPCA(const CartesianVector &centroid, const CartesianVector &primaryAxis, const CartesianVector &secondaryAxis,
    const CartesianVector &tertiaryAxis, const CartesianVector &eigenValues) :
    m_centroid(centroid),
    m_primaryAxis(primaryAxis),
    m_secondaryAxis(secondaryAxis),
    m_tertiaryAxis(tertiaryAxis),
    m_eigenValues(eigenValues),
    m_axisLengths(((eigenValues.GetX() > std::numeric_limits<float>::epsilon()) ? 6.f * std::sqrt(eigenValues.GetX()) : 0.f),
        ((eigenValues.GetY() > std::numeric_limits<float>::epsilon()) ? 6.f * std::sqrt(eigenValues.GetY()) : 0.f),
        ((eigenValues.GetZ() > std::numeric_limits<float>::epsilon()) ? 6.f * std::sqrt(eigenValues.GetZ()) : 0.f))
{
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetCentroid() const
{
    return m_centroid;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetPrimaryAxis() const
{
    return m_primaryAxis;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetSecondaryAxis() const
{
    return m_secondaryAxis;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetTertiaryAxis() const
{
    return m_tertiaryAxis;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetEigenValues() const
{
    return m_eigenValues;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
const CartesianVector &LArShowerPCA::GetAxisLengths() const
{
    return m_axisLengths;
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
float LArShowerPCA::GetPrimaryLength() const
{
    return m_axisLengths.GetX();
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
float LArShowerPCA::GetSecondaryLength() const
{
    return m_axisLengths.GetY();
}
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
float LArShowerPCA::GetTertiaryLength() const
{
    return m_axisLengths.GetZ();
}
 
} // namespace lar_content