CosmicRayExtensionAlgorithm.h

Go to the documentation of this file.

/**
 *  @file   larpandoracontent/LArTwoDReco/LArCosmicRay/CosmicRayExtensionAlgorithm.h
 *
 *  @brief  Header file for the cosmic-ray extension algorithm class.
 *
 *  $Log: $
 */
#ifndef LAR_COSMIC_RAY_EXTENSION_ALGORITHM_H
#define LAR_COSMIC_RAY_EXTENSION_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

#include "larpandoracontent/LArObjects/LArPointingCluster.h"

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

namespace lar_content
{

/**
 *  @brief  CosmicRayExtensionAlgorithm class
 */
class CosmicRayExtensionAlgorithm : public ClusterExtensionAlgorithm
{
public:
    /**
     *  @brief  Default constructor
     */
    CosmicRayExtensionAlgorithm();

private:
    void GetListOfCleanClusters(const pandora::ClusterList *const pClusterList, pandora::ClusterVector &clusterVector) const;
    void FillClusterAssociationMatrix(const pandora::ClusterVector &clusterVector, ClusterAssociationMatrix &clusterAssociationMatrix) const;
    void FillClusterMergeMap(const ClusterAssociationMatrix &clusterAssociationMatrix, ClusterMergeMap &clusterMergeMap) const;

    /**
     *  @brief  Form association between two pointing clusters
     *
     *  @param  clusterI the first pointing cluster
     *  @param  clusterJ the second pointing cluster
     *  @param  clusterAssociationMatrix the matrix of cluster associations
     */
    void FillClusterAssociationMatrix(
        const LArPointingCluster &clusterI, const LArPointingCluster &clusterJ, ClusterAssociationMatrix &clusterAssociationMatrix) const;

    /**
     *  @brief  Calculate RMS deviation of a cluster with respect to the reference line
     *
     *  @param  pCluster the input cluster
     *  @param  position the intercept of the reference line
     *  @param  direction the direction of the reference line
     */
    float CalculateRms(const pandora::Cluster *const pCluster, const pandora::CartesianVector &position, const pandora::CartesianVector &direction) const;

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

    float m_minClusterLength;            ///<
    float m_minSeedClusterLength;        ///<
    float m_maxLongitudinalDisplacement; ///<
    float m_maxTransverseDisplacement;   ///<
    float m_minCosRelativeAngle;         ///<
    float m_maxAverageRms;               ///<
};

} // namespace lar_content

#endif // #ifndef LAR_COSMIC_RAY_EXTENSION_ALGORITHM_H