Gaussian.cc

Go to the documentation of this file.

// -*- C++ -*-
// $Id: Gaussian.cc,v 1.8 2010/06/16 18:22:01 garren Exp $
#include "CLHEP/GenericFunctions/defs.h"
#include "CLHEP/GenericFunctions/Gaussian.hh"
#include "CLHEP/GenericFunctions/Variable.hh"
#include <assert.h>
#include <cmath>      // for exp()
#include <iostream>
 
#if (defined __STRICT_ANSI__) || (defined _WIN32)
#ifndef M_PI
#define M_PI            3.14159265358979323846
#endif // M_PI
#endif // __STRICT_ANSI__
 
namespace Genfun {
FUNCTION_OBJECT_IMP(Gaussian)
 
Gaussian::Gaussian():
  _mean("Mean", 0.0,-10,10),
  _sigma("Sigma",1.0,0, 10)
{}
 
Gaussian::~Gaussian() {
}
 
Gaussian::Gaussian(const Gaussian & right):
AbsFunction(right),
_mean(right._mean),
_sigma(right._sigma)
{
}
 
double Gaussian::operator() (double x) const {
  double s   = _sigma.getValue();
  double x0  = _mean.getValue();
  return (1.0/(sqrt(2*M_PI)*s))*
      exp(-(x-x0)*(x-x0)/(2.0*s*s));
}
 
Parameter & Gaussian::mean() {
  return _mean;
}
 
Parameter & Gaussian::sigma() {
  return _sigma;
}
 
const Parameter & Gaussian::mean() const {
  return _mean;
}
 
const Parameter & Gaussian::sigma() const {
  return _sigma;
}
 
// don't generate warnings about unused parameter inside assert
#if defined __GNUC__ 
  #if __GNUC__ > 3 && __GNUC_MINOR__ > 6
    #pragma GCC diagnostic push
    #pragma GCC diagnostic ignored "-Wunused-parameter"
  #endif
#endif
#ifdef __clang__
  #pragma clang diagnostic push
  #pragma clang diagnostic ignored "-Wunused-parameter"
#endif
Derivative Gaussian::partial(unsigned int index) const {
  assert(index==0);
  Variable x;
  const AbsFunction & fPrime  = (*this)*(_mean-x)/_sigma/_sigma;
  return Derivative(&fPrime);
}
#if defined __GNUC__ 
  #if __GNUC__ > 3 && __GNUC_MINOR__ > 6
    #pragma GCC diagnostic pop
  #endif
#endif
#ifdef __clang__
  #pragma clang diagnostic pop
#endif
 
} // namespace Genfun