ClusterizingHistogram Class Reference

#include <ClusterizingHistogram.h>

Public Member Functions
int	bin (float x) const
int	bin (double x) const
float	bin_pos (int i) const
std::vector< float >	clusterize (float resolution)
	ClusterizingHistogram (int nb, float xmi, float xma)
void	dump () const
void	dump (int i1, int i2) const
void	dump (float x1, float x2) const
void	dump (double x1, double x2) const
void	dump (float x1, double x2) const
void	dump (double x1, float x2) const
int	entries () const
void	fill (float x)
float	max_x () const
float	min_x () const
int	nbins () const
int	overflows () const
void	reset ()
int	underflows () const
	~ClusterizingHistogram ()

Private Member Functions
	ClusterizingHistogram ()

Private Attributes
int *	bin_entries
float *	bin_means
float	binsiz
int	my_entries
int	my_nbins
int	my_overflows
int	my_underflows
float	xmax
float	xmin

Detailed Description

A very simple 1D equidistant bin histogram that has the ability to clusterize it's contents. The bin entries are averaged in X, giving more accurate indication of where the bin contents are than the center of the bin.

Definition at line 13 of file ClusterizingHistogram.h.

Constructor & Destructor Documentation

ClusterizingHistogram::ClusterizingHistogram	(	int	nb,
		float	xmi,
		float	xma
	)

Definition at line 6 of file ClusterizingHistogram.cc.

References bin_entries, bin_means, binsiz, mps_fire::i, my_nbins, xmax, and xmin.

    : my_nbins(nb), xmin(xmi), xmax(xma), my_entries(0), my_underflows(0), my_overflows(0) {
  bin_entries = new int[my_nbins];
  bin_means = new float[my_nbins];
  binsiz = (xmax - xmin) / my_nbins;
  for (int i = 0; i < my_nbins; i++) {
    bin_entries[i] = 0;
    bin_means[i] = 0.;
  }
}

ClusterizingHistogram::~ClusterizingHistogram

(

)

Definition at line 16 of file ClusterizingHistogram.cc.

References bin_entries, and bin_means.

                                              {
  delete[] bin_entries;
  delete[] bin_means;
}

ClusterizingHistogram::ClusterizingHistogram ( )

inlineprivate

Definition at line 39 of file ClusterizingHistogram.h.

{}  // Prohibit

Member Function Documentation

int ClusterizingHistogram::bin

(

float

x

)

const

Definition at line 21 of file ClusterizingHistogram.cc.

References binsiz, my_nbins, xmax, and xmin.

Referenced by dump(), and fill().

                                            {
  if (x < xmin)
    return -1;
  else if (x > xmax)
    return my_nbins;
  else
    return int((x - xmin) / binsiz);
}

int ClusterizingHistogram::bin

(

double

x

)

const

Definition at line 29 of file ClusterizingHistogram.cc.

References binsiz, my_nbins, xmax, and xmin.

                                             {
  if (x < xmin)
    return -1;
  else if (x > xmax)
    return my_nbins;
  else
    return int((x - xmin) / binsiz);
}

float ClusterizingHistogram::bin_pos ( int  i ) const

inline

Definition at line 25 of file ClusterizingHistogram.h.

References bin_entries, bin_means, and mps_fire::i.

Referenced by clusterize(), and dump().

{ return (bin_entries[i] != 0) ? bin_means[i] / bin_entries[i] : 0; }

vector< float > ClusterizingHistogram::clusterize

(

float

resolution

)

Definition at line 54 of file ClusterizingHistogram.cc.

References bin_entries, bin_means, bin_pos(), mps_fire::i, my_nbins, DeadROCCounter::nclust, l1ParticleFlow_cff::resol, and xmin.

Referenced by magneticfield::bRod::bRod(), magneticfield::bSector::bSector(), magneticfield::MagGeoBuilder::build(), MagGeoBuilderFromDDD::build(), and magneticfield::eSector::eSector().

                                                           {
  vector<float> clust;
  int nclust = 0;
  bool inclust = false;
  float last_pos = xmin - 1000. * resol;
  int sum = 0;
  float sumx = 0;
  for (int i = 0; i < my_nbins; i++) {
    if (bin_entries[i] != 0) {
      if (fabs(bin_pos(i) - last_pos) > resol) {
        inclust = false;
        if (nclust != 0)
          clust.push_back(sumx / sum);  // create cluster
      }
      if (!inclust) {
        nclust++;
        sumx = 0.;
        sum = 0;
      }
      sum += bin_entries[i];
      sumx += bin_means[i];
      last_pos = bin_pos(i);
      inclust = true;
    }
  }
  if (nclust != 0)
    clust.push_back(sumx / sum);  // create last cluster
  return clust;
}

void ClusterizingHistogram::dump

(

void

)

const

Definition at line 84 of file ClusterizingHistogram.cc.

References my_nbins.

Referenced by dump().

{ dump(0, my_nbins); }

void ClusterizingHistogram::dump	(	int	i1,
		int	i2
	)		const

Definition at line 86 of file ClusterizingHistogram.cc.

References bin_entries, bin_pos(), gather_cfg::cout, mps_fire::i, SiStripPI::max, SiStripPI::min, my_entries, my_nbins, my_overflows, and my_underflows.

                                                     {
  cout << "Dumping ClusterizingHistogram contents:" << endl;
  for (int i = max(i1, 0); i < min(i2, my_nbins); i++) {
    cout << i << "  " << bin_entries[i] << "   " << bin_pos(i) << endl;
  }
  cout << "Underflows: " << my_underflows << endl;
  cout << "Overflows:  " << my_overflows << endl;
  cout << "Total number of entries: " << my_entries << endl;
}

void ClusterizingHistogram::dump	(	float	x1,
		float	x2
	)		const

Definition at line 96 of file ClusterizingHistogram.cc.

References bin(), and dump().

{ dump(bin(x1), bin(x2)); }

void ClusterizingHistogram::dump	(	double	x1,
		double	x2
	)		const

Definition at line 97 of file ClusterizingHistogram.cc.

References bin(), and dump().

{ dump(bin(x1), bin(x2)); }

void ClusterizingHistogram::dump	(	float	x1,
		double	x2
	)		const

Definition at line 98 of file ClusterizingHistogram.cc.

References bin(), and dump().

{ dump(bin(x1), bin(x2)); }

void ClusterizingHistogram::dump	(	double	x1,
		float	x2
	)		const

Definition at line 99 of file ClusterizingHistogram.cc.

References bin(), and dump().

{ dump(bin(x1), bin(x2)); }

int ClusterizingHistogram::entries ( ) const

inline

Definition at line 22 of file ClusterizingHistogram.h.

References my_entries.

{ return my_entries; }

void ClusterizingHistogram::fill

(

float

x

)

Definition at line 38 of file ClusterizingHistogram.cc.

References bin(), bin_entries, bin_means, binsiz, my_entries, my_nbins, my_overflows, my_underflows, x, xmax, and xmin.

Referenced by magneticfield::bRod::bRod(), magneticfield::bSector::bSector(), magneticfield::MagGeoBuilder::build(), MagGeoBuilderFromDDD::build(), and magneticfield::eSector::eSector().

                                        {
  if (x < xmin)
    my_underflows++;
  else if (x > xmax)
    my_overflows++;
  else {
    int bin = int((x - xmin) / binsiz);
    if (bin > my_nbins - 1)
      bin = my_nbins - 1;
    ++bin_entries[bin];
    bin_means[bin] += x;
    my_entries++;
    // may be problematic for negative x; check!
  }
}

float ClusterizingHistogram::max_x ( ) const

inline

Definition at line 21 of file ClusterizingHistogram.h.

References xmax.

{ return xmax; }

float ClusterizingHistogram::min_x ( ) const

inline

Definition at line 20 of file ClusterizingHistogram.h.

References xmin.

{ return xmin; }

int ClusterizingHistogram::nbins ( ) const

inline

Definition at line 19 of file ClusterizingHistogram.h.

References my_nbins.

{ return my_nbins; }

int ClusterizingHistogram::overflows ( ) const

inline

Definition at line 24 of file ClusterizingHistogram.h.

References my_overflows.

{ return my_overflows; }

void ClusterizingHistogram::reset

(

void

)

Definition at line 101 of file ClusterizingHistogram.cc.

References bin_entries, bin_means, mps_fire::i, my_entries, my_nbins, my_overflows, and my_underflows.

                                  {
  my_entries = 0;
  my_underflows = 0;
  my_overflows = 0;
  for (int i = 0; i < my_nbins; i++) {
    bin_entries[i] = 0;
    bin_means[i] = 0.;
  }
}

int ClusterizingHistogram::underflows ( ) const

inline

Definition at line 23 of file ClusterizingHistogram.h.

References my_underflows.

{ return my_underflows; }

Member Data Documentation

int* ClusterizingHistogram::bin_entries

private

Definition at line 46 of file ClusterizingHistogram.h.

Referenced by bin_pos(), clusterize(), ClusterizingHistogram(), dump(), fill(), reset(), and ~ClusterizingHistogram().

float* ClusterizingHistogram::bin_means

private

Definition at line 47 of file ClusterizingHistogram.h.

Referenced by bin_pos(), clusterize(), ClusterizingHistogram(), fill(), reset(), and ~ClusterizingHistogram().

float ClusterizingHistogram::binsiz

private

Definition at line 48 of file ClusterizingHistogram.h.

Referenced by bin(), ClusterizingHistogram(), and fill().

int ClusterizingHistogram::my_entries

private

Definition at line 43 of file ClusterizingHistogram.h.

Referenced by dump(), entries(), fill(), and reset().

int ClusterizingHistogram::my_nbins

private

Definition at line 40 of file ClusterizingHistogram.h.

Referenced by bin(), clusterize(), ClusterizingHistogram(), dump(), fill(), nbins(), and reset().

int ClusterizingHistogram::my_overflows

private

Definition at line 45 of file ClusterizingHistogram.h.

Referenced by dump(), fill(), overflows(), and reset().

int ClusterizingHistogram::my_underflows

private

Definition at line 44 of file ClusterizingHistogram.h.

Referenced by dump(), fill(), reset(), and underflows().

float ClusterizingHistogram::xmax

private

Definition at line 42 of file ClusterizingHistogram.h.

Referenced by svgfig.XAxis::__repr__(), bin(), ClusterizingHistogram(), fill(), and max_x().

float ClusterizingHistogram::xmin

private

Definition at line 41 of file ClusterizingHistogram.h.

Referenced by svgfig.XAxis::__repr__(), svgfig.Axes::__repr__(), svgfig.HGrid::__repr__(), svgfig.Grid::__repr__(), bin(), clusterize(), ClusterizingHistogram(), fill(), min_x(), and svgfig.Axes::SVG().