TwoViewCosmicRayRemovalTool.h

Go to the documentation of this file.

#ifndef TWO_VIEW_COSMIC_RAY_REMOVAL_TOOL_H
#define TWO_VIEW_COSMIC_RAY_REMOVAL_TOOL_H 1

#include "larpandoracontent/LArThreeDReco/LArCosmicRay/TwoViewDeltaRayMatchingAlgorithm.h"

namespace lar_content
{
class TwoViewCosmicRayRemovalTool : public DeltaRayMatrixTool
{
public:
    typedef std::vector<pandora::HitType> HitTypeVector;
    TwoViewCosmicRayRemovalTool();

private:
    bool Run(TwoViewDeltaRayMatchingAlgorithm *const pAlgorithm, MatrixType &overlapMatrix);
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);

    bool RemoveCosmicRayHits(const MatrixType::ElementList &elementList) const;

    bool PassElementChecks(const MatrixType::Element &element, const pandora::HitType hitType) const;

    bool IsMuonEndpoint(const MatrixType::Element &element, const pandora::HitType hitType) const;

    bool IsContaminated(const MatrixType::Element &element, const pandora::HitType hitType) const;

    bool IsCloseToLine(const pandora::CartesianVector &hitPosition, const pandora::CartesianVector &lineStart,
        const pandora::CartesianVector &lineEnd, const float distanceToLine) const;

    bool IsInLineSegment(const pandora::CartesianVector &lowerBoundary, const pandora::CartesianVector &upperBoundary,
        const pandora::CartesianVector &point) const;

    bool IsBestElement(const MatrixType::Element &element, const pandora::HitType hitType, const MatrixType::ElementList &elementList) const;

    void CreateSeed(const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits) const;

    pandora::StatusCode GrowSeed(const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits,
        pandora::CaloHitList &deltaRayRemantHits) const;

    void CollectHitsFromDeltaRay(const pandora::CartesianVector &positionOnMuon, const pandora::CartesianVector &muonDirection,
        const pandora::Cluster *const pDeltaRayCluster, const float &minDistanceFromMuon, const bool demandCloserToCollected,
        const pandora::CaloHitList &protectedHits, pandora::CaloHitList &collectedHits) const;

    void SplitDeltaRayCluster(const MatrixType::Element &element, const pandora::HitType hitType, pandora::CaloHitList &collectedHits,
        pandora::CaloHitList &deltaRayRemnantHits) const;

    void ReclusterRemnant(const pandora::HitType hitType, const pandora::Cluster *const pMuonCluster,
        const pandora::Cluster *const pDeltaRayRemnant, pandora::ClusterVector &clusterVector, pandora::PfoVector &pfoVector) const;

    pandora::StatusCode ProjectDeltaRayPositions(
        const MatrixType::Element &element, const pandora::HitType hitType, pandora::CartesianPointVector &projectedPositions) const;

    float m_minSeparation;           
    unsigned int m_slidingFitWindow; 
    float m_distanceToLine;          
    float m_minContaminationLength; 
    float m_maxDistanceToHit;             
    unsigned int m_minRemnantClusterSize; 
    float m_maxDistanceToTrack; 
};

} // namespace lar_content

#endif // #ifndef TWO_VIEW_COSMIC_RAY_REMOVAL_TOOL_H