#include <SigmaPtDiff.h>

Public Member Functions
double	operator() (const double &eta)
void	setParErr (const std::vector< double > &parameters, const std::vector< double > &errors)
double	sigma (const double &eta)
	SigmaPt (const std::vector< double > &parameters_, const std::vector< double > &errors_)
	SigmaPt ()
Protected Attributes
double	b_0
double	b_1
double	b_2
double	b_3
double	b_4
double	c
double	sb_0
double	sb_1
double	sb_2
double	sb_3
double	sb_4

Detailed Description

Definition at line 4 of file SigmaPtDiff.h.

Constructor & Destructor Documentation

SigmaPt::SigmaPt	(	const std::vector< double > &	parameters_,
		const std::vector< double > &	errors_
	)		`[inline]`

Definition at line 7 of file SigmaPtDiff.h.

References setParErr().

    {
      setParErr(parameters_, errors_);
    }

SigmaPt::SigmaPt ( ) [inline]

Definition at line 13 of file SigmaPtDiff.h.

{};

Member Function Documentation

double SigmaPt::operator() ( const double & eta ) [inline]

Definition at line 31 of file SigmaPtDiff.h.

References b_0, b_1, b_2, b_3, b_4, and c.

  {
    if( fabs(eta) <= b_0 ) {
      return( c + b_1*eta*eta );
    }
    return( b_2 + b_3*(fabs(eta) - b_4)*(fabs(eta) - b_4) );
  }

void SigmaPt::setParErr	(	const std::vector< double > &	parameters,
		const std::vector< double > &	errors
	)		`[inline]`

Definition at line 15 of file SigmaPtDiff.h.

References b_0, b_1, b_2, b_3, b_4, c, sb_0, sb_1, sb_2, sb_3, and sb_4.

Referenced by SigmaPt(), and SigmaPtDiff::SigmaPtDiff().

  {
    b_0 = parameters[0];
    b_1 = parameters[5];
    b_2 = parameters[1];
    b_3 = parameters[7];
    b_4 = parameters[8];
    sb_0 = errors[0];
    sb_1 = errors[5];
    sb_2 = errors[1];
    sb_3 = errors[7];
    sb_4 = errors[8];
    c = b_2 + b_3*(b_0 - b_4)*(b_0 - b_4) - b_1*b_0*b_0;
  }

double SigmaPt::sigma ( const double & eta ) [inline]

Definition at line 38 of file SigmaPtDiff.h.

References b_0, b_1, b_3, b_4, funct::pow(), sb_0, sb_1, sb_2, sb_3, sb_4, and mathSSE::sqrt().

Referenced by SigmaPtDiff::squaredDiff().

  {
    if( fabs(eta) <= b_0 ) {
      return sqrt( (eta*eta - b_0*b_0)*(eta*eta - b_0*b_0)*sb_1*sb_1 +
                   sb_2*sb_2 +
                   pow(b_0 - b_4, 4)*sb_3*sb_3 +
                   pow(-2*b_3*pow(b_0-b_4,2), 2)*sb_4*sb_4 +
                   pow(2*b_3*(b_0 - b_4) - 2*b_1*b_0, 2)*sb_0*sb_0 );
    }
    return sqrt( sb_2*sb_2 +
                 pow(fabs(eta) - b_4, 4)*sb_3*sb_3 +
                 pow(-2*b_3*pow(fabs(eta)-b_4,2), 2)*sb_4*sb_4 );
  }