lar_content::TwoDSlidingFitResult Class Reference

TwoDSlidingFitResult class. More...

#include "LArTwoDSlidingFitResult.h"

Public Member Functions
template<typename T >
	TwoDSlidingFitResult (const T *const pT, const unsigned int layerFitHalfWindow, const float layerPitch)
	Constructor using internal definition of primary axis. More...
template<typename T >
	TwoDSlidingFitResult (const T *const pT, const unsigned int layerFitHalfWindow, const float layerPitch, const pandora::CartesianVector &axisIntercept, const pandora::CartesianVector &axisDirection, const pandora::CartesianVector &orthoDirection)
	Constructor using specified primary axis. The orthogonal axis must be perpendicular to the primary axis. More...
	TwoDSlidingFitResult (const unsigned int layerFitHalfWindow, const float layerPitch, const pandora::CartesianVector &axisIntercept, const pandora::CartesianVector &axisDirection, const pandora::CartesianVector &orthoDirection, const LayerFitContributionMap &layerFitContributionMap)
	Constructor using specified primary axis and layer fit contribution map. User is responsible for ensuring that z-pitch, axis intercept and axis direction agree with calculations used to fill the layer fit contribution map. The orthogonal axis must be perpendicular to the primary axis. More...
const pandora::Cluster *	GetCluster () const
	Get the address of the cluster, if originally provided. More...
unsigned int	GetLayerFitHalfWindow () const
	Get the layer fit half window. More...
float	GetLayerPitch () const
	Get the layer pitch, units cm. More...
const pandora::CartesianVector &	GetAxisIntercept () const
	Get the axis intercept position. More...
const pandora::CartesianVector &	GetAxisDirection () const
	Get the axis direction vector. More...
const pandora::CartesianVector &	GetOrthoDirection () const
	Get the orthogonal direction vector. More...
const LayerFitResultMap &	GetLayerFitResultMap () const
	Get the layer fit result map. More...
const LayerFitContributionMap &	GetLayerFitContributionMap () const
	Get the layer fit contribution map. More...
const FitSegmentList &	GetFitSegmentList () const
	Get the fit segment list. More...
float	GetLayerFitHalfWindowLength () const
	Get the layer fit half window length. More...
int	GetMinLayer () const
	Get the minimum occupied layer in the sliding fit. More...
int	GetMaxLayer () const
	Get the maximum occupied layer in the sliding fit. More...
void	GetMinAndMaxX (float &minX, float &maxX) const
	Get the minimum and maximum x coordinates associated with the sliding fit. More...
void	GetMinAndMaxZ (float &minZ, float &maxZ) const
	Get the minimum and maximum z coordinates associated with the sliding fit. More...
int	GetLayer (const float rL) const
	Get layer number for given sliding linear fit longitudinal coordinate. More...
float	GetL (const int layer) const
	Get longitudinal coordinate for a given sliding linear fit layer number. More...
void	GetLocalPosition (const pandora::CartesianVector &position, float &rL, float &rT) const
	Get local sliding fit coordinates for a given global position. More...
void	GetLocalDirection (const pandora::CartesianVector &direction, float &dTdL) const
	Get local sliding fit gradient for a given global direction. More...
void	GetGlobalPosition (const float rL, const float rT, pandora::CartesianVector &position) const
	Get global coordinates for given sliding linear fit coordinates. More...
void	GetGlobalDirection (const float dTdL, pandora::CartesianVector &direction) const
	Get global direction coordinates for given sliding linear fit gradient. More...
pandora::CartesianVector	GetGlobalMinLayerPosition () const
	Get global position corresponding to the fit result in minimum fit layer. More...
pandora::CartesianVector	GetGlobalMaxLayerPosition () const
	Get global position corresponding to the fit result in maximum fit layer. More...
pandora::CartesianVector	GetGlobalMinLayerDirection () const
	Get global direction corresponding to the fit result in minimum fit layer. More...
pandora::CartesianVector	GetGlobalMaxLayerDirection () const
	Get global direction corresponding to the fit result in maximum fit layer. More...
float	GetMinLayerRms () const
	Get rms at minimum layer. More...
float	GetMaxLayerRms () const
	Get rms at maximum layer. More...
float	GetFitRms (const float rL) const
	Get fit rms for a given longitudinal coordinate. More...
float	GetCosScatteringAngle (const float rL) const
	Get scattering angle for a given longitudinal coordinate. More...
pandora::StatusCode	GetGlobalFitPosition (const float rL, pandora::CartesianVector &position) const
	Get global fit position for a given longitudinal coordinate. More...
pandora::StatusCode	GetGlobalFitDirection (const float rL, pandora::CartesianVector &direction) const
	Get global fit direction for a given longitudinal coordinate. More...
pandora::StatusCode	GetGlobalFitPositionAtX (const float x, pandora::CartesianVector &position) const
	Get global fit position for a given input x coordinate. More...
pandora::StatusCode	GetGlobalFitDirectionAtX (const float x, pandora::CartesianVector &direction) const
	Get global fit direction for a given input x coordinate. More...
pandora::StatusCode	GetGlobalFitProjection (const pandora::CartesianVector &inputPosition, pandora::CartesianVector &projectedPosition) const
	Get projected position on global fit for a given position vector. More...
pandora::StatusCode	GetGlobalFitPositionListAtX (const float x, pandora::CartesianPointVector &positionList) const
	Get a list of projected positions for a given input x coordinate. More...
pandora::StatusCode	GetTransverseProjection (const float x, const FitSegment &fitSegment, pandora::CartesianVector &position) const
	Get projected position for a given input x coordinate and fit segment. More...
pandora::StatusCode	GetTransverseProjection (const float x, const FitSegment &fitSegment, pandora::CartesianVector &position, pandora::CartesianVector &direction) const
	Get projected position and direction for a given input x coordinate and fit segment. More...
pandora::StatusCode	GetExtrapolatedPosition (const float rL, pandora::CartesianVector &position) const
	Get extrapolated position (beyond span) for a given input coordinate. More...
pandora::StatusCode	GetExtrapolatedDirection (const float rL, pandora::CartesianVector &direction) const
	Get extrapolated direction (beyond span) for a given input coordinate. More...
pandora::StatusCode	GetExtrapolatedPositionAtX (const float x, pandora::CartesianVector &position) const
	Get extrapolated position (beyond span) for a given input x coordinate. More...
const FitSegment &	GetFitSegment (const float rL) const
	Get fit segment for a given longitudinal coordinate. More...
template<>
	TwoDSlidingFitResult (const Cluster *const pCluster, const unsigned int layerFitHalfWindow, const float layerPitch)
template<>
	TwoDSlidingFitResult (const CartesianPointVector *const pPointVector, const unsigned int layerFitHalfWindow, const float layerPitch)
template<>
	TwoDSlidingFitResult (const Cluster *const pCluster, const unsigned int layerFitHalfWindow, const float layerPitch, const CartesianVector &axisIntercept, const CartesianVector &axisDirection, const CartesianVector &orthoDirection)
template<>
	TwoDSlidingFitResult (const CartesianPointVector *const pPointVector, const unsigned int layerFitHalfWindow, const float layerPitch, const CartesianVector &axisIntercept, const CartesianVector &axisDirection, const CartesianVector &orthoDirection)

Private Member Functions
void	CalculateAxes (const pandora::CartesianPointVector &coordinateVector, const float layerPitch)
	Calculate the longitudinal and transverse axes. More...
void	FillLayerFitContributionMap (const pandora::CartesianPointVector &coordinateVector)
	Fill the layer fit contribution map. More...
void	PerformSlidingLinearFit ()
	Perform the sliding linear fit. More...
void	FindSlidingFitSegments ()
	Find sliding fit segments; sections with tramsverse direction. More...
void	GetMinAndMaxCoordinate (const bool isX, float &min, float &max) const
	Get the minimum and maximum x or z coordinates associated with the sliding fit. More...
pandora::CartesianVector	GetGlobalFitPosition (const LayerInterpolation &layerInterpolation) const
	Interpolate a position between two layers. More...
pandora::CartesianVector	GetGlobalFitDirection (const LayerInterpolation &layerInterpolation) const
	Interpolate a direction between two layers. More...
float	GetFitRms (const LayerInterpolation &layerInterpolation) const
	Interpolate a rms between two layers. More...
pandora::StatusCode	LongitudinalInterpolation (const float rL, LayerInterpolation &layerInterpolation) const
	Get the pair of layers surrounding a specified longitudinal position. More...
pandora::StatusCode	TransverseInterpolation (const float x, const FitSegment &fitSegment, LayerInterpolation &layerInterpolation) const
	Get the surrounding pair of layers for a specified transverse position and fit segment. More...
pandora::StatusCode	TransverseInterpolation (const float x, LayerInterpolationList &layerInterpolationList) const
	Get the a list of surrounding layer pairs for a specified transverse position. More...
pandora::StatusCode	GetLongitudinalSurroundingLayers (const float rL, LayerFitResultMap::const_iterator &firstLayerIter, LayerFitResultMap::const_iterator &secondLayerIter) const
	Get iterators for layers surrounding the specified longitudinal position. More...
pandora::StatusCode	GetTransverseSurroundingLayers (const float x, const int minLayer, const int maxLayer, LayerFitResultMap::const_iterator &firstLayerIter, LayerFitResultMap::const_iterator &secondLayerIter) const
	Get iterators for layers surrounding a specified transverse position. More...
void	GetLongitudinalInterpolationWeights (const float rL, const LayerFitResultMap::const_iterator &firstLayerIter, const LayerFitResultMap::const_iterator &secondLayerIter, double &firstWeight, double &secondWeight) const
	Get interpolation weights for layers surrounding a specified longitudinal position. More...
void	GetTransverseInterpolationWeights (const float x, const LayerFitResultMap::const_iterator &firstLayerIter, const LayerFitResultMap::const_iterator &secondLayerIter, double &firstWeight, double &secondWeight) const
	Get interpolation weights for layers surrounding a specified transverse position. More...

Private Attributes
const pandora::Cluster *	m_pCluster
	The address of the cluster. More...
unsigned int	m_layerFitHalfWindow
	The layer fit half window. More...
float	m_layerPitch
	The layer pitch, units cm. More...
pandora::CartesianVector	m_axisIntercept
	The axis intercept position. More...
pandora::CartesianVector	m_axisDirection
	The axis direction vector. More...
pandora::CartesianVector	m_orthoDirection
	The orthogonal direction vector. More...
LayerFitResultMap	m_layerFitResultMap
	The layer fit result map. More...
LayerFitContributionMap	m_layerFitContributionMap
	The layer fit contribution map. More...
FitSegmentList	m_fitSegmentList
	The fit segment list. More...

Detailed Description

TwoDSlidingFitResult class.

Definition at line 23 of file LArTwoDSlidingFitResult.h.

Constructor & Destructor Documentation

template<typename T >

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const T *const	pT,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch
	)

Constructor using internal definition of primary axis.

Parameters

pT	describing the positions to be fitted
layerFitHalfWindow	the layer fit half window
layerPitch	the layer pitch, units cm

Referenced by TwoDSlidingFitResult().

template<typename T >

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const T *const	pT,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch,
		const pandora::CartesianVector &	axisIntercept,
		const pandora::CartesianVector &	axisDirection,
		const pandora::CartesianVector &	orthoDirection
	)

Constructor using specified primary axis. The orthogonal axis must be perpendicular to the primary axis.

Parameters

pT	describing the positions to be fitted
layerFitHalfWindow	the layer fit half window
layerPitch	the layer pitch, units cm
axisIntercept	the axis intercept position
axisDirection	the axis direction vector
orthoDirection	the orthogonal direction vector

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const unsigned int	layerFitHalfWindow,
		const float	layerPitch,
		const pandora::CartesianVector &	axisIntercept,
		const pandora::CartesianVector &	axisDirection,
		const pandora::CartesianVector &	orthoDirection,
		const LayerFitContributionMap &	layerFitContributionMap
	)

Constructor using specified primary axis and layer fit contribution map. User is responsible for ensuring that z-pitch, axis intercept and axis direction agree with calculations used to fill the layer fit contribution map. The orthogonal axis must be perpendicular to the primary axis.

Parameters

layerFitHalfWindow	the layer fit half window
layerPitch	the layer pitch, units cm
axisIntercept	the axis intercept position
axisDirection	the axis direction vector
orthoDirection	the orthogonal direction vector
layerFitContributionMap	the layer fit contribution map

template<>

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const Cluster *const	pCluster,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch
	)

Definition at line 25 of file LArTwoDSlidingFitResult.cc.

References CalculateAxes(), FillLayerFitContributionMap(), FindSlidingFitSegments(), lar_content::LArClusterHelper::GetCoordinateVector(), and PerformSlidingLinearFit().

                                                                                                                                       :
    m_pCluster(pCluster),
    m_layerFitHalfWindow(layerFitHalfWindow),
    m_layerPitch(layerPitch),
    m_axisIntercept(0.f, 0.f, 0.f),
    m_axisDirection(0.f, 0.f, 0.f),
    m_orthoDirection(0.f, 0.f, 0.f)
{
    CartesianPointVector pointVector;
    LArClusterHelper::GetCoordinateVector(pCluster, pointVector);
    this->CalculateAxes(pointVector, layerPitch);
    this->FillLayerFitContributionMap(pointVector);
    this->PerformSlidingLinearFit();
    this->FindSlidingFitSegments();
}

template<>

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const CartesianPointVector *const	pPointVector,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch
	)

Definition at line 42 of file LArTwoDSlidingFitResult.cc.

References CalculateAxes(), FillLayerFitContributionMap(), FindSlidingFitSegments(), and PerformSlidingLinearFit().

                                                                                                                                                        :
    m_pCluster(nullptr),
    m_layerFitHalfWindow(layerFitHalfWindow),
    m_layerPitch(layerPitch),
    m_axisIntercept(0.f, 0.f, 0.f),
    m_axisDirection(0.f, 0.f, 0.f),
    m_orthoDirection(0.f, 0.f, 0.f)
{
    this->CalculateAxes(*pPointVector, layerPitch);
    this->FillLayerFitContributionMap(*pPointVector);
    this->PerformSlidingLinearFit();
    this->FindSlidingFitSegments();
}

template<>

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const Cluster *const	pCluster,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch,
		const CartesianVector &	axisIntercept,
		const CartesianVector &	axisDirection,
		const CartesianVector &	orthoDirection
	)

Definition at line 59 of file LArTwoDSlidingFitResult.cc.

References FillLayerFitContributionMap(), FindSlidingFitSegments(), lar_content::LArClusterHelper::GetCoordinateVector(), and PerformSlidingLinearFit().

                                                                                                                           :
    m_pCluster(pCluster),
    m_layerFitHalfWindow(layerFitHalfWindow),
    m_layerPitch(layerPitch),
    m_axisIntercept(axisIntercept),
    m_axisDirection(axisDirection),
    m_orthoDirection(orthoDirection)
{
    CartesianPointVector pointVector;
    LArClusterHelper::GetCoordinateVector(pCluster, pointVector);
    this->FillLayerFitContributionMap(pointVector);
    this->PerformSlidingLinearFit();
    this->FindSlidingFitSegments();
}

template<>

lar_content::TwoDSlidingFitResult::TwoDSlidingFitResult	(	const CartesianPointVector *const	pPointVector,
		const unsigned int	layerFitHalfWindow,
		const float	layerPitch,
		const CartesianVector &	axisIntercept,
		const CartesianVector &	axisDirection,
		const CartesianVector &	orthoDirection
	)

Definition at line 76 of file LArTwoDSlidingFitResult.cc.

References FillLayerFitContributionMap(), FindSlidingFitSegments(), m_axisDirection, m_axisIntercept, m_layerFitContributionMap, m_layerFitHalfWindow, m_layerPitch, m_orthoDirection, m_pCluster, PerformSlidingLinearFit(), and TwoDSlidingFitResult().

                                                                                                                           :
    m_pCluster(nullptr),
    m_layerFitHalfWindow(layerFitHalfWindow),
    m_layerPitch(layerPitch),
    m_axisIntercept(axisIntercept),
    m_axisDirection(axisDirection),
    m_orthoDirection(orthoDirection)
{
    this->FillLayerFitContributionMap(*pPointVector);
    this->PerformSlidingLinearFit();
    this->FindSlidingFitSegments();
}

Member Function Documentation

void lar_content::TwoDSlidingFitResult::CalculateAxes ( const pandora::CartesianPointVector &  coordinateVector, const float  layerPitch  )

private

Calculate the longitudinal and transverse axes.

Definition at line 555 of file LArTwoDSlidingFitResult.cc.

References f, m_axisDirection, m_axisIntercept, m_orthoDirection, max, min, and lar_content::LArPcaHelper::RunPca().

Referenced by TwoDSlidingFitResult().

{
    if (coordinateVector.size() < 2)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    CartesianVector centroid(0.f, 0.f, 0.f);
    LArPcaHelper::EigenVectors eigenVecs;
    LArPcaHelper::EigenValues eigenValues(0.f, 0.f, 0.f);
    LArPcaHelper::RunPca(coordinateVector, centroid, eigenValues, eigenVecs);

    float minProjection(std::numeric_limits<float>::max());
    CartesianVector fitDirection(eigenVecs.at(0));

    // ATTN TwoDSlidingFitResult convention is to point in direction of increasing z
    if (fitDirection.GetZ() < 0.f)
        fitDirection *= -1.f;

    for (const CartesianVector &coordinate : coordinateVector)
        minProjection = std::min(minProjection, fitDirection.GetDotProduct(coordinate - centroid));

    // Define layers based on centroid rather than individual extremal hits
    const float fitProjection(layerPitch * std::floor(minProjection / layerPitch));

    m_axisIntercept = centroid + (fitDirection * fitProjection);
    m_axisDirection = fitDirection;

    // Use y-axis to generate an orthogonal axis (assuming that cluster occupies X-Z plane)
    CartesianVector yAxis(0.f, 1.f, 0.f);
    m_orthoDirection = yAxis.GetCrossProduct(m_axisDirection).GetUnitVector();
}

void lar_content::TwoDSlidingFitResult::FillLayerFitContributionMap ( const pandora::CartesianPointVector &  coordinateVector )

private

Fill the layer fit contribution map.

Definition at line 588 of file LArTwoDSlidingFitResult.cc.

References f, GetLayer(), GetLocalPosition(), m_axisDirection, m_layerFitContributionMap, m_layerPitch, and m_orthoDirection.

Referenced by TwoDSlidingFitResult().

{
    if (m_layerPitch < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    if ((m_axisDirection.GetMagnitudeSquared() < std::numeric_limits<float>::epsilon()) ||
        (m_orthoDirection.GetMagnitudeSquared() < std::numeric_limits<float>::epsilon()))
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    if (!m_layerFitContributionMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    for (CartesianPointVector::const_iterator iter = coordinateVector.begin(), iterEnd = coordinateVector.end(); iter != iterEnd; ++iter)
    {
        float rL(0.f), rT(0.f);
        this->GetLocalPosition(*iter, rL, rT);
        m_layerFitContributionMap[this->GetLayer(rL)].AddPoint(rL, rT);
    }
}

void lar_content::TwoDSlidingFitResult::FindSlidingFitSegments ( )

private

Find sliding fit segments; sections with tramsverse direction.

Definition at line 707 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), GetLayerFitResultMap(), m_fitSegmentList, lar_content::NEGATIVE_IN_X, lar_content::POSITIVE_IN_X, and lar_content::UNKNOWN.

Referenced by TwoDSlidingFitResult().

{
    if (!m_fitSegmentList.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    unsigned int nSustainedSteps(0);
    float sustainedDirectionStartX(0.f), sustainedDirectionEndX(0.f);

    CartesianVector previousPosition(0.f, 0.f, 0.f);
    TransverseDirection previousDirection(UNKNOWN), sustainedDirection(UNKNOWN);

    LayerFitResultMap::const_iterator sustainedDirectionStartIter, sustainedDirectionEndIter;
    const LayerFitResultMap &layerFitResultMap(this->GetLayerFitResultMap());

    for (LayerFitResultMap::const_iterator iter = layerFitResultMap.begin(), iterEnd = layerFitResultMap.end(); iter != iterEnd; ++iter)
    {
        CartesianVector position(0.f, 0.f, 0.f);
        this->GetGlobalPosition(iter->second.GetL(), iter->second.GetFitT(), position);

        // TODO currentDirection could also be UNCHANGED_IN_X
        const TransverseDirection currentDirection(((position - previousPosition).GetX() > 0.f) ? POSITIVE_IN_X : NEGATIVE_IN_X);

        if (previousDirection == currentDirection)
        {
            ++nSustainedSteps;

            if (nSustainedSteps > 2)
            {
                sustainedDirection = currentDirection;
                sustainedDirectionEndIter = iter;
                sustainedDirectionEndX = position.GetX();
            }
        }
        else
        {
            if ((POSITIVE_IN_X == sustainedDirection) || (NEGATIVE_IN_X == sustainedDirection))
                m_fitSegmentList.push_back(FitSegment(sustainedDirectionStartIter->first, sustainedDirectionEndIter->first, sustainedDirectionStartX, sustainedDirectionEndX));

            nSustainedSteps = 0;
            sustainedDirection = UNKNOWN;
            sustainedDirectionStartIter = iter;
            sustainedDirectionStartX = position.GetX();
        }

        previousPosition = position;
        previousDirection = currentDirection;
    }

    if ((POSITIVE_IN_X == sustainedDirection) || (NEGATIVE_IN_X == sustainedDirection))
        m_fitSegmentList.push_back(FitSegment(sustainedDirectionStartIter->first, sustainedDirectionEndIter->first, sustainedDirectionStartX, sustainedDirectionEndX));
}

const pandora::CartesianVector & lar_content::TwoDSlidingFitResult::GetAxisDirection ( ) const

inline

Get the axis direction vector.

Returns

the axis direction vector

Definition at line 570 of file LArTwoDSlidingFitResult.h.

References m_axisDirection.

Referenced by lar_content::TwoDSlidingShowerFitResult::GetShowerEdges().

{
    return m_axisDirection;
}

const pandora::CartesianVector & lar_content::TwoDSlidingFitResult::GetAxisIntercept ( ) const

inline

Get the axis intercept position.

Returns

the axis intercept position

Definition at line 563 of file LArTwoDSlidingFitResult.h.

References m_axisIntercept.

Referenced by lar_content::TwoDSlidingShowerFitResult::GetShowerEdges().

{
    return m_axisIntercept;
}

const pandora::Cluster * lar_content::TwoDSlidingFitResult::GetCluster

(

)

const

Get the address of the cluster, if originally provided.

Returns

the address of the cluster

Exceptions

StatusCodeException

Definition at line 109 of file LArTwoDSlidingFitResult.cc.

References m_pCluster.

Referenced by lar_content::CrossGapsAssociationAlgorithm::AreClustersAssociated(), lar_content::LArPointingCluster::BuildPointingCluster(), lar_content::TracksCrossingGapsTool::CheckXPositionInGap(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointVertex(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::VertexSplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::ThreeDSlidingFitResult::GetCluster(), lar_content::LArParticleIdPlugins::LArMuonId::GetMuonTrackWidth(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::CrossGapsAssociationAlgorithm::IsAssociated(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), lar_content::ThreeDTrackFragmentsAlgorithm::PerformMainLoop(), lar_content::TwoDSlidingFitSplittingAndSwitchingAlgorithm::Run(), lar_content::CosmicRaySplittingAlgorithm::Run(), lar_content::TwoDSlidingFitMultiSplitAlgorithm::SplitCluster(), lar_content::TrackConsolidationAlgorithm::TrackConsolidationAlgorithm(), and lar_content::TwoDSlidingFitSplittingAlgorithm::TwoDSlidingFitSplittingAlgorithm().

{
    if (!m_pCluster)
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    return m_pCluster;
}

float lar_content::TwoDSlidingFitResult::GetCosScatteringAngle

(

const float

rL

)

const

Get scattering angle for a given longitudinal coordinate.

Parameters

rL	the longitudinal coordinate

Definition at line 292 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalFitDirection(), GetGlobalFitPosition(), GetLayerFitHalfWindowLength(), and LongitudinalInterpolation().

Referenced by lar_content::BranchSplittingAlgorithm::FindBestSplitPosition().

{
    const float halfWindowLength(this->GetLayerFitHalfWindowLength());

    CartesianVector firstDirection(0.f,0.f,0.f);
    CartesianVector secondDirection(0.f,0.f,0.f);

    const StatusCode firstStatusCode(this->GetGlobalFitDirection(rL - halfWindowLength, firstDirection));
    const StatusCode secondStatusCode(this->GetGlobalFitDirection(rL + halfWindowLength, secondDirection));

    if (STATUS_CODE_SUCCESS != firstStatusCode)
        throw StatusCodeException(firstStatusCode);

    if (STATUS_CODE_SUCCESS != secondStatusCode)
        throw StatusCodeException(secondStatusCode);

    return firstDirection.GetDotProduct(secondDirection);
}

StatusCode lar_content::TwoDSlidingFitResult::GetExtrapolatedDirection	(	const float	rL,
		pandora::CartesianVector &	direction
	)		const

Get extrapolated direction (beyond span) for a given input coordinate.

Parameters

rL	the input coordinate
position	the extrapolated direction at these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 462 of file LArTwoDSlidingFitResult.cc.

References GetGlobalFitDirection(), GetGlobalMaxLayerDirection(), GetGlobalMinLayerDirection(), GetLayer(), GetMaxLayer(), and GetMinLayer().

{
    const StatusCode statusCode(this->GetGlobalFitDirection(rL, direction));

    if (STATUS_CODE_NOT_FOUND != statusCode)
        return statusCode;

    const int thisLayer(this->GetLayer(rL));
    const int minLayer(this->GetMinLayer());
    const int maxLayer(this->GetMaxLayer());

    if (thisLayer <= minLayer)
    {
        direction = this->GetGlobalMinLayerDirection();
    }
    else if (thisLayer >= maxLayer)
    {
        direction = this->GetGlobalMaxLayerDirection();
    }
    else
    {
        return STATUS_CODE_FAILURE;
    }

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::GetExtrapolatedPosition	(	const float	rL,
		pandora::CartesianVector &	position
	)		const

Get extrapolated position (beyond span) for a given input coordinate.

Parameters

rL	the input coordinate
position	the extrapolated position at these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 433 of file LArTwoDSlidingFitResult.cc.

References GetGlobalFitPosition(), GetGlobalMaxLayerDirection(), GetGlobalMaxLayerPosition(), GetGlobalMinLayerDirection(), GetGlobalMinLayerPosition(), GetL(), GetLayer(), GetMaxLayer(), and GetMinLayer().

Referenced by lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints().

{
    const StatusCode statusCode(this->GetGlobalFitPosition(rL, position));

    if (STATUS_CODE_NOT_FOUND != statusCode)
        return statusCode;

    const int thisLayer(this->GetLayer(rL));
    const int minLayer(this->GetMinLayer());
    const int maxLayer(this->GetMaxLayer());

    if (thisLayer <= minLayer)
    {
        position = (this->GetGlobalMinLayerPosition() + this->GetGlobalMinLayerDirection() * (rL - this->GetL(minLayer)));
    }
    else if (thisLayer >= maxLayer)
    {
        position = (this->GetGlobalMaxLayerPosition() + this->GetGlobalMaxLayerDirection() * (rL - this->GetL(maxLayer)));
    }
    else
    {
        return STATUS_CODE_FAILURE;
    }

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::GetExtrapolatedPositionAtX	(	const float	x,
		pandora::CartesianVector &	position
	)		const

Get extrapolated position (beyond span) for a given input x coordinate.

Parameters

x	the input coordinate
position	the extrapolated position at these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 491 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalFitPositionAtX(), GetGlobalMaxLayerDirection(), GetGlobalMaxLayerPosition(), GetGlobalMinLayerDirection(), GetGlobalMinLayerPosition(), GetMaxLayer(), and GetMinLayer().

Referenced by lar_content::CandidateVertexCreationAlgorithm::CreateEndpointVertex(), and lar_content::ClearTransverseTrackHitsTool::GetTransverseTrackHit3D().

{
    const StatusCode statusCode(this->GetGlobalFitPositionAtX(x, position));

    if (STATUS_CODE_NOT_FOUND != statusCode)
        return statusCode;

    const float minLayerX(this->GetGlobalMinLayerPosition().GetX());
    const float maxLayerX(this->GetGlobalMaxLayerPosition().GetX());

    const int minLayer(this->GetMinLayer());
    const int maxLayer(this->GetMaxLayer());

    const int innerLayer((minLayerX < maxLayerX) ? minLayer : maxLayer);
    const int outerLayer((minLayerX < maxLayerX) ? maxLayer : minLayer);

    const CartesianVector innerVertex((innerLayer == minLayer) ? this->GetGlobalMinLayerPosition() : this->GetGlobalMaxLayerPosition());
    const CartesianVector innerDirection((innerLayer == minLayer) ? this->GetGlobalMinLayerDirection() * -1.f : this->GetGlobalMaxLayerDirection());

    const CartesianVector outerVertex((outerLayer == minLayer) ? this->GetGlobalMinLayerPosition() : this->GetGlobalMaxLayerPosition());
    const CartesianVector outerDirection((outerLayer == minLayer) ? this->GetGlobalMinLayerDirection() * -1.f : this->GetGlobalMaxLayerDirection());

    if (innerDirection.GetX() > -std::numeric_limits<float>::epsilon() || outerDirection.GetX() < +std::numeric_limits<float>::epsilon() ||
        outerVertex.GetX() - innerVertex.GetX() < +std::numeric_limits<float>::epsilon())
    {
        return STATUS_CODE_NOT_FOUND;
    }
    else if (x >= outerVertex.GetX())
    {
        position = outerVertex + outerDirection * ((x - outerVertex.GetX()) / outerDirection.GetX());
    }
    else if (x <= innerVertex.GetX())
    {
        position = innerVertex + innerDirection * ((x - innerVertex.GetX()) / innerDirection.GetX());
    }
    else
    {
        return STATUS_CODE_NOT_FOUND;
    }

    // TODO How to assign an uncertainty on the extrapolated position?
    return STATUS_CODE_SUCCESS;
}

float lar_content::TwoDSlidingFitResult::GetFitRms

(

const float

rL

)

const

Get fit rms for a given longitudinal coordinate.

Parameters

rL	the longitudinal coordinate

Returns

the fit rms

Definition at line 279 of file LArTwoDSlidingFitResult.cc.

References LongitudinalInterpolation().

Referenced by lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), and lar_content::ThreeDSlidingFitResult::GetFitRms().

{
    LayerInterpolation layerInterpolation;
    const StatusCode statusCode(this->LongitudinalInterpolation(rL, layerInterpolation));

    if (STATUS_CODE_SUCCESS != statusCode)
        throw StatusCodeException(statusCode);

    return this->GetFitRms(layerInterpolation);
}

float lar_content::TwoDSlidingFitResult::GetFitRms ( const LayerInterpolation &  layerInterpolation ) const

private

Interpolate a rms between two layers.

Parameters

layerInterpolation

the pair of surrounding layers

Returns

the interpolated rms

Definition at line 837 of file LArTwoDSlidingFitResult.cc.

References lar_content::LayerInterpolation::GetEndLayerIter(), lar_content::LayerInterpolation::GetEndLayerWeight(), lar_content::LayerInterpolation::GetStartLayerIter(), lar_content::LayerInterpolation::GetStartLayerWeight(), and m_layerFitResultMap.

{
    const LayerFitResultMap::const_iterator firstLayerIter(layerInterpolation.GetStartLayerIter());
    const LayerFitResultMap::const_iterator secondLayerIter(layerInterpolation.GetEndLayerIter());

    const float firstWeight(layerInterpolation.GetStartLayerWeight());
    const float secondWeight(layerInterpolation.GetEndLayerWeight());

    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        throw StatusCodeException(STATUS_CODE_FAILURE);

    const float firstLayerRms(firstLayerIter->second.GetRms());

    if (firstLayerIter == secondLayerIter)
        return firstLayerRms;

    const float secondLayerRms(secondLayerIter->second.GetRms());

    if (firstWeight + secondWeight < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    return ((firstLayerRms * firstWeight + secondLayerRms * secondWeight) / (firstWeight + secondWeight));
}

const FitSegment & lar_content::TwoDSlidingFitResult::GetFitSegment

(

const float

rL

)

const

Get fit segment for a given longitudinal coordinate.

Parameters

rL	the longitudinal coordinate

Definition at line 537 of file LArTwoDSlidingFitResult.cc.

References lar_content::FitSegment::GetEndLayer(), GetLayer(), lar_content::FitSegment::GetStartLayer(), and m_fitSegmentList.

Referenced by lar_content::ClearLongitudinalTrackHitsTool::GetLongitudinalTrackHit3D(), and lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions().

{
    int layer(this->GetLayer(rL));

    for (FitSegmentList::const_iterator iter = m_fitSegmentList.begin(), iterEnd = m_fitSegmentList.end(); iter != iterEnd; ++iter)
    {
        const FitSegment &fitSegment = *iter;

        if (layer >= fitSegment.GetStartLayer() && layer <= fitSegment.GetEndLayer())
            return fitSegment;
    }

    throw StatusCodeException(STATUS_CODE_NOT_FOUND);
}

const FitSegmentList & lar_content::TwoDSlidingFitResult::GetFitSegmentList ( ) const

inline

Get the fit segment list.

Returns

the fit segment list

Definition at line 598 of file LArTwoDSlidingFitResult.h.

References m_fitSegmentList.

Referenced by lar_content::ThreeDTransverseTracksAlgorithm::GetFitSegmentTensor().

{
    return m_fitSegmentList;
}

void lar_content::TwoDSlidingFitResult::GetGlobalDirection	(	const float	dTdL,
		pandora::CartesianVector &	direction
	)		const

Get global direction coordinates for given sliding linear fit gradient.

Parameters

dTdL	the transverse coordinate
direction	to receive the direction cartesian vector

Definition at line 193 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), and m_axisIntercept.

Referenced by GetGlobalFitDirection(), GetGlobalMaxLayerDirection(), and GetGlobalMinLayerDirection().

{
    const float pL(1.f / std::sqrt(1.f + dTdL * dTdL));
    const float pT(dTdL / std::sqrt(1.f + dTdL * dTdL));

    CartesianVector globalCoordinates(0.f, 0.f, 0.f);
    this->GetGlobalPosition(pL, pT, globalCoordinates);
    direction = (globalCoordinates - m_axisIntercept).GetUnitVector();
}

pandora::StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitDirection	(	const float	rL,
		pandora::CartesianVector &	direction
	)		const

Get global fit direction for a given longitudinal coordinate.

Parameters

rL	the longitudinal coordinate
direction	the fitted direction at these coordinates

Returns

status code, faster than throwing in regular use-cases

Referenced by lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), GetCosScatteringAngle(), GetExtrapolatedDirection(), lar_content::ThreeDSlidingFitResult::GetGlobalFitDirection(), GetGlobalFitDirectionAtX(), GetGlobalFitPositionListAtX(), and lar_content::ShowerAsymmetryFeatureTool::GetShowerAsymmetryForView().

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalFitDirection ( const LayerInterpolation &  layerInterpolation ) const

private

Interpolate a direction between two layers.

Parameters

layerInterpolation

the pair of surrounding layers

Returns

the interpolated direction

Definition at line 809 of file LArTwoDSlidingFitResult.cc.

References f, lar_content::LayerInterpolation::GetEndLayerIter(), lar_content::LayerInterpolation::GetEndLayerWeight(), GetGlobalDirection(), lar_content::LayerInterpolation::GetStartLayerIter(), lar_content::LayerInterpolation::GetStartLayerWeight(), and m_layerFitResultMap.

{
    const LayerFitResultMap::const_iterator firstLayerIter(layerInterpolation.GetStartLayerIter());
    const LayerFitResultMap::const_iterator secondLayerIter(layerInterpolation.GetEndLayerIter());

    const float firstWeight(layerInterpolation.GetStartLayerWeight());
    const float secondWeight(layerInterpolation.GetEndLayerWeight());

    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        throw StatusCodeException(STATUS_CODE_FAILURE);

    CartesianVector firstLayerDirection(0.f,0.f,0.f);
    this->GetGlobalDirection(firstLayerIter->second.GetGradient(),firstLayerDirection);

    if (firstLayerIter == secondLayerIter)
        return firstLayerDirection;

    CartesianVector secondLayerDirection(0.f,0.f,0.f);
    this->GetGlobalDirection(secondLayerIter->second.GetGradient(),secondLayerDirection);

    if (firstWeight + secondWeight < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    return ((firstLayerDirection * firstWeight + secondLayerDirection * secondWeight).GetUnitVector());
}

StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitDirectionAtX	(	const float	x,
		pandora::CartesianVector &	direction
	)		const

Get global fit direction for a given input x coordinate.

Parameters

x	the input coordinate
direction	the fitted direction at these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 358 of file LArTwoDSlidingFitResult.cc.

References GetGlobalFitDirection(), and TransverseInterpolation().

{
    LayerInterpolationList layerInterpolationList;
    const StatusCode statusCode(this->TransverseInterpolation(x, layerInterpolationList));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    if (layerInterpolationList.size() != 1)
        return STATUS_CODE_NOT_FOUND;

    direction = this->GetGlobalFitDirection(layerInterpolationList.at(0));
    return STATUS_CODE_SUCCESS;
}

pandora::StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitPosition	(	const float	rL,
		pandora::CartesianVector &	position
	)		const

Get global fit position for a given longitudinal coordinate.

Parameters

rL	the longitudinal coordinate
position	the fitted position at these coordinates

Returns

status code, faster than throwing in regular use-cases

Referenced by lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), GetCosScatteringAngle(), GetExtrapolatedPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalFitPosition(), GetGlobalFitPositionAtX(), GetGlobalFitPositionListAtX(), GetGlobalFitProjection(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), and lar_content::CrossGapsAssociationAlgorithm::IsNearCluster().

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalFitPosition ( const LayerInterpolation &  layerInterpolation ) const

private

Interpolate a position between two layers.

Parameters

layerInterpolation

the pair of surrounding layers

Returns

the interpolated position

Definition at line 781 of file LArTwoDSlidingFitResult.cc.

References f, lar_content::LayerInterpolation::GetEndLayerIter(), lar_content::LayerInterpolation::GetEndLayerWeight(), GetGlobalPosition(), lar_content::LayerInterpolation::GetStartLayerIter(), lar_content::LayerInterpolation::GetStartLayerWeight(), and m_layerFitResultMap.

{
    const LayerFitResultMap::const_iterator firstLayerIter(layerInterpolation.GetStartLayerIter());
    const LayerFitResultMap::const_iterator secondLayerIter(layerInterpolation.GetEndLayerIter());

    const float firstWeight(layerInterpolation.GetStartLayerWeight());
    const float secondWeight(layerInterpolation.GetEndLayerWeight());

    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        throw StatusCodeException(STATUS_CODE_FAILURE);

    CartesianVector firstLayerPosition(0.f, 0.f, 0.f);
    this->GetGlobalPosition(firstLayerIter->second.GetL(), firstLayerIter->second.GetFitT(), firstLayerPosition);

    if (firstLayerIter == secondLayerIter)
        return firstLayerPosition;

    CartesianVector secondLayerPosition(0.f, 0.f, 0.f);
    this->GetGlobalPosition(secondLayerIter->second.GetL(), secondLayerIter->second.GetFitT(), secondLayerPosition);

    if (firstWeight + secondWeight < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    return ((firstLayerPosition * firstWeight + secondLayerPosition * secondWeight) * (1.f / (firstWeight + secondWeight)));
}

StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitPositionAtX	(	const float	x,
		pandora::CartesianVector &	position
	)		const

Get global fit position for a given input x coordinate.

Parameters

x	the input coordinate
position	the fitted position at these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 341 of file LArTwoDSlidingFitResult.cc.

References GetGlobalFitPosition(), and TransverseInterpolation().

Referenced by lar_content::TracksCrossingGapsTool::CheckXPositionInGap(), GetExtrapolatedPositionAtX(), lar_content::CrossGapsAssociationAlgorithm::IsNearCluster(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), and lar_content::TracksCrossingGapsTool::PassesGapChecks().

{
    LayerInterpolationList layerInterpolationList;
    const StatusCode statusCode(this->TransverseInterpolation(x, layerInterpolationList));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    if (layerInterpolationList.size() != 1)
        return STATUS_CODE_NOT_FOUND;

    position = this->GetGlobalFitPosition(layerInterpolationList.at(0));
    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitPositionListAtX	(	const float	x,
		pandora::CartesianPointVector &	positionList
	)		const

Get a list of projected positions for a given input x coordinate.

Parameters

x	the input x coordinate
positionList	the output list of positions

Returns

status code, faster than throwing in regular use-cases

Definition at line 384 of file LArTwoDSlidingFitResult.cc.

References GetGlobalFitDirection(), GetGlobalFitPosition(), GetTransverseProjection(), TransverseInterpolation(), and x.

Referenced by lar_content::TrackSplittingTool::CheckSplitPosition(), and lar_content::MultiValuedTransverseTrackHitsTool::GetTransverseTrackHit3D().

{
    LayerInterpolationList layerInterpolationList;
    const StatusCode statusCode(this->TransverseInterpolation(x, layerInterpolationList));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    for (LayerInterpolationList::const_iterator iter = layerInterpolationList.begin(), iterEnd = layerInterpolationList.end();
        iter != iterEnd; ++iter)
    {
        positionList.push_back(this->GetGlobalFitPosition(*iter));
    }

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::GetGlobalFitProjection	(	const pandora::CartesianVector &	inputPosition,
		pandora::CartesianVector &	projectedPosition
	)		const

Get projected position on global fit for a given position vector.

Parameters

inputPosition	the input coordinate
projectedPosition	the projected position on the global fit for these coordinates

Returns

status code, faster than throwing in regular use-cases

Definition at line 375 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalFitPosition(), and GetLocalPosition().

Referenced by lar_content::VertexSplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::MultiValuedLongitudinalTrackHitsTool::GetLongitudinalTrackHit3D(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::ThreeDLongitudinalTracksAlgorithm::ThreeDLongitudinalTracksAlgorithm(), and lar_content::TrackConsolidationAlgorithm::TrackConsolidationAlgorithm().

{
    float rL(0.f), rT(0.f);
    this->GetLocalPosition(inputPosition, rL, rT);
    return this->GetGlobalFitPosition(rL, projectedPosition);
}

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerDirection

(

)

const

Get global direction corresponding to the fit result in maximum fit layer.

Returns

the position

Definition at line 244 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalDirection(), and m_layerFitResultMap.

Referenced by lar_content::CrossGapsAssociationAlgorithm::AreClustersAssociated(), lar_content::LArPointingCluster::BuildPointingCluster(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), GetExtrapolatedDirection(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::TransverseTrackHitsBaseTool::GetTransverseChi2(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), and lar_content::VertexSelectionBaseAlgorithm::SlidingFitData::SlidingFitData().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_reverse_iterator iter = m_layerFitResultMap.rbegin();
    CartesianVector direction(0.f, 0.f, 0.f);
    this->GetGlobalDirection(iter->second.GetGradient(), direction);
    return direction;
}

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalMaxLayerPosition

(

)

const

Get global position corresponding to the fit result in maximum fit layer.

Returns

the position

Definition at line 218 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), and m_layerFitResultMap.

Referenced by lar_content::CrossGapsAssociationAlgorithm::AreClustersAssociated(), lar_content::LArPointingCluster::BuildPointingCluster(), lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointCandidates(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointVertex(), lar_content::CutPfoCharacterisationAlgorithm::CutPfoCharacterisationAlgorithm(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::VertexSplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::TransverseTrackHitsBaseTool::GetTransverseChi2(), lar_content::LongitudinalTrackHitsBaseTool::GetVertexAndEndPositions(), lar_content::CrossGapsAssociationAlgorithm::IsAssociated(), lar_content::CutClusterCharacterisationAlgorithm::IsClearTrack(), lar_content::TracksCrossingGapsTool::IsEndOfCluster(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), lar_content::CosmicRayTrackRecoveryAlgorithm::MatchClusters(), lar_content::TwoDShowerFitFeatureTool::Run(), lar_content::TwoDVertexDistanceFeatureTool::Run(), lar_content::VertexSelectionBaseAlgorithm::SlidingFitData::SlidingFitData(), and lar_content::ThreeDLongitudinalTracksAlgorithm::ThreeDLongitudinalTracksAlgorithm().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_reverse_iterator iter = m_layerFitResultMap.rbegin();
    CartesianVector position(0.f, 0.f, 0.f);
    this->GetGlobalPosition(iter->second.GetL(), iter->second.GetFitT(), position);
    return position;
}

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalMinLayerDirection

(

)

const

Get global direction corresponding to the fit result in minimum fit layer.

Returns

the position

Definition at line 231 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalDirection(), and m_layerFitResultMap.

Referenced by lar_content::CrossGapsAssociationAlgorithm::AreClustersAssociated(), lar_content::LArPointingCluster::BuildPointingCluster(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), GetExtrapolatedDirection(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::TransverseTrackHitsBaseTool::GetTransverseChi2(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), and lar_content::VertexSelectionBaseAlgorithm::SlidingFitData::SlidingFitData().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_iterator iter = m_layerFitResultMap.begin();
    CartesianVector direction(0.f, 0.f, 0.f);
    this->GetGlobalDirection(iter->second.GetGradient(), direction);
    return direction;
}

CartesianVector lar_content::TwoDSlidingFitResult::GetGlobalMinLayerPosition

(

)

const

Get global position corresponding to the fit result in minimum fit layer.

Returns

the position

Definition at line 205 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), and m_layerFitResultMap.

Referenced by lar_content::CrossGapsAssociationAlgorithm::AreClustersAssociated(), lar_content::LArPointingCluster::BuildPointingCluster(), lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::CosmicRaySplittingAlgorithm::ConfirmSplitPosition(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointCandidates(), lar_content::CandidateVertexCreationAlgorithm::CreateEndpointVertex(), lar_content::CutPfoCharacterisationAlgorithm::CutPfoCharacterisationAlgorithm(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::VertexSplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::TransverseTrackHitsBaseTool::GetTransverseChi2(), lar_content::LongitudinalTrackHitsBaseTool::GetVertexAndEndPositions(), lar_content::CrossGapsAssociationAlgorithm::IsAssociated(), lar_content::CutClusterCharacterisationAlgorithm::IsClearTrack(), lar_content::TracksCrossingGapsTool::IsEndOfCluster(), lar_content::LArGeometryHelper::IsXSamplingPointInGap(), lar_content::CosmicRayTrackRecoveryAlgorithm::MatchClusters(), lar_content::TwoDShowerFitFeatureTool::Run(), lar_content::TwoDVertexDistanceFeatureTool::Run(), lar_content::VertexSelectionBaseAlgorithm::SlidingFitData::SlidingFitData(), and lar_content::ThreeDLongitudinalTracksAlgorithm::ThreeDLongitudinalTracksAlgorithm().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_iterator iter = m_layerFitResultMap.begin();
    CartesianVector position(0.f, 0.f, 0.f);
    this->GetGlobalPosition(iter->second.GetL(), iter->second.GetFitT(), position);
    return position;
}

void lar_content::TwoDSlidingFitResult::GetGlobalPosition	(	const float	rL,
		const float	rT,
		pandora::CartesianVector &	position
	)		const

Get global coordinates for given sliding linear fit coordinates.

Parameters

rL	the longitudinal coordinate
rT	the transverse coordinate
position	to receive the position cartesian vector

Definition at line 186 of file LArTwoDSlidingFitResult.cc.

References m_axisDirection, m_axisIntercept, and m_orthoDirection.

Referenced by lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), FindSlidingFitSegments(), GetGlobalDirection(), GetGlobalFitPosition(), GetGlobalMaxLayerPosition(), GetGlobalMinLayerPosition(), GetMinAndMaxCoordinate(), GetTransverseInterpolationWeights(), GetTransverseSurroundingLayers(), and lar_content::CutClusterCharacterisationAlgorithm::IsClearTrack().

{
    position = m_axisIntercept + (m_axisDirection * rL) + (m_orthoDirection * rT);
}

float lar_content::TwoDSlidingFitResult::GetL

(

const int

layer

)

const

Get longitudinal coordinate for a given sliding linear fit layer number.

Parameters

layer

the layer number

Definition at line 156 of file LArTwoDSlidingFitResult.cc.

References f, and m_layerPitch.

Referenced by lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), GetExtrapolatedPosition(), lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), PerformSlidingLinearFit(), lar_content::TwoDSlidingFitMultiSplitAlgorithm::SplitCluster(), and lar_content::ThreeDSlidingFitResult::ThreeDSlidingFitResult().

{
    return (static_cast<float>(layer) + 0.5f) * m_layerPitch;
}

int lar_content::TwoDSlidingFitResult::GetLayer

(

const float

rL

)

const

Get layer number for given sliding linear fit longitudinal coordinate.

Parameters

rL	the longitudinal coordinate

Definition at line 146 of file LArTwoDSlidingFitResult.cc.

References m_layerPitch.

Referenced by lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), FillLayerFitContributionMap(), GetExtrapolatedDirection(), GetExtrapolatedPosition(), GetFitSegment(), lar_content::ThreeDSlidingFitResult::GetGlobalFitDirection(), lar_content::ThreeDSlidingFitResult::GetGlobalFitPosition(), GetLongitudinalSurroundingLayers(), lar_content::LArParticleIdPlugins::LArMuonId::GetMuonTrackWidth(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), GetTransverseSurroundingLayers(), and lar_content::ThreeDKinkBaseTool::GetXSamplingPoint().

{
    if (m_layerPitch < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    return static_cast<int>(std::floor(rL / m_layerPitch));
}

const LayerFitContributionMap & lar_content::TwoDSlidingFitResult::GetLayerFitContributionMap ( ) const

inline

Get the layer fit contribution map.

Returns

the layer fit contribution map

Definition at line 591 of file LArTwoDSlidingFitResult.h.

References m_layerFitContributionMap.

Referenced by PerformSlidingLinearFit().

{
    return m_layerFitContributionMap;
}

unsigned int lar_content::TwoDSlidingFitResult::GetLayerFitHalfWindow ( ) const

inline

Get the layer fit half window.

Returns

the layer fit half window

Definition at line 549 of file LArTwoDSlidingFitResult.h.

References m_layerFitHalfWindow.

Referenced by lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::NeutrinoHierarchyAlgorithm::PfoInfo::operator=(), PerformSlidingLinearFit(), and lar_content::NeutrinoHierarchyAlgorithm::PfoInfo::PfoInfo().

{
    return m_layerFitHalfWindow;
}

float lar_content::TwoDSlidingFitResult::GetLayerFitHalfWindowLength

(

)

const

Get the layer fit half window length.

Returns

the layer fit half window length

Definition at line 119 of file LArTwoDSlidingFitResult.cc.

References m_layerFitHalfWindow, and m_layerPitch.

Referenced by lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), and GetCosScatteringAngle().

{
    return (static_cast<float>(m_layerFitHalfWindow)) * m_layerPitch;
}

const LayerFitResultMap & lar_content::TwoDSlidingFitResult::GetLayerFitResultMap ( ) const

inline

Get the layer fit result map.

Returns

the layer fit result map

Definition at line 584 of file LArTwoDSlidingFitResult.h.

References m_layerFitResultMap.

Referenced by lar_content::TwoDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ThreeDLinearFitFeatureTool::CalculateVariablesSlidingLinearFit(), lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), lar_content::CutPfoCharacterisationAlgorithm::CutPfoCharacterisationAlgorithm(), lar_content::KinkSplittingAlgorithm::FindBestSplitPosition(), FindSlidingFitSegments(), lar_content::LArParticleIdPlugins::LArMuonId::GetMuonTrackWidth(), lar_content::CutClusterCharacterisationAlgorithm::GetShowerFitWidth(), and lar_content::CutClusterCharacterisationAlgorithm::IsClearTrack().

{
    return m_layerFitResultMap;
}

float lar_content::TwoDSlidingFitResult::GetLayerPitch ( ) const

inline

Get the layer pitch, units cm.

Returns

the layer pitch

Definition at line 556 of file LArTwoDSlidingFitResult.h.

References m_layerPitch.

Referenced by lar_content::CosmicRayTaggingTool::CRCandidate::CalculateFitVariables(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::TracksCrossingGapsTool::IsEndOfCluster(), lar_content::NeutrinoHierarchyAlgorithm::PfoInfo::operator=(), and lar_content::NeutrinoHierarchyAlgorithm::PfoInfo::PfoInfo().

{
    return m_layerPitch;
}

void lar_content::TwoDSlidingFitResult::GetLocalDirection	(	const pandora::CartesianVector &	direction,
		float &	dTdL
	)		const

Get local sliding fit gradient for a given global direction.

Parameters

direction	the direction cartesian vector
dTdL	to receive the local gradient

Definition at line 173 of file LArTwoDSlidingFitResult.cc.

References f, GetLocalPosition(), and m_axisIntercept.

Referenced by lar_content::ThreeDSlidingFitResult::GetGlobalFitDirection().

{
    float pL(0.f), pT(0.f);
    this->GetLocalPosition((direction + m_axisIntercept), pL, pT);

    if (std::fabs(pL) < std::numeric_limits<float>::epsilon())
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    dTdL = pT / pL;
}

void lar_content::TwoDSlidingFitResult::GetLocalPosition	(	const pandora::CartesianVector &	position,
		float &	rL,
		float &	rT
	)		const

Get local sliding fit coordinates for a given global position.

Parameters

position	the position cartesian vector
rL	to receive the longitudinal coordinate
rT	to receive the transverse coordinate

Definition at line 163 of file LArTwoDSlidingFitResult.cc.

References m_axisDirection, m_axisIntercept, and m_orthoDirection.

Referenced by lar_content::OvershootSplittingAlgorithm::BuildSortedIntersectionMap(), lar_content::ConeClusterMopUpAlgorithm::ClusterMopUp(), FillLayerFitContributionMap(), lar_content::BranchSplittingAlgorithm::FindBestSplitPosition(), lar_content::DeltaRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::CrossedTrackSplittingAlgorithm::FindBestSplitPosition(), lar_content::CosmicRaySplittingAlgorithm::FindBestSplitPosition(), lar_content::ThreeDSlidingFitResult::GetGlobalFitPosition(), GetGlobalFitProjection(), GetLocalDirection(), lar_content::ClearLongitudinalTrackHitsTool::GetLongitudinalTrackHit3D(), lar_content::LArParticleIdPlugins::LArMuonId::GetMuonTrackWidth(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), lar_content::ShowerAsymmetryFeatureTool::GetShowerAsymmetryForView(), lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), lar_content::CrossGapsAssociationAlgorithm::IsNearCluster(), lar_content::TwoDSlidingFitMultiSplitAlgorithm::SplitCluster(), lar_content::TrackConsolidationAlgorithm::TrackConsolidationAlgorithm(), and lar_content::TwoDSlidingFitSplittingAlgorithm::TwoDSlidingFitSplittingAlgorithm().

{
    const CartesianVector displacement(position - m_axisIntercept);

    rL = displacement.GetDotProduct(m_axisDirection);
    rT = displacement.GetDotProduct(m_orthoDirection);
}

void lar_content::TwoDSlidingFitResult::GetLongitudinalInterpolationWeights ( const float  rL, const LayerFitResultMap::const_iterator &  firstLayerIter, const LayerFitResultMap::const_iterator &  secondLayerIter, double &  firstWeight, double &  secondWeight  ) const

private

Get interpolation weights for layers surrounding a specified longitudinal position.

Parameters

rL	the longitudinal coordinate
firstLayerIter	the iterator for the layer below the input coordinate
secondLayerIter	the iterator for the layer above the input coordinate
firstWeight	the weight assigned to the layer below the input coordinate
secondWeight	the weight assigned to the layer above the input coordinate

Definition at line 1062 of file LArTwoDSlidingFitResult.cc.

References m_layerFitResultMap.

Referenced by LongitudinalInterpolation().

{
    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    const double deltaL(rL - firstLayerIter->second.GetL());
    const double deltaLLayers(secondLayerIter->second.GetL() - firstLayerIter->second.GetL());

    if (std::fabs(deltaLLayers) > std::numeric_limits<float>::epsilon())
    {
        firstWeight = 1. - deltaL / deltaLLayers;
        secondWeight = deltaL / deltaLLayers;
    }
    else
    {
        firstWeight = 0.5;
        secondWeight = 0.5;
    }
}

StatusCode lar_content::TwoDSlidingFitResult::GetLongitudinalSurroundingLayers ( const float  rL, LayerFitResultMap::const_iterator &  firstLayerIter, LayerFitResultMap::const_iterator &  secondLayerIter  ) const

private

Get iterators for layers surrounding the specified longitudinal position.

Parameters

rL	the longitudinal coordinate
firstLayerIter	to receive the iterator for the layer just below the input coordinate
secondLayerIter	to receive the iterator for the layer just above the input coordinate

Returns

status code, faster than throwing in regular use-cases

Definition at line 921 of file LArTwoDSlidingFitResult.cc.

References GetLayer(), and m_layerFitResultMap.

Referenced by LongitudinalInterpolation().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    // Get minimum, maximum and input layers
    const int minLayer(m_layerFitResultMap.begin()->first), maxLayer(m_layerFitResultMap.rbegin()->first);
    const int thisLayer(this->GetLayer(rL));

    // Allow special case of single-layer sliding fit result
    if (minLayer == thisLayer && thisLayer == maxLayer)
    {
        firstLayerIter = m_layerFitResultMap.find(minLayer);
        secondLayerIter = m_layerFitResultMap.find(maxLayer);
        return STATUS_CODE_SUCCESS;
    }

    // For multi-layer sliding fit result, set start layer and bail out if out of range
    const int startLayer((thisLayer >= maxLayer) ? thisLayer - 1 : thisLayer);

    if ((startLayer < minLayer) || (startLayer >= maxLayer))
        return STATUS_CODE_NOT_FOUND;

    // First layer iterator
    firstLayerIter = m_layerFitResultMap.end();

    for (int iLayer = startLayer; iLayer >= minLayer; --iLayer)
    {
        firstLayerIter = m_layerFitResultMap.find(iLayer);

        if (m_layerFitResultMap.end() != firstLayerIter)
            break;
    }

    if (m_layerFitResultMap.end() == firstLayerIter)
        return STATUS_CODE_NOT_FOUND;

    // Second layer iterator
    secondLayerIter = m_layerFitResultMap.end();

    for (int iLayer = startLayer + 1; iLayer <= maxLayer; ++iLayer)
    {
        secondLayerIter = m_layerFitResultMap.find(iLayer);

        if (m_layerFitResultMap.end() != secondLayerIter)
            break;
    }

    if (m_layerFitResultMap.end() == secondLayerIter)
        return STATUS_CODE_NOT_FOUND;

    return STATUS_CODE_SUCCESS;
}

int lar_content::TwoDSlidingFitResult::GetMaxLayer

(

)

const

Get the maximum occupied layer in the sliding fit.

Parameters

the	maximum occupied layer in the sliding fit

Definition at line 136 of file LArTwoDSlidingFitResult.cc.

References m_layerFitResultMap.

Referenced by lar_content::LArPointingCluster::BuildPointingCluster(), GetExtrapolatedDirection(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::ThreeDSlidingFitResult::GetMaxLayer(), lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), lar_content::TwoDSlidingFitMultiSplitAlgorithm::SplitCluster(), lar_content::ThreeDTransverseTracksAlgorithm::ThreeDTransverseTracksAlgorithm(), lar_content::BoundedClusterMopUpAlgorithm::XSampling::XSampling(), and lar_content::ThreeDShowersAlgorithm::XSampling::XSampling().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    return m_layerFitResultMap.rbegin()->first;
}

float lar_content::TwoDSlidingFitResult::GetMaxLayerRms

(

)

const

Get rms at maximum layer.

Returns

the rms

Definition at line 268 of file LArTwoDSlidingFitResult.cc.

References m_layerFitResultMap.

Referenced by lar_content::LArPointingCluster::BuildPointingCluster(), and lar_content::ThreeDSlidingFitResult::GetMaxLayerRms().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_reverse_iterator iter = m_layerFitResultMap.rbegin();
    return iter->second.GetRms();
}

void lar_content::TwoDSlidingFitResult::GetMinAndMaxCoordinate ( const bool  isX, float &  min, float &  max  ) const

private

Get the minimum and maximum x or z coordinates associated with the sliding fit.

Parameters

isX	whether to provide extremal x or z coordinates
to	receive the min coordinate value
to	receive the max coordinate value

Definition at line 761 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), m_layerFitResultMap, max, and min.

Referenced by GetMinAndMaxX(), and GetMinAndMaxZ().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    min = std::numeric_limits<float>::max();
    max = -std::numeric_limits<float>::max();

    for (LayerFitResultMap::const_iterator iter = m_layerFitResultMap.begin(), iterEnd = m_layerFitResultMap.end(); iter != iterEnd; ++iter)
    {
        CartesianVector globalPosition(0.f, 0.f, 0.f);
        this->GetGlobalPosition(iter->second.GetL(), iter->second.GetFitT(), globalPosition);
        const float coordinate(isX ? globalPosition.GetX() : globalPosition.GetZ());
        min = std::min(min, coordinate);
        max = std::max(max, coordinate);
    }
}

void lar_content::TwoDSlidingFitResult::GetMinAndMaxX ( float &  minX, float &  maxX  ) const

inline

Get the minimum and maximum x coordinates associated with the sliding fit.

Parameters

to	receive the min x value
to	receive the max x value

Definition at line 605 of file LArTwoDSlidingFitResult.h.

References GetMinAndMaxCoordinate().

Referenced by lar_content::TwoDSlidingShowerFitResult::GetShowerEdges(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), lar_content::BoundedClusterMopUpAlgorithm::XSampling::XSampling(), and lar_content::ThreeDShowersAlgorithm::XSampling::XSampling().

{
    return this->GetMinAndMaxCoordinate(true, minX, maxX);
}

void lar_content::TwoDSlidingFitResult::GetMinAndMaxZ ( float &  minZ, float &  maxZ  ) const

inline

Get the minimum and maximum z coordinates associated with the sliding fit.

Parameters

to	receive the min z value
to	receive the max z value

Definition at line 612 of file LArTwoDSlidingFitResult.h.

References GetMinAndMaxCoordinate().

Referenced by lar_content::TwoDSlidingShowerFitResult::GetShowerEdges().

{
    return this->GetMinAndMaxCoordinate(false, minZ, maxZ);
}

int lar_content::TwoDSlidingFitResult::GetMinLayer

(

)

const

Get the minimum occupied layer in the sliding fit.

Parameters

the	minimum occupied layer in the sliding fit

Definition at line 126 of file LArTwoDSlidingFitResult.cc.

References m_layerFitResultMap.

Referenced by lar_content::LArPointingCluster::BuildPointingCluster(), GetExtrapolatedDirection(), GetExtrapolatedPosition(), GetExtrapolatedPositionAtX(), lar_content::ThreeDSlidingFitResult::GetMinLayer(), lar_content::CandidateVertexCreationAlgorithm::GetSpacepoints(), lar_content::ThreeDKinkBaseTool::GetXSamplingPoint(), lar_content::ThreeDTransverseTracksAlgorithm::ThreeDTransverseTracksAlgorithm(), lar_content::BoundedClusterMopUpAlgorithm::XSampling::XSampling(), and lar_content::ThreeDShowersAlgorithm::XSampling::XSampling().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    return m_layerFitResultMap.begin()->first;
}

float lar_content::TwoDSlidingFitResult::GetMinLayerRms

(

)

const

Get rms at minimum layer.

Returns

the rms

Definition at line 257 of file LArTwoDSlidingFitResult.cc.

References m_layerFitResultMap.

Referenced by lar_content::LArPointingCluster::BuildPointingCluster(), and lar_content::ThreeDSlidingFitResult::GetMinLayerRms().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_iterator iter = m_layerFitResultMap.begin();
    return iter->second.GetRms();
}

const pandora::CartesianVector & lar_content::TwoDSlidingFitResult::GetOrthoDirection ( ) const

inline

Get the orthogonal direction vector.

Returns

the orthogonal direction vector

Definition at line 577 of file LArTwoDSlidingFitResult.h.

References m_orthoDirection.

Referenced by lar_content::TwoDSlidingShowerFitResult::GetShowerEdges().

{
    return m_orthoDirection;
}

void lar_content::TwoDSlidingFitResult::GetTransverseInterpolationWeights ( const float  x, const LayerFitResultMap::const_iterator &  firstLayerIter, const LayerFitResultMap::const_iterator &  secondLayerIter, double &  firstWeight, double &  secondWeight  ) const

private

Get interpolation weights for layers surrounding a specified transverse position.

Parameters

x	the transverse coordinate
firstLayerIter	the iterator for the layer below the input coordinate
secondLayerIter	the iterator for the layer above the input coordinate
firstWeight	the weight assigned to the layer below the input coordinate
firstWeight	the weight assigned to the layer above the input coordinate

Definition at line 1085 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), and m_layerFitResultMap.

Referenced by TransverseInterpolation().

{
    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    CartesianVector firstLayerPosition(0.f, 0.f, 0.f);
    CartesianVector secondLayerPosition(0.f, 0.f, 0.f);

    this->GetGlobalPosition(firstLayerIter->second.GetL(), firstLayerIter->second.GetFitT(), firstLayerPosition);
    this->GetGlobalPosition(secondLayerIter->second.GetL(), secondLayerIter->second.GetFitT(), secondLayerPosition);

    const double deltaP(x - firstLayerPosition.GetX());
    const double deltaPLayers(secondLayerPosition.GetX() - firstLayerPosition.GetX());

    if (std::fabs(deltaPLayers) > std::numeric_limits<float>::epsilon())
    {
        firstWeight = 1. - deltaP / deltaPLayers;
        secondWeight = deltaP / deltaPLayers;
    }
    else
    {
        firstWeight = 0.5;
        secondWeight = 0.5;
    }
}

pandora::StatusCode lar_content::TwoDSlidingFitResult::GetTransverseProjection	(	const float	x,
		const FitSegment &	fitSegment,
		pandora::CartesianVector &	position
	)		const

Get projected position for a given input x coordinate and fit segment.

Parameters

x	the input x coordinate
fitSegment	the portion of sliding linear fit
position	the output position

Returns

status code, faster than throwing in regular use-cases

Referenced by GetGlobalFitPositionListAtX(), lar_content::ClearLongitudinalTrackHitsTool::GetLongitudinalTrackHit3D(), lar_content::ThreeDTrackFragmentsAlgorithm::GetProjectedPositions(), and lar_content::ThreeDTransverseTracksAlgorithm::GetSegmentOverlap().

pandora::StatusCode lar_content::TwoDSlidingFitResult::GetTransverseProjection	(	const float	x,
		const FitSegment &	fitSegment,
		pandora::CartesianVector &	position,
		pandora::CartesianVector &	direction
	)		const

Get projected position and direction for a given input x coordinate and fit segment.

Parameters

x	the input x coordinate
fitSegment	the portion of sliding linear fit
position	the output position
position	the output direction

Returns

status code, faster than throwing in regular use-cases

StatusCode lar_content::TwoDSlidingFitResult::GetTransverseSurroundingLayers ( const float  x, const int  minLayer, const int  maxLayer, LayerFitResultMap::const_iterator &  firstLayerIter, LayerFitResultMap::const_iterator &  secondLayerIter  ) const

private

Get iterators for layers surrounding a specified transverse position.

Parameters

x	the transverse coordinate
minLayer	the minimum allowed layer
maxLayer	the maximum allowed layer
firstLayerIter	to receive the iterator for the layer just below the input coordinate
secondLayerIter	to receive the iterator for the layer just above the input coordinate

Returns

status code, faster than throwing in regular use-cases

Definition at line 978 of file LArTwoDSlidingFitResult.cc.

References f, GetGlobalPosition(), GetLayer(), m_layerFitResultMap, max, min, and x.

Referenced by TransverseInterpolation().

{
    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);

    LayerFitResultMap::const_iterator minLayerIter = m_layerFitResultMap.find(minLayer);
    if (m_layerFitResultMap.end() == minLayerIter)
        throw StatusCodeException(STATUS_CODE_FAILURE);

    LayerFitResultMap::const_iterator maxLayerIter = m_layerFitResultMap.find(maxLayer);
    if (m_layerFitResultMap.end() == maxLayerIter)
        throw StatusCodeException(STATUS_CODE_FAILURE);

    CartesianVector minPosition(0.f, 0.f, 0.f), maxPosition(0.f, 0.f, 0.f);
    this->GetGlobalPosition(minLayerIter->second.GetL(), minLayerIter->second.GetFitT(), minPosition);
    this->GetGlobalPosition(maxLayerIter->second.GetL(), maxLayerIter->second.GetFitT(), maxPosition);

    if ((std::fabs(maxPosition.GetX() - minPosition.GetX()) < std::numeric_limits<float>::epsilon()))
        return STATUS_CODE_NOT_FOUND;

    // Find start layer
    const float minL(minLayerIter->second.GetL());
    const float maxL(maxLayerIter->second.GetL());
    const float startL(minL + (maxL - minL) * (x - minPosition.GetX()) / (maxPosition.GetX() - minPosition.GetX()));
    const int startLayer(std::max(minLayer, std::min(maxLayer, this->GetLayer(startL))));

    // Find nearest layer iterator to start layer
    LayerFitResultMap::const_iterator startLayerIter = m_layerFitResultMap.end();
    CartesianVector startLayerPosition(0.f, 0.f, 0.f);

    for (int iLayer = startLayer; iLayer <= maxLayer; ++iLayer)
    {
        startLayerIter = m_layerFitResultMap.find(iLayer);

        if (m_layerFitResultMap.end() != startLayerIter)
            break;
    }

    if (m_layerFitResultMap.end() == startLayerIter)
        return STATUS_CODE_NOT_FOUND;

    this->GetGlobalPosition(startLayerIter->second.GetL(), startLayerIter->second.GetFitT(), startLayerPosition);

    const bool startIsAhead((startLayerPosition.GetX() - x) > std::numeric_limits<float>::epsilon());
    const bool increasesWithLayers(maxPosition.GetX() > minPosition.GetX());
    const int increment = ((startIsAhead == increasesWithLayers) ? -1 : +1);

    // Find surrounding layer iterators
    // (Second layer iterator comes immediately after the fit has crossed the target X coordinate
    //  and first layer iterator comes immediately before the second layer iterator).
    firstLayerIter = m_layerFitResultMap.end();
    secondLayerIter = m_layerFitResultMap.end();

    CartesianVector firstLayerPosition(0.f, 0.f, 0.f);
    CartesianVector secondLayerPosition(0.f, 0.f, 0.f);

    for (int iLayer = startLayerIter->first; (iLayer >= minLayer) && (iLayer <= maxLayer); iLayer += increment)
    {
        LayerFitResultMap::const_iterator tempIter = m_layerFitResultMap.find(iLayer);
        if (m_layerFitResultMap.end() == tempIter)
            continue;

        firstLayerIter = secondLayerIter;
        firstLayerPosition = secondLayerPosition;
        secondLayerIter = tempIter;

        this->GetGlobalPosition(secondLayerIter->second.GetL(), secondLayerIter->second.GetFitT(), secondLayerPosition);
        const bool isAhead(secondLayerPosition.GetX() > x);

        if (startIsAhead != isAhead)
            break;

        firstLayerIter = m_layerFitResultMap.end();
    }

    if (m_layerFitResultMap.end() == firstLayerIter || m_layerFitResultMap.end() == secondLayerIter)
        return STATUS_CODE_NOT_FOUND;

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::LongitudinalInterpolation ( const float  rL, LayerInterpolation &  layerInterpolation  ) const

private

Get the pair of layers surrounding a specified longitudinal position.

Parameters

rL	the longitudinal coordinate
layerInterpolation	to receive the populated layer interpolation object

Returns

status code, faster than throwing in regular use-cases

Definition at line 863 of file LArTwoDSlidingFitResult.cc.

References GetLongitudinalInterpolationWeights(), and GetLongitudinalSurroundingLayers().

Referenced by GetCosScatteringAngle(), and GetFitRms().

{
    double firstWeight(0.), secondWeight(0.);
    LayerFitResultMap::const_iterator firstLayerIter, secondLayerIter;

    const StatusCode statusCode(this->GetLongitudinalSurroundingLayers(rL, firstLayerIter, secondLayerIter));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    this->GetLongitudinalInterpolationWeights(rL, firstLayerIter, secondLayerIter, firstWeight, secondWeight);
    layerInterpolation = LayerInterpolation(firstLayerIter, secondLayerIter, firstWeight, secondWeight);

    return STATUS_CODE_SUCCESS;
}

void lar_content::TwoDSlidingFitResult::PerformSlidingLinearFit ( )

private

Perform the sliding linear fit.

Definition at line 610 of file LArTwoDSlidingFitResult.cc.

References GetL(), GetLayerFitContributionMap(), GetLayerFitHalfWindow(), m_layerFitContributionMap, m_layerFitResultMap, and m_layerPitch.

Referenced by TwoDSlidingFitResult().

{
    if (!m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_FAILURE);

    if ((m_layerPitch < std::numeric_limits<float>::epsilon()) || (m_layerFitContributionMap.empty()))
        throw StatusCodeException(STATUS_CODE_INVALID_PARAMETER);

    unsigned int slidingNPoints(0);
    double slidingSumT(0.), slidingSumL(0.), slidingSumTT(0.), slidingSumLT(0.), slidingSumLL(0.);

    const LayerFitContributionMap &layerFitContributionMap(this->GetLayerFitContributionMap());
    const int innerLayer(layerFitContributionMap.begin()->first);
    const int layerFitHalfWindow(static_cast<int>(this->GetLayerFitHalfWindow()));

    for (int iLayer = innerLayer; iLayer < innerLayer + layerFitHalfWindow; ++iLayer)
    {
        LayerFitContributionMap::const_iterator lyrIter = layerFitContributionMap.find(iLayer);

        if (layerFitContributionMap.end() != lyrIter)
        {
            slidingSumT += lyrIter->second.GetSumT();
            slidingSumL += lyrIter->second.GetSumL();
            slidingSumTT += lyrIter->second.GetSumTT();
            slidingSumLT += lyrIter->second.GetSumLT();
            slidingSumLL += lyrIter->second.GetSumLL();
            slidingNPoints += lyrIter->second.GetNPoints();
        }
    }

    const int outerLayer(layerFitContributionMap.rbegin()->first);

    for (int iLayer = innerLayer; iLayer <= outerLayer; ++iLayer)
    {
        const int fwdLayer(iLayer + layerFitHalfWindow);
        LayerFitContributionMap::const_iterator fwdIter = layerFitContributionMap.find(fwdLayer);

        if (layerFitContributionMap.end() != fwdIter)
        {
            slidingSumT += fwdIter->second.GetSumT();
            slidingSumL += fwdIter->second.GetSumL();
            slidingSumTT += fwdIter->second.GetSumTT();
            slidingSumLT += fwdIter->second.GetSumLT();
            slidingSumLL += fwdIter->second.GetSumLL();
            slidingNPoints += fwdIter->second.GetNPoints();
        }

        const int bwdLayer(iLayer - layerFitHalfWindow - 1);
        LayerFitContributionMap::const_iterator bwdIter = layerFitContributionMap.find(bwdLayer);

        if (layerFitContributionMap.end() != bwdIter)
        {
            slidingSumT -= bwdIter->second.GetSumT();
            slidingSumL -= bwdIter->second.GetSumL();
            slidingSumTT -= bwdIter->second.GetSumTT();
            slidingSumLT -= bwdIter->second.GetSumLT();
            slidingSumLL -= bwdIter->second.GetSumLL();
            slidingNPoints -= bwdIter->second.GetNPoints();
        }

        // require three points for meaningful results
        if (slidingNPoints <= 2)
            continue;

        // only fill the result map if there is an entry in the contribution map
        if (layerFitContributionMap.end() == layerFitContributionMap.find(iLayer))
            continue;

        const double denominator(slidingSumLL - slidingSumL * slidingSumL / static_cast<double>(slidingNPoints));

        if (std::fabs(denominator) < std::numeric_limits<float>::epsilon())
            continue;

        const double gradient((slidingSumLT - slidingSumL * slidingSumT / static_cast<double>(slidingNPoints)) / denominator);
        const double intercept((slidingSumLL * slidingSumT / static_cast<double>(slidingNPoints) - slidingSumL * slidingSumLT / static_cast<double>(slidingNPoints)) / denominator);
        double variance((slidingSumTT - 2. * intercept * slidingSumT - 2. * gradient * slidingSumLT + intercept * intercept * static_cast<double>(slidingNPoints) + 2. * gradient * intercept * slidingSumL + gradient * gradient * slidingSumLL) / (1. + gradient * gradient));

        if (variance < -std::numeric_limits<float>::epsilon())
            continue;

        if (variance < std::numeric_limits<float>::epsilon())
            variance = 0.;

        const double rms(std::sqrt(variance / static_cast<double>(slidingNPoints)));
        const double l(this->GetL(iLayer));
        const double fitT(intercept + gradient * l);

        const LayerFitResult layerFitResult(l, fitT, gradient, rms);
        (void) m_layerFitResultMap.insert(LayerFitResultMap::value_type(iLayer, layerFitResult));
    }

    if (m_layerFitResultMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_INITIALIZED);
}

StatusCode lar_content::TwoDSlidingFitResult::TransverseInterpolation ( const float  x, const FitSegment &  fitSegment, LayerInterpolation &  layerInterpolation  ) const

private

Get the surrounding pair of layers for a specified transverse position and fit segment.

Parameters

x	the input coordinate
fitSegment	the fit segment
layerInterpolation	to receive the populated layer interpolation object

Returns

status code, faster than throwing in regular use-cases

Definition at line 881 of file LArTwoDSlidingFitResult.cc.

References lar_content::FitSegment::GetEndLayer(), lar_content::FitSegment::GetStartLayer(), GetTransverseInterpolationWeights(), and GetTransverseSurroundingLayers().

Referenced by GetGlobalFitDirectionAtX(), GetGlobalFitPositionAtX(), GetGlobalFitPositionListAtX(), and TransverseInterpolation().

{
    double firstWeight(0.), secondWeight(0.);
    LayerFitResultMap::const_iterator firstLayerIter, secondLayerIter;

    const StatusCode statusCode(this->GetTransverseSurroundingLayers(x, fitSegment.GetStartLayer(), fitSegment.GetEndLayer(), firstLayerIter, secondLayerIter));

    if (STATUS_CODE_SUCCESS != statusCode)
        return statusCode;

    this->GetTransverseInterpolationWeights(x, firstLayerIter, secondLayerIter, firstWeight, secondWeight);
    layerInterpolation = LayerInterpolation(firstLayerIter, secondLayerIter, firstWeight, secondWeight);

    return STATUS_CODE_SUCCESS;
}

StatusCode lar_content::TwoDSlidingFitResult::TransverseInterpolation ( const float  x, LayerInterpolationList &  layerInterpolationList  ) const

private

Get the a list of surrounding layer pairs for a specified transverse position.

Parameters

x	the input coordinate
layerInterpolationList	the output list of layer interpolation objects

Returns

status code, faster than throwing in regular use-cases

Definition at line 899 of file LArTwoDSlidingFitResult.cc.

References m_fitSegmentList, and TransverseInterpolation().

{
    for (FitSegmentList::const_iterator iter = m_fitSegmentList.begin(), iterEnd = m_fitSegmentList.end(); iter != iterEnd; ++iter)
    {
        LayerInterpolation layerInterpolation;
        const StatusCode statusCode(this->TransverseInterpolation(x, *iter, layerInterpolation));

        if (STATUS_CODE_SUCCESS == statusCode)
        {
            layerInterpolationList.push_back(layerInterpolation);
        }
        else if (STATUS_CODE_NOT_FOUND != statusCode)
        {
            return statusCode;
        }
    }

    return STATUS_CODE_SUCCESS;
}

Member Data Documentation

pandora::CartesianVector lar_content::TwoDSlidingFitResult::m_axisDirection

private

The axis direction vector.

Definition at line 536 of file LArTwoDSlidingFitResult.h.

Referenced by CalculateAxes(), FillLayerFitContributionMap(), GetAxisDirection(), GetGlobalPosition(), GetLocalPosition(), and TwoDSlidingFitResult().

pandora::CartesianVector lar_content::TwoDSlidingFitResult::m_axisIntercept

private

The axis intercept position.

Definition at line 535 of file LArTwoDSlidingFitResult.h.

Referenced by CalculateAxes(), GetAxisIntercept(), GetGlobalDirection(), GetGlobalPosition(), GetLocalDirection(), GetLocalPosition(), and TwoDSlidingFitResult().

FitSegmentList lar_content::TwoDSlidingFitResult::m_fitSegmentList

private

The fit segment list.

Definition at line 540 of file LArTwoDSlidingFitResult.h.

Referenced by FindSlidingFitSegments(), GetFitSegment(), GetFitSegmentList(), and TransverseInterpolation().

LayerFitContributionMap lar_content::TwoDSlidingFitResult::m_layerFitContributionMap

private

The layer fit contribution map.

Definition at line 539 of file LArTwoDSlidingFitResult.h.

Referenced by FillLayerFitContributionMap(), GetLayerFitContributionMap(), PerformSlidingLinearFit(), and TwoDSlidingFitResult().

unsigned int lar_content::TwoDSlidingFitResult::m_layerFitHalfWindow

private

The layer fit half window.

Definition at line 533 of file LArTwoDSlidingFitResult.h.

Referenced by GetLayerFitHalfWindow(), GetLayerFitHalfWindowLength(), and TwoDSlidingFitResult().

LayerFitResultMap lar_content::TwoDSlidingFitResult::m_layerFitResultMap

private

The layer fit result map.

Definition at line 538 of file LArTwoDSlidingFitResult.h.

Referenced by GetFitRms(), GetGlobalFitDirection(), GetGlobalFitPosition(), GetGlobalMaxLayerDirection(), GetGlobalMaxLayerPosition(), GetGlobalMinLayerDirection(), GetGlobalMinLayerPosition(), GetLayerFitResultMap(), GetLongitudinalInterpolationWeights(), GetLongitudinalSurroundingLayers(), GetMaxLayer(), GetMaxLayerRms(), GetMinAndMaxCoordinate(), GetMinLayer(), GetMinLayerRms(), GetTransverseInterpolationWeights(), GetTransverseSurroundingLayers(), and PerformSlidingLinearFit().

float lar_content::TwoDSlidingFitResult::m_layerPitch

private

The layer pitch, units cm.

Definition at line 534 of file LArTwoDSlidingFitResult.h.

Referenced by FillLayerFitContributionMap(), GetL(), GetLayer(), GetLayerFitHalfWindowLength(), GetLayerPitch(), PerformSlidingLinearFit(), and TwoDSlidingFitResult().

pandora::CartesianVector lar_content::TwoDSlidingFitResult::m_orthoDirection

private

The orthogonal direction vector.

Definition at line 537 of file LArTwoDSlidingFitResult.h.

Referenced by CalculateAxes(), FillLayerFitContributionMap(), GetGlobalPosition(), GetLocalPosition(), GetOrthoDirection(), and TwoDSlidingFitResult().

const pandora::Cluster* lar_content::TwoDSlidingFitResult::m_pCluster

private

The address of the cluster.

Definition at line 532 of file LArTwoDSlidingFitResult.h.

Referenced by GetCluster(), and TwoDSlidingFitResult().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArObjects/LArTwoDSlidingFitResult.h
• /cvmfs/larsoft.opensciencegrid.org/products/larpandoracontent/v03_14_03/source/larpandoracontent/LArObjects/LArTwoDSlidingFitResult.cc