/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/JetMETCorrections/InterpolationTables/interface/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:

        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:
        inline GridAxis() : npt_(0), useLogSpace_(false) {}

        void initialize();

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

    };
}

#endif // NPSTAT_GRIDAXIS_HH_