#include <HistoPdf.h>

Inheritance diagram for funct::HistoPdf:

Public Member Functions
void	dump ()

template<typename Iterator >
	HistoPdf (double xMin, double xMax, const Iterator &begin, const Iterator &end)

	HistoPdf ()

template<typename Iterator >
void	init (double xMin, double xMax, const Iterator &begin, const Iterator &end)

double	operator() (double x) const

void	rebin (unsigned int r)

Private Attributes
double	binSize_

double	delta_

double	xMax_

double	xMin_

std::vector< double >	y_

Detailed Description

Definition at line 7 of file HistoPdf.h.

Constructor & Destructor Documentation

template<typename Iterator >

funct::HistoPdf::HistoPdf	(	double	xMin,
		double	xMax,
		const Iterator &	begin,
		const Iterator &	end
	)

inline

Definition at line 10 of file HistoPdf.h.

References end, i, asciidump::s, and y_.

                                                        :
       xMin_(xMin), xMax_(xMax), delta_(xMax - xMin), binSize_(delta_ / (end - begin)), y_(end - begin) {
       double s = 0;
       unsigned int i = 0;
       for(Iterator it = begin; it != end; ++it)
     s += (y_[i++] = *it);
       for(std::vector<double>::iterator i = y_.begin(); i != y_.end(); ++i)
     *i /= s;
     }

funct::HistoPdf::HistoPdf ( )

inline

Definition at line 20 of file HistoPdf.h.

20 { }

Member Function Documentation

void funct::HistoPdf::dump ( void )

inline

Definition at line 54 of file HistoPdf.h.

References binSize_, gather_cfg::cout, delta_, i, operator()(), asciidump::s, ExpressReco_HICollisions_FallBack::x, xMax_, xMin_, ExpressReco_HICollisions_FallBack::y, and y_.

                 {
       std::cout << ">>> range: [" << xMin_ << ", " << xMax_ << "], bin size: " 
        << delta_ << "/" << y_.size() << " = " << binSize_ << std::endl;
        double s = 0;
       for(unsigned int i = 0; i != y_.size(); ++i) {
     double x = xMin_ + (0.5 + i)*binSize_;
     double y = operator()(x);
     std::cout << ">>> pdf(" << x << ") = " << y << std::endl;
     s+= y*binSize_;
       }
      std::cout << ">>>: PDF normalization is " << s << std::endl;
     }

template<typename Iterator >

void funct::HistoPdf::init	(	double	xMin,
		double	xMax,
		const Iterator &	begin,
		const Iterator &	end
	)

inline

Definition at line 22 of file HistoPdf.h.

References begin, binSize_, delta_, end, i, n, asciidump::s, xMax_, xMin_, and y_.

Referenced by funct::RootHistoPdf::RootHistoPdf().

                                                         {
       xMin_ = xMin;
       xMax_ = xMax;
       delta_ = xMax - xMin;
       unsigned int n = end - begin;
       binSize_ = delta_ / n;
       y_.resize(n);
       double s = 0;
       unsigned int i = 0;
       for(Iterator it = begin; it != end; ++it)
     s += (y_[i++] = *it);
       for(std::vector<double>::iterator i = y_.begin(); i != y_.end(); ++i)
     *i /= s;
     }

double funct::HistoPdf::operator() ( double x ) const

inline

Definition at line 37 of file HistoPdf.h.

References binSize_, delta_, xMax_, xMin_, and y_.

Referenced by dump().

                                       {
       if (x < xMin_ || x > xMax_) return 0;
       double pdf = y_[static_cast<unsigned int>(((x -xMin_)/delta_)*y_.size())] / binSize_;
       return pdf;
     }

void funct::HistoPdf::rebin ( unsigned int r )

inline

Definition at line 42 of file HistoPdf.h.

References binSize_, edm::errors::Configuration, i, j, gen::k, n, csvReporter::r, ExpressReco_HICollisions_FallBack::y, and y_.

                                {
       if(y_.size() % r != 0) 
     throw edm::Exception(edm::errors::Configuration) <<
       "HistoPdf: can't rebin histogram of " << y_.size() << " entries by " << r << "\n";
       unsigned int n = y_.size() / r;
       std::vector<double> y(n, 0);
       for(unsigned int i = 0, j = 0; i < n; ++i)
     for(unsigned int k = 0; k < r; ++k) 
       y[i] += y_[j++];
       y_ = y;
       binSize_ *= r;
     }