DeltaRayIdentificationAlgorithm.h

Go to the documentation of this file.

 
#ifndef LAR_DELTA_RAY_IDENTIFICATION_ALGORITHM_H
#define LAR_DELTA_RAY_IDENTIFICATION_ALGORITHM_H 1
 
#include "Pandora/Algorithm.h"
 
#include <unordered_map>
 
namespace lar_content
{
 
class DeltaRayIdentificationAlgorithm : public pandora::Algorithm
{
public:
    DeltaRayIdentificationAlgorithm();
 
private:
    pandora::StatusCode Run();
 
    typedef std::unordered_map<const pandora::ParticleFlowObject *, const pandora::ParticleFlowObject *> PfoAssociationMap;
 
    void GetPfos(const std::string &inputPfoListName, pandora::PfoVector &outputPfoVector) const;
 
    void BuildAssociationMap(const pandora::PfoVector &inputPfos, const pandora::PfoVector &outputPfos, PfoAssociationMap &pfoAssociationMap) const;
 
    bool IsAssociated(const pandora::ParticleFlowObject *const pDaughterPfo, const pandora::ParticleFlowObject *const pParentPfo, float &displacement) const;
 
    float GetTwoDSeparation(const pandora::ParticleFlowObject *const pDaughterPfo, const pandora::ParticleFlowObject *const pParentPfo) const;
 
    void GetTwoDVertices(const pandora::ParticleFlowObject *const pPfo, const pandora::HitType &hitType, pandora::CartesianPointVector &vertexVector) const;
 
    float GetClosestDistance(const pandora::CartesianPointVector &vertexVector, const pandora::ClusterList &clusterList) const;
 
    void BuildParentDaughterLinks(const PfoAssociationMap &pfoAssociationMap, pandora::PfoList &outputPfoList) const;
 
    const pandora::ParticleFlowObject *GetParent(const PfoAssociationMap &pfoAssociationMap, const pandora::ParticleFlowObject *const pPfo) const;
 
    pandora::StatusCode ReadSettings(const pandora::TiXmlHandle xmlHandle);
 
    std::string m_parentPfoListName;   
    std::string m_daughterPfoListName; 
 
    float m_distanceForMatching;      
    float m_minParentLengthSquared;   
    float m_maxDaughterLengthSquared; 
};
 
} // namespace lar_content
 
#endif // #ifndef LAR_DELTA_RAY_IDENTIFICATION_ALGORITHM_H