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