HyperPointSet Class Reference

#include <HyperPointSet.h>

Public Member Functions
const int &	getDimension () const
bool	compatible (const HyperPoint &other, bool printError=true) const
const HyperPoint &	at (int i) const
HyperPoint &	at (int i)
unsigned int	size () const
void	push_back (const HyperPoint &point)
void	addHyperPointSet (const HyperPointSet &other)
void	removeDuplicates ()
void	sort ()
double	getSumW () const
double	getSumW2 () const
double	getCorrelation (int i, int j) const
TMatrixD	getCorrelationMatrix () const
double	getCovarience (int i, int j) const
TMatrixD	getCovarienceMatrix () const
HyperPoint	getMin () const
	Get a HyperCuboid that surrounds the points. More...
HyperPoint	getMax () const
	Get a HyperPoint that gives the maximum in each dim. More...
HyperPoint	mean () const
HyperPoint	geometricMean () const
HyperPoint	harmonicMean () const
void	save (TString path)
void	load (TString path)
void	print (std::ostream &os=std::cout) const
	HyperPointSet (int dimension)
	HyperPointSet (const HyperPoint &point)
	HyperPointSet (const HyperPoint &point1, const HyperPoint &point2)
	HyperPointSet (const HyperPoint &point1, const HyperPoint &point2, const HyperPoint &point3)
	HyperPointSet (int npoints, const HyperPoint &point)
	HyperPointSet (TString path)
virtual	~HyperPointSet ()
bool	linearlyIndependant () const
HyperPointSet	gramSchmidtProcess () const
void	normalise ()

Private Member Functions
void	save ()

Private Attributes
int	_dimension
std::vector< HyperPoint >	_points

Detailed Description

HyperPlot, Author: Sam Harnew, sam.h.nosp@m.arne.nosp@m.w@gma.nosp@m.il.c.nosp@m.om , Date: Dec 2015

A vector of HyperPoints

Definition at line 31 of file HyperPointSet.h.

Constructor & Destructor Documentation

HyperPointSet() [1/6]

HyperPointSet::HyperPointSet

(

int

dimension

)

Standard constuctor where only the dimensionality is given

Definition at line 6 of file HyperPointSet.cpp.

                                          : 
  _dimension(dimension)
{
  WELCOME_LOG << "Hello from the HyperPointSet() Constructor";
}

HyperPointSet() [2/6]

HyperPointSet::HyperPointSet

(

const HyperPoint &

point

)

Constuctor for a HyperPointSet that takes dimensionality from a given HyperPoint and adds that point to the HyperPointSet

Definition at line 35 of file HyperPointSet.cpp.

                                                    :
  _dimension(point.getDimension())
{
  push_back(point);
}

HyperPointSet() [3/6]

HyperPointSet::HyperPointSet	(	const HyperPoint &	point1,
		const HyperPoint &	point2
	)

Constuctor for a HyperPointSet that takes dimensionality from a given HyperPoint and adds that point and another to the HyperPointSet

Definition at line 43 of file HyperPointSet.cpp.

                                                                               :
  _dimension(point1.getDimension())
{
  push_back(point1);
  push_back(point2);
}

HyperPointSet() [4/6]

HyperPointSet::HyperPointSet	(	const HyperPoint &	point1,
		const HyperPoint &	point2,
		const HyperPoint &	point3
	)

Constuctor for a HyperPointSet that takes dimensionality from a given HyperPoint and adds that point and another 2 to the HyperPointSet

Definition at line 52 of file HyperPointSet.cpp.

                                                                                                         :
  _dimension(point1.getDimension())
{
  push_back(point1);
  push_back(point2);
  push_back(point3);
}

HyperPointSet() [5/6]

HyperPointSet::HyperPointSet	(	int	npoints,
		const HyperPoint &	point
	)

Constuctor for a HyperPointSet that repeats the same point a total of npoints times.

Definition at line 24 of file HyperPointSet.cpp.

                                                                 :
  _dimension(point.getDimension())
{

  for (int i = 0; i < npoints; i++) push_back(point);

}

HyperPointSet() [6/6]

HyperPointSet::HyperPointSet

(

TString

path

)

Constuctor for loading exisiting HyperPointSet from a file

Definition at line 14 of file HyperPointSet.cpp.

                                         :
  _dimension(-1)
{

  load(path);
 
}

~HyperPointSet()

HyperPointSet::~HyperPointSet ( )

virtual

Destructor

Definition at line 506 of file HyperPointSet.cpp.

                             {
  GOODBYE_LOG << "Goodbye from the HyperPointSet() Destructor";
}

Member Function Documentation

addHyperPointSet()

void HyperPointSet::addHyperPointSet

(

const HyperPointSet &

other

)

Add all the points from another HyperPointSet to this one.

Definition at line 202 of file HyperPointSet.cpp.

                                                              {
  
  if (_dimension != other._dimension){
    ERROR_LOG << "Removing all duplicates from the HyperPointSet - did you remember to call sort() first?!?";

  }

  for (unsigned int i = 0; i < other.size(); i++){
    this->push_back(other.at(i));
  }

}

at() [1/2]

const HyperPoint & HyperPointSet::at

(

int

i

)

const

Get a point from the HyperPointSet

Definition at line 72 of file HyperPointSet.cpp.

                                                    { 
  return _points.at(i);
}

at() [2/2]

HyperPoint & HyperPointSet::at

(

int

i

)

Get a point from the HyperPointSet

Definition at line 78 of file HyperPointSet.cpp.

                                                    { 
  return _points.at(i);
}

compatible()

bool HyperPointSet::compatible	(	const HyperPoint &	other,
		bool	printError = true
	)		const

Is this HyperPoint compatible with this HyperPointSet? (same dimension)

Definition at line 62 of file HyperPointSet.cpp.

                                                                            {

  if (_dimension == other.getDimension()) return true;

  if (printError) ERROR_LOG << "This HyperPoint is NOT compatible with this HyperPointSet";
  return false;
}

geometricMean()

HyperPoint HyperPointSet::geometricMean

(

)

const

Get the geometric mean of all points. This does not use the weights.

Definition at line 460 of file HyperPointSet.cpp.

                                             {

  HyperPoint geoMean(getDimension());
  for (unsigned i = 0; i < size(); i++){
    HyperPoint temp(getDimension());
    for (int j= 0; j < getDimension(); j++){
      temp.at(j) = log( this->at(i).at(j) );
    }
    geoMean = geoMean + temp;
  }
  
  geoMean = geoMean/double(size());

  for (int j= 0; j < getDimension(); j++){
    geoMean.at(j) = exp( geoMean.at(j) );
  }
  
  return geoMean;

}

getCorrelation()

double HyperPointSet::getCorrelation	(	int	i,
		int	j
	)		const

Get the correlation between two variables.

Definition at line 217 of file HyperPointSet.cpp.

                                                      {

  WidthFinder statsi;
  WidthFinder statsj;

  for (unsigned evt = 0; evt < size(); evt++){
    statsi.add(at(evt).at(i));
    statsj.add(at(evt).at(j));
  }
  
  double meani  = statsi.mean();
  double meanj  = statsj.mean();
  double widthi = statsi.width();
  double widthj = statsj.width();

  double correlation = 0.0;

  for (unsigned evt = 0; evt < size(); evt++){
    correlation += (at(evt).at(i) - meani)*(at(evt).at(j) - meanj);
  } 
  
  correlation = correlation / (widthi * widthj * size());
  
  return correlation;
}

getCorrelationMatrix()

TMatrixD HyperPointSet::getCorrelationMatrix

(

)

const

Get the correlation matrix associated with the HyperPointSet.

Definition at line 286 of file HyperPointSet.cpp.

                                                  {

  TMatrixD matrix(getDimension(), getDimension());

  for (int i = 0; i < getDimension(); i++){
    for (int j = 0; j < getDimension(); j++){
      matrix(i,j) = getCorrelation(i,j);
    }    
  }
  
  return matrix;
}

getCovarience()

double HyperPointSet::getCovarience	(	int	i,
		int	j
	)		const

Get the covarience between two variables.

Definition at line 245 of file HyperPointSet.cpp.

                                                     {

  MeanFinder statsi;
  MeanFinder statsj;

  for (unsigned evt = 0; evt < size(); evt++){
    statsi.add(at(evt).at(i));
    statsj.add(at(evt).at(j));
  }
  
  double meani  = statsi.mean();
  double meanj  = statsj.mean();

  double correlation = 0.0;

  for (unsigned evt = 0; evt < size(); evt++){
    correlation += (at(evt).at(i) - meani)*(at(evt).at(j) - meanj);
  } 
  
  correlation = correlation / size();
  
  return correlation;
}

getCovarienceMatrix()

TMatrixD HyperPointSet::getCovarienceMatrix

(

)

const

Get the covarience matrix associated with the HyperPointSet.

Definition at line 271 of file HyperPointSet.cpp.

                                                 {

  TMatrixD matrix(getDimension(), getDimension());

  for (int i = 0; i < getDimension(); i++){
    for (int j = 0; j < getDimension(); j++){
      matrix(i,j) = getCovarience(i,j);
    }    
  }
  
  return matrix;
}

getDimension()

const int& HyperPointSet::getDimension ( ) const

inline

Get the dimensionality of HyperPoints in the HyperPointSet

Definition at line 40 of file HyperPointSet.h.

getMax()

HyperPoint HyperPointSet::getMax

(

)

const

Get a HyperPoint that gives the maximum in each dim.

Definition at line 434 of file HyperPointSet.cpp.

                                      {
  HyperPoint val(getDimension());

  for (int i = 0; i < getDimension(); i++){
    MinMaxFinder minMax;
    for (unsigned j = 0; j < size(); j++){
      minMax.add( this->at(j).at(i) );
    }
    val.at(i) = minMax.getMax();
  }
  return val;
}

getMin()

HyperPoint HyperPointSet::getMin

(

)

const

Get a HyperCuboid that surrounds the points.

Get a HyperPoint that gives the minimum in each dim

Definition at line 420 of file HyperPointSet.cpp.

                                      {
  HyperPoint val(getDimension());

  for (int i = 0; i < getDimension(); i++){
    MinMaxFinder minMax;
    for (unsigned j = 0; j < size(); j++){
      minMax.add( this->at(j).at(i) );
    }
    val.at(i) = minMax.getMin();
  }
  return val;
}

getSumW()

double HyperPointSet::getSumW

(

)

const

Find the Sum of weights for all HyperPoints in the HyperPointSet

Definition at line 99 of file HyperPointSet.cpp.

                                   {
  
  double sumW = 0.0;

  for(unsigned i = 0; i < size(); i++){
    sumW += at(i).getWeight();
  }
  return sumW;
}

getSumW2()

double HyperPointSet::getSumW2

(

)

const

Find the Sum of weights squared for all HyperPoints in the HyperPointSet

Definition at line 111 of file HyperPointSet.cpp.

                                    {

  double sumW2 = 0.0;
  double w = 0.0;
  for(unsigned i = 0; i < size(); i++){
    w = at(i).getWeight();
    sumW2 += w*w;
  }
  return sumW2;

}

gramSchmidtProcess()

HyperPointSet HyperPointSet::gramSchmidtProcess

(

)

const

Use the gramSchmidtProcess to make the first nDim vectors orthanormal. If there are less than this, only size() vectors will be returned (but they will be orthanormal).

Definition at line 127 of file HyperPointSet.cpp.

                                                     {

  int nPoints = size();
  if ( nPoints > getDimension() ) nPoints = getDimension();
  
  if (linearlyIndependant() == 0) {
    ERROR_LOG << "Trying to perform the Gram Schmidt Process on a set of L.D. vectors";
  }

  HyperPointSet orthogonalPoints(at(0));

  for (int i = 1; i < nPoints; i++){
    HyperPoint nextPoint = at(i);
    for (int j = 0; j < i; j++){
      nextPoint = nextPoint - orthogonalPoints.at(j).project(at(i));
    }
    orthogonalPoints.push_back(nextPoint);
  }
  
  orthogonalPoints.normalise();
  return orthogonalPoints;

}

harmonicMean()

HyperPoint HyperPointSet::harmonicMean

(

)

const

Get the harmonic mean of all points. This does not use the weights.

Definition at line 483 of file HyperPointSet.cpp.

                                            {

  HyperPoint geoMean(getDimension());
  for (unsigned i = 0; i < size(); i++){
    HyperPoint temp(getDimension());
    for (int j= 0; j < getDimension(); j++){
      temp.at(j) = 1.0 / this->at(i).at(j);
    }
    geoMean = geoMean + temp;
  }
  
  geoMean = geoMean/double(size());

  for (int j= 0; j < getDimension(); j++){
    geoMean.at(j) = 1.0 / geoMean.at(j);
  }
  
  return geoMean;

}

linearlyIndependant()

bool HyperPointSet::linearlyIndependant

(

)

const

Check if the HyperPoints in the HyperPointSet are linearly independent (this isn't implemented yet!!!)

Definition at line 153 of file HyperPointSet.cpp.

                                             {
  //need to implement this...
  return 1;
}

load()

void HyperPointSet::load

(

TString

path

)

Load exisiting HyperPointSet from a file. HyperPoints are appended to the end of the

Definition at line 362 of file HyperPointSet.cpp.

                                    {
  
  TFile* file = new TFile(path, "READ");
  if (file == 0) {
    ERROR_LOG << "Cannot open file at " << path;
    return;
  }

  TTree* tree = (TTree*)file->Get("HyperPointSet");
  if (tree == 0) {
    ERROR_LOG << "Cannot open tree";
    return;
  }
  
  std::vector<double>*  values  = 0;
  std::vector<double>*  weights = 0;

  tree->SetBranchAddress("values" , &values );
  tree->SetBranchAddress("weights", &weights);
  
  int nEntries = tree->GetEntries();

  for(int i = 0; i < nEntries; i++ ){
    tree->GetEntry(i);
    
    if (_dimension == -1) _dimension = values->size();
    if ( int(values->size()) != _dimension) ERROR_LOG << "You are loading HyperPoints that do not have the correct dimensionality for this HyperPointSet";

    HyperPoint point(*values);
    for(unsigned w = 0; w < weights->size(); w++){
      point.addWeight(weights->at(w));
    }
    this->push_back(point);
  }
  
  tree->ResetBranchAddresses();
  file->Close();

}

mean()

HyperPoint HyperPointSet::mean

(

)

const

Get the arithmatic mean of all points. This does not use the weights.

Definition at line 450 of file HyperPointSet.cpp.

                                    {
  HyperPoint temp(getDimension());
  for (unsigned i = 0; i < size(); i++){
    temp = temp + at(i);
  }
  return temp/double(size());
}

normalise()

void HyperPointSet::normalise

(

)

Normalise all points in the HyperPointSet.

Definition at line 160 of file HyperPointSet.cpp.

                             {

  for(unsigned i = 0; i < size(); i++){
    at(i) = at(i)/at(i).norm();
  }

}

print()

void HyperPointSet::print

(

std::ostream &

os = std::cout

)

const

Print out the entire HyperPointSet.

Definition at line 301 of file HyperPointSet.cpp.

                                             {
  
  os << "---------- HyperPointSet ------------" << std::endl;
  for (unsigned i = 0; i < size(); i++){
    os << "HyperPoint " << i << ":   ";
    at(i).print(os, 0);
    os << std::endl;
  }

}

push_back()

void HyperPointSet::push_back

(

const HyperPoint &

point

)

Add new HyperPoint to the HyperPointSet

Definition at line 91 of file HyperPointSet.cpp.

                                                    { 
  if (compatible(point) == false) return;
  _points.push_back(point); 
} 

removeDuplicates()

void HyperPointSet::removeDuplicates

(

)

Remove all duplicate events. IMPORTANT must be followed by sort()

Definition at line 179 of file HyperPointSet.cpp.

                                    {

  //INFO_LOG << "Removing all duplicates from the HyperPointSet - did you remember to call sort() first?!?";

  HyperPointSet hyperPointSet(getDimension());
  
  HyperPoint currentHyperPoint = _points.at(0); 
  hyperPointSet.push_back(currentHyperPoint);

  for(unsigned int i = 1; i < size(); i++){
    if ( _points.at(i) != currentHyperPoint ) {
      currentHyperPoint = _points.at(i); 
      hyperPointSet.push_back(currentHyperPoint);
    }
  }
  
  _points.clear();
  *this = hyperPointSet;

}

save() [1/2]

void HyperPointSet::save ( )

private

Save the HyperPointSet to the file currently in scope.

Definition at line 332 of file HyperPointSet.cpp.

                        {
 
  TTree* tree = new TTree("HyperPointSet", "HyperPointSet");
  
  if (tree == 0) {
    ERROR_LOG << "Cannot create tree";
    return;
  }

  std::vector<double>  values;
  std::vector<double>  weights;

  tree->Branch("values" , &values );
  tree->Branch("weights", &weights);
  
  for(unsigned i = 0; i < size(); i++ ){
    weights = at(i).getWeightsVector();
    values  = at(i).getVector();
    tree->Fill();
  }
  
  tree->ResetBranchAddresses();

  tree->Write();

}

save() [2/2]

void HyperPointSet::save

(

TString

path

)

Save the HyperPointSet to a file (opened with RECREATE).

Definition at line 314 of file HyperPointSet.cpp.

                                    {

  TFile* file = new TFile(path, "RECREATE");
  
  if (file == 0) {
    ERROR_LOG << "Cannot open file at " << path;
    return;
  }

  save();

  file->Write();
  file->Close();  

}

size()

unsigned int HyperPointSet::size

(

)

const

Get the number of points in the HyperPointSet

Definition at line 84 of file HyperPointSet.cpp.

                                                    { 
  return _points.size();
}

sort()

void HyperPointSet::sort

(

)

Sort all the HyperPoints in acending order. First sort by the first element, then by the second etc.

Definition at line 171 of file HyperPointSet.cpp.

                        {

  std::sort(_points.begin(), _points.end());

}

Member Data Documentation

_dimension

int HyperPointSet::_dimension

private

The dimensionality of HyperPoints in the HyperPointSet

Definition at line 33 of file HyperPointSet.h.

_points

std::vector<HyperPoint> HyperPointSet::_points

private

std::vector containing the HyperPoints

Definition at line 34 of file HyperPointSet.h.

---

The documentation for this class was generated from the following files:

• Mint/HyperPointSet.h
• src/Mojito/Chi2Test/HyperPointSet.cpp