lar_content::TwoViewAmbiguousDeltaRayTool Class Reference

TwoViewAmbiguousDeltaRayTool class. More...

#include <TwoViewAmbiguousDeltaRayTool.h>

Inheritance diagram for lar_content::TwoViewAmbiguousDeltaRayTool:

Public Types
typedef std::vector< pandora::HitType >	HitTypeVector
Public Types inherited from lar_content::DeltaRayMatrixTool
typedef TwoViewDeltaRayMatchingAlgorithm::MatchingType::MatrixType	MatrixType
typedef std::vector< MatrixType::ElementList::const_iterator >	IteratorList

Public Member Functions
	TwoViewAmbiguousDeltaRayTool ()
	Default constructor. More...

Private Member Functions
bool	Run (TwoViewDeltaRayMatchingAlgorithm *const pAlgorithm, MatrixType &overlapMatrix)
	Run the algorithm tool. More...
pandora::StatusCode	ReadSettings (const pandora::TiXmlHandle xmlHandle)
void	ExamineConnectedElements (MatrixType &overlapMatrix) const
	Identify ambiguous matches (e.g. 3:2:1) and, if possible, create pfos out of the best 1:1:1 cluster match. More...
bool	PickOutGoodMatches (const MatrixType::ElementList &elementList) const
	Identify the best 1:1:1 match in a group of connected elements and from it create a pfo. More...

Additional Inherited Members
Public Attributes inherited from lar_content::DeltaRayMatrixTool
TwoViewDeltaRayMatchingAlgorithm *	m_pParentAlgorithm
	Address of the parent matching algorithm. More...

Detailed Description

TwoViewAmbiguousDeltaRayTool class.

Definition at line 19 of file TwoViewAmbiguousDeltaRayTool.h.

Member Typedef Documentation

typedef std::vector<pandora::HitType> lar_content::TwoViewAmbiguousDeltaRayTool::HitTypeVector

Definition at line 22 of file TwoViewAmbiguousDeltaRayTool.h.

Constructor & Destructor Documentation

lar_content::TwoViewAmbiguousDeltaRayTool::TwoViewAmbiguousDeltaRayTool

(

)

Default constructor.

Definition at line 18 of file TwoViewAmbiguousDeltaRayTool.cc.

{
}

Member Function Documentation

void lar_content::TwoViewAmbiguousDeltaRayTool::ExamineConnectedElements ( MatrixType &  overlapMatrix ) const

private

Identify ambiguous matches (e.g. 3:2:1) and, if possible, create pfos out of the best 1:1:1 cluster match.

Parameters

overlapMatrix

the overlap matrix

Definition at line 39 of file TwoViewAmbiguousDeltaRayTool.cc.

{
    bool particleCreated(false);

    do
    {
        particleCreated = false;

        ClusterVector sortedKeyClusters;
        overlapMatrix.GetSortedKeyClusters(sortedKeyClusters);

        ClusterSet usedKeyClusters;

        for (const Cluster *const pKeyCluster : sortedKeyClusters)
        {
            if (usedKeyClusters.count(pKeyCluster))
                continue;

            MatrixType::ElementList elementList;
            overlapMatrix.GetConnectedElements(pKeyCluster, true, elementList);

            for (const MatrixType::Element &element : elementList)
                usedKeyClusters.insert(element.GetCluster1());

            if (this->PickOutGoodMatches(elementList))
            {
                particleCreated = true;
                break;
            }
        }
    } while (particleCreated);
}

bool lar_content::TwoViewAmbiguousDeltaRayTool::PickOutGoodMatches ( const MatrixType::ElementList &  elementList ) const

private

Identify the best 1:1:1 match in a group of connected elements and from it create a pfo.

Parameters

elementList

the tensor element list

Returns

whether any particles were created

Definition at line 74 of file TwoViewAmbiguousDeltaRayTool.cc.

{
    unsigned int highestHitCount(0);
    float bestChiSquared(std::numeric_limits<float>::max());
    auto bestElementIter(elementList.end());

    for (auto iter = elementList.begin(); iter != elementList.end(); ++iter)
    {
        const MatrixType::Element &element(*iter);
        const float chiSquared(element.GetOverlapResult().GetReducedChiSquared());
        const Cluster *const pCluster1(element.GetCluster1()), *const pCluster2(element.GetCluster2());
        const unsigned int hitSum(pCluster1->GetNCaloHits() + pCluster2->GetNCaloHits());

        if ((hitSum > highestHitCount) || ((hitSum == highestHitCount) && (chiSquared < bestChiSquared)))
        {
            bestChiSquared = chiSquared;
            highestHitCount = hitSum;
            bestElementIter = iter;
        }
    }

    if (bestElementIter != elementList.end())
    {
        m_pParentAlgorithm->CreatePfo(*bestElementIter);
        return true;
    }

    return false;
}

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

private

Definition at line 106 of file TwoViewAmbiguousDeltaRayTool.cc.

{
    return STATUS_CODE_SUCCESS;
}

bool lar_content::TwoViewAmbiguousDeltaRayTool::Run ( TwoViewDeltaRayMatchingAlgorithm *const  pAlgorithm, MatrixType &  matrixTensor  )

privatevirtual

Run the algorithm tool.

Parameters

pAlgorithm	address of the calling algorithm
overlapTensor	the overlap tensor

Returns

whether changes have been made by the tool

Implements lar_content::DeltaRayMatrixTool.

Definition at line 24 of file TwoViewAmbiguousDeltaRayTool.cc.

{
    m_pParentAlgorithm = pAlgorithm;

    if (PandoraContentApi::GetSettings(*m_pParentAlgorithm)->ShouldDisplayAlgorithmInfo())
        std::cout << "----> Running Algorithm Tool: " << this->GetInstanceName() << ", " << this->GetType() << std::endl;

    this->ExamineConnectedElements(overlapMatrix);

    // ATTN: Prevent tensor tool loop running again
    return false;
}

---

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

• larpandoracontent/larpandoracontent/LArThreeDReco/LArCosmicRay/TwoViewAmbiguousDeltaRayTool.h
• larpandoracontent/larpandoracontent/LArThreeDReco/LArCosmicRay/TwoViewAmbiguousDeltaRayTool.cc