HitWidthClusterMergingAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_HIT_WIDTH_CLUSTER_MERGING_ALGORITHM_H
#define LAR_HIT_WIDTH_CLUSTER_MERGING_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

#include "larpandoracontent/LArTwoDReco/LArClusterAssociation/ClusterAssociationAlgorithm.h"

#include "larpandoracontent/LArHelpers/LArHitWidthHelper.h"

namespace lar_content
{

class HitWidthClusterMergingAlgorithm : public ClusterAssociationAlgorithm
{
public:
    HitWidthClusterMergingAlgorithm();

private:
    void GetListOfCleanClusters(const pandora::ClusterList *const pClusterList, pandora::ClusterVector &clusterVector) const;
    void PopulateClusterAssociationMap(const pandora::ClusterVector &clusterVector, ClusterAssociationMap &clusterAssociationMap) const;
    bool IsExtremalCluster(const bool isForward, const pandora::Cluster *const pCurrentCluster, const pandora::Cluster *const pTestCluster) const;

    bool AreClustersAssociated(const LArHitWidthHelper::ClusterParameters &currentClusterParameters,
        const LArHitWidthHelper::ClusterParameters &testClusterParameters) const;

    void FindClosestPointToPosition(const pandora::CartesianVector &position,
        const LArHitWidthHelper::ConstituentHitVector &constituentHitVector, pandora::CartesianVector &closestPoint) const;

    void GetClusterDirection(const LArHitWidthHelper::ConstituentHitVector &constituentHitVector, pandora::CartesianVector &direction,
        const pandora::CartesianVector &fitReferencePoint, const float fittingWeight) const;

    void GetConstituentHitSubsetVector(const LArHitWidthHelper::ConstituentHitVector &constituentHitVector, const pandora::CartesianVector &fitReferencePoint,
        const float fittingWeight, LArHitWidthHelper::ConstituentHitVector &constituentHitSubsetVector) const;

    void GetFittingAxes(const LArHitWidthHelper::ConstituentHitVector &constituentHitSubsetVector, pandora::CartesianVector &axisDirection,
        pandora::CartesianVector &orthoDirection) const;

    void GetFittingCoordinates(
        const pandora::CartesianVector &axisDirection, const pandora::CartesianVector &constituentHitPosition, float &rL, float &rT) const;

    void GetGlobalDirection(const pandora::CartesianVector &axisDirection, const float gradient, pandora::CartesianVector &globalDirection) const;

    void RemoveShortcutAssociations(const pandora::ClusterVector &clusterVector, ClusterAssociationMap &clusterAssociationMap) const;

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

    float m_maxConstituentHitWidth;  
    float m_hitWidthScalingFactor;   
    float m_fittingWeight;           
    float m_minClusterWeight;        
    float m_maxXMergeDistance;       
    float m_maxZMergeDistance;       
    float m_minMergeCosOpeningAngle; 
    float m_minDirectionDeviationCosAngle; 
    float m_minClusterSparseness;          

    // ATTN Dangling pointers emerge during cluster merging, here explicitly not dereferenced
    mutable LArHitWidthHelper::ClusterToParametersMap m_clusterToParametersMap; 
};

} //namespace lar_content

#endif //LAR_HIT_WIDTH_CLUSTER_MERGING_ALGORITHM_H