Gaussian.h

Go to the documentation of this file.

#ifndef PhysicsTools_Utilities_Gaussian_h
#define PhysicsTools_Utilities_Gaussian_h
#include "PhysicsTools/Utilities/interface/Parameter.h"
#include <boost/shared_ptr.hpp>
#include <cmath>

namespace funct {

  const double oneOverSqrtTwoPi = 1/sqrt(2*M_PI);
  
  struct Gaussian {
    Gaussian(const Parameter & m, const Parameter & s) : 
      mean(m.ptr()), sigma(s.ptr()) { }
    Gaussian(boost::shared_ptr<double> m, boost::shared_ptr<double> s): 
      mean(m), sigma(s) { }
    Gaussian(double m, double s): 
      mean(new double(m)), sigma(new double(s)){}
    double operator()(double x) const {
      double z = (x - *mean)/ *sigma;
      if(fabs(z)>8) return 0;
      return oneOverSqrtTwoPi/ *sigma * exp(-z*z/2);
    }
    boost::shared_ptr<double> mean, sigma;
  };

}

#endif