DeltaRayParentAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_DELTA_RAY_PARENT_ALGORITHM_H
#define LAR_DELTA_RAY_PARENT_ALGORITHM_H 1

#include "Pandora/Algorithm.h"

namespace lar_content
{

class DeltaRayParentAlgorithm : public pandora::Algorithm
{
public:
    typedef std::map<const pandora::ParticleFlowObject *, float> PfoLengthMap;

    DeltaRayParentAlgorithm();

    pandora::StatusCode Run();

private:
    void InitialisePfoLengthMap(const pandora::PfoList *const muonPfoList, const pandora::PfoList *const deltaRayPfoList, PfoLengthMap &pfoLengthMap) const;

    void FindParentPfo(const PfoLengthMap &pfoLengthMap, const pandora::ParticleFlowObject *const pPfo, const pandora::ParticleFlowObject *&pParentPfo) const;

    pandora::StatusCode GetTwoDSeparation(
        const pandora::ParticleFlowObject *const pPfo1, const pandora::ParticleFlowObject *const pPfo2, float &separation) const;

    void AssignToParentPfo(const pandora::PfoList *const muonPfoList, const pandora::PfoList *const deltaRayPfoList,
        const pandora::ParticleFlowObject *const pPfo, const pandora::ParticleFlowObject *const pParentPfo, PfoLengthMap &pfoLengthMap) const;

    void UpdatePfoLengthMap(const pandora::PfoList &pfosToRemove, const pandora::ParticleFlowObject *const pPfoToAdd, PfoLengthMap &pfoLengthMap) const;

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

    std::string m_muonPfoListName;     
    std::string m_deltaRayPfoListName; 
    float m_distanceForMatching;       
};

} // namespace lar_content

#endif // #ifndef LAR_DELTA_RAY_PARENT_ALGORITHM_H