GridAxis.h

Go to the documentation of this file.

#ifndef NPSTAT_GRIDAXIS_HH_
#define NPSTAT_GRIDAXIS_HH_

#include <vector>
#include <utility>
#include <string>
#include <iostream>

#include "Alignment/Geners/interface/ClassId.hh"

namespace npstat {
  class GridAxis {
  public:

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


    inline const std::vector<double>& coords() const { return coords_; }
    inline const std::string& label() const { return label_; }
    inline bool usesLogSpace() const { return useLogSpace_; }

    std::pair<unsigned, double> getInterval(double coordinate) const;

    std::pair<unsigned, double> linearInterval(double coordinate) const;


    inline unsigned nCoords() const { return npt_; }
    inline double coordinate(const unsigned i) const { return coords_.at(i); }
    inline double min() const { return coords_.front(); }
    inline double max() const { return coords_.back(); }
    inline double length() const { return coords_.back() - coords_.front(); }
    inline bool isUniform() const { return false; }
    inline unsigned nIntervals() const { return coords_.size() - 1; }
    inline double intervalWidth(const unsigned i = 0) const { return coords_.at(i + 1) - coords_.at(i); }

    bool operator==(const GridAxis& r) const;

    inline bool operator!=(const GridAxis& r) const { return !(*this == r); }

    bool isClose(const GridAxis& r, double tol) const;

    inline void setLabel(const char* newlabel) { label_ = newlabel ? newlabel : ""; }


    inline gs::ClassId classId() const { return gs::ClassId(*this); }
    bool write(std::ostream& of) const;

    static inline const char* classname() { return "npstat::GridAxis"; }
    static inline unsigned version() { return 2; }
    static GridAxis* read(const gs::ClassId& id, std::istream& in);

  private:
    void initialize();

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

    inline GridAxis() : npt_(0), useLogSpace_(false) {}
  };
}  // namespace npstat

#endif  // NPSTAT_GRIDAXIS_HH_