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