lar_content::LongitudinalTrackHitsBaseTool Class Reference

LongitudinalTrackHitsBaseTool class. More...

#include <LongitudinalTrackHitsBaseTool.h>

Inheritance diagram for lar_content::LongitudinalTrackHitsBaseTool:

Public Member Functions
	LongitudinalTrackHitsBaseTool ()
	Default constructor. More...
Public Member Functions inherited from lar_content::TrackHitsBaseTool
	TrackHitsBaseTool ()
	Default constructor. More...
virtual void	Run (ThreeDHitCreationAlgorithm *const pAlgorithm, const pandora::ParticleFlowObject *const pPfo, const pandora::CaloHitVector &inputTwoDHits, ProtoHitVector &protoHitVector)
	Run the algorithm tool. More...
Public Member Functions inherited from lar_content::HitCreationBaseTool
	HitCreationBaseTool ()
	Default constructor. More...
virtual	~HitCreationBaseTool ()
	Destructor. More...

Protected Member Functions
virtual void	GetLongitudinalTrackHit3D (const MatchedSlidingFitMap &matchedSlidingFitMap, const pandora::CartesianVector &vtx3D, const pandora::CartesianVector &end3D, ProtoHit &protoHit) const =0
	Get the three dimensional position using a provided two dimensional calo hit and sliding linear fits in the other two views. More...
virtual void	GetTrackHits3D (const pandora::CaloHitVector &inputTwoDHits, const MatchedSlidingFitMap &matchedSlidingFitMap, ProtoHitVector &protoHitVector) const
	Calculate 3D hits from an input list of 2D hits. More...
virtual pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)
Protected Member Functions inherited from lar_content::TrackHitsBaseTool
virtual void	BuildSlidingFitMap (const pandora::ParticleFlowObject *const pPfo, MatchedSlidingFitMap &matchedSlidingFitMap) const
	Calculate sliding fit results for clusters from each view. More...
Protected Member Functions inherited from lar_content::HitCreationBaseTool
virtual void	GetBestPosition3D (const pandora::HitType hitType1, const pandora::HitType hitType2, const pandora::CartesianPointVector &fitPositionList1, const pandora::CartesianPointVector &fitPositionList2, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views. More...
virtual void	GetBestPosition3D (const pandora::HitType hitType1, const pandora::HitType hitType2, const pandora::CartesianVector &fitPosition1, const pandora::CartesianVector &fitPosition2, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and candidate fit positions from the other two views. More...
virtual void	GetBestPosition3D (const pandora::HitType hitType, const pandora::CartesianVector &fitPosition, ProtoHit &protoHit) const
	Get the three dimensional position using a provided two dimensional calo hit and a candidate fit position from another view. More...

Private Member Functions
void	GetVertexAndEndPositions (const MatchedSlidingFitMap &inputSlidingFitMap, MatchedSlidingFitMap &outputSlidingFitMap, pandora::CartesianVector &outputVtx3D, pandora::CartesianVector &outputEnd3D) const
	Get reconstructed vertex and end positions for this 3D track. More...
void	UpdateBestPosition (const pandora::HitType hitType1, const pandora::HitType hitType2, const pandora::CartesianVector &vtx1, const pandora::CartesianVector &vtx2, pandora::CartesianVector &bestVtx, float &bestChi2) const
	Combine two 2D coordinates to give a 3D coordinate. More...

Private Attributes
float	m_vtxDisplacementCutSquared
float	m_minTrackLengthSquared

Additional Inherited Members
Public Types inherited from lar_content::HitCreationBaseTool
typedef ThreeDHitCreationAlgorithm::ProtoHit	ProtoHit
typedef ThreeDHitCreationAlgorithm::ProtoHitVector	ProtoHitVector
typedef ThreeDHitCreationAlgorithm::TrajectorySample	TrajectorySample
Protected Types inherited from lar_content::TrackHitsBaseTool
typedef std::map < pandora::HitType, TwoDSlidingFitResult >	MatchedSlidingFitMap
Protected Attributes inherited from lar_content::TrackHitsBaseTool
unsigned int	m_minViews
	The minimum number of views required for building hits. More...
unsigned int	m_slidingFitWindow
	The layer window for the sliding linear fits. More...
Protected Attributes inherited from lar_content::HitCreationBaseTool
double	m_sigmaX2
	The sigmaX squared value, for calculation of chi2 deltaX term. More...
double	m_chiSquaredCut
	The chi squared cut (accept only values below the cut value) More...

Detailed Description

LongitudinalTrackHitsBaseTool class.

Definition at line 19 of file LongitudinalTrackHitsBaseTool.h.

Constructor & Destructor Documentation

lar_content::LongitudinalTrackHitsBaseTool::LongitudinalTrackHitsBaseTool

(

)

Default constructor.

Definition at line 21 of file LongitudinalTrackHitsBaseTool.cc.

                                                             :
    m_vtxDisplacementCutSquared(5.f * 5.f),
    m_minTrackLengthSquared(7.5f * 7.5f)
{
}

Member Function Documentation

virtual void lar_content::LongitudinalTrackHitsBaseTool::GetLongitudinalTrackHit3D ( const MatchedSlidingFitMap &  matchedSlidingFitMap, const pandora::CartesianVector &  vtx3D, const pandora::CartesianVector &  end3D, ProtoHit &  protoHit  ) const

protectedpure virtual

Get the three dimensional position using a provided two dimensional calo hit and sliding linear fits in the other two views.

Parameters

matchedSlidingFitMap	map of sliding fit results from each view
vtx3D	the 3D vertex position
end3D	the 3D end position
protoHit	to receive the populated proto hit

Implemented in lar_content::ClearLongitudinalTrackHitsTool, and lar_content::MultiValuedLongitudinalTrackHitsTool.

void lar_content::LongitudinalTrackHitsBaseTool::GetTrackHits3D ( const pandora::CaloHitVector &  inputTwoDHits, const MatchedSlidingFitMap &  matchedSlidingFitMap, ProtoHitVector &  protoHitVector  ) const

protectedvirtual

Calculate 3D hits from an input list of 2D hits.

Parameters

pAlgorithm	the hit creation algorithm
inputTwoDHits	the input vector of 2D hits
matchedSlidingFitMap	the group of sliding fit results
protoHitVector	to receive the new three dimensional proto hits

Implements lar_content::TrackHitsBaseTool.

Definition at line 29 of file LongitudinalTrackHitsBaseTool.cc.

{
    MatchedSlidingFitMap matchedSlidingFitMap;
    CartesianVector vtx3D(0.f, 0.f, 0.f), end3D(0.f, 0.f, 0.f);
    this->GetVertexAndEndPositions(inputSlidingFitMap, matchedSlidingFitMap, vtx3D, end3D);

    for (const CaloHit *const pCaloHit2D : inputTwoDHits)
    {
        try
        {
            ProtoHit protoHit(pCaloHit2D);
            this->GetLongitudinalTrackHit3D(matchedSlidingFitMap, vtx3D, end3D, protoHit);

            if (protoHit.IsPositionSet() && (protoHit.GetChi2() < m_chiSquaredCut))
                protoHitVector.push_back(protoHit);
        }
        catch (StatusCodeException &)
        {
        }
    }
}

void lar_content::LongitudinalTrackHitsBaseTool::GetVertexAndEndPositions ( const MatchedSlidingFitMap &  inputSlidingFitMap, MatchedSlidingFitMap &  outputSlidingFitMap, pandora::CartesianVector &  outputVtx3D, pandora::CartesianVector &  outputEnd3D  ) const

private

Get reconstructed vertex and end positions for this 3D track.

Parameters

inputSlidingFitMap	input map of sliding fit results from each view
outputSlidingFitMap	output map of clean sliding fit results from each view
outputVtx3D	reconstructed start position of 3D track
outputEnd3D	reconstructed end position of 3D track

Definition at line 54 of file LongitudinalTrackHitsBaseTool.cc.

{
    // TODO Tidy up: The code below is quite repetitive...
    MatchedSlidingFitMap::const_iterator iterU = inputSlidingFitMap.find(TPC_VIEW_U);
    const bool foundU(inputSlidingFitMap.end() != iterU);

    MatchedSlidingFitMap::const_iterator iterV = inputSlidingFitMap.find(TPC_VIEW_V);
    const bool foundV(inputSlidingFitMap.end() != iterV);

    MatchedSlidingFitMap::const_iterator iterW = inputSlidingFitMap.find(TPC_VIEW_W);
    const bool foundW(inputSlidingFitMap.end() != iterW);

    bool useU(false), useV(false), useW(false);
    float bestChi2(std::numeric_limits<float>::max());

    for (unsigned int iPermutation = 0; iPermutation < 4; ++iPermutation)
    {
        const bool isForwardU((1 == iPermutation) ? false : true);
        const bool isForwardV((2 == iPermutation) ? false : true);
        const bool isForwardW((3 == iPermutation) ? false : true);

        CartesianVector vtxU(0.f, 0.f, 0.f), endU(0.f, 0.f, 0.f);
        CartesianVector vtxV(0.f, 0.f, 0.f), endV(0.f, 0.f, 0.f);
        CartesianVector vtxW(0.f, 0.f, 0.f), endW(0.f, 0.f, 0.f);

        if (foundU)
        {
            const TwoDSlidingFitResult &slidingFitResultU = iterU->second;
            vtxU = (isForwardU ? slidingFitResultU.GetGlobalMinLayerPosition() : slidingFitResultU.GetGlobalMaxLayerPosition());
            endU = (isForwardU ? slidingFitResultU.GetGlobalMaxLayerPosition() : slidingFitResultU.GetGlobalMinLayerPosition());
        }

        if (foundV)
        {
            const TwoDSlidingFitResult &slidingFitResultV = iterV->second;
            vtxV = (isForwardV ? slidingFitResultV.GetGlobalMinLayerPosition() : slidingFitResultV.GetGlobalMaxLayerPosition());
            endV = (isForwardV ? slidingFitResultV.GetGlobalMaxLayerPosition() : slidingFitResultV.GetGlobalMinLayerPosition());
        }

        if (foundW)
        {
            const TwoDSlidingFitResult &slidingFitResultW = iterW->second;
            vtxW = (isForwardW ? slidingFitResultW.GetGlobalMinLayerPosition() : slidingFitResultW.GetGlobalMaxLayerPosition());
            endW = (isForwardW ? slidingFitResultW.GetGlobalMaxLayerPosition() : slidingFitResultW.GetGlobalMinLayerPosition());
        }

        CartesianVector vtx3D(0.f, 0.f, 0.f), end3D(0.f, 0.f, 0.f);
        float vtxChi2(std::numeric_limits<float>::max()), endChi2(std::numeric_limits<float>::max());

        if (foundU && foundV)
        {
            this->UpdateBestPosition(TPC_VIEW_U, TPC_VIEW_V, vtxU, vtxV, vtx3D, vtxChi2);
            this->UpdateBestPosition(TPC_VIEW_U, TPC_VIEW_V, endU, endV, end3D, endChi2);
        }

        if (foundV && foundW)
        {
            this->UpdateBestPosition(TPC_VIEW_V, TPC_VIEW_W, vtxV, vtxW, vtx3D, vtxChi2);
            this->UpdateBestPosition(TPC_VIEW_V, TPC_VIEW_W, endV, endW, end3D, endChi2);
        }

        if (foundW && foundU)
        {
            this->UpdateBestPosition(TPC_VIEW_W, TPC_VIEW_U, vtxW, vtxU, vtx3D, vtxChi2);
            this->UpdateBestPosition(TPC_VIEW_W, TPC_VIEW_U, endW, endU, end3D, endChi2);
        }

        bool matchedU(false), matchedV(false), matchedW(false);
        unsigned int matchedViews(0);

        if (foundU)
        {
            const CartesianVector projVtxU(LArGeometryHelper::ProjectPosition(this->GetPandora(), vtx3D, TPC_VIEW_U));
            const CartesianVector projEndU(LArGeometryHelper::ProjectPosition(this->GetPandora(), end3D, TPC_VIEW_U));

            if ((endU - vtxU).GetMagnitudeSquared() > m_minTrackLengthSquared &&
                (projVtxU - vtxU).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projVtxU - endU).GetMagnitudeSquared()) &&
                (projEndU - endU).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projEndU - vtxU).GetMagnitudeSquared()))
            {
                matchedU = true;
                ++matchedViews;
            }
        }

        if (foundV)
        {
            const CartesianVector projVtxV(LArGeometryHelper::ProjectPosition(this->GetPandora(), vtx3D, TPC_VIEW_V));
            const CartesianVector projEndV(LArGeometryHelper::ProjectPosition(this->GetPandora(), end3D, TPC_VIEW_V));

            if ((endV - vtxV).GetMagnitudeSquared() > m_minTrackLengthSquared &&
                (projVtxV - vtxV).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projVtxV - endV).GetMagnitudeSquared()) &&
                (projEndV - endV).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projEndV - vtxV).GetMagnitudeSquared()))
            {
                matchedV = true;
                ++matchedViews;
            }
        }

        if (foundW)
        {
            const CartesianVector projVtxW(LArGeometryHelper::ProjectPosition(this->GetPandora(), vtx3D, TPC_VIEW_W));
            const CartesianVector projEndW(LArGeometryHelper::ProjectPosition(this->GetPandora(), end3D, TPC_VIEW_W));

            if ((endW - vtxW).GetMagnitudeSquared() > m_minTrackLengthSquared &&
                (projVtxW - vtxW).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projVtxW - endW).GetMagnitudeSquared()) &&
                (projEndW - endW).GetMagnitudeSquared() < std::min(m_vtxDisplacementCutSquared, (projEndW - vtxW).GetMagnitudeSquared()))
            {
                matchedW = true;
                ++matchedViews;
            }
        }

        if (matchedViews < 2)
            continue;

        if (vtxChi2 + endChi2 < bestChi2)
        {
            useU = matchedU;
            useV = matchedV;
            useW = matchedW;

            bestVtx3D = vtx3D;
            bestEnd3D = end3D;
            bestChi2 = vtxChi2 + endChi2;
        }
    }

    if (useU)
        outputSlidingFitMap.insert(MatchedSlidingFitMap::value_type(iterU->first, iterU->second));

    if (useV)
        outputSlidingFitMap.insert(MatchedSlidingFitMap::value_type(iterV->first, iterV->second));

    if (useW)
        outputSlidingFitMap.insert(MatchedSlidingFitMap::value_type(iterW->first, iterW->second));

    if (outputSlidingFitMap.empty())
        throw StatusCodeException(STATUS_CODE_NOT_FOUND);
}

StatusCode lar_content::LongitudinalTrackHitsBaseTool::ReadSettings ( const pandora::TiXmlHandle  xmlHandle )

protectedvirtual

Reimplemented from lar_content::TrackHitsBaseTool.

Definition at line 214 of file LongitudinalTrackHitsBaseTool.cc.

{
    float vtxDisplacementCut = std::sqrt(m_vtxDisplacementCutSquared);
    PANDORA_RETURN_RESULT_IF_AND_IF(
        STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "VertexDisplacementCut", vtxDisplacementCut));
    m_vtxDisplacementCutSquared = vtxDisplacementCut * vtxDisplacementCut;

    float minTrackLength = std::sqrt(m_minTrackLengthSquared);
    PANDORA_RETURN_RESULT_IF_AND_IF(STATUS_CODE_SUCCESS, STATUS_CODE_NOT_FOUND, !=, XmlHelper::ReadValue(xmlHandle, "MinTrackLength", minTrackLength));
    m_minTrackLengthSquared = minTrackLength * minTrackLength;

    return TrackHitsBaseTool::ReadSettings(xmlHandle);
}

void lar_content::LongitudinalTrackHitsBaseTool::UpdateBestPosition ( const pandora::HitType  hitType1, const pandora::HitType  hitType2, const pandora::CartesianVector &  vtx1, const pandora::CartesianVector &  vtx2, pandora::CartesianVector &  bestVtx, float &  bestChi2  ) const

private

Combine two 2D coordinates to give a 3D coordinate.

Parameters

hitType1	the view corresponding to the first position
hitType2	the view corresponding to the second position
vtx1	the first position
vtx2	the second position
bestVtx	the combined vertex position
bestChi2	the chi-squared from the combination

Definition at line 197 of file LongitudinalTrackHitsBaseTool.cc.

{
    CartesianVector mergedVtx(0.f, 0.f, 0.f);
    float mergedChi2(std::numeric_limits<float>::max());

    LArGeometryHelper::MergeTwoPositions3D(this->GetPandora(), hitType1, hitType2, vtx1, vtx2, mergedVtx, mergedChi2);

    if (mergedChi2 < bestChi2)
    {
        bestVtx = mergedVtx;
        bestChi2 = mergedChi2;
    }
}

Member Data Documentation

float lar_content::LongitudinalTrackHitsBaseTool::m_minTrackLengthSquared

private

Definition at line 70 of file LongitudinalTrackHitsBaseTool.h.

float lar_content::LongitudinalTrackHitsBaseTool::m_vtxDisplacementCutSquared

private

Definition at line 69 of file LongitudinalTrackHitsBaseTool.h.

---

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

• LongitudinalTrackHitsBaseTool.h
• LongitudinalTrackHitsBaseTool.cc