TwoViewDeltaRayMatchingAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_TWO_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H
#define LAR_TWO_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H 1

#include "Pandora/Algorithm.h"
#include "Pandora/AlgorithmTool.h"

#include "larpandoracontent/LArObjects/LArTrackTwoViewOverlapResult.h"

#include "larpandoracontent/LArThreeDReco/LArCosmicRay/NViewDeltaRayMatchingAlgorithm.h"
#include "larpandoracontent/LArThreeDReco/LArThreeDBase/TwoViewMatchingControl.h"

#include "larpandoracontent/LArUtility/KDTreeLinkerAlgoT.h"

namespace lar_content
{

class DeltaRayMatrixTool;

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

class TwoViewDeltaRayMatchingAlgorithm : public NViewDeltaRayMatchingAlgorithm<TwoViewMatchingControl<TwoViewDeltaRayOverlapResult>>
{
public:
    typedef NViewDeltaRayMatchingAlgorithm<TwoViewMatchingControl<TwoViewDeltaRayOverlapResult>> BaseAlgorithm;
    typedef TwoViewDeltaRayMatchingAlgorithm::MatchingType::MatrixType MatrixType;

    TwoViewDeltaRayMatchingAlgorithm();

    const std::string &GetThirdViewClusterListName() const;

    const std::string &GetClusteringAlgName() const;

    HitTypeVector GetHitTypeVector();

    const pandora::Cluster *GetCluster(const MatrixType::Element &element, const pandora::HitType hitType);

    bool CreatePfo(const MatrixType::Element &element);

    void UpdateForThirdViewClusterModification(const pandora::Cluster *const pModifiedCluster, const bool isMuon);

private:
    void CalculateOverlapResult(const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, const pandora::Cluster *const pCluster3);

    virtual bool DoesClusterPassTensorThreshold(const pandora::Cluster *const pCluster) const;

    pandora::StatusCode CalculateOverlapResult(
        const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, TwoViewDeltaRayOverlapResult &overlapResult) const;

    void FindCommonMuonParents(const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2, pandora::PfoList &commonMuonPfoList) const;

    void CollectThirdViewClusters(const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2,
        const pandora::CartesianPointVector &projectedPositions, pandora::ClusterList &matchedClusters) const;

    const pandora::Cluster *GetBestMatchedCluster(const pandora::Cluster *const pCluster1, const pandora::Cluster *const pCluster2,
        const pandora::PfoList &commonMuonPfoList, const pandora::ClusterList &matchedClusters, float &reducedChiSquared) const;

    void FormThirdViewCluster(const MatrixType::Element &element, ProtoParticle &protoParticle);

    void MergeThirdView(const MatrixType::Element &element, const pandora::Cluster *const pSeedCluster);

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

    typedef std::vector<DeltaRayMatrixTool *> MatrixToolVector;

    std::string m_inputClusterListName;      
    std::string m_reclusteringAlgorithmName; 
    MatrixToolVector m_algorithmToolVector;  
    unsigned int m_nMaxMatrixToolRepeats;    
    unsigned int m_minClusterCaloHits;       
    float m_maxDistanceFromPrediction;       
    float m_maxGoodMatchReducedChiSquared;   
    float m_minDistanceFromMuon;             
    float m_maxDistanceToCollected;          
};

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

inline const std::string &TwoViewDeltaRayMatchingAlgorithm::GetThirdViewClusterListName() const
{
    return m_inputClusterListName;
}

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

inline const std::string &TwoViewDeltaRayMatchingAlgorithm::GetClusteringAlgName() const
{
    return m_reclusteringAlgorithmName;
}

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

class DeltaRayMatrixTool : public pandora::AlgorithmTool
{
public:
    typedef TwoViewDeltaRayMatchingAlgorithm::MatchingType::MatrixType MatrixType;
    typedef std::vector<MatrixType::ElementList::const_iterator> IteratorList;

    virtual bool Run(TwoViewDeltaRayMatchingAlgorithm *const pAlgorithm, MatrixType &matrixTensor) = 0;

    TwoViewDeltaRayMatchingAlgorithm *m_pParentAlgorithm; 
};

} // namespace lar_content

#endif // #ifndef LAR_TWO_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H