npstat::HistoAxis Class Reference

#include <HistoAxis.h>

Public Member Functions
double	binCenter (const int binNum) const
Interval< double >	binInterval (const int binNum) const
int	binNumber (double x) const
LinearMapper1d	binNumberMapper (bool mapLeftEdgeTo0=true) const
double	binWidth (const int=0) const
gs::ClassId	classId () const
unsigned	closestValidBin (double x) const
double	fltBinNumber (const double x, const bool mapLeftEdgeTo0=true) const
	HistoAxis (unsigned nBins, double min, double max, const char *label=nullptr)
Interval< double >	interval () const
bool	isClose (const HistoAxis &, double tol) const
bool	isUniform () const
CircularMapper1d	kernelScanMapper (bool doubleRange) const
const std::string &	label () const
double	leftBinEdge (const int binNum) const
double	length () const
double	max () const
double	min () const
unsigned	nBins () const
bool	operator!= (const HistoAxis &) const
bool	operator== (const HistoAxis &) const
HistoAxis	rebin (unsigned newBins) const
double	rightBinEdge (const int binNum) const
void	setLabel (const char *newlabel)
bool	write (std::ostream &of) const

Static Public Member Functions
static const char *	classname ()
static HistoAxis *	read (const gs::ClassId &id, std::istream &in)
static unsigned	version ()

Private Member Functions
	HistoAxis ()
unsigned	overflowIndex (const double x, unsigned *binNumber) const
unsigned	overflowIndexWeighted (double x, unsigned *binNumber, double *weight) const

Private Attributes
double	bw_
std::string	label_
double	max_
double	min_
unsigned	nBins_

Friends
class	DualHistoAxis
template<typename Numeric , class Axis >
class	HistoND

Detailed Description

Class which contain the information needed to define a histogram axis. All bins will have the same width. See NUHistoAxis and DualHistoAxis classes for non-uniform binning.

Definition at line 31 of file HistoAxis.h.

Constructor & Destructor Documentation

HistoAxis() [1/2]

npstat::HistoAxis::HistoAxis	(	unsigned	nBins,
		double	min,
		double	max,
		const char *	label = nullptr
	)

Minimum and maximum will be internally swapped if the minimum parameter is larger than the maximum

Definition at line 13 of file HistoAxis.cc.

      : min_(min), max_(max), label_(label ? label : ""), nBins_(nbins) {
    if (!(nBins_ && nBins_ < UINT_MAX / 2U - 1U))
      throw npstat::NpstatInvalidArgument(
          "In npstat::HistoAxis constructor: "
          "number of bins is out of range");
    if (min_ > max_)
      std::swap(min_, max_);
    bw_ = (max_ - min_) / nBins_;
  }

References bw_, max_, min_, nBins_, std::swap(), and mitigatedMETSequence_cff::U.

HistoAxis() [2/2]

npstat::HistoAxis::HistoAxis ( )

inlineprivate

Definition at line 149 of file HistoAxis.h.

: min_(0.0), max_(0.0), bw_(0.0), nBins_(0) {}

Referenced by read(), and rebin().

Member Function Documentation

binCenter()

double npstat::HistoAxis::binCenter ( const int  binNum ) const

inline

Return the coordinate of the given bin center

Definition at line 52 of file HistoAxis.h.

{ return min_ + (binNum + 0.5) * bw_; }

References bw_, and min_.

Referenced by npstat::DualHistoAxis::binCenter(), and npstat::convertToGridAxis().

binInterval()

Interval<double> npstat::HistoAxis::binInterval ( const int  binNum ) const

inline

Return the coordinate interval occupied by the given bin

Definition at line 61 of file HistoAxis.h.

                                                                {
      return Interval<double>(min_ + binNum * bw_, min_ + (binNum + 1) * bw_);
    }

References bw_, and min_.

Referenced by npstat::DualHistoAxis::binInterval().

binNumber()

int npstat::HistoAxis::binNumber

(

double

x

)

const

This method returns arbitrary integer bin number, including negative numbers and numbers which can exceed nBins()-1

Definition at line 37 of file HistoAxis.cc.

                                               {
    if (bw_) {
      int binnum = static_cast<int>(floor((x - min_) / bw_));
      if (x >= max_) {
        if (binnum < static_cast<int>(nBins_))
          binnum = nBins_;
      } else {
        if (binnum >= static_cast<int>(nBins_))
          binnum = nBins_ - 1U;
      }
      return binnum;
    } else {
      if (x < min_)
        return -1;
      else
        return nBins_;
    }
  }

References bw_, max_, min_, nBins_, and mitigatedMETSequence_cff::U.

Referenced by npstat::DualHistoAxis::binNumber(), overflowIndex(), and overflowIndexWeighted().

binNumberMapper()

LinearMapper1d npstat::HistoAxis::binNumberMapper

(

bool

mapLeftEdgeTo0 = true

)

const

Return the mapper which calculates floating point bin number given the coordinate. The resulting bin number can go above and below the axis range. If "mapLeftEdgeTo0" is specified as "false", it is the center of the first bin which gets mapped to 0.

Definition at line 67 of file HistoAxis.cc.

                                                                           {
    if (!bw_)
      throw npstat::NpstatDomainError(
          "In npstat::HistoAxis::binNumberMapper: "
          "bin width is zero. Mapper can not be constructed.");
    const double base = mapLeftEdgeTo0 ? min_ / bw_ : min_ / bw_ + 0.5;
    return LinearMapper1d(1.0 / bw_, -base);
  }

References newFWLiteAna::base, bw_, and min_.

binWidth()

double npstat::HistoAxis::binWidth ( const int  = 0 ) const

inline

Definition at line 46 of file HistoAxis.h.

{ return bw_; }

References bw_.

Referenced by npstat::DualHistoAxis::binWidth().

classId()

gs::ClassId npstat::HistoAxis::classId ( ) const

inline

Method related to "geners" I/O

Definition at line 117 of file HistoAxis.h.

{ return gs::ClassId(*this); }

Referenced by npstat::DualHistoAxis::write().

classname()

static const char* npstat::HistoAxis::classname ( )

inlinestatic

Definition at line 121 of file HistoAxis.h.

{ return "npstat::HistoAxis"; }

closestValidBin()

unsigned npstat::HistoAxis::closestValidBin

(

double

x

)

const

This method returns the closest valid bin number (above 0 and below nBins() )

Definition at line 56 of file HistoAxis.cc.

                                                          {
    if (x <= min_)
      return 0U;
    else if (bw_ && x < max_) {
      const unsigned binnum = static_cast<unsigned>(floor((x - min_) / bw_));
      if (binnum < nBins_)
        return binnum;
    }
    return nBins_ - 1U;
  }

References bw_, max_, min_, nBins_, and mitigatedMETSequence_cff::U.

Referenced by npstat::DualHistoAxis::closestValidBin().

fltBinNumber()

double npstat::HistoAxis::fltBinNumber ( const double  x, const bool  mapLeftEdgeTo0 = true  ) const

inline

Floating point bin number given the coordinate (no bin number truncation of any kind is performed). Works in exactly the same way as the mapper returned by the previous method.

Definition at line 94 of file HistoAxis.h.

                                                                                       {
      return (x - min_) / bw_ - (mapLeftEdgeTo0 ? 0.0 : 0.5);
    }

References bw_, and min_.

Referenced by npstat::DualHistoAxis::fltBinNumber().

interval()

Interval<double> npstat::HistoAxis::interval ( ) const

inline

Definition at line 43 of file HistoAxis.h.

{ return Interval<double>(min_, max_); }

References max_, and min_.

Referenced by npstat::DualHistoAxis::interval().

isClose()

bool npstat::HistoAxis::isClose	(	const HistoAxis &	r,
		double	tol
	)		const

Comparison of axis coordinates within given tolerance

Definition at line 26 of file HistoAxis.cc.

                                                                    {
    return closeWithinTolerance(min_, r.min_, tol) && closeWithinTolerance(max_, r.max_, tol) && label_ == r.label_ &&
           nBins_ == r.nBins_;
  }

References npstat::closeWithinTolerance(), label_, max_, min_, nBins_, and alignCSCRings::r.

Referenced by npstat::DualHistoAxis::isClose().

isUniform()

bool npstat::HistoAxis::isUniform ( ) const

inline

Definition at line 48 of file HistoAxis.h.

{ return true; }

kernelScanMapper()

CircularMapper1d npstat::HistoAxis::kernelScanMapper

(

bool

doubleRange

)

const

The following function returns a mapper that can be helpful in scanning a kernel (a density function) for subsequent convolution with the histogram which contains this axis.

Definition at line 76 of file HistoAxis.cc.

                                                                           {
    if (!bw_)
      throw npstat::NpstatDomainError(
          "In npstat::HistoAxis::kernelScanMapper: "
          "bin width is zero. Mapper can not be constructed.");
    double range = max_ - min_;
    if (doubleRange)
      range *= 2.0;
    return CircularMapper1d(bw_, 0.0, range);
  }

References bw_, max_, min_, and FastTimerService_cff::range.

label()

const std::string& npstat::HistoAxis::label ( ) const

inline

Definition at line 47 of file HistoAxis.h.

{ return label_; }

References label_.

Referenced by npstat::convertToGridAxis(), npstat::DualHistoAxis::label(), and read().

leftBinEdge()

double npstat::HistoAxis::leftBinEdge ( const int  binNum ) const

inline

Return the coordinate of the given bin left edge

Definition at line 55 of file HistoAxis.h.

{ return min_ + binNum * bw_; }

References bw_, and min_.

Referenced by npstat::DualHistoAxis::leftBinEdge().

length()

double npstat::HistoAxis::length ( ) const

inline

Definition at line 44 of file HistoAxis.h.

{ return max_ - min_; }

References max_, and min_.

Referenced by npstat::DualHistoAxis::length().

max()

double npstat::HistoAxis::max ( ) const

inline

Definition at line 42 of file HistoAxis.h.

{ return max_; }

References max_.

Referenced by npstat::DualHistoAxis::max(), and read().

min()

double npstat::HistoAxis::min ( ) const

inline

Examine axis properties

Definition at line 41 of file HistoAxis.h.

{ return min_; }

References min_.

Referenced by npstat::DualHistoAxis::min(), and read().

nBins()

unsigned npstat::HistoAxis::nBins ( ) const

inline

Definition at line 45 of file HistoAxis.h.

{ return nBins_; }

References nBins_.

Referenced by npstat::convertToGridAxis(), npstat::DualHistoAxis::nBins(), and read().

operator!=()

bool npstat::HistoAxis::operator!=

(

const HistoAxis &

r

)

const

Definition at line 35 of file HistoAxis.cc.

{ return !(*this == r); }

References alignCSCRings::r.

operator==()

bool npstat::HistoAxis::operator==

(

const HistoAxis &

r

)

const

Definition at line 31 of file HistoAxis.cc.

                                                     {
    return min_ == r.min_ && max_ == r.max_ && label_ == r.label_ && nBins_ == r.nBins_;
  }

References label_, max_, min_, nBins_, and alignCSCRings::r.

overflowIndex()

unsigned npstat::HistoAxis::overflowIndex ( const double  x, unsigned *  binNumber  ) const

inlineprivate

Definition at line 136 of file HistoAxis.h.

                                                                             {
      if (x < min_)
        return 0U;
      else if (x >= max_)
        return 2U;
      else {
        const unsigned bin = static_cast<unsigned>((x - min_) / bw_);
        *binNumber = bin >= nBins_ ? nBins_ - 1U : bin;
        return 1U;
      }
    }

References newFWLiteAna::bin, binNumber(), bw_, max_, min_, nBins_, mitigatedMETSequence_cff::U, and x.

Referenced by npstat::DualHistoAxis::overflowIndex().

overflowIndexWeighted()

unsigned npstat::HistoAxis::overflowIndexWeighted ( double  x, unsigned *  binNumber, double *  weight  ) const

private

Definition at line 87 of file HistoAxis.cc.

                                                                                                     {
    if (x < min_)
      return 0U;
    else if (x >= max_)
      return 2U;
    else {
      if (nBins_ <= 1U)
        throw npstat::NpstatInvalidArgument(
            "In npstat::HistoAxis::overflowIndexWeighted: "
            "must have more than one bin");
      const double dbin = (x - min_) / bw_;
      if (dbin <= 0.5) {
        *binNumber = 0;
        *weight = 1.0;
      } else if (dbin >= nBins_ - 0.5) {
        *binNumber = nBins_ - 2;
        *weight = 0.0;
      } else {
        const unsigned bin = static_cast<unsigned>(dbin - 0.5);
        *binNumber = bin >= nBins_ - 1U ? nBins_ - 2U : bin;
        *weight = 1.0 - (dbin - 0.5 - *binNumber);
      }
      return 1U;
    }
  }

References newFWLiteAna::bin, binNumber(), bw_, max_, min_, nBins_, and mitigatedMETSequence_cff::U.

read()

HistoAxis * npstat::HistoAxis::read ( const gs::ClassId &  id, std::istream &  in  )

static

Definition at line 121 of file HistoAxis.cc.

                                                                {
    static const gs::ClassId current(gs::ClassId::makeId<HistoAxis>());
    current.ensureSameId(id);
 
    double min = 0.0, max = 0.0;
    std::string label;
    unsigned nBins = 0;
 
    gs::read_pod(in, &min);
    gs::read_pod(in, &max);
    gs::read_pod(in, &label);
    gs::read_pod(in, &nBins);
 
    if (!in.fail())
      return new HistoAxis(nBins, min, max, label.c_str());
    else
      throw gs::IOReadFailure(
          "In npstat::HistoAxis::read: "
          "input stream failure");
  }

References HistoAxis(), recoMuon::in, label(), max(), min(), nBins(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by edmIntegrityCheck.PublishToFileSystem::get(), and npstat::DualHistoAxis::read().

rebin()

HistoAxis npstat::HistoAxis::rebin

(

unsigned

newBins

)

const

Return rebinned axis

Definition at line 24 of file HistoAxis.cc.

{ return HistoAxis(nbins, min_, max_, label_.c_str()); }

References HistoAxis(), label_, max_, min_, and LaserClient_cfi::nbins.

rightBinEdge()

double npstat::HistoAxis::rightBinEdge ( const int  binNum ) const

inline

Return the coordinate of the given bin right edge

Definition at line 58 of file HistoAxis.h.

{ return min_ + (binNum + 1) * bw_; }

References bw_, and min_.

Referenced by npstat::DualHistoAxis::rightBinEdge().

setLabel()

void npstat::HistoAxis::setLabel ( const char *  newlabel )

inline

Change the axis label

Definition at line 66 of file HistoAxis.h.

{ label_ = newlabel ? newlabel : ""; }

References label_.

Referenced by npstat::DualHistoAxis::setLabel().

version()

static unsigned npstat::HistoAxis::version ( )

inlinestatic

Definition at line 122 of file HistoAxis.h.

{ return 1; }

Referenced by validation.Sample::datasetpattern(), and validation.Sample::filename().

write()

bool npstat::HistoAxis::write

(

std::ostream &

of

)

const

Definition at line 113 of file HistoAxis.cc.

                                            {
    gs::write_pod(of, min_);
    gs::write_pod(of, max_);
    gs::write_pod(of, label_);
    gs::write_pod(of, nBins_);
    return !of.fail();
  }

References label_, max_, min_, and nBins_.

Referenced by npstat::DualHistoAxis::write().

Friends And Related Function Documentation

DualHistoAxis

friend class DualHistoAxis

friend

Definition at line 134 of file HistoAxis.h.

HistoND

template<typename Numeric , class Axis >

friend class HistoND

friend

Definition at line 133 of file HistoAxis.h.

Member Data Documentation

bw_

double npstat::HistoAxis::bw_

private

Definition at line 128 of file HistoAxis.h.

Referenced by binCenter(), binInterval(), binNumber(), binNumberMapper(), binWidth(), closestValidBin(), fltBinNumber(), HistoAxis(), kernelScanMapper(), leftBinEdge(), overflowIndex(), overflowIndexWeighted(), and rightBinEdge().

label_

std::string npstat::HistoAxis::label_

private

Definition at line 129 of file HistoAxis.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), isClose(), label(), operator==(), rebin(), setLabel(), Mixins._Labelable::setLabel(), and write().

max_

double npstat::HistoAxis::max_

private

Definition at line 127 of file HistoAxis.h.

Referenced by binNumber(), closestValidBin(), HistoAxis(), interval(), isClose(), kernelScanMapper(), length(), max(), operator==(), overflowIndex(), overflowIndexWeighted(), rebin(), and write().

min_

double npstat::HistoAxis::min_

private

Definition at line 126 of file HistoAxis.h.

Referenced by binCenter(), binInterval(), binNumber(), binNumberMapper(), closestValidBin(), fltBinNumber(), HistoAxis(), interval(), isClose(), kernelScanMapper(), leftBinEdge(), length(), min(), operator==(), overflowIndex(), overflowIndexWeighted(), rebin(), rightBinEdge(), and write().

nBins_

unsigned npstat::HistoAxis::nBins_

private

Definition at line 130 of file HistoAxis.h.

Referenced by binNumber(), closestValidBin(), HistoAxis(), isClose(), nBins(), operator==(), overflowIndex(), overflowIndexWeighted(), and write().