TwoViewTransverseTracksAlgorithm.h

Go to the documentation of this file.

 
#ifndef LAR_TWO_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H
#define LAR_TWO_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H 1
 
#include "Pandora/Algorithm.h"
#include "Pandora/AlgorithmTool.h"
 
#include "larpandoracontent/LArObjects/LArDiscreteProbabilityVector.h"
#include "larpandoracontent/LArObjects/LArTrackTwoViewOverlapResult.h"
 
#include "larpandoracontent/LArThreeDReco/LArThreeDBase/NViewTrackMatchingAlgorithm.h"
#include "larpandoracontent/LArThreeDReco/LArThreeDBase/TwoViewMatchingControl.h"
 
#include <random>
 
namespace lar_content
{
 
class TransverseMatrixTool;
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
class TwoViewTransverseTracksAlgorithm : public NViewTrackMatchingAlgorithm<TwoViewMatchingControl<TwoViewTransverseOverlapResult>>
{
public:
    typedef NViewTrackMatchingAlgorithm<TwoViewMatchingControl<TwoViewTransverseOverlapResult>> BaseAlgorithm;
    typedef std::set<unsigned int> UIntSet;
 
    TwoViewTransverseTracksAlgorithm();
 
private:
    void CalculateOverlapResult(const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, const pandora::Cluster *const);
 
    pandora::StatusCode CalculateOverlapResult(
        const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, TwoViewTransverseOverlapResult &overlapResult);
 
    unsigned int CalculateNumberOfLocallyMatchingSamplingPoints(const DiscreteProbabilityVector &discreteProbabilityVector1,
        const DiscreteProbabilityVector &discreteProbabilityVector2, std::mt19937 &randomNumberGenerator);
 
    float GetPrimaryAxisDotDriftAxis(const pandora::Cluster *const pCluster);
 
    void ExamineOverlapContainer();
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);
 
    typedef std::vector<TransverseMatrixTool *> MatrixToolVector;
    MatrixToolVector m_algorithmToolVector; 
 
    unsigned int m_nMaxMatrixToolRepeats;     
    unsigned int m_downsampleFactor;          
    unsigned int m_minSamples;                
    unsigned int m_nPermutations;             
    float m_localMatchingScoreThreshold;      
    float m_maxDotProduct;                    
    float m_minOverallMatchingScore;          
    float m_minOverallLocallyMatchedFraction; 
    std::mt19937 m_randomNumberGenerator;     
};
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
class TransverseMatrixTool : public pandora::AlgorithmTool
{
public:
    typedef TwoViewTransverseTracksAlgorithm::MatchingType::MatrixType MatrixType;
    typedef std::vector<MatrixType::ElementList::const_iterator> IteratorList;
 
    virtual bool Run(TwoViewTransverseTracksAlgorithm *const pAlgorithm, MatrixType &overlapMatrix) = 0;
};
 
} // namespace lar_content
 
#endif // #ifndef LAR_TWO_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H