OneViewDeltaRayMatchingAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_ONE_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H
#define LAR_ONE_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

#include "larpandoracontent/LArThreeDReco/LArCosmicRay/DeltaRayMatchingContainers.h"

#include "larpandoracontent/LArUtility/KDTreeLinkerAlgoT.h"

namespace lar_content
{

class OneViewDeltaRayMatchingAlgorithm : public pandora::Algorithm
{
public:
    OneViewDeltaRayMatchingAlgorithm();

private:
    pandora::StatusCode Run();
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    const pandora::ClusterList GetInputClusterList(const pandora::HitType hitType);

    const pandora::PfoList GetMuonPfoList();

    const pandora::PfoList GetDeltaRayPfoList();

    void PerformOneViewMatching(const pandora::HitType hitType);

    bool IsMuonPfo(const pandora::Cluster *const pCluster);

    bool AddIntoExistingDeltaRay(const pandora::Cluster *const pAvailableCluster, const pandora::PfoVector &nearbyMuonPfoVector);

    const pandora::Cluster *GetBestProjectedCluster(const pandora::ClusterList &deltaRayClusterGroup,
        const pandora::ParticleFlowObject *const pNearbyMuonPfo, const pandora::HitType hitType, const bool findAvailable);

    bool IsDeltaRayPfo(const pandora::Cluster *const pCluster);

    void GetClusterSpanX(const pandora::ClusterList &clusterList, float &spanMinX, float &spanMaxX);

    void CreateDeltaRay(const pandora::Cluster *const pAvailableCluster, const pandora::PfoVector &nearbyMuonPfoVector, pandora::ClusterSet &modifiedClusters);

    void GetNearbyAvailableClusters(const pandora::Cluster *const pCluster, pandora::ClusterList &consideredClusters, pandora::ClusterList &foundClusters);

    const pandora::Cluster *MergeClusterGroup(const pandora::ClusterList &clusterGroup);

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

    void PerformRecovery(const pandora::HitType hitType);

    std::string m_muonPfoListName;                           
    std::string m_deltaRayPfoListName;                       
    std::string m_inputClusterListNameU;                     
    std::string m_inputClusterListNameV;                     
    std::string m_inputClusterListNameW;                     
    std::string m_outputPfoListName;                         
    DeltaRayMatchingContainers m_deltaRayMatchingContainers; 
    float m_overlapExtension;      
    unsigned int m_minClusterHits; 
};

} // namespace lar_content

#endif // #ifndef LAR_ONE_VIEW_DELTA_RAY_MATCHING_ALGORITHM_H