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
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)
bool	isClose (const HistoAxis &, double tol) const
CircularMapper1d	kernelScanMapper (bool doubleRange) const
double	leftBinEdge (const int binNum) 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)
double	min () const
double	max () const
Interval< double >	interval () const
double	length () const
unsigned	nBins () const
double	binWidth (const int=0) const
const std::string &	label () const
bool	isUniform () const
gs::ClassId	classId () const
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.

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

      : 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_;
  }

HistoAxis() [2/2]

npstat::HistoAxis::HistoAxis ( )

inlineprivate

Definition at line 149 of file HistoAxis.h.

Referenced by read(), and rebin().

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

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.

References bw_, and min_.

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

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

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.

References bw_, and min_.

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

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

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.

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

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

                                               {
    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_;
    }
  }

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.

References edmMakeDummyCfis::base, bw_, and min_.

                                                                           {
    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);
  }

binWidth()

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

inline

Definition at line 46 of file HistoAxis.h.

References bw_.

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

{ return bw_; }

classId()

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

inline

Method related to "geners" I/O

Definition at line 117 of file HistoAxis.h.

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

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

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.

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

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

                                                          {
    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;
  }

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.

References bw_, min_, and x.

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

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

interval()

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

inline

Definition at line 43 of file HistoAxis.h.

References max_, and min_.

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

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

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.

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

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

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

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.

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

                                                                           {
    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);
  }

label()

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

inline

Definition at line 47 of file HistoAxis.h.

References label_.

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

{ return label_; }

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.

References bw_, and min_.

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

{ return min_ + binNum * bw_; }

length()

double npstat::HistoAxis::length ( ) const

inline

Definition at line 44 of file HistoAxis.h.

References max_, and min_.

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

{ return max_ - min_; }

max()

double npstat::HistoAxis::max ( ) const

inline

Definition at line 42 of file HistoAxis.h.

References max_.

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

{ return max_; }

min()

double npstat::HistoAxis::min ( ) const

inline

Examine axis properties

Definition at line 41 of file HistoAxis.h.

References min_.

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

{ return min_; }

nBins()

unsigned npstat::HistoAxis::nBins ( ) const

inline

Definition at line 45 of file HistoAxis.h.

References nBins_.

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

{ return nBins_; }

operator!=()

bool npstat::HistoAxis::operator!=

(

const HistoAxis &

r

)

const

Definition at line 35 of file HistoAxis.cc.

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

operator==()

bool npstat::HistoAxis::operator==

(

const HistoAxis &

r

)

const

Definition at line 31 of file HistoAxis.cc.

References label_, max_, min_, and nBins_.

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

overflowIndex()

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

inlineprivate

Definition at line 136 of file HistoAxis.h.

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

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

                                                                             {
      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;
      }
    }

overflowIndexWeighted()

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

private

Definition at line 87 of file HistoAxis.cc.

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

                                                                                                     {
    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;
    }
  }

read()

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

static

Definition at line 121 of file HistoAxis.cc.

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

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

                                                                {
    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");
  }

rebin()

HistoAxis npstat::HistoAxis::rebin

(

unsigned

newBins

)

const

Return rebinned axis

Definition at line 24 of file HistoAxis.cc.

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

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

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.

References bw_, and min_.

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

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

setLabel()

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

inline

Change the axis label

Definition at line 66 of file HistoAxis.h.

References label_.

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

{ label_ = newlabel ? newlabel : ""; }

version()

static unsigned npstat::HistoAxis::version ( )

inlinestatic

Definition at line 122 of file HistoAxis.h.

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

{ return 1; }

write()

bool npstat::HistoAxis::write

(

std::ostream &

of

)

const

Definition at line 113 of file HistoAxis.cc.

References label_, max_, min_, nBins_, and haddnano::of.

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

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

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 Modules.SwitchProducer::__addParameter(), SequenceTypes.DummyModule::__repr__(), Modules.SwitchProducer::__setattr__(), 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().