PeakDirectionFinderTool.h

Go to the documentation of this file.

#ifndef LAR_PEAK_DIRECTION_FINDER_TOOL_H
#define LAR_PEAK_DIRECTION_FINDER_TOOL_H 1

#include "Pandora/AlgorithmHeaders.h"
#include "Pandora/AlgorithmTool.h"

namespace lar_content
{

class PeakDirectionFinderTool : public pandora::AlgorithmTool
{
public:
    PeakDirectionFinderTool();

    pandora::StatusCode Run(const pandora::ParticleFlowObject *const pShowerPfo, const pandora::CartesianVector &nuVertex3D,
        const pandora::CaloHitList *const pViewHitList, const pandora::HitType hitType, pandora::CartesianPointVector &peakDirectionVector);

private:
    typedef std::map<int, float> AngularDecompositionMap;

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

    void CollectHitsWithinROI(const pandora::CaloHitList &showerHitList, const pandora::CaloHitList *const pViewHitList,
        const pandora::CartesianVector &nuVertex2D, pandora::CaloHitList &viewROIHits) const;

    void GetAngularExtrema(const pandora::CaloHitList &showerHitList, const pandora::CartesianVector &nuVertex2D, float &lowestTheta,
        float &highestTheta) const;

    void CollectHitsWithinExtrema(const pandora::CaloHitList *const pViewHitList, const pandora::CartesianVector &nuVertex2D,
        const float lowestTheta, const float highestTheta, pandora::CaloHitList &viewROIHits) const;

    void FillAngularDecompositionMap(const pandora::CaloHitList &viewShowerHitList, const pandora::CartesianVector &nuVertex2D,
        AngularDecompositionMap &angularDecompositionMap) const;

    void SmoothAngularDecompositionMap(AngularDecompositionMap &angularDecompositionMap) const;

    void RetrievePeakDirections(const AngularDecompositionMap &angularDecompositionMap, pandora::CartesianPointVector &peakDirectionVector) const;

    float m_pathwaySearchRegion;  
    float m_theta0XZBinSize;      
    int m_smoothingWindow;        
    bool m_ambiguousParticleMode; 
};

//------------------------------------------------------------------------------------------------------------------------------------------
} // namespace lar_content

#endif // #ifndef LAR_PEAK_DIRECTION_FINDER_TOOL_H