CosmicRaySplittingAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_COSMIC_RAY_SPLITTING_ALGORITHM_H
#define LAR_COSMIC_RAY_SPLITTING_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

#include "larpandoracontent/LArObjects/LArTwoDSlidingFitResult.h"

namespace lar_content
{

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

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

    void GetListOfCleanClusters(const pandora::ClusterList *const pClusterList, pandora::ClusterVector &clusterVector) const;

    void BuildSlidingFitResultMap(const pandora::ClusterVector &clusterVector, TwoDSlidingFitResultMap &slidingFitResultMap) const;

    pandora::StatusCode FindBestSplitPosition(const TwoDSlidingFitResult &slidingFitResult, pandora::CartesianVector &splitPosition,
        pandora::CartesianVector &splitDirection1, pandora::CartesianVector &splitDirection2) const;

    pandora::StatusCode ConfirmSplitPosition(const TwoDSlidingFitResult &branchSlidingFitResult,
        const TwoDSlidingFitResult &replacementSlidingFitResult, const pandora::CartesianVector &splitPosition,
        const pandora::CartesianVector &splitDirection1, const pandora::CartesianVector &splitDirection2,
        float &lengthSquared1, float &lengthSquared2) const;

    pandora::StatusCode PerformSingleSplit(const pandora::Cluster *const pBranchCluster, const pandora::Cluster *const pReplacementCluster,
        const pandora::CartesianVector &splitPosition, const  pandora::CartesianVector &forwardDirection,
        const pandora::CartesianVector &backwardDirection) const;

    pandora::StatusCode PerformDoubleSplit(const pandora::Cluster *const pBranchCluster, const pandora::Cluster *const pReplacementCluster1,
        const pandora::Cluster *const pReplacementCluster2, const pandora::CartesianVector &splitPosition,
        const pandora::CartesianVector &splitDirection1, const pandora::CartesianVector &splitDirection2) const;

    void GetCaloHitListToMove(const pandora::Cluster *const pBranchCluster, const pandora::Cluster *const pReplacementCluster,
        const pandora::CartesianVector &splitPosition, const pandora::CartesianVector &forwardDirection,
        const pandora::CartesianVector &backwardDirection, pandora::CaloHitList &caloHitList) const;

    void GetCaloHitListsToMove(const pandora::Cluster *const pBranchCluster, const pandora::CartesianVector &splitPosition,
        const pandora::CartesianVector &splitDirection1, const pandora::CartesianVector &splitDirection2,
        pandora::CaloHitList &caloHitList1, pandora::CaloHitList &caloHitList2) const;

    bool IdentifyCrossedTracks(const pandora::Cluster *const pBranchCluster, const pandora::Cluster *const pReplacementCluster1,
        const pandora::Cluster *const pReplacementCluster2, const pandora::CartesianVector &splitPosition) const;

    pandora::StatusCode GetCaloHitListToKeep(const pandora::Cluster *const pBranchCluster, const pandora::CaloHitList &caloHitListToMove,
        pandora::CaloHitList &caloHitListToKeep) const;

    pandora::StatusCode SplitCluster(const pandora::Cluster *const pBranchCluster, const pandora::Cluster *const pReplacementCluster,
        const pandora::CaloHitList &caloHitListToMove) const;

    float          m_clusterMinLength;                     
    unsigned int   m_halfWindowLayers;                     
    float          m_samplingPitch;                        
    float          m_maxCosSplittingAngle;                 
    float          m_minCosMergingAngle;                   
    float          m_maxTransverseDisplacement;            
    float          m_maxLongitudinalDisplacement;          
    float          m_maxLongitudinalDisplacementSquared;   
};

} // namespace lar_content

#endif // #ifndef LAR_COSMIC_RAY_SPLITTING_ALGORITHM_H