#include <TrivariateGaussian.hh>

Inheritance diagram for Genfun::TrivariateGaussian:

Public Member Functions
	TrivariateGaussian ()

	TrivariateGaussian (const TrivariateGaussian &right)

virtual	~TrivariateGaussian ()

virtual double	operator() (double argument) const override

virtual double	operator() (const Argument &a) const override

virtual unsigned int	dimensionality () const override

Parameter &	mean0 ()

const Parameter &	mean0 () const

Parameter &	mean1 ()

const Parameter &	mean1 () const

Parameter &	mean2 ()

const Parameter &	mean2 () const

Parameter &	sigma0 ()

const Parameter &	sigma0 () const

Parameter &	sigma1 ()

const Parameter &	sigma1 () const

Parameter &	sigma2 ()

const Parameter &	sigma2 () const

Parameter &	corr01 ()

const Parameter &	corr01 () const

Parameter &	corr02 ()

const Parameter &	corr02 () const

Parameter &	corr12 ()

const Parameter &	corr12 () const

Public Member Functions inherited from Genfun::AbsFunction
	AbsFunction ()

	AbsFunction (const AbsFunction &right)

virtual	~AbsFunction ()

virtual unsigned int	dimensionality () const

virtual double	operator() (double argument) const =0

virtual double	operator() (const Argument &argument) const =0

virtual AbsFunction *	clone () const =0

virtual FunctionComposition	operator() (const AbsFunction &f) const

virtual ParameterComposition	operator() (const AbsParameter &p) const

Derivative	derivative (const Variable &v) const

Derivative	prime () const

virtual bool	hasAnalyticDerivative () const

virtual Derivative	partial (unsigned int) const

Detailed Description

Author

Definition at line 20 of file TrivariateGaussian.hh.

Constructor & Destructor Documentation

◆ TrivariateGaussian() [1/2]

Genfun::TrivariateGaussian::TrivariateGaussian ( )

Definition at line 20 of file TrivariateGaussian.cc.

                                      :
  _mean0("Mean0", 0.0,-10,10),
  _mean1("Mean1", 0.0,-10,10),
  _mean2("Mean2", 0.0,-10,10),
  _sigma0("Sigma0",1.0,0, 10),
  _sigma1("Sigma1",1.0,0, 10),
  _sigma2("Sigma2",1.0,0, 10),
  _corr01("Corr01", 0.0, -1.0, 1.0),
  _corr02("Corr02", 0.0, -1.0, 1.0),
  _corr12("Corr12", 0.0, -1.0, 1.0)
{}

◆ TrivariateGaussian() [2/2]

Genfun::TrivariateGaussian::TrivariateGaussian ( const TrivariateGaussian & right )

Definition at line 35 of file TrivariateGaussian.cc.

                                                                      :
  AbsFunction(right),
  _mean0(right._mean0),
  _mean1(right._mean1),
  _mean2(right._mean2),
  _sigma0(right._sigma0),
  _sigma1(right._sigma1),
  _sigma2(right._sigma2),
  _corr01(right._corr01),
  _corr02(right._corr02),
  _corr12(right._corr12)
{
}

◆ ~TrivariateGaussian()

Genfun::TrivariateGaussian::~TrivariateGaussian ( )

virtual

Definition at line 32 of file TrivariateGaussian.cc.

32 {

33}

Member Function Documentation

◆ corr01() [1/2]

Parameter & Genfun::TrivariateGaussian::corr01 ( )

Definition at line 142 of file TrivariateGaussian.cc.

                                       {
  return _corr01;
}

◆ corr01() [2/2]

const Parameter & Genfun::TrivariateGaussian::corr01 ( ) const

Definition at line 146 of file TrivariateGaussian.cc.

                                                   {
  return _corr01;
}

◆ corr02() [1/2]

Parameter & Genfun::TrivariateGaussian::corr02 ( )

Definition at line 150 of file TrivariateGaussian.cc.

                                       {
  return _corr02;
}

◆ corr02() [2/2]

const Parameter & Genfun::TrivariateGaussian::corr02 ( ) const

Definition at line 154 of file TrivariateGaussian.cc.

                                                   {
  return _corr02;
}

◆ corr12() [1/2]

Parameter & Genfun::TrivariateGaussian::corr12 ( )

Definition at line 158 of file TrivariateGaussian.cc.

                                       {
  return _corr12;
}

◆ corr12() [2/2]

const Parameter & Genfun::TrivariateGaussian::corr12 ( ) const

Definition at line 162 of file TrivariateGaussian.cc.

                                                   {
  return _corr12;
}

◆ dimensionality()

unsigned int Genfun::TrivariateGaussian::dimensionality ( ) const

overridevirtual

Reimplemented from Genfun::AbsFunction.

Definition at line 167 of file TrivariateGaussian.cc.

                                                      {
  return 3;
}

◆ mean0() [1/2]

Parameter & Genfun::TrivariateGaussian::mean0 ( )

Definition at line 91 of file TrivariateGaussian.cc.

                                      {
  return _mean0;
}

◆ mean0() [2/2]

const Parameter & Genfun::TrivariateGaussian::mean0 ( ) const

Definition at line 99 of file TrivariateGaussian.cc.

                                                  {
  return _mean0;
}

◆ mean1() [1/2]

Parameter & Genfun::TrivariateGaussian::mean1 ( )

Definition at line 107 of file TrivariateGaussian.cc.

                                      {
  return _mean1;
}

◆ mean1() [2/2]

const Parameter & Genfun::TrivariateGaussian::mean1 ( ) const

Definition at line 115 of file TrivariateGaussian.cc.

                                                  {
  return _mean1;
}

◆ mean2() [1/2]

Parameter & Genfun::TrivariateGaussian::mean2 ( )

Definition at line 123 of file TrivariateGaussian.cc.

                                      {
  return _mean2;
}

◆ mean2() [2/2]

const Parameter & Genfun::TrivariateGaussian::mean2 ( ) const

Definition at line 132 of file TrivariateGaussian.cc.

                                                  {
  return _mean2;
}

◆ operator()() [1/2]

double Genfun::TrivariateGaussian::operator() ( const Argument & a ) const

overridevirtual

Implements Genfun::AbsFunction.

Definition at line 50 of file TrivariateGaussian.cc.

                                                               {
  assert (a.dimension()==3);
  double x = a[0];
  double y = a[1];
  double z = a[2];
 
 
  double x0   = _mean0.getValue();
  double y0   = _mean1.getValue();
  double z0   = _mean2.getValue();
 
  double dx   = x-x0;
  double dy   = y-y0;
  double dz   = z-z0;
 
  double sx     = _sigma0.getValue();
  double sy     = _sigma1.getValue();
  double sz     = _sigma2.getValue();
 
 
  double sxs    = sx*sx;
  double sys    = sy*sy;
  double szs    = sz*sz;
 
 
  double rho1    = _corr01.getValue();
  double rho2    = _corr12.getValue();
  double rho3    = _corr02.getValue();
 
  double dt     = (1.0+rho1*rho2*rho3-rho1*rho1-rho2*rho2-rho3*rho3);
  double tmp1 ,tmp2;
 
  tmp1= 1.0/((2*M_PI)*sqrt(2*M_PI)*sx*sy*sz*sqrt(dt));
  tmp2= exp(-0.5/dt*(dx*dx*(1.0-rho2*rho2)/sxs+dy*dy*(1.0-rho3*rho3)/sys+dz*dz*(1.0-rho1*rho1)/szs+2.0*dx*dy*(rho2*rho3-rho1)/sx/sy+2.0*dy*dz*(rho1*rho3-rho2)/sy/sz+2.0*dx*dz*(rho1*rho2-rho3)/sx/sz));
  
 
  return tmp1*tmp2; 
  
 
}

◆ operator()() [2/2]

double Genfun::TrivariateGaussian::operator() ( double argument ) const

overridevirtual

Implements Genfun::AbsFunction.

Definition at line 171 of file TrivariateGaussian.cc.

{
  std::cerr
    << "Warning.  trivariate Gaussian called with scalar argument"
    << std::endl;
  assert(0);
  return 0;
}

◆ sigma0() [1/2]

Parameter & Genfun::TrivariateGaussian::sigma0 ( )

Definition at line 95 of file TrivariateGaussian.cc.

                                       {
  return _sigma0;
}

◆ sigma0() [2/2]

const Parameter & Genfun::TrivariateGaussian::sigma0 ( ) const

Definition at line 103 of file TrivariateGaussian.cc.

                                                   {
  return _sigma0;
}

◆ sigma1() [1/2]

Parameter & Genfun::TrivariateGaussian::sigma1 ( )

Definition at line 111 of file TrivariateGaussian.cc.

                                       {
  return _sigma1;
}

◆ sigma1() [2/2]

const Parameter & Genfun::TrivariateGaussian::sigma1 ( ) const

Definition at line 119 of file TrivariateGaussian.cc.

                                                   {
  return _sigma1;
}

◆ sigma2() [1/2]

Parameter & Genfun::TrivariateGaussian::sigma2 ( )

Definition at line 128 of file TrivariateGaussian.cc.

                                       {
  return _sigma2;
}

◆ sigma2() [2/2]

const Parameter & Genfun::TrivariateGaussian::sigma2 ( ) const

Definition at line 136 of file TrivariateGaussian.cc.

                                                   {
  return _sigma2;
}

The documentation for this class was generated from the following files:

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ TrivariateGaussian() [1/2]

◆ TrivariateGaussian() [2/2]

◆ ~TrivariateGaussian()

Member Function Documentation

◆ corr01() [1/2]

◆ corr01() [2/2]

◆ corr02() [1/2]

◆ corr02() [2/2]

◆ corr12() [1/2]

◆ corr12() [2/2]

◆ dimensionality()

◆ mean0() [1/2]

◆ mean0() [2/2]

◆ mean1() [1/2]

◆ mean1() [2/2]

◆ mean2() [1/2]

◆ mean2() [2/2]

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ sigma0() [1/2]

◆ sigma0() [2/2]

◆ sigma1() [1/2]

◆ sigma1() [2/2]

◆ sigma2() [1/2]

◆ sigma2() [2/2]