npstat::GridAxis Class Reference

#include <GridAxis.h>

Public Member Functions
std::pair< unsigned, double >	getInterval (double coordinate) const
bool	isClose (const GridAxis &r, double tol) const
std::pair< unsigned, double >	linearInterval (double coordinate) const
bool	operator!= (const GridAxis &r) const
bool	operator== (const GridAxis &r) const
void	setLabel (const char *newlabel)
	GridAxis (const std::vector< double > &coords, bool useLogSpace=false)
	GridAxis (const std::vector< double > &coords, const char *label, bool useLogSpace=false)
const std::vector< double > &	coords () const
const std::string &	label () const
bool	usesLogSpace () const
unsigned	nCoords () const
double	coordinate (const unsigned i) const
double	min () const
double	max () const
double	length () const
bool	isUniform () const
unsigned	nIntervals () const
double	intervalWidth (const unsigned i=0) const
gs::ClassId	classId () const
bool	write (std::ostream &of) const

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

Private Member Functions
	GridAxis ()
void	initialize ()

Private Attributes
std::vector< double >	coords_
std::string	label_
std::vector< double >	logs_
unsigned	npt_
bool	useLogSpace_

Detailed Description

Information needed to define an axis of a rectangular grid. The distance between grid points can change from point to point.

The UniformAxis class will be more efficient in representing equidistant grids.

Definition at line 30 of file GridAxis.h.

Constructor & Destructor Documentation

GridAxis() [1/3]

npstat::GridAxis::GridAxis ( const std::vector< double > &  coords, bool  useLogSpace = false  )

explicit

The number of grid coordinates provided must be at least 2. Coordinates will be sorted internally in the increasing order.

Definition at line 38 of file GridAxis.cc.

References initialize().

      : coords_(coords), npt_(coords_.size()), useLogSpace_(useLogSpace) {
    initialize();
  }

GridAxis() [2/3]

npstat::GridAxis::GridAxis	(	const std::vector< double > &	coords,
		const char *	label,
		bool	useLogSpace = false
	)

Definition at line 43 of file GridAxis.cc.

References initialize().

      : coords_(coords), label_(label ? label : ""), npt_(coords_.size()), useLogSpace_(useLogSpace) {
    initialize();
  }

GridAxis() [3/3]

npstat::GridAxis::GridAxis ( )

inlineprivate

Definition at line 125 of file GridAxis.h.

Referenced by read().

: npt_(0), useLogSpace_(false) {}

Member Function Documentation

classId()

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

inline

Method related to "geners" I/O

Definition at line 108 of file GridAxis.h.

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

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

classname()

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

inlinestatic

Definition at line 112 of file GridAxis.h.

{ return "npstat::GridAxis"; }

coordinate()

double npstat::GridAxis::coordinate ( const unsigned  i ) const

inline

Definition at line 82 of file GridAxis.h.

References coords_, and mps_fire::i.

Referenced by npstat::convertToHistoAxis(), and npstat::DualAxis::coordinate().

{ return coords_.at(i); }

coords()

const std::vector<double>& npstat::GridAxis::coords ( ) const

inline

Basic accessor returning a parameter provided in the constructor

Definition at line 43 of file GridAxis.h.

References coords_.

Referenced by npstat::DualAxis::coords(), and read().

{ return coords_; }

getInterval()

std::pair< unsigned, double > npstat::GridAxis::getInterval

(

double

coordinate

)

const

This method returns the grid interval number and the weight of the point at the left side of the interval. The weight will be set to 1 if the given coordinate coincides with the grid point and will decay to 0 linearly as the coordinate moves towards the next point on the right.

The coordinates below the leftmost grid point are mapped into the 0th interval with weight 1. The coordinates above the rightmost grid point are mapped into the last interval with weight 0 for the left point (it is expected that weight 1 will then be assigned to the right point).

Definition at line 48 of file GridAxis.cc.

References c, coords_, cmsLHEtoEOSManager::l, dqm-mbProfile::log, logs_, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, npt_, mitigatedMETSequence_cff::U, useLogSpace_, w(), and x.

Referenced by npstat::DualAxis::getInterval().

                                                                      {
    if (useLogSpace_)
      if (x <= 0.0)
        throw npstat::NpstatDomainError("In GridAxis::getInterval: argument must be positive");
    const double* c = &coords_[0];
    if (x <= c[0])
      return std::pair<unsigned, double>(0U, 1.0);
    else if (x >= c[npt_ - 1U])
      return std::pair<unsigned, double>(npt_ - 2U, 0.0);
    else {
      unsigned ibnd = lower_bound(coords_.begin(), coords_.end(), x) - coords_.begin();
      --ibnd;
      if (useLogSpace_) {
        const double* l = &logs_[0];
        const double w = 1.0 - (log(x) - l[ibnd]) / (l[ibnd + 1U] - l[ibnd]);
        return std::pair<unsigned, double>(ibnd, w);
      } else {
        const double w = 1.0 - (x - c[ibnd]) / (c[ibnd + 1U] - c[ibnd]);
        return std::pair<unsigned, double>(ibnd, w);
      }
    }
  }

initialize()

void npstat::GridAxis::initialize ( void  )

private

Definition at line 12 of file GridAxis.cc.

References c, coords_, mps_fire::i, dqm-mbProfile::log, logs_, npt_, jetsAK4_CHS_cff::sort, mitigatedMETSequence_cff::U, and useLogSpace_.

Referenced by GridAxis().

                            {
    if (npt_ <= 1U)
      throw npstat::NpstatInvalidArgument(
          "In npstat::GridAxis::initialize: insufficient number "
          "of coordinates (at least 2 are required)");
    std::sort(coords_.begin(), coords_.end());
    const double* c = &coords_[0];
    if (useLogSpace_) {
      logs_.reserve(npt_);
      for (unsigned i = 0; i < npt_; ++i) {
        if (c[i] <= 0.0)
          throw npstat::NpstatInvalidArgument(
              "In npstat::GridAxis::initialize: in log space"
              "all coordinates must be positive");
        logs_.push_back(log(c[i]));
      }
    }

    // Can't have duplicate coordinates
    for (unsigned i = 1; i < npt_; ++i)
      if (c[i] <= c[i - 1U])
        throw npstat::NpstatInvalidArgument(
            "In npstat::GridAxis::initialize: "
            "all coordinates must be distinct");
  }

intervalWidth()

double npstat::GridAxis::intervalWidth ( const unsigned  i = 0 ) const

inline

Definition at line 88 of file GridAxis.h.

References coords_, and mps_fire::i.

Referenced by npstat::DualAxis::intervalWidth().

{ return coords_.at(i + 1) - coords_.at(i); }

isClose()

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

Check for closeness of coordinates with another axis within the given relative tolerance

Definition at line 141 of file GridAxis.cc.

References npstat::closeWithinTolerance(), coords_, mps_fire::i, label_, dqmiodumpmetadata::n, and useLogSpace_.

                                                                  {
    if (!(useLogSpace_ == r.useLogSpace_ && label_ == r.label_))
      return false;
    const unsigned long n = coords_.size();
    if (n != r.coords_.size())
      return false;
    for (unsigned long i = 0; i < n; ++i)
      if (!closeWithinTolerance(coords_[i], r.coords_[i], tol))
        return false;
    return true;
  }

isUniform()

bool npstat::GridAxis::isUniform ( ) const

inline

Definition at line 86 of file GridAxis.h.

{ return false; }

label()

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

inline

Definition at line 44 of file GridAxis.h.

References label_.

Referenced by npstat::convertToHistoAxis(), npstat::DualAxis::label(), and read().

{ return label_; }

length()

double npstat::GridAxis::length ( ) const

inline

Definition at line 85 of file GridAxis.h.

References coords_.

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

{ return coords_.back() - coords_.front(); }

linearInterval()

std::pair< unsigned, double > npstat::GridAxis::linearInterval

(

double

coordinate

)

const

This method returns the grid interval number and the weight of the point at the left side of the interval. The weight will be set to 1 if the given coordinate coincides with the grid point and will decay to 0 linearly as the coordinate moves towards the next point on the right. The weight for the point on the right should be set to one minus the weight on the left.

The coordinates outside of grid boundaries will result in weights which are less than zero or more than one. They will be calculated by linear extrapolation from the closest interval in the grid (i.e., leftmost or rightmost).

Definition at line 71 of file GridAxis.cc.

References c, coords_, cmsLHEtoEOSManager::l, dqm-mbProfile::log, logs_, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, npt_, mitigatedMETSequence_cff::U, useLogSpace_, w(), and x.

Referenced by npstat::DualAxis::linearInterval().

                                                                         {
    if (useLogSpace_)
      if (x <= 0.0)
        throw npstat::NpstatDomainError("In GridAxis::linearInterval: argument must be positive");
    const double* c = &coords_[0];
    if (x <= c[0]) {
      if (useLogSpace_) {
        const double* l = &logs_[0];
        const double bw = l[1] - l[0];
        return std::pair<unsigned, double>(0U, 1.0 - (log(x) - l[0]) / bw);
      } else {
        const double bw = c[1] - c[0];
        return std::pair<unsigned, double>(0U, 1.0 - (x - c[0]) / bw);
      }
    } else if (x >= c[npt_ - 1U]) {
      if (useLogSpace_) {
        const double* l = &logs_[0];
        const double bw = l[npt_ - 1U] - l[npt_ - 2U];
        return std::pair<unsigned, double>(npt_ - 2U, (l[npt_ - 1U] - log(x)) / bw);
      } else {
        const double bw = c[npt_ - 1U] - c[npt_ - 2U];
        return std::pair<unsigned, double>(npt_ - 2U, (c[npt_ - 1U] - x) / bw);
      }
    } else {
      unsigned ibnd = lower_bound(coords_.begin(), coords_.end(), x) - coords_.begin();
      --ibnd;
      if (useLogSpace_) {
        const double* l = &logs_[0];
        const double w = 1.0 - (log(x) - l[ibnd]) / (l[ibnd + 1U] - l[ibnd]);
        return std::pair<unsigned, double>(ibnd, w);
      } else {
        const double w = 1.0 - (x - c[ibnd]) / (c[ibnd + 1U] - c[ibnd]);
        return std::pair<unsigned, double>(ibnd, w);
      }
    }
  }

max()

double npstat::GridAxis::max ( ) const

inline

Definition at line 84 of file GridAxis.h.

References coords_.

Referenced by npstat::DualAxis::max().

{ return coords_.back(); }

min()

double npstat::GridAxis::min ( ) const

inline

Definition at line 83 of file GridAxis.h.

References coords_.

Referenced by npstat::DualAxis::min().

{ return coords_.front(); }

nCoords()

unsigned npstat::GridAxis::nCoords ( ) const

inline

Convenience accessor

Definition at line 81 of file GridAxis.h.

References npt_.

Referenced by npstat::convertToHistoAxis(), and npstat::DualAxis::nCoords().

{ return npt_; }

nIntervals()

unsigned npstat::GridAxis::nIntervals ( ) const

inline

Definition at line 87 of file GridAxis.h.

References coords_.

Referenced by npstat::DualAxis::nIntervals().

{ return coords_.size() - 1; }

operator!=()

bool npstat::GridAxis::operator!= ( const GridAxis &  r ) const

inline

Logical negation of operator==

Definition at line 95 of file GridAxis.h.

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

operator==()

bool npstat::GridAxis::operator==

(

const GridAxis &

r

)

const

Compare two grids for equality

Definition at line 137 of file GridAxis.cc.

References coords_, label_, and useLogSpace_.

                                                   {
    return useLogSpace_ == r.useLogSpace_ && coords_ == r.coords_ && label_ == r.label_;
  }

read()

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

static

Definition at line 117 of file GridAxis.cc.

References coords(), GridAxis(), recoMuon::in, label(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

                                                              {
    static const gs::ClassId current(gs::ClassId::makeId<GridAxis>());
    current.ensureSameId(id);

    std::vector<double> coords;
    gs::read_pod_vector(in, &coords);

    std::string label;
    gs::read_pod(in, &label);

    bool useLogSpace;
    gs::read_pod(in, &useLogSpace);

    if (in.fail())
      throw gs::IOReadFailure(
          "In npstat::GridAxis::read: "
          "input stream failure");
    return new GridAxis(coords, label.c_str(), useLogSpace);
  }

setLabel()

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

inline

Modify the axis label

Definition at line 104 of file GridAxis.h.

References label_.

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

{ label_ = newlabel ? newlabel : ""; }

usesLogSpace()

bool npstat::GridAxis::usesLogSpace ( ) const

inline

Definition at line 45 of file GridAxis.h.

References useLogSpace_.

Referenced by npstat::DualAxis::usesLogSpace().

{ return useLogSpace_; }

version()

static unsigned npstat::GridAxis::version ( )

inlinestatic

Definition at line 113 of file GridAxis.h.

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

{ return 2; }

write()

bool npstat::GridAxis::write

(

std::ostream &

of

)

const

Definition at line 108 of file GridAxis.cc.

References coords_, label_, haddnano::of, and useLogSpace_.

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

                                           {
    // It is unlikely that this object will be written in isolation.
    // So, do not bother with too many checks.
    gs::write_pod_vector(of, coords_);
    gs::write_pod(of, label_);
    gs::write_pod(of, useLogSpace_);
    return !of.fail();
  }

Member Data Documentation

coords_

std::vector<double> npstat::GridAxis::coords_

private

Definition at line 119 of file GridAxis.h.

Referenced by coordinate(), coords(), getInterval(), initialize(), intervalWidth(), isClose(), length(), linearInterval(), max(), min(), nIntervals(), operator==(), and write().

label_

std::string npstat::GridAxis::label_

private

Definition at line 121 of file GridAxis.h.

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

logs_

std::vector<double> npstat::GridAxis::logs_

private

Definition at line 120 of file GridAxis.h.

Referenced by getInterval(), initialize(), and linearInterval().

npt_

unsigned npstat::GridAxis::npt_

private

Definition at line 122 of file GridAxis.h.

Referenced by getInterval(), initialize(), linearInterval(), and nCoords().

useLogSpace_

bool npstat::GridAxis::useLogSpace_

private

Definition at line 123 of file GridAxis.h.

Referenced by getInterval(), initialize(), isClose(), linearInterval(), operator==(), usesLogSpace(), and write().