TwoViewThreeDKinkTool.h

Go to the documentation of this file.

#ifndef TWO_VIEW_THREE_D_KINK_TOOL_H
#define TWO_VIEW_THREE_D_KINK_TOOL_H 1

#include "larpandoracontent/LArObjects/LArPointingCluster.h"

#include "larpandoracontent/LArThreeDReco/LArTwoViewMatching/TwoViewTransverseTracksAlgorithm.h"

namespace lar_content
{

class TwoViewThreeDKinkTool : public TransverseMatrixTool
{
public:
    TwoViewThreeDKinkTool();

    virtual ~TwoViewThreeDKinkTool();

    bool Run(TwoViewTransverseTracksAlgorithm *const pAlgorithm, MatrixType &overlapMatrix);

private:
    class Modification
    {
    public:
        SplitPositionMap m_splitPositionMap;     
        ClusterMergeMap m_clusterMergeMap;       
        pandora::ClusterList m_affectedClusters; 
    };

    typedef std::vector<Modification> ModificationList;

    class Particle
    {
    public:
        Particle(const MatrixType::Element &elementA, const MatrixType::Element &elementB);

        const pandora::Cluster *m_pClusterA;      
        const pandora::Cluster *m_pClusterB;      
        const pandora::Cluster *m_pCommonCluster; 
    };

    void GetModifications(TwoViewTransverseTracksAlgorithm *const pAlgorithm, const MatrixType &overlapMatrix, ModificationList &modificationList) const;

    bool PassesElementCuts(MatrixType::ElementList::const_iterator eIter, const pandora::ClusterSet &usedClusters) const;

    void SelectMatrixElements(MatrixType::ElementList::const_iterator eIter, const MatrixType::ElementList &elementList,
        const pandora::ClusterSet &usedClusters, IteratorList &iteratorList) const;

    void GetIteratorListModifications(
        TwoViewTransverseTracksAlgorithm *const pAlgorithm, const IteratorList &iteratorList, ModificationList &modificationList) const;

    static bool IsALowestInX(const LArPointingCluster &pointingClusterA, const LArPointingCluster &pointingClusterB);

    bool IsThreeDKink(TwoViewTransverseTracksAlgorithm *const pAlgorithm, const Particle &particle,
        const pandora::CartesianVector &splitPosition, const bool isALowestInX) const;

    float GetXSamplingPoint(const pandora::CartesianVector &splitPosition1, const bool isForwardInX, const TwoDSlidingFitResult &fitResult1,
        const TwoDSlidingFitResult &fitResult2) const;

    bool ApplyChanges(TwoViewTransverseTracksAlgorithm *const pAlgorithm, const ModificationList &modificationList) const;

    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    float m_minXOverlapFraction;            
    float m_minMatchingScore;               
    float m_minLocallyMatchedFraction;      
    float m_minLongitudinalImpactParameter; 
    int m_nLayersForKinkSearch;             
    float m_additionalXStepForKinkSearch;   
    float m_maxTransverseImpactParameter;   
    float m_minImpactParameterCosTheta;     
    float m_cosThetaCutForKinkSearch;       
};

} // namespace lar_content

#endif // #ifndef TWO_VIEW_THREE_D_KINK_TOOL_H