ThreeDRemnantsAlgorithm.h

Go to the documentation of this file.

#ifndef LAR_THREE_D_REMNANTS_ALGORITHM_H
#define LAR_THREE_D_REMNANTS_ALGORITHM_H 1

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

#include "larpandoracontent/LArObjects/LArOverlapTensor.h"

#include "larpandoracontent/LArThreeDReco/LArThreeDBase/ThreeDBaseAlgorithm.h"

namespace lar_content
{

class RemnantTensorTool;

//------------------------------------------------------------------------------------------------------------------------------------------

class ThreeDRemnantsAlgorithm : public ThreeDBaseAlgorithm<float>
{
public:
    ThreeDRemnantsAlgorithm();

    void SelectInputClusters(const pandora::ClusterList *const pInputClusterList, pandora::ClusterList &selectedClusterList) const;
    void SetPfoParameters(const ProtoParticle &protoParticle, PandoraContentApi::ParticleFlowObject::Parameters &pfoParameters) const;

private:
    void CalculateOverlapResult(const pandora::Cluster *const pClusterU, const pandora::Cluster *const pClusterV, const pandora::Cluster *const pClusterW);
    void ExamineTensor();

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

    typedef std::vector<RemnantTensorTool*> RemnantTensorToolVector;
    RemnantTensorToolVector     m_algorithmToolVector;      

    unsigned int                m_nMaxTensorToolRepeats;    
    unsigned int                m_minClusterCaloHits;       
    float                       m_xOverlapWindow;           
    float                       m_pseudoChi2Cut;            
};

//------------------------------------------------------------------------------------------------------------------------------------------

class RemnantTensorTool : public pandora::AlgorithmTool
{
public:
    typedef ThreeDRemnantsAlgorithm::TensorType TensorType;
    typedef std::vector<TensorType::ElementList::const_iterator> IteratorList;

    virtual bool Run(ThreeDRemnantsAlgorithm *const pAlgorithm, TensorType &overlapTensor) = 0;
};

} // namespace lar_content

#endif // #ifndef LAR_THREE_D_REMNANTS_ALGORITHM_H