Histogram2D Class Reference

Two dimensional histogram. More...

Inheritance diagram for Histogram2D:

Collaboration diagram for Histogram2D:

Public Types
enum	DataType { INTEGRAL , AVERAGE , STANDARD_ERROR , NENTRIES }

Public Member Functions
	Histogram2D (const IAxis &axis_x, const IAxis &axis_y)
	Constructor for 2D histogram with custom axes. More...
	Histogram2D (const OutputData< double > &data)
	Constructor for 2D histograms from basic OutputData object. More...
	Histogram2D (int nbinsx, const std::vector< double > &xbins, int nbinsy, const std::vector< double > &ybins)
	Constructor for variable bin size histograms. More...
	Histogram2D (int nbinsx, double xlow, double xup, int nbinsy, double ylow, double yup)
	Constructor for fix bin size histograms. More...
	Histogram2D (std::vector< std::vector< double >> data)
	Constructor for 2D histograms from numpy array (thanks to swig) More...
void	addBinContent (size_t i, double value)
	Add the value to the bin. More...
void	addContent (const std::vector< std::vector< double >> &data)
	Add to values in histograms channels from numpy array,. More...
PyObject *	array (DataType dataType=DataType::INTEGRAL) const
	Returns numpy array with bin content (accumulated values). More...
double	binAverage (size_t binx, size_t biny) const
	Returns average value of the bin with given indices (for 2D histograms). More...
double	binAverage (size_t i) const
	Returns average value in the bin with given index. More...
double	binContent (size_t binx, size_t biny) const
	Returns content (accumulated value) of the 2D histogram bin. More...
double	binContent (size_t i) const
	Returns content (accumulated value) of bin i. More...
double	binError (size_t binx, size_t biny) const
	Returns error of the bin with given indices (for 2D histograms). More...
double	binError (size_t i) const
	Returns error of the bin with given index. More...
int	binNumberOfEntries (size_t binx, size_t biny) const
	Returns number of entries in the bin with given indices (for 2D histograms). More...
int	binNumberOfEntries (size_t i) const
	Returns number of entries in the bin with given index. More...
Histogram2D *	clone () const
	Returns clone of other histogram. More...
OutputData< double > *	createOutputData (DataType dataType=DataType::INTEGRAL) const
	creates new OutputData with histogram's shape and values corresponding to DataType More...
Histogram2D *	crop (double xmin, double ymin, double xmax, double ymax)
	Creates new histogram by applying rectangular clip. More...
int	fill (double x, double y, double weight=1.0)
	Increment bin with abscissa x and ordinate y with a weight. More...
size_t	findGlobalBin (double x, double y) const
	Returns closest global bin index for given axes coordinates. More...
PyObject *	getArray (DataType dataType=DataType::INTEGRAL) const
OutputData< CumulativeValue > &	getData ()
const OutputData< CumulativeValue > &	getData () const
size_t	getGlobalBin (size_t binx, size_t biny=0) const
	Returns global bin index for given axes indices. For 1D, just returns binx. More...
double	getMaximum () const
	Returns histogram maximum value (maximum of binContent() over all bins) More...
size_t	getMaximumBinIndex () const
	Returns globalbin index with maximum content. More...
double	getMinimum () const
	Returns histogram minimum value (minimum of binContent() over all bins) More...
size_t	getMinimumBinIndex () const
	Returns globalbin index with minimum content. More...
size_t	getNbinsX () const
	Returns number of bins on x-axis. More...
size_t	getNbinsY () const
	Returns number of bins on y-axis. More...
size_t	getTotalNumberOfBins () const
	Returns total number of histogram bins. More...
double	getXmax () const
	Returns x-axis max (upper edge of last bin). More...
double	getXmin () const
	Returns x-axis min (lower edge of first bin). More...
double	getYmax () const
	Returns y-axis max (upper edge of last bin) for 2D histograms. More...
double	getYmin () const
	Returns y-axis min (lower edge of first bin) for 2D histograms. More...
bool	hasSameDimensions (const IHistogram &other) const
	Returns true if object have same rank and number of axes bins. More...
bool	hasSameShape (const IHistogram &other) const
	Returns true if objects a) have same dimensions b) bin boundaries of axes coincide. More...
double	integral () const
	Returns integral of bins content (computed as a sum of all bin content). More...
void	load (const std::string &filename)
	Loads histogram from file, the shape of array in file should match Following formats are available: *.txt, *.tif, *.int (*.txt.gz, *.tif.gz, *.int.gz) Only bin content will be loaded, histogram axes remain the same. More...
const IHistogram &	operator+= (const IHistogram &right)
	addition-assignment operator for two histograms More...
Histogram1D *	projectionX ()
	Project a 2D histogram into 1D histogram along X. More...
Histogram1D *	projectionX (double ylow, double yup)
	Project a 2D histogram into 1D histogram along X. More...
Histogram1D *	projectionX (double yvalue)
	Project a 2D histogram into 1D histogram along X. More...
Histogram1D *	projectionY ()
	Project a 2D histogram into 1D histogram along Y. More...
Histogram1D *	projectionY (double xlow, double xup)
	Project a 2D histogram into 1D histogram along Y. More...
Histogram1D *	projectionY (double xvalue)
	Project a 2D histogram into 1D histogram along Y. More...
size_t	rank () const
	Returns the number of histogram dimensions. More...
IHistogram *	relativeDifferenceHistogram (const IHistogram &rhs)
	returns histogram representing relative difference of two histograms. More...
void	reset ()
	Reset histogram content (axes remains) More...
void	save (const std::string &filename)
	Saves histogram in file Following formats are available: *.txt, *.tif, *.int (*.txt.gz, *.tif.gz, *.int.gz) More...
void	scale (double value)
	Multiply this histogram (every bin content value) by a constant. More...
void	setBinContent (size_t i, double value)
	Sets content of the bin corresponding to the globalbin number. More...
void	setContent (const std::vector< std::vector< double >> &data)
	Sets the values in histograms channels from numpy array,. More...
const IAxis &	xAxis () const
	returns x-axis More...
size_t	xAxisIndex (size_t i) const
	Returns x-axis index for global bin index i. More...
double	xAxisValue (size_t i)
	Returns the center of bin i of the x axis. More...
const IAxis &	yAxis () const
	returns y-axis for 2D histograms More...
size_t	yAxisIndex (size_t i) const
	Returns y-axis index for global bin index i. More...
double	yAxisValue (size_t i)
	Returns the center of bin i of the y axis. More...

Static Public Member Functions
static IHistogram *	createFrom (const std::string &filename)
	create new histogram from file content More...
static IHistogram *	createFrom (const std::vector< std::vector< double >> &data)
	create new histogram from numpy array More...
static IHistogram *	createHistogram (const OutputData< double > &source)

Protected Member Functions
double	binData (size_t i, DataType dataType) const
	returns data of requested type for globalbin number More...
void	check_x_axis () const
void	check_y_axis () const
void	copyContentFrom (const IHistogram &other)
	Copy content (but not the axes) from other histogram. Dimensions should be the same. More...
Histogram1D *	create_projectionX (int ybinlow, int ybinup)
	Creates projection along X. More...
Histogram1D *	create_projectionY (int xbinlow, int xbinup)
	Creates projection along Y. More...
std::vector< double >	getDataVector (DataType dataType) const
	returns vector of values of requested DataType More...
void	init_from_data (const OutputData< double > &source)
template<typename T >
void	initFromShape (const T &data)

Protected Attributes
OutputData< CumulativeValue >	m_data

Detailed Description

Two dimensional histogram.

Definition at line 24 of file Histogram2D.h.

Member Enumeration Documentation

DataType

enum IHistogram::DataType

inherited

Enumerator
INTEGRAL
AVERAGE
STANDARD_ERROR
NENTRIES

Definition at line 29 of file IHistogram.h.

{ INTEGRAL, AVERAGE, STANDARD_ERROR, NENTRIES };

Constructor & Destructor Documentation

Histogram2D() [1/5]

Histogram2D::Histogram2D	(	int	nbinsx,
		double	xlow,
		double	xup,
		int	nbinsy,
		double	ylow,
		double	yup
	)

Constructor for fix bin size histograms.

Parameters

nbinsx	number of bins on X-axis
xlow	low edge of the first bin of X-axis
xup	upper edge of the last bin of X-axis
nbinsy	number of bins on Y axis
ylow	low edge of the first bin of Y-axis
yup	upper edge of the last bin of Y-axis

Definition at line 20 of file Histogram2D.cpp.

{
    m_data.addAxis(FixedBinAxis("x-axis", nbinsx, xlow, xup));
    m_data.addAxis(FixedBinAxis("y-axis", nbinsy, ylow, yup));
}

References OutputData< T >::addAxis(), and IHistogram::m_data.

Referenced by clone(), and crop().

Here is the call graph for this function:

Histogram2D() [2/5]

Histogram2D::Histogram2D	(	int	nbinsx,
		const std::vector< double > &	xbins,
		int	nbinsy,
		const std::vector< double > &	ybins
	)

Constructor for variable bin size histograms.

Parameters

nbinsx	number of bins on X-axis
xbins	Array of size nbins+1 containing low-edges for each bin and upper edge of last bin.
nbinsy	number of bins on Y-axis
ybins	Array of size nbins+1 containing low-edges for each bin and upper edge of last bin.

Definition at line 26 of file Histogram2D.cpp.

{
    m_data.addAxis(VariableBinAxis("x-axis", nbinsx, xbins));
    m_data.addAxis(VariableBinAxis("y-axis", nbinsy, ybins));
}

References OutputData< T >::addAxis(), and IHistogram::m_data.

Here is the call graph for this function:

Histogram2D() [3/5]

Histogram2D::Histogram2D	(	const IAxis &	axis_x,
		const IAxis &	axis_y
	)

Constructor for 2D histogram with custom axes.

Definition at line 33 of file Histogram2D.cpp.

: IHistogram(axis_x, axis_y) {}

Histogram2D() [4/5]

Histogram2D::Histogram2D

(

const OutputData< double > &

data

)

Constructor for 2D histograms from basic OutputData object.

Definition at line 35 of file Histogram2D.cpp.

{
    init_from_data(data);
}

References IHistogram::init_from_data().

Here is the call graph for this function:

Histogram2D() [5/5]

Histogram2D::Histogram2D

(

std::vector< std::vector< double >>

data

)

Constructor for 2D histograms from numpy array (thanks to swig)

Definition at line 41 of file Histogram2D.cpp.

{
    initFromShape(data);
    this->setContent(data);
}

References initFromShape(), and setContent().

Here is the call graph for this function:

Member Function Documentation

addBinContent()

void IHistogram::addBinContent ( size_t  i, double  value  )

inherited

Add the value to the bin.

Definition at line 151 of file IHistogram.cpp.

{
    m_data[i].add(value);
}

References IHistogram::m_data.

Referenced by IHistogram::operator+=().

addContent()

void Histogram2D::addContent

(

const std::vector< std::vector< double >> &

data

)

Add to values in histograms channels from numpy array,.

Definition at line 125 of file Histogram2D.cpp.

{
    auto shape = ArrayUtils::getShape(data);
    const size_t nrows = shape.first;
    const size_t ncols = shape.second;
 
    if (nrows != m_data.axis(1).size() || ncols != m_data.axis(0).size()) {
        std::ostringstream ostr;
        ostr << "Histogram2D::addContent() -> Shape of input array [" << nrows << ", " << ncols
             << "] doesn't mach histogram axes. "
             << "X-axis size: " << m_data.axis(0).size()
             << "Y-axis size: " << m_data.axis(1).size();
        throw std::runtime_error(ostr.str());
    }
 
    for (size_t row = 0; row < nrows; ++row) {
        for (size_t col = 0; col < ncols; ++col) {
            size_t globalbin = nrows - row - 1 + col * nrows;
            m_data[globalbin].add(data[row][col]);
        }
    }
}

References OutputData< T >::axis(), ArrayUtils::getShape(), IHistogram::m_data, and IAxis::size().

Referenced by setContent().

Here is the call graph for this function:

array()

PyObject * IHistogram::array ( DataType  dataType = DataType::INTEGRAL ) const

inherited

Returns numpy array with bin content (accumulated values).

Definition at line 226 of file IHistogram.cpp.

{
    const std::unique_ptr<OutputData<double>> data(createOutputData(dataType));
    return data->getArray();
}

References IHistogram::createOutputData().

Referenced by IHistogram::getArray().

Here is the call graph for this function:

binAverage() [1/2]

double IHistogram::binAverage ( size_t  binx, size_t  biny  ) const

inherited

Returns average value of the bin with given indices (for 2D histograms).

Definition at line 171 of file IHistogram.cpp.

{
    return binAverage(getGlobalBin(binx, biny));
}

References IHistogram::binAverage(), and IHistogram::getGlobalBin().

Here is the call graph for this function:

binAverage() [2/2]

double IHistogram::binAverage ( size_t  i ) const

inherited

Returns average value in the bin with given index.

Definition at line 166 of file IHistogram.cpp.

{
    return m_data[i].getAverage();
}

References IHistogram::m_data.

Referenced by IHistogram::binAverage(), and IHistogram::binData().

binContent() [1/2]

double IHistogram::binContent ( size_t  binx, size_t  biny  ) const

inherited

Returns content (accumulated value) of the 2D histogram bin.

Definition at line 141 of file IHistogram.cpp.

{
    return binContent(getGlobalBin(binx, biny));
}

References IHistogram::binContent(), and IHistogram::getGlobalBin().

Here is the call graph for this function:

binContent() [2/2]

double IHistogram::binContent ( size_t  i ) const

inherited

Returns content (accumulated value) of bin i.

Definition at line 136 of file IHistogram.cpp.

{
    return m_data[i].getContent();
}

References IHistogram::m_data.

Referenced by IHistogram::binContent(), IHistogram::binData(), create_projectionX(), create_projectionY(), IHistogram::operator+=(), and IHistogram::relativeDifferenceHistogram().

binData()

double IHistogram::binData ( size_t  i, IHistogram::DataType  dataType  ) const

protectedinherited

returns data of requested type for globalbin number

Definition at line 305 of file IHistogram.cpp.

{
    if (dataType == DataType::INTEGRAL) {
        return binContent(i);
    } else if (dataType == DataType::AVERAGE) {
        return binAverage(i);
    } else if (dataType == DataType::STANDARD_ERROR) {
        return binError(i);
    } else if (dataType == DataType::NENTRIES) {
        return binNumberOfEntries(i);
    } else
        throw std::runtime_error("IHistogram::binData() -> Error. Unknown data type.");
}

References IHistogram::binAverage(), IHistogram::binContent(), IHistogram::binError(), and IHistogram::binNumberOfEntries().

Referenced by IHistogram::createOutputData(), and IHistogram::getDataVector().

Here is the call graph for this function:

binError() [1/2]

double IHistogram::binError ( size_t  binx, size_t  biny  ) const

inherited

Returns error of the bin with given indices (for 2D histograms).

Definition at line 161 of file IHistogram.cpp.

{
    return binError(getGlobalBin(binx, biny));
}

References IHistogram::binError(), and IHistogram::getGlobalBin().

Here is the call graph for this function:

binError() [2/2]

double IHistogram::binError ( size_t  i ) const

inherited

Returns error of the bin with given index.

Definition at line 156 of file IHistogram.cpp.

{
    return m_data[i].getRMS();
}

References IHistogram::m_data.

Referenced by IHistogram::binData(), and IHistogram::binError().

binNumberOfEntries() [1/2]

int IHistogram::binNumberOfEntries ( size_t  binx, size_t  biny  ) const

inherited

Returns number of entries in the bin with given indices (for 2D histograms).

Definition at line 181 of file IHistogram.cpp.

{
    return binNumberOfEntries(getGlobalBin(binx, biny));
}

References IHistogram::binNumberOfEntries(), and IHistogram::getGlobalBin().

Here is the call graph for this function:

binNumberOfEntries() [2/2]

int IHistogram::binNumberOfEntries ( size_t  i ) const

inherited

Returns number of entries in the bin with given index.

Definition at line 176 of file IHistogram.cpp.

{
    return m_data[i].getNumberOfEntries();
}

References IHistogram::m_data.

Referenced by IHistogram::binData(), and IHistogram::binNumberOfEntries().

check_x_axis()

void IHistogram::check_x_axis ( ) const

protectedinherited

Definition at line 268 of file IHistogram.cpp.

{
    if (rank() < 1) {
        std::ostringstream message;
        message << "IHistogram::check_x_axis() -> Error. X-xis does not exist. ";
        message << "Rank of histogram " << rank() << "." << std::endl;
        throw std::runtime_error(message.str());
    }
}

References IHistogram::rank().

Referenced by IHistogram::xAxis(), and IHistogram::xAxisValue().

Here is the call graph for this function:

check_y_axis()

void IHistogram::check_y_axis ( ) const

protectedinherited

Definition at line 278 of file IHistogram.cpp.

{
    if (rank() < 2) {
        std::ostringstream message;
        message << "IHistogram::check_y_axis() -> Error. Y-axis does not exist. ";
        message << "Rank of histogram " << rank() << "." << std::endl;
        throw std::runtime_error(message.str());
    }
}

References IHistogram::rank().

Referenced by IHistogram::yAxis(), and IHistogram::yAxisValue().

Here is the call graph for this function:

clone()

Histogram2D * Histogram2D::clone ( ) const

virtual

Returns clone of other histogram.

Implements IHistogram.

Definition at line 47 of file Histogram2D.cpp.

{
    return new Histogram2D(*this);
}

References Histogram2D().

Here is the call graph for this function:

copyContentFrom()

void IHistogram::copyContentFrom ( const IHistogram &  other )

protectedinherited

Copy content (but not the axes) from other histogram. Dimensions should be the same.

Definition at line 331 of file IHistogram.cpp.

{
    if (!hasSameDimensions(other))
        throw std::runtime_error(
            "IHistogram::copyContentFrom() -> Error. Can't copy the data of different shape.");
    reset();
    for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
        m_data[i] = other.m_data[i];
    }
}

References IHistogram::getTotalNumberOfBins(), IHistogram::hasSameDimensions(), IHistogram::m_data, and IHistogram::reset().

Referenced by IHistogram::load().

Here is the call graph for this function:

create_projectionX()

Histogram1D * Histogram2D::create_projectionX ( int  ybinlow, int  ybinup  )

protected

Creates projection along X.

The projections is made by collecting the data in the range between [ybinlow, ybinup].

Definition at line 148 of file Histogram2D.cpp.

{
    Histogram1D* result = new Histogram1D(this->xAxis());
 
    for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
 
        int ybin = static_cast<int>(yAxisIndex(index));
 
        if (ybin >= ybinlow && ybin <= ybinup) {
            result->fill(xAxisValue(index), binContent(index));
        }
    }
    return result;
}

References IHistogram::binContent(), Histogram1D::fill(), IHistogram::getTotalNumberOfBins(), IHistogram::xAxis(), IHistogram::xAxisValue(), and IHistogram::yAxisIndex().

Referenced by projectionX().

Here is the call graph for this function:

create_projectionY()

Histogram1D * Histogram2D::create_projectionY ( int  xbinlow, int  xbinup  )

protected

Creates projection along Y.

The projections is made by collecting the data in the range between [xbinlow, xbinup].

Definition at line 163 of file Histogram2D.cpp.

{
    Histogram1D* result = new Histogram1D(this->yAxis());
 
    for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
 
        int xbin = static_cast<int>(xAxisIndex(index));
 
        if (xbin >= xbinlow && xbin <= xbinup) {
            result->fill(yAxisValue(index), binContent(index));
        }
    }
    return result;
}

References IHistogram::binContent(), Histogram1D::fill(), IHistogram::getTotalNumberOfBins(), IHistogram::xAxisIndex(), IHistogram::yAxis(), and IHistogram::yAxisValue().

Referenced by projectionY().

Here is the call graph for this function:

createFrom() [1/2]

IHistogram * IHistogram::createFrom ( const std::string &  filename )

staticinherited

create new histogram from file content

Definition at line 258 of file IHistogram.cpp.

{
    return IntensityDataIOFactory::readIntensityData(filename);
}

References FileSystemUtils::filename(), and IntensityDataIOFactory::readIntensityData().

Here is the call graph for this function:

createFrom() [2/2]

IHistogram * IHistogram::createFrom ( const std::vector< std::vector< double >> &  data )

staticinherited

create new histogram from numpy array

Definition at line 263 of file IHistogram.cpp.

{
    return new Histogram2D(data);
}

createHistogram()

IHistogram * IHistogram::createHistogram ( const OutputData< double > &  source )

staticinherited

Definition at line 243 of file IHistogram.cpp.

{
    if (source.rank() == 1) {
        return new Histogram1D(source);
    } else if (source.rank() == 2) {
        return new Histogram2D(source);
    } else {
        std::ostringstream message;
        message << "IHistogram::createHistogram(const OutputData<double>& source) -> Error. ";
        message << "The rank of source " << source.rank() << " ";
        message << "is not suitable for creation neither 1-dim nor 2-dim histograms.";
        throw std::runtime_error(message.str());
    }
}

References OutputData< T >::rank().

Referenced by DetectorMask::createHistogram(), and IntensityDataIOFactory::readIntensityData().

Here is the call graph for this function:

createOutputData()

OutputData< double > * IHistogram::createOutputData ( IHistogram::DataType  dataType = DataType::INTEGRAL ) const

inherited

creates new OutputData with histogram's shape and values corresponding to DataType

creates new OutputData with histogram's shape and put there values corresponding to DataType

Definition at line 343 of file IHistogram.cpp.

{
    OutputData<double>* result = new OutputData<double>;
    result->copyShapeFrom(m_data);
    for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
        (*result)[i] = binData(i, dataType);
    }
    return result;
}

References IHistogram::binData(), OutputData< T >::copyShapeFrom(), IHistogram::getTotalNumberOfBins(), and IHistogram::m_data.

Referenced by IHistogram::array(), HistoUtils::FindPeaks(), and IntensityDataIOFactory::writeIntensityData().

Here is the call graph for this function:

crop()

Histogram2D * Histogram2D::crop	(	double	xmin,
		double	ymin,
		double	xmax,
		double	ymax
	)

Creates new histogram by applying rectangular clip.

Definition at line 99 of file Histogram2D.cpp.

{
    const std::unique_ptr<IAxis> xaxis(xAxis().createClippedAxis(xmin, xmax));
    const std::unique_ptr<IAxis> yaxis(yAxis().createClippedAxis(ymin, ymax));
 
    Histogram2D* result = new Histogram2D(*xaxis, *yaxis);
    OutputData<CumulativeValue>::const_iterator it_origin = m_data.begin();
    OutputData<CumulativeValue>::iterator it_result = result->m_data.begin();
    while (it_origin != m_data.end()) {
        double x = m_data.getAxisValue(it_origin.getIndex(), 0);
        double y = m_data.getAxisValue(it_origin.getIndex(), 1);
        if (result->xAxis().contains(x) && result->yAxis().contains(y)) {
            *it_result = *it_origin;
            ++it_result;
        }
        ++it_origin;
    }
    return result;
}

References Histogram2D(), OutputData< T >::begin(), IAxis::contains(), OutputData< T >::end(), OutputData< T >::getAxisValue(), IHistogram::m_data, IHistogram::xAxis(), and IHistogram::yAxis().

Here is the call graph for this function:

fill()

int Histogram2D::fill	(	double	x,
		double	y,
		double	weight = 1.0
	)

Increment bin with abscissa x and ordinate y with a weight.

Definition at line 52 of file Histogram2D.cpp.

{
    if (!xAxis().contains(x))
        return -1;
    if (!yAxis().contains(y))
        return -1;
    size_t index = m_data.findGlobalIndex({x, y});
    m_data[index].add(weight);
    return (int)index;
}

References OutputData< T >::findGlobalIndex(), IHistogram::m_data, IHistogram::xAxis(), and IHistogram::yAxis().

Here is the call graph for this function:

findGlobalBin()

size_t IHistogram::findGlobalBin ( double  x, double  y  ) const

inherited

Returns closest global bin index for given axes coordinates.

Definition at line 95 of file IHistogram.cpp.

{
    std::vector<double> coordinates;
    coordinates.push_back(x);
    if (rank() == 2)
        coordinates.push_back(y);
    return m_data.findGlobalIndex(coordinates);
}

References OutputData< T >::findGlobalIndex(), IHistogram::m_data, and IHistogram::rank().

Here is the call graph for this function:

getArray()

PyObject * IHistogram::getArray ( DataType  dataType = DataType::INTEGRAL ) const

inherited

Deprecated:

Use array() instead.

Definition at line 232 of file IHistogram.cpp.

{
    return array(dataType);
}

References IHistogram::array().

Here is the call graph for this function:

getData() [1/2]

OutputData< CumulativeValue > & IHistogram::getData ( )

inherited

Definition at line 131 of file IHistogram.cpp.

{
    return m_data;
}

References IHistogram::m_data.

getData() [2/2]

const OutputData< CumulativeValue > & IHistogram::getData ( ) const

inherited

Definition at line 126 of file IHistogram.cpp.

{
    return m_data;
}

References IHistogram::m_data.

Referenced by HistoUtils::getRelativeDifference().

getDataVector()

std::vector< double > IHistogram::getDataVector ( IHistogram::DataType  dataType ) const

protectedinherited

returns vector of values of requested DataType

Definition at line 320 of file IHistogram.cpp.

{
    std::vector<double> result;
    result.resize(getTotalNumberOfBins(), 0.0);
    for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
        result[index] = binData(index, dataType);
    }
    return result;
}

References IHistogram::binData(), and IHistogram::getTotalNumberOfBins().

Referenced by Histogram1D::binErrors(), and Histogram1D::binValues().

Here is the call graph for this function:

getGlobalBin()

size_t IHistogram::getGlobalBin ( size_t  binx, size_t  biny = 0  ) const

inherited

Returns global bin index for given axes indices. For 1D, just returns binx.

Definition at line 86 of file IHistogram.cpp.

{
    std::vector<unsigned> axes_indices;
    axes_indices.push_back(static_cast<unsigned>(binx));
    if (rank() == 2)
        axes_indices.push_back(static_cast<unsigned>(biny));
    return m_data.toGlobalIndex(axes_indices);
}

References IHistogram::m_data, IHistogram::rank(), and OutputData< T >::toGlobalIndex().

Referenced by IHistogram::binAverage(), IHistogram::binContent(), IHistogram::binError(), and IHistogram::binNumberOfEntries().

Here is the call graph for this function:

getMaximum()

double IHistogram::getMaximum ( ) const

inherited

Returns histogram maximum value (maximum of binContent() over all bins)

Definition at line 186 of file IHistogram.cpp.

{
    OutputData<CumulativeValue>::const_iterator it = std::max_element(m_data.begin(), m_data.end());
    return it->getContent();
}

References OutputData< T >::begin(), OutputData< T >::end(), and IHistogram::m_data.

Here is the call graph for this function:

getMaximumBinIndex()

size_t IHistogram::getMaximumBinIndex ( ) const

inherited

Returns globalbin index with maximum content.

Definition at line 192 of file IHistogram.cpp.

{
    OutputData<CumulativeValue>::const_iterator it = std::max_element(m_data.begin(), m_data.end());
    return std::distance(m_data.begin(), it);
}

References OutputData< T >::begin(), OutputData< T >::end(), and IHistogram::m_data.

Here is the call graph for this function:

getMinimum()

double IHistogram::getMinimum ( ) const

inherited

Returns histogram minimum value (minimum of binContent() over all bins)

Definition at line 198 of file IHistogram.cpp.

{
    OutputData<CumulativeValue>::const_iterator it = std::min_element(m_data.begin(), m_data.end());
    return it->getContent();
}

References OutputData< T >::begin(), OutputData< T >::end(), and IHistogram::m_data.

Here is the call graph for this function:

getMinimumBinIndex()

size_t IHistogram::getMinimumBinIndex ( ) const

inherited

Returns globalbin index with minimum content.

Definition at line 204 of file IHistogram.cpp.

{
    return std::distance(m_data.begin(), std::min_element(m_data.begin(), m_data.end()));
}

References OutputData< T >::begin(), OutputData< T >::end(), and IHistogram::m_data.

Here is the call graph for this function:

getNbinsX()

size_t IHistogram::getNbinsX ( ) const

inherited

Returns number of bins on x-axis.

Definition at line 66 of file IHistogram.cpp.

{
    return xAxis().size();
}

References IAxis::size(), and IHistogram::xAxis().

Here is the call graph for this function:

getNbinsY()

size_t IHistogram::getNbinsY ( ) const

inherited

Returns number of bins on y-axis.

Definition at line 81 of file IHistogram.cpp.

{
    return yAxis().size();
}

References IAxis::size(), and IHistogram::yAxis().

Here is the call graph for this function:

getTotalNumberOfBins()

size_t IHistogram::getTotalNumberOfBins ( ) const

inherited

Returns total number of histogram bins.

For 2D histograms the result will be the product of bin numbers along X and Y axes.

Definition at line 39 of file IHistogram.cpp.

{
    return m_data.getAllocatedSize();
}

References OutputData< T >::getAllocatedSize(), and IHistogram::m_data.

Referenced by IHistogram::copyContentFrom(), create_projectionX(), create_projectionY(), IHistogram::createOutputData(), IHistogram::getDataVector(), IHistogram::integral(), IHistogram::operator+=(), IHistogram::relativeDifferenceHistogram(), and IHistogram::scale().

Here is the call graph for this function:

getXmax()

double IHistogram::getXmax ( ) const

inherited

Returns x-axis max (upper edge of last bin).

Definition at line 61 of file IHistogram.cpp.

{
    return xAxis().upperBound();
}

References IAxis::upperBound(), and IHistogram::xAxis().

Here is the call graph for this function:

getXmin()

double IHistogram::getXmin ( ) const

inherited

Returns x-axis min (lower edge of first bin).

Definition at line 56 of file IHistogram.cpp.

{
    return xAxis().lowerBound();
}

References IAxis::lowerBound(), and IHistogram::xAxis().

Here is the call graph for this function:

getYmax()

double IHistogram::getYmax ( ) const

inherited

Returns y-axis max (upper edge of last bin) for 2D histograms.

Definition at line 76 of file IHistogram.cpp.

{
    return yAxis().upperBound();
}

References IAxis::upperBound(), and IHistogram::yAxis().

Here is the call graph for this function:

getYmin()

double IHistogram::getYmin ( ) const

inherited

Returns y-axis min (lower edge of first bin) for 2D histograms.

Definition at line 71 of file IHistogram.cpp.

{
    return yAxis().lowerBound();
}

References IAxis::lowerBound(), and IHistogram::yAxis().

Here is the call graph for this function:

hasSameDimensions()

bool IHistogram::hasSameDimensions ( const IHistogram &  other ) const

inherited

Returns true if object have same rank and number of axes bins.

Definition at line 358 of file IHistogram.cpp.

{
    return m_data.hasSameDimensions(other.m_data);
}

References OutputData< T >::hasSameDimensions(), and IHistogram::m_data.

Referenced by IHistogram::copyContentFrom(), IHistogram::operator+=(), and IHistogram::relativeDifferenceHistogram().

Here is the call graph for this function:

hasSameShape()

bool IHistogram::hasSameShape ( const IHistogram &  other ) const

inherited

Returns true if objects a) have same dimensions b) bin boundaries of axes coincide.

Definition at line 353 of file IHistogram.cpp.

{
    return m_data.hasSameShape(other.m_data);
}

References OutputData< T >::hasSameShape(), and IHistogram::m_data.

Here is the call graph for this function:

init_from_data()

void IHistogram::init_from_data ( const OutputData< double > &  source )

protectedinherited

Definition at line 288 of file IHistogram.cpp.

{
    if (rank() != source.rank()) {
        std::ostringstream message;
        message << "IHistogram::IHistogram(const OutputData<double>& data) -> Error. ";
        message << "The dimension of this histogram " << rank() << " ";
        message << "is differ from the dimension of source " << m_data.rank() << std::endl;
        throw std::runtime_error(message.str());
    }
 
    m_data.copyShapeFrom(source);
    for (size_t i = 0; i < source.getAllocatedSize(); ++i) {
        m_data[i].add(source[i]);
    }
}

References OutputData< T >::copyShapeFrom(), OutputData< T >::getAllocatedSize(), IHistogram::m_data, OutputData< T >::rank(), and IHistogram::rank().

Referenced by Histogram1D::Histogram1D(), and Histogram2D().

Here is the call graph for this function:

initFromShape()

template<typename T >

void Histogram2D::initFromShape ( const T &  data )

protected

Definition at line 114 of file Histogram2D.h.

{
    auto shape = ArrayUtils::getShape(data);
    const size_t nrows = shape.first;
    const size_t ncols = shape.second;
 
    if (nrows == 0 || ncols == 0)
        throw std::runtime_error("Histogram2D::Histogram2D() -> Error. "
                                 "Not a two-dimensional numpy array");
 
    m_data.addAxis(FixedBinAxis("x-axis", ncols, 0.0, static_cast<double>(ncols)));
    m_data.addAxis(FixedBinAxis("y-axis", nrows, 0.0, static_cast<double>(nrows)));
}

References OutputData< T >::addAxis(), ArrayUtils::getShape(), and IHistogram::m_data.

Referenced by Histogram2D().

Here is the call graph for this function:

integral()

double IHistogram::integral ( ) const

inherited

Returns integral of bins content (computed as a sum of all bin content).

Definition at line 216 of file IHistogram.cpp.

{
    double result(0.0);
    for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
        result += m_data[index].getContent();
    }
    return result;
}

References IHistogram::getTotalNumberOfBins(), and IHistogram::m_data.

Here is the call graph for this function:

load()

void IHistogram::load ( const std::string &  filename )

inherited

Loads histogram from file, the shape of array in file should match Following formats are available: *.txt, *.tif, *.int (*.txt.gz, *.tif.gz, *.int.gz) Only bin content will be loaded, histogram axes remain the same.

Definition at line 394 of file IHistogram.cpp.

{
    const std::unique_ptr<IHistogram> hist(IntensityDataIOFactory::readIntensityData(filename));
    copyContentFrom(*hist);
}

References IHistogram::copyContentFrom(), FileSystemUtils::filename(), and IntensityDataIOFactory::readIntensityData().

Here is the call graph for this function:

operator+=()

const IHistogram & IHistogram::operator+= ( const IHistogram &  right )

inherited

addition-assignment operator for two histograms

Definition at line 363 of file IHistogram.cpp.

{
    if (!hasSameDimensions(right))
        throw std::runtime_error(
            "IHistogram::operator+=() -> Error. Histograms have different dimension");
    for (size_t i = 0; i < getTotalNumberOfBins(); ++i)
        addBinContent(i, right.binContent(i));
    return *this;
}

References IHistogram::addBinContent(), IHistogram::binContent(), IHistogram::getTotalNumberOfBins(), and IHistogram::hasSameDimensions().

Here is the call graph for this function:

projectionX() [1/3]

Histogram1D * Histogram2D::projectionX

(

)

Project a 2D histogram into 1D histogram along X.

The projection is made from all bins along y-axis.

Definition at line 63 of file Histogram2D.cpp.

{
    return create_projectionX(0, static_cast<int>(xAxis().size()) - 1);
}

References create_projectionX(), and IHistogram::xAxis().

Here is the call graph for this function:

projectionX() [2/3]

Histogram1D * Histogram2D::projectionX	(	double	ylow,
		double	yup
	)

Project a 2D histogram into 1D histogram along X.

The projection is made from all y-bins corresponding to ordinate between ylow and yup.

Parameters

ylow	lower edje on y-axis
yup	upper edje on y-axis

Definition at line 74 of file Histogram2D.cpp.

{
    int ybinlow = static_cast<int>(yAxis().findClosestIndex(ylow));
    int ybinup = static_cast<int>(yAxis().findClosestIndex(yup));
    return create_projectionX(ybinlow, ybinup);
}

References create_projectionX(), IAxis::findClosestIndex(), and IHistogram::yAxis().

Here is the call graph for this function:

projectionX() [3/3]

Histogram1D * Histogram2D::projectionX

(

double

yvalue

)

Project a 2D histogram into 1D histogram along X.

The projection is made from the y-bin closest to given ordinate yvalue.

Parameters

yvalue

the value on y-axis at which projection is taken

Definition at line 68 of file Histogram2D.cpp.

{
    int ybin_selected = static_cast<int>(yAxis().findClosestIndex(yvalue));
    return create_projectionX(ybin_selected, ybin_selected);
}

References create_projectionX(), IAxis::findClosestIndex(), and IHistogram::yAxis().

Here is the call graph for this function:

projectionY() [1/3]

Histogram1D * Histogram2D::projectionY

(

)

Project a 2D histogram into 1D histogram along Y.

The projection is made from all bins along x-axis.

Definition at line 81 of file Histogram2D.cpp.

{
    return create_projectionY(0, static_cast<int>(xAxis().size()) - 1);
}

References create_projectionY(), and IHistogram::xAxis().

Here is the call graph for this function:

projectionY() [2/3]

Histogram1D * Histogram2D::projectionY	(	double	xlow,
		double	xup
	)

Project a 2D histogram into 1D histogram along Y.

The projection is made from all x-bins corresponding to abscissa between xlow and xup.

Parameters

xlow	lower edje on x-axis
xup	upper edje on x-axis

Definition at line 92 of file Histogram2D.cpp.

{
    int xbinlow = static_cast<int>(xAxis().findClosestIndex(xlow));
    int xbinup = static_cast<int>(xAxis().findClosestIndex(xup));
    return create_projectionY(xbinlow, xbinup);
}

References create_projectionY(), IAxis::findClosestIndex(), and IHistogram::xAxis().

Here is the call graph for this function:

projectionY() [3/3]

Histogram1D * Histogram2D::projectionY

(

double

xvalue

)

Project a 2D histogram into 1D histogram along Y.

The projection is made from the x-bin closest to given abscissa xvalue.

Parameters

xvalue

the value on x-axis at which projection is taken

Definition at line 86 of file Histogram2D.cpp.

{
    int xbin_selected = static_cast<int>(xAxis().findClosestIndex(xvalue));
    return create_projectionY(xbin_selected, xbin_selected);
}

References create_projectionY(), IAxis::findClosestIndex(), and IHistogram::xAxis().

Here is the call graph for this function:

rank()

size_t Histogram2D::rank ( ) const

inlinevirtual

Returns the number of histogram dimensions.

Implements IHistogram.

Definition at line 58 of file Histogram2D.h.

{ return 2; }

relativeDifferenceHistogram()

IHistogram * IHistogram::relativeDifferenceHistogram ( const IHistogram &  rhs )

inherited

returns histogram representing relative difference of two histograms.

Definition at line 373 of file IHistogram.cpp.

{
    if (!hasSameDimensions(rhs))
        throw std::runtime_error("IHistogram::relativeDifferenceHistogram() -> Error. "
                                 "Histograms have different dimensions");
 
    IHistogram* result = this->clone();
    result->reset();
 
    for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
        double diff = Numeric::GetRelativeDifference(binContent(i), rhs.binContent(i));
        result->setBinContent(i, diff);
    }
    return result;
}

References IHistogram::binContent(), IHistogram::clone(), Numeric::GetRelativeDifference(), IHistogram::getTotalNumberOfBins(), IHistogram::hasSameDimensions(), IHistogram::reset(), and IHistogram::setBinContent().

Here is the call graph for this function:

reset()

void IHistogram::reset ( )

inherited

Reset histogram content (axes remains)

Definition at line 238 of file IHistogram.cpp.

{
    m_data.setAllTo(CumulativeValue());
}

References IHistogram::m_data, and OutputData< T >::setAllTo().

Referenced by IHistogram::copyContentFrom(), IHistogram::relativeDifferenceHistogram(), and setContent().

Here is the call graph for this function:

save()

void IHistogram::save ( const std::string &  filename )

inherited

Saves histogram in file Following formats are available: *.txt, *.tif, *.int (*.txt.gz, *.tif.gz, *.int.gz)

Definition at line 389 of file IHistogram.cpp.

{
    IntensityDataIOFactory::writeIntensityData(*this, filename);
}

References FileSystemUtils::filename(), and IntensityDataIOFactory::writeIntensityData().

Here is the call graph for this function:

scale()

void IHistogram::scale ( double  value )

inherited

Multiply this histogram (every bin content value) by a constant.

Definition at line 209 of file IHistogram.cpp.

{
    for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
        m_data[index].setContent(value * m_data[index].getContent());
    }
}

References IHistogram::getTotalNumberOfBins(), and IHistogram::m_data.

Here is the call graph for this function:

setBinContent()

void IHistogram::setBinContent ( size_t  i, double  value  )

inherited

Sets content of the bin corresponding to the globalbin number.

Definition at line 146 of file IHistogram.cpp.

{
    m_data[i].setContent(value);
}

References IHistogram::m_data.

Referenced by IHistogram::relativeDifferenceHistogram().

setContent()

void Histogram2D::setContent

(

const std::vector< std::vector< double >> &

data

)

Sets the values in histograms channels from numpy array,.

Definition at line 119 of file Histogram2D.cpp.

{
    reset();
    addContent(data);
}

References addContent(), and IHistogram::reset().

Referenced by Histogram2D().

Here is the call graph for this function:

xAxis()

const IAxis & IHistogram::xAxis ( ) const

inherited

returns x-axis

Definition at line 44 of file IHistogram.cpp.

{
    check_x_axis();
    return m_data.axis(0);
}

References OutputData< T >::axis(), IHistogram::check_x_axis(), and IHistogram::m_data.

Referenced by Histogram1D::binCenters(), create_projectionX(), Histogram1D::crop(), crop(), Histogram1D::fill(), fill(), HistoUtils::FindPeaks(), IHistogram::getNbinsX(), IHistogram::getXmax(), IHistogram::getXmin(), projectionX(), and projectionY().

Here is the call graph for this function:

xAxisIndex()

size_t IHistogram::xAxisIndex ( size_t  i ) const

inherited

Returns x-axis index for global bin index i.

Definition at line 104 of file IHistogram.cpp.

{
    return m_data.getAxisBinIndex(i, 0);
}

References OutputData< T >::getAxisBinIndex(), and IHistogram::m_data.

Referenced by create_projectionY().

Here is the call graph for this function:

xAxisValue()

double IHistogram::xAxisValue ( size_t  i )

inherited

Returns the center of bin i of the x axis.

Definition at line 114 of file IHistogram.cpp.

{
    check_x_axis();
    return m_data.getAxisValue(i, 0);
}

References IHistogram::check_x_axis(), OutputData< T >::getAxisValue(), and IHistogram::m_data.

Referenced by create_projectionX().

Here is the call graph for this function:

yAxis()

const IAxis & IHistogram::yAxis ( ) const

inherited

returns y-axis for 2D histograms

Definition at line 50 of file IHistogram.cpp.

{
    check_y_axis();
    return m_data.axis(1);
}

References OutputData< T >::axis(), IHistogram::check_y_axis(), and IHistogram::m_data.

Referenced by create_projectionY(), crop(), fill(), HistoUtils::FindPeaks(), IHistogram::getNbinsY(), IHistogram::getYmax(), IHistogram::getYmin(), and projectionX().

Here is the call graph for this function:

yAxisIndex()

size_t IHistogram::yAxisIndex ( size_t  i ) const

inherited

Returns y-axis index for global bin index i.

Definition at line 109 of file IHistogram.cpp.

{
    return m_data.getAxisBinIndex(i, 1);
}

References OutputData< T >::getAxisBinIndex(), and IHistogram::m_data.

Referenced by create_projectionX().

Here is the call graph for this function:

yAxisValue()

double IHistogram::yAxisValue ( size_t  i )

inherited

Returns the center of bin i of the y axis.

Definition at line 120 of file IHistogram.cpp.

{
    check_y_axis();
    return m_data.getAxisValue(i, 1);
}

References IHistogram::check_y_axis(), OutputData< T >::getAxisValue(), and IHistogram::m_data.

Referenced by create_projectionY().

Here is the call graph for this function:

Member Data Documentation

m_data

OutputData<CumulativeValue> IHistogram::m_data

protectedinherited

Definition at line 193 of file IHistogram.h.

Referenced by Histogram1D::Histogram1D(), Histogram2D(), IHistogram::IHistogram(), IHistogram::addBinContent(), addContent(), IHistogram::binAverage(), IHistogram::binContent(), IHistogram::binError(), IHistogram::binNumberOfEntries(), IHistogram::copyContentFrom(), IHistogram::createOutputData(), Histogram1D::crop(), crop(), Histogram1D::fill(), fill(), IHistogram::findGlobalBin(), IHistogram::getData(), IHistogram::getGlobalBin(), IHistogram::getMaximum(), IHistogram::getMaximumBinIndex(), IHistogram::getMinimum(), IHistogram::getMinimumBinIndex(), IHistogram::getTotalNumberOfBins(), IHistogram::hasSameDimensions(), IHistogram::hasSameShape(), IHistogram::init_from_data(), initFromShape(), IHistogram::integral(), IHistogram::reset(), IHistogram::scale(), IHistogram::setBinContent(), IHistogram::xAxis(), IHistogram::xAxisIndex(), IHistogram::xAxisValue(), IHistogram::yAxis(), IHistogram::yAxisIndex(), and IHistogram::yAxisValue().

---

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

• /home/www/ba/Device/Histo/Histogram2D.h
• /home/www/ba/Device/Histo/Histogram2D.cpp