ThreeViewTransverseTracksAlgorithm.h

Go to the documentation of this file.

 
#ifndef LAR_THREE_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H
#define LAR_THREE_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H 1
 
#include "Pandora/Algorithm.h"
#include "Pandora/AlgorithmTool.h"
 
#include "larpandoracontent/LArObjects/LArTrackOverlapResult.h"
 
#include "larpandoracontent/LArThreeDReco/LArThreeDBase/NViewTrackMatchingAlgorithm.h"
#include "larpandoracontent/LArThreeDReco/LArThreeDBase/ThreeViewMatchingControl.h"
 
namespace lar_content
{
 
class TransverseTensorTool;
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
class ThreeViewTransverseTracksAlgorithm : public NViewTrackMatchingAlgorithm<ThreeViewMatchingControl<TransverseOverlapResult>>
{
public:
    typedef NViewTrackMatchingAlgorithm<ThreeViewMatchingControl<TransverseOverlapResult>> BaseAlgorithm;
 
    ThreeViewTransverseTracksAlgorithm();
 
private:
    typedef std::map<unsigned int, TransverseOverlapResult> FitSegmentToOverlapResultMap;
    typedef std::map<unsigned int, FitSegmentToOverlapResultMap> FitSegmentMatrix;
    typedef std::map<unsigned int, FitSegmentMatrix> FitSegmentTensor;
 
    void CalculateOverlapResult(const pandora::Cluster *const pClusterU, const pandora::Cluster *const pClusterV, const pandora::Cluster *const pClusterW);
 
    pandora::StatusCode CalculateOverlapResult(const pandora::Cluster *const pClusterU, const pandora::Cluster *const pClusterV,
        const pandora::Cluster *const pClusterW, TransverseOverlapResult &overlapResult);
 
    void GetFitSegmentTensor(const TwoDSlidingFitResult &slidingFitResultU, const TwoDSlidingFitResult &slidingFitResultV,
        const TwoDSlidingFitResult &slidingFitResultW, FitSegmentTensor &fitSegmentTensor) const;
 
    pandora::StatusCode GetSegmentOverlap(const FitSegment &fitSegmentU, const FitSegment &fitSegmentV, const FitSegment &fitSegmentW,
        const TwoDSlidingFitResult &slidingFitResultU, const TwoDSlidingFitResult &slidingFitResultV,
        const TwoDSlidingFitResult &slidingFitResultW, TransverseOverlapResult &transverseOverlapResult) const;
 
    void GetBestOverlapResult(const FitSegmentTensor &fitSegmentTensor, TransverseOverlapResult &bestTransverseOverlapResult) const;
 
    void GetPreviousOverlapResults(const unsigned int indexU, const unsigned int indexV, const unsigned int indexW,
        FitSegmentTensor &fitSegmentSumTensor, TransverseOverlapResultVector &transverseOverlapResultVector) const;
 
    void ExamineOverlapContainer();
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);
 
    typedef std::vector<TransverseTensorTool *> TensorToolVector;
    TensorToolVector m_algorithmToolVector; 
 
    unsigned int m_nMaxTensorToolRepeats;   
    unsigned int m_maxFitSegmentIndex;      
    float m_pseudoChi2Cut;                  
    float m_minSegmentMatchedFraction;      
    unsigned int m_minSegmentMatchedPoints; 
    float m_minOverallMatchedFraction;      
    unsigned int m_minOverallMatchedPoints; 
    float m_minSamplingPointsPerLayer;      
};
 
//------------------------------------------------------------------------------------------------------------------------------------------
 
class TransverseTensorTool : public pandora::AlgorithmTool
{
public:
    typedef ThreeViewTransverseTracksAlgorithm::MatchingType::TensorType TensorType;
    typedef std::vector<TensorType::ElementList::const_iterator> IteratorList;
 
    virtual bool Run(ThreeViewTransverseTracksAlgorithm *const pAlgorithm, TensorType &overlapTensor) = 0;
};
 
} // namespace lar_content
 
#endif // #ifndef LAR_THREE_VIEW_TRANSVERSE_TRACKS_ALGORITHM_H