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