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