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