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