lar_content::LArEigenHelper Class Reference

LArEigenHelper class. More...

#include "LArEigenHelper.h"

Static Public Member Functions
template<class T >
static void	Vectorize (const T &caloHitContainer, Eigen::MatrixXf &hitMatrix)
	Convert a container of calo hits into an Eigen matrix. More...
template<class T >
static void	Vectorize (const T &caloHitContainer, Eigen::MatrixXf &centre, Eigen::MatrixXf &low, Eigen::MatrixXf &high)
	Convert a container of calo hits into a collection of Eigen matrices representing the centre, low and high coordinates of hits. More...
static void	GetAngles (const Eigen::MatrixXf &hitMatrix, const Eigen::RowVectorXf &origin, Eigen::RowVectorXf &phis)
	Retrieve the angle, coutner-clockwise relative to the x axis, between all hits in a matrix and a specified origin. More...

Detailed Description

LArEigenHelper class.

Definition at line 23 of file LArEigenHelper.h.

Member Function Documentation

void lar_content::LArEigenHelper::GetAngles ( const Eigen::MatrixXf &  hitMatrix, const Eigen::RowVectorXf &  origin, Eigen::RowVectorXf &  phis  )

static

Retrieve the angle, coutner-clockwise relative to the x axis, between all hits in a matrix and a specified origin.

Parameters

hitMatrix	the input collection of hits
origin	the origin from which to measure hits
phis	the output vector of angles in the range [0, 2pi]

Definition at line 58 of file LArEigenHelper.cc.

References geo::origin(), and util::pi().

{
    const float pi{static_cast<float>(M_PI)};
    Eigen::RowVectorXf piVec{Eigen::RowVectorXf::Constant(hitMatrix.rows(), pi)};
    Eigen::RowVectorXf zeroVec{Eigen::RowVectorXf::Zero(hitMatrix.rows())};
    Eigen::MatrixXf deltas(hitMatrix.rowwise() - origin);
    for (int i = 0; i < deltas.rows(); ++i)
        phis(i) = std::atan2(deltas(i, 1), deltas(i, 0));
    // Move from [-pi, +pi] to [0, 2pi)
    phis = (phis.array() < 0).select(piVec * 2 + phis, phis);
    phis = (phis.array() < 2 * pi).select(phis, zeroVec);
}

template<class T >

void lar_content::LArEigenHelper::Vectorize ( const T &  caloHitContainer, Eigen::MatrixXf &  hitMatrix  )

static

Convert a container of calo hits into an Eigen matrix.

Parameters

caloHitContainer	the calo hit list containing the hits from which to construct a maxtrix
hitMatrix	the output Eigen matrix

Definition at line 26 of file LArEigenHelper.cc.

{
    int i{0};
    for (const CaloHit *const pCaloHit : caloHitContainer)
    {
        const CartesianVector &pos{pCaloHit->GetPositionVector()};
        hitMatrix(i, 0) = pos.GetX();
        hitMatrix(i, 1) = pos.GetZ();
        ++i;
    }
}

template<class T >

void lar_content::LArEigenHelper::Vectorize ( const T &  caloHitContainer, Eigen::MatrixXf &  centre, Eigen::MatrixXf &  low, Eigen::MatrixXf &  high  )

static

Convert a container of calo hits into a collection of Eigen matrices representing the centre, low and high coordinates of hits.

Parameters

caloHitContainer	the calo hit list containing the hits from which to construct matrices
centre	the output Eigen matrix of hit centres
low	the output Eigen matrix of hit low edges (centre - width / 2)
high	the output Eigen matrix of hit high edges (centre + width / 2)

Definition at line 41 of file LArEigenHelper.cc.

{
    LArEigenHelper::Vectorize(caloHitContainer, centre);
    int i{0};
    for (const CaloHit *const pCaloHit : caloHitContainer)
    {
        const CartesianVector &pos{pCaloHit->GetPositionVector()};
        low(i, 0) = pos.GetX() - pCaloHit->GetCellSize1() * 0.5f;
        low(i, 1) = pos.GetZ();
        high(i, 0) = pos.GetX() + pCaloHit->GetCellSize1() * 0.5f;
        high(i, 1) = pos.GetZ();
        ++i;
    }
}

---

The documentation for this class was generated from the following files:

• larpandoracontent/v04_14_00/source/larpandoracontent/LArHelpers/LArEigenHelper.h
• larpandoracontent/v04_14_00/source/larpandoracontent/LArHelpers/LArEigenHelper.cc