#ifndef TDensityCluster_seen
#define TDensityCluster_seen

#include <iostream>
#include <vector>
#include <functional>

#include "IHitSelection.hxx"

namespace COMET {
    class IDensityCluster;
    typedef COMET::IHitSelection Cluster;
};

class COMET::IDensityCluster {
public:

    IDensityCluster(unsigned int minPts, double epsilon);
    virtual ~IDensityCluster() { }
    
    void SetMinClusterSize(unsigned int size) { fMinClusterSize = size;}
    void Cluster(const COMET::IHitSelection&);
    unsigned int GetClusterCount() { return fClusters.size(); }
    const COMET::IHitSelection& GetCluster(unsigned int i) { return fClusters.at(i); }

    unsigned int ls();

protected:
    COMET::IHitSelection FindSeeds(const COMET::IHitSelection&);
    COMET::IHitSelection GetNeighbors(COMET::IHandle<COMET::IHit>, const COMET::IHitSelection&);
    double CalculateDistance(COMET::IHandle<COMET::IHit>, COMET::IHandle<COMET::IHit>);
    
private:
    unsigned int fMinClusterSize;
    unsigned int fMinPoints;
    double fMinDistance;
    std::vector<COMET::IHitSelection> fClusters;
    COMET::IHitSelection fRemainingHits;


};
#endif