SplitShowersTool.h

Go to the documentation of this file.

#ifndef SPLIT_SHOWERS_TOOL_H
#define SPLIT_SHOWERS_TOOL_H 1

#include "larpandoracontent/LArThreeDReco/LArShowerMatching/ThreeDShowersAlgorithm.h"

namespace lar_content
{

class SplitShowersTool : public ShowerTensorTool
{
public:
    SplitShowersTool();

    bool Run(ThreeDShowersAlgorithm *const pAlgorithm, TensorType &overlapTensor);

private:
    void FindSplitShowers(ThreeDShowersAlgorithm *const pAlgorithm, const TensorType &overlapTensor, ClusterMergeMap &clusterMergeMap) const;

    bool PassesElementCuts(TensorType::ElementList::const_iterator eIter, const pandora::ClusterSet &usedClusters) const;

    void SelectTensorElements(TensorType::ElementList::const_iterator eIter, const TensorType::ElementList &elementList,
        const pandora::ClusterSet &usedClusters, IteratorList &iteratorList) const;

    void FindShowerMerges(ThreeDShowersAlgorithm *const pAlgorithm, const IteratorList &iteratorList, pandora::ClusterSet &usedClusters,
        ClusterMergeMap &clusterMergeMap) const;

    bool CheckClusterProximities(ThreeDShowersAlgorithm *const pAlgorithm, const pandora::ClusterList &clusterList) const;

    bool CheckClusterVertexRelations(ThreeDShowersAlgorithm *const pAlgorithm, const pandora::ClusterList &clusterList) const;

    bool CheckClusterSplitPositions(ThreeDShowersAlgorithm *const pAlgorithm, const pandora::ClusterList &clusterListU,
        const pandora::ClusterList &clusterListV, const pandora::ClusterList &clusterListW) const;

    void GetSplitXDetails(ThreeDShowersAlgorithm *const pAlgorithm, const pandora::Cluster *const pClusterA, const pandora::Cluster *const pClusterB,
        float &splitXPosition, float &overlapX) const;

    void SpecifyClusterMerges(ThreeDShowersAlgorithm *const pAlgorithm, const pandora::ClusterList &clusterList, ClusterMergeMap &clusterMergeMap) const;

    bool ApplyChanges(ThreeDShowersAlgorithm *const pAlgorithm, const ClusterMergeMap &clusterMergeMap) const;

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

    unsigned int    m_nCommonClusters;                  
    float           m_minMatchedFraction;               
    unsigned int    m_minMatchedSamplingPoints;         

    bool            m_checkClusterProximities;          
    float           m_maxClusterSeparation;             

    bool            m_checkClusterVertexRelations;      
    float           m_minVertexLongitudinalDistance;    
    float           m_maxVertexLongitudinalDistance;    
    float           m_maxVertexTransverseDistance;      
    float           m_vertexAngularAllowance;           
    unsigned int    m_maxVertexAssociations;            

    bool            m_checkClusterSplitPositions;       
    float           m_vetoMergeXDifference;             
    float           m_vetoMergeXOverlap;                
};

} // namespace lar_content

#endif // #ifndef SPLIT_SHOWERS_TOOL_H