PCASeedFinderAlg.h

Go to the documentation of this file.

#ifndef PCASeedFinderAlg_h
#define PCASeedFinderAlg_h

// Framework Includes
#include "fhiclcpp/ParameterSet.h"

#include "larreco/RecoAlg/Cluster3DAlgs/SeedFinderAlgBase.h"
#include "larreco/RecoAlg/Cluster3DAlgs/PrincipalComponentsAlg.h"

// LArSoft includes
#include "larcore/Geometry/Geometry.h"
#include "lardataobj/RecoBase/Seed.h"

// Algorithm includes
#include "larreco/RecoAlg/Cluster3DAlgs/Cluster3D.h"

// ROOT includes
#include "TCanvas.h"
#include "TFrame.h"
#include "TH2D.h"

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

namespace lar_cluster3d
{

class PCASeedFinderAlg : virtual public SeedFinderAlgBase
{
public:
    PCASeedFinderAlg(fhicl::ParameterSet const &pset);

    virtual ~PCASeedFinderAlg();
    
    virtual void reconfigure(fhicl::ParameterSet const &pset);

    virtual bool findTrackSeeds(reco::HitPairListPtr&      hitPairListPtr,
                                reco::PrincipalComponents& inputPCA,
                                SeedHitPairListPairVec&    seedHitMap) const;

private:
    
    bool getHitsAtEnd(reco::HitPairListPtr& hit3DList, reco::PrincipalComponents& seedPca) const;
    
    void LineFit2DHits(const reco::HitPairListPtr& hitList, double XOrigin, TVector3& Pos, TVector3& Dir, double& ChiDOF) const;

    geo::Geometry*                         m_geometry;         // pointer to the Geometry service
    //    const detinfo::DetectorProperties*    m_detector;         // Pointer to the detector properties

    double                                 m_gapDistance;      
    size_t                                 m_numSeed2DHits;    
    double                                 m_minAllowedCosAng; 
    
    PrincipalComponentsAlg                 m_pcaAlg;           // For running Principal Components Analysis
};

} // namespace lar_cluster3d
#endif