ClusteringComparisonCP.h

#ifndef ClusteringComparisonCP_h
#define ClusteringComparisonCP_h

#include "ClusteringComparisonIndexCP.h"

namespace stab {

class TransversePairwiseIndexCP: public ClusteringComparisonIndexCP {

public:
  inline virtual ~TransversePairwiseIndexCP() { }

  double index( const PseudoJet * inputParticles,
        const unsigned * clustering1,
        const unsigned * clustering2,
        unsigned size,
        double * constitStab ) const;

private:
  virtual double combineCounts(double c1, double c2,
                               double c3, double c4) const = 0;

}; // end TransversepairwiseIndexCP class 



class TransverseRandIndexCP : public TransversePairwiseIndexCP {

  inline double combineCounts(const double c1, const double c2,
                  const double c3, const double c4) const
  {
    return (c1 + c2)/(c1 + c2 + c3 + c4);
  }

};  // end TransRandIndexCP class


class TransverseJaccardIndexCP : public TransversePairwiseIndexCP {

  inline double combineCounts(const double c1, const double c2,
                  const double c3, const double c4) const
  {
    return c1/(c1 + c3 + c4);
  }

};  // end TransJaccardIndexCP class

} // end stab namespace

#endif