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