OvershootSplittingAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_OVERSHOOT_SPLITTING_ALGORITHM_H
#define LAR_OVERSHOOT_SPLITTING_ALGORITHM_H 1

#include "larpandoracontent/LArTwoDReco/LArClusterSplitting/TwoDSlidingFitMultiSplitAlgorithm.h"

namespace lar_content
{

class OvershootSplittingAlgorithm : public TwoDSlidingFitMultiSplitAlgorithm
{
public:
    OvershootSplittingAlgorithm();

private:
    void GetListOfCleanClusters(const pandora::ClusterList *const pClusterList, pandora::ClusterVector &clusterVector) const;
    void FindBestSplitPositions(const TwoDSlidingFitResultMap &slidingFitResultMap, ClusterPositionMap &clusterSplittingMap) const;

    typedef std::pair<float, pandora::CartesianVector> MyTrajectoryPoint;
    typedef std::vector<MyTrajectoryPoint> MyTrajectoryPointList;

    void BuildIntersectionMap(const TwoDSlidingFitResultMap &slidingFitResultMap, ClusterPositionMap &clusterIntersectionMap) const;

    void BuildSortedIntersectionMap(const TwoDSlidingFitResultMap &slidingFitResultMap, const ClusterPositionMap &clusterIntersectionMap,
        ClusterPositionMap &sortedIntersectionMap) const;

    void PopulateSplitPositionMap(const ClusterPositionMap &sortedIntersectionMap, ClusterPositionMap &clusterSplittingMap) const;

    static bool SortByHitProjection(const MyTrajectoryPoint &lhs, const MyTrajectoryPoint &rhs);

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

    float m_minClusterLength;      
    float m_maxClusterSeparation;  
    float m_minVertexDisplacement; 
    float m_maxIntersectDisplacement;
    float m_minSplitDisplacement;
};

} // namespace lar_content

#endif // #ifndef LAR_OVERSHOOT_SPLITTING_ALGORITHM_H