BornAgain  1.19.0
Simulate and fit neutron and x-ray scattering at grazing incidence
IHistogram Class Referenceabstract

Base class for 1D and 2D histograms holding values of double type. More...

Inheritance diagram for IHistogram:
[legend]
Collaboration diagram for IHistogram:
[legend]

Public Types

enum  DataType { INTEGRAL , AVERAGE , STANDARD_ERROR , NENTRIES }
 

Public Member Functions

 IHistogram ()
 
 IHistogram (const IAxis &axis_x)
 
 IHistogram (const IAxis &axis_x, const IAxis &axis_y)
 
 IHistogram (const IHistogram &other)
 
virtual ~IHistogram ()=default
 
void addBinContent (size_t i, double value)
 Add the value to the bin. More...
 
PyObjectarray (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...
 
virtual IHistogramclone () const =0
 
OutputData< double > * createOutputData (DataType dataType=DataType::INTEGRAL) const
 creates new OutputData with histogram's shape and values corresponding to DataType More...
 
size_t findGlobalBin (double x, double y) const
 Returns closest global bin index for given axes coordinates. More...
 
PyObjectgetArray (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 IHistogramoperator+= (const IHistogram &right)
 addition-assignment operator for two histograms More...
 
virtual size_t rank () const =0
 Returns number of histogram dimensions. More...
 
IHistogramrelativeDifferenceHistogram (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...
 
const IAxisxAxis () 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 IAxisyAxis () 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 IHistogramcreateFrom (const std::string &filename)
 create new histogram from file content More...
 
static IHistogramcreateFrom (const std::vector< std::vector< double >> &data)
 create new histogram from numpy array More...
 
static IHistogramcreateHistogram (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...
 
std::vector< double > getDataVector (DataType dataType) const
 returns vector of values of requested DataType More...
 
void init_from_data (const OutputData< double > &source)
 

Protected Attributes

OutputData< CumulativeValuem_data
 

Detailed Description

Base class for 1D and 2D histograms holding values of double type.

Definition at line 27 of file IHistogram.h.

Member Enumeration Documentation

◆ DataType

Enumerator
INTEGRAL 
AVERAGE 
STANDARD_ERROR 
NENTRIES 

Definition at line 29 of file IHistogram.h.

Constructor & Destructor Documentation

◆ IHistogram() [1/4]

IHistogram::IHistogram ( )
default

◆ IHistogram() [2/4]

IHistogram::IHistogram ( const IHistogram other)

Definition at line 23 of file IHistogram.cpp.

24 {
25  m_data.copyFrom(other.m_data);
26 }
OutputData< CumulativeValue > m_data
Definition: IHistogram.h:193
void copyFrom(const OutputData< T > &x)
Definition: OutputData.h:267

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

Here is the call graph for this function:

◆ ~IHistogram()

virtual IHistogram::~IHistogram ( )
virtualdefault

◆ IHistogram() [3/4]

IHistogram::IHistogram ( const IAxis axis_x)

Definition at line 28 of file IHistogram.cpp.

29 {
30  m_data.addAxis(axis_x);
31 }
void addAxis(const IAxis &new_axis)
Definition: OutputData.h:295

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

Here is the call graph for this function:

◆ IHistogram() [4/4]

IHistogram::IHistogram ( const IAxis axis_x,
const IAxis axis_y 
)

Definition at line 33 of file IHistogram.cpp.

34 {
35  m_data.addAxis(axis_x);
36  m_data.addAxis(axis_y);
37 }

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

Here is the call graph for this function:

Member Function Documentation

◆ addBinContent()

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

Add the value to the bin.

Definition at line 151 of file IHistogram.cpp.

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

References m_data.

Referenced by operator+=().

◆ array()

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

Returns numpy array with bin content (accumulated values).

Definition at line 226 of file IHistogram.cpp.

227 {
228  const std::unique_ptr<OutputData<double>> data(createOutputData(dataType));
229  return data->getArray();
230 }
OutputData< double > * createOutputData(DataType dataType=DataType::INTEGRAL) const
creates new OutputData with histogram's shape and values corresponding to DataType
Definition: IHistogram.cpp:343

References createOutputData().

Referenced by getArray().

Here is the call graph for this function:

◆ binAverage() [1/2]

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

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

Definition at line 171 of file IHistogram.cpp.

172 {
173  return binAverage(getGlobalBin(binx, biny));
174 }
size_t getGlobalBin(size_t binx, size_t biny=0) const
Returns global bin index for given axes indices. For 1D, just returns binx.
Definition: IHistogram.cpp:86
double binAverage(size_t i) const
Returns average value in the bin with given index.
Definition: IHistogram.cpp:166

References binAverage(), and getGlobalBin().

Here is the call graph for this function:

◆ binAverage() [2/2]

double IHistogram::binAverage ( size_t  i) const

Returns average value in the bin with given index.

Definition at line 166 of file IHistogram.cpp.

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

References m_data.

Referenced by binAverage(), and binData().

◆ binContent() [1/2]

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

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

Definition at line 141 of file IHistogram.cpp.

142 {
143  return binContent(getGlobalBin(binx, biny));
144 }
double binContent(size_t i) const
Returns content (accumulated value) of bin i.
Definition: IHistogram.cpp:136

References binContent(), and getGlobalBin().

Here is the call graph for this function:

◆ binContent() [2/2]

double IHistogram::binContent ( size_t  i) const

Returns content (accumulated value) of bin i.

Definition at line 136 of file IHistogram.cpp.

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

References m_data.

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

◆ binData()

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

returns data of requested type for globalbin number

Definition at line 305 of file IHistogram.cpp.

306 {
307  if (dataType == DataType::INTEGRAL) {
308  return binContent(i);
309  } else if (dataType == DataType::AVERAGE) {
310  return binAverage(i);
311  } else if (dataType == DataType::STANDARD_ERROR) {
312  return binError(i);
313  } else if (dataType == DataType::NENTRIES) {
314  return binNumberOfEntries(i);
315  } else
316  throw std::runtime_error("IHistogram::binData() -> Error. Unknown data type.");
317 }
double binError(size_t i) const
Returns error of the bin with given index.
Definition: IHistogram.cpp:156
int binNumberOfEntries(size_t i) const
Returns number of entries in the bin with given index.
Definition: IHistogram.cpp:176

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

Referenced by createOutputData(), and getDataVector().

Here is the call graph for this function:

◆ binError() [1/2]

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

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

Definition at line 161 of file IHistogram.cpp.

162 {
163  return binError(getGlobalBin(binx, biny));
164 }

References binError(), and getGlobalBin().

Here is the call graph for this function:

◆ binError() [2/2]

double IHistogram::binError ( size_t  i) const

Returns error of the bin with given index.

Definition at line 156 of file IHistogram.cpp.

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

References m_data.

Referenced by binData(), and binError().

◆ binNumberOfEntries() [1/2]

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

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

Definition at line 181 of file IHistogram.cpp.

182 {
183  return binNumberOfEntries(getGlobalBin(binx, biny));
184 }

References binNumberOfEntries(), and getGlobalBin().

Here is the call graph for this function:

◆ binNumberOfEntries() [2/2]

int IHistogram::binNumberOfEntries ( size_t  i) const

Returns number of entries in the bin with given index.

Definition at line 176 of file IHistogram.cpp.

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

References m_data.

Referenced by binData(), and binNumberOfEntries().

◆ check_x_axis()

void IHistogram::check_x_axis ( ) const
protected

Definition at line 268 of file IHistogram.cpp.

269 {
270  if (rank() < 1) {
271  std::ostringstream message;
272  message << "IHistogram::check_x_axis() -> Error. X-xis does not exist. ";
273  message << "Rank of histogram " << rank() << "." << std::endl;
274  throw std::runtime_error(message.str());
275  }
276 }
virtual size_t rank() const =0
Returns number of histogram dimensions.

References rank().

Referenced by xAxis(), and xAxisValue().

Here is the call graph for this function:

◆ check_y_axis()

void IHistogram::check_y_axis ( ) const
protected

Definition at line 278 of file IHistogram.cpp.

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

References rank().

Referenced by yAxis(), and yAxisValue().

Here is the call graph for this function:

◆ clone()

virtual IHistogram* IHistogram::clone ( ) const
pure virtual

Implemented in Histogram2D, and Histogram1D.

Referenced by relativeDifferenceHistogram().

◆ copyContentFrom()

void IHistogram::copyContentFrom ( const IHistogram other)
protected

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

Definition at line 331 of file IHistogram.cpp.

332 {
333  if (!hasSameDimensions(other))
334  throw std::runtime_error(
335  "IHistogram::copyContentFrom() -> Error. Can't copy the data of different shape.");
336  reset();
337  for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
338  m_data[i] = other.m_data[i];
339  }
340 }
size_t getTotalNumberOfBins() const
Returns total number of histogram bins.
Definition: IHistogram.cpp:39
void reset()
Reset histogram content (axes remains)
Definition: IHistogram.cpp:238
bool hasSameDimensions(const IHistogram &other) const
Returns true if object have same rank and number of axes bins.
Definition: IHistogram.cpp:358

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

Referenced by load().

Here is the call graph for this function:

◆ createFrom() [1/2]

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

create new histogram from file content

Definition at line 258 of file IHistogram.cpp.

259 {
261 }
static IHistogram * readIntensityData(const std::string &file_name)
Reads file and returns newly created Histogram object.
std::string filename(const std::string &path)
Returns path without directory part ("Foo/Bar/Doz.int.gz" -> "Doz.int.gz")

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)
static

create new histogram from numpy array

Definition at line 263 of file IHistogram.cpp.

264 {
265  return new Histogram2D(data);
266 }
Two dimensional histogram.
Definition: Histogram2D.h:24

◆ createHistogram()

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

Definition at line 243 of file IHistogram.cpp.

244 {
245  if (source.rank() == 1) {
246  return new Histogram1D(source);
247  } else if (source.rank() == 2) {
248  return new Histogram2D(source);
249  } else {
250  std::ostringstream message;
251  message << "IHistogram::createHistogram(const OutputData<double>& source) -> Error. ";
252  message << "The rank of source " << source.rank() << " ";
253  message << "is not suitable for creation neither 1-dim nor 2-dim histograms.";
254  throw std::runtime_error(message.str());
255  }
256 }
One dimensional histogram.
Definition: Histogram1D.h:23
size_t rank() const
Returns number of dimensions.
Definition: OutputData.h:56

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

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.

344 {
346  result->copyShapeFrom(m_data);
347  for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
348  (*result)[i] = binData(i, dataType);
349  }
350  return result;
351 }
double binData(size_t i, DataType dataType) const
returns data of requested type for globalbin number
Definition: IHistogram.cpp:305
void copyShapeFrom(const OutputData< U > &other)
Definition: OutputData.h:277

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

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

Here is the call graph for this function:

◆ findGlobalBin()

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

Returns closest global bin index for given axes coordinates.

Definition at line 95 of file IHistogram.cpp.

96 {
97  std::vector<double> coordinates;
98  coordinates.push_back(x);
99  if (rank() == 2)
100  coordinates.push_back(y);
101  return m_data.findGlobalIndex(coordinates);
102 }
size_t findGlobalIndex(const std::vector< double > &coordinates) const
Returns global index for specified axes values.
Definition: OutputData.h:416

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

Here is the call graph for this function:

◆ getArray()

PyObject * IHistogram::getArray ( DataType  dataType = DataType::INTEGRAL) const
Deprecated:
Use array() instead.

Definition at line 232 of file IHistogram.cpp.

233 {
234  return array(dataType);
235 }
PyObject * array(DataType dataType=DataType::INTEGRAL) const
Returns numpy array with bin content (accumulated values).
Definition: IHistogram.cpp:226

References array().

Here is the call graph for this function:

◆ getData() [1/2]

OutputData< CumulativeValue > & IHistogram::getData ( )

Definition at line 131 of file IHistogram.cpp.

132 {
133  return m_data;
134 }

References m_data.

◆ getData() [2/2]

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

Definition at line 126 of file IHistogram.cpp.

127 {
128  return m_data;
129 }

References m_data.

Referenced by HistoUtils::getRelativeDifference().

◆ getDataVector()

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

returns vector of values of requested DataType

Definition at line 320 of file IHistogram.cpp.

321 {
322  std::vector<double> result;
323  result.resize(getTotalNumberOfBins(), 0.0);
324  for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
325  result[index] = binData(index, dataType);
326  }
327  return result;
328 }

References binData(), and 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

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

Definition at line 86 of file IHistogram.cpp.

87 {
88  std::vector<unsigned> axes_indices;
89  axes_indices.push_back(static_cast<unsigned>(binx));
90  if (rank() == 2)
91  axes_indices.push_back(static_cast<unsigned>(biny));
92  return m_data.toGlobalIndex(axes_indices);
93 }
size_t toGlobalIndex(const std::vector< unsigned > &axes_indices) const
Returns global index for specified indices of axes.
Definition: OutputData.h:392

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

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

Here is the call graph for this function:

◆ getMaximum()

double IHistogram::getMaximum ( ) const

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

Definition at line 186 of file IHistogram.cpp.

187 {
189  return it->getContent();
190 }
iterator end()
Returns read/write iterator that points to the one past last element.
Definition: OutputData.h:93
iterator begin()
Returns read/write iterator that points to the first element.
Definition: OutputData.h:343

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

Here is the call graph for this function:

◆ getMaximumBinIndex()

size_t IHistogram::getMaximumBinIndex ( ) const

Returns globalbin index with maximum content.

Definition at line 192 of file IHistogram.cpp.

193 {
195  return std::distance(m_data.begin(), it);
196 }

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

Here is the call graph for this function:

◆ getMinimum()

double IHistogram::getMinimum ( ) const

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

Definition at line 198 of file IHistogram.cpp.

199 {
201  return it->getContent();
202 }

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

Here is the call graph for this function:

◆ getMinimumBinIndex()

size_t IHistogram::getMinimumBinIndex ( ) const

Returns globalbin index with minimum content.

Definition at line 204 of file IHistogram.cpp.

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

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

Here is the call graph for this function:

◆ getNbinsX()

size_t IHistogram::getNbinsX ( ) const

Returns number of bins on x-axis.

Definition at line 66 of file IHistogram.cpp.

67 {
68  return xAxis().size();
69 }
virtual size_t size() const =0
retrieve the number of bins
const IAxis & xAxis() const
returns x-axis
Definition: IHistogram.cpp:44

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

Here is the call graph for this function:

◆ getNbinsY()

size_t IHistogram::getNbinsY ( ) const

Returns number of bins on y-axis.

Definition at line 81 of file IHistogram.cpp.

82 {
83  return yAxis().size();
84 }
const IAxis & yAxis() const
returns y-axis for 2D histograms
Definition: IHistogram.cpp:50

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

Here is the call graph for this function:

◆ getTotalNumberOfBins()

size_t IHistogram::getTotalNumberOfBins ( ) const

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.

40 {
41  return m_data.getAllocatedSize();
42 }
size_t getAllocatedSize() const
Returns total size of data buffer (product of bin number in every dimension).
Definition: OutputData.h:59

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

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

Here is the call graph for this function:

◆ getXmax()

double IHistogram::getXmax ( ) const

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

Definition at line 61 of file IHistogram.cpp.

62 {
63  return xAxis().upperBound();
64 }
virtual double upperBound() const =0
Returns value of last point of axis.

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

Here is the call graph for this function:

◆ getXmin()

double IHistogram::getXmin ( ) const

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

Definition at line 56 of file IHistogram.cpp.

57 {
58  return xAxis().lowerBound();
59 }
virtual double lowerBound() const =0
Returns value of first point of axis.

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

Here is the call graph for this function:

◆ getYmax()

double IHistogram::getYmax ( ) const

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

Definition at line 76 of file IHistogram.cpp.

77 {
78  return yAxis().upperBound();
79 }

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

Here is the call graph for this function:

◆ getYmin()

double IHistogram::getYmin ( ) const

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

Definition at line 71 of file IHistogram.cpp.

72 {
73  return yAxis().lowerBound();
74 }

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

Here is the call graph for this function:

◆ hasSameDimensions()

bool IHistogram::hasSameDimensions ( const IHistogram other) const

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

Definition at line 358 of file IHistogram.cpp.

359 {
360  return m_data.hasSameDimensions(other.m_data);
361 }
bool hasSameDimensions(const OutputData< U > &right) const
Returns true if object have same dimensions and number of axes bins.
Definition: OutputData.h:575

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

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

Here is the call graph for this function:

◆ hasSameShape()

bool IHistogram::hasSameShape ( const IHistogram other) const

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

Definition at line 353 of file IHistogram.cpp.

354 {
355  return m_data.hasSameShape(other.m_data);
356 }
bool hasSameShape(const OutputData< U > &right) const
Returns true if objects a) have same dimensions b) bin boundaries of axes coincide.
Definition: OutputData.h:592

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

Here is the call graph for this function:

◆ init_from_data()

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

Definition at line 288 of file IHistogram.cpp.

289 {
290  if (rank() != source.rank()) {
291  std::ostringstream message;
292  message << "IHistogram::IHistogram(const OutputData<double>& data) -> Error. ";
293  message << "The dimension of this histogram " << rank() << " ";
294  message << "is differ from the dimension of source " << m_data.rank() << std::endl;
295  throw std::runtime_error(message.str());
296  }
297 
298  m_data.copyShapeFrom(source);
299  for (size_t i = 0; i < source.getAllocatedSize(); ++i) {
300  m_data[i].add(source[i]);
301  }
302 }

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

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

Here is the call graph for this function:

◆ integral()

double IHistogram::integral ( ) const

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

Definition at line 216 of file IHistogram.cpp.

217 {
218  double result(0.0);
219  for (size_t index = 0; index < getTotalNumberOfBins(); ++index) {
220  result += m_data[index].getContent();
221  }
222  return result;
223 }

References getTotalNumberOfBins(), and m_data.

Here is the call graph for this function:

◆ load()

void IHistogram::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.

Definition at line 394 of file IHistogram.cpp.

395 {
396  const std::unique_ptr<IHistogram> hist(IntensityDataIOFactory::readIntensityData(filename));
397  copyContentFrom(*hist);
398 }
void copyContentFrom(const IHistogram &other)
Copy content (but not the axes) from other histogram. Dimensions should be the same.
Definition: IHistogram.cpp:331

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

Here is the call graph for this function:

◆ operator+=()

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

addition-assignment operator for two histograms

Definition at line 363 of file IHistogram.cpp.

364 {
365  if (!hasSameDimensions(right))
366  throw std::runtime_error(
367  "IHistogram::operator+=() -> Error. Histograms have different dimension");
368  for (size_t i = 0; i < getTotalNumberOfBins(); ++i)
369  addBinContent(i, right.binContent(i));
370  return *this;
371 }
void addBinContent(size_t i, double value)
Add the value to the bin.
Definition: IHistogram.cpp:151

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

Here is the call graph for this function:

◆ rank()

virtual size_t IHistogram::rank ( ) const
pure virtual

Returns number of histogram dimensions.

Implemented in Histogram2D, and Histogram1D.

Referenced by check_x_axis(), check_y_axis(), findGlobalBin(), getGlobalBin(), and init_from_data().

◆ relativeDifferenceHistogram()

IHistogram * IHistogram::relativeDifferenceHistogram ( const IHistogram rhs)

returns histogram representing relative difference of two histograms.

Definition at line 373 of file IHistogram.cpp.

374 {
375  if (!hasSameDimensions(rhs))
376  throw std::runtime_error("IHistogram::relativeDifferenceHistogram() -> Error. "
377  "Histograms have different dimensions");
378 
379  IHistogram* result = this->clone();
380  result->reset();
381 
382  for (size_t i = 0; i < getTotalNumberOfBins(); ++i) {
383  double diff = Numeric::GetRelativeDifference(binContent(i), rhs.binContent(i));
384  result->setBinContent(i, diff);
385  }
386  return result;
387 }
Base class for 1D and 2D histograms holding values of double type.
Definition: IHistogram.h:27
void setBinContent(size_t i, double value)
Sets content of the bin corresponding to the globalbin number.
Definition: IHistogram.cpp:146
virtual IHistogram * clone() const =0
double GetRelativeDifference(double a, double b)
Returns the safe relative difference, which is 2(|a-b|)/(|a|+|b|) except in special cases.
Definition: Numeric.cpp:29

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

Here is the call graph for this function:

◆ reset()

void IHistogram::reset ( )

Reset histogram content (axes remains)

Definition at line 238 of file IHistogram.cpp.

239 {
241 }
The cumulative value with average and rms on-the-flight calculations.
void setAllTo(const T &value)
Sets content of output data to specific value.
Definition: OutputData.h:476

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

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

Here is the call graph for this function:

◆ save()

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

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.

390 {
392 }
static void writeIntensityData(const IHistogram &histogram, const std::string &file_name)
Writes histogram in file.

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

Here is the call graph for this function:

◆ scale()

void IHistogram::scale ( double  value)

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

Definition at line 209 of file IHistogram.cpp.

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

References getTotalNumberOfBins(), and m_data.

Here is the call graph for this function:

◆ setBinContent()

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

Sets content of the bin corresponding to the globalbin number.

Definition at line 146 of file IHistogram.cpp.

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

References m_data.

Referenced by relativeDifferenceHistogram().

◆ xAxis()

const IAxis & IHistogram::xAxis ( ) const

returns x-axis

Definition at line 44 of file IHistogram.cpp.

45 {
46  check_x_axis();
47  return m_data.axis(0);
48 }
void check_x_axis() const
Definition: IHistogram.cpp:268
const IAxis & axis(size_t serial_number) const
returns axis with given serial number
Definition: OutputData.h:318

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

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

Here is the call graph for this function:

◆ xAxisIndex()

size_t IHistogram::xAxisIndex ( size_t  i) const

Returns x-axis index for global bin index i.

Definition at line 104 of file IHistogram.cpp.

105 {
106  return m_data.getAxisBinIndex(i, 0);
107 }
size_t getAxisBinIndex(size_t global_index, size_t i_selected_axis) const
Returns axis bin index for given global index.
Definition: OutputData.h:370

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

Referenced by Histogram2D::create_projectionY().

Here is the call graph for this function:

◆ xAxisValue()

double IHistogram::xAxisValue ( size_t  i)

Returns the center of bin i of the x axis.

Definition at line 114 of file IHistogram.cpp.

115 {
116  check_x_axis();
117  return m_data.getAxisValue(i, 0);
118 }
double getAxisValue(size_t global_index, size_t i_selected_axis) const
Returns the value of selected axis for given global_index.
Definition: OutputData.h:430

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

Referenced by Histogram2D::create_projectionX().

Here is the call graph for this function:

◆ yAxis()

const IAxis & IHistogram::yAxis ( ) const

returns y-axis for 2D histograms

Definition at line 50 of file IHistogram.cpp.

51 {
52  check_y_axis();
53  return m_data.axis(1);
54 }
void check_y_axis() const
Definition: IHistogram.cpp:278

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

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

Here is the call graph for this function:

◆ yAxisIndex()

size_t IHistogram::yAxisIndex ( size_t  i) const

Returns y-axis index for global bin index i.

Definition at line 109 of file IHistogram.cpp.

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

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

Referenced by Histogram2D::create_projectionX().

Here is the call graph for this function:

◆ yAxisValue()

double IHistogram::yAxisValue ( size_t  i)

Returns the center of bin i of the y axis.

Definition at line 120 of file IHistogram.cpp.

121 {
122  check_y_axis();
123  return m_data.getAxisValue(i, 1);
124 }

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

Referenced by Histogram2D::create_projectionY().

Here is the call graph for this function:

Member Data Documentation

◆ m_data


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