RecursivePfoMopUpAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_RECURSIVE_PFO_MOP_UP_ALGORITHM_H
#define LAR_RECURSIVE_PFO_MOP_UP_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

namespace lar_content
{

class RecursivePfoMopUpAlgorithm : public pandora::Algorithm
{
public:
    RecursivePfoMopUpAlgorithm();

private:
    typedef std::vector<unsigned int> ClusterNumHitsList;

    struct PfoMergeStats
    {
    public:
        PfoMergeStats(const ClusterNumHitsList &numClusterHits, const float trackScore);

        const ClusterNumHitsList m_numClusterHits; 
        const float m_trackScore;                  
    };

    typedef std::vector<PfoMergeStats> PfoMergeStatsList;

    static bool PfoMergeStatsComp(const PfoMergeStats &lhs, const PfoMergeStats &rhs);

    PfoMergeStatsList GetPfoMergeStats() const;

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

    unsigned int m_maxIterations;            
    pandora::StringVector m_pfoListNames;    
    pandora::StringVector m_mopUpAlgorithms; 
};

//------------------------------------------------------------------------------------------------------------------------------------------

inline RecursivePfoMopUpAlgorithm::RecursivePfoMopUpAlgorithm() :
    m_maxIterations(10)
{
}

//------------------------------------------------------------------------------------------------------------------------------------------

inline RecursivePfoMopUpAlgorithm::PfoMergeStats::PfoMergeStats(const ClusterNumHitsList &numClusterHits, const float trackScore) :
    m_numClusterHits(numClusterHits),
    m_trackScore(trackScore)
{
}

//------------------------------------------------------------------------------------------------------------------------------------------

inline bool RecursivePfoMopUpAlgorithm::PfoMergeStatsComp(
    const RecursivePfoMopUpAlgorithm::PfoMergeStats &lhs, const RecursivePfoMopUpAlgorithm::PfoMergeStats &rhs)
{
    return ((lhs.m_numClusterHits == rhs.m_numClusterHits) && (std::abs(lhs.m_trackScore - rhs.m_trackScore) < std::numeric_limits<float>::epsilon()));
}

} // namespace lar_content

#endif // #ifndef LAR_RECURSIVE_PFO_MOP_UP_ALGORITHM_H