#include <ChisqConsistency.h>

Inheritance diagram for ChisqConsistency:

Public Member Functions
	ChisqConsistency ()

	ChisqConsistency (double chisq, double nDof)

	ChisqConsistency (unsigned nDof, double consistency)

	ChisqConsistency (const ChisqConsistency &)

ChisqConsistency &	operator= (const ChisqConsistency &)

virtual	~ChisqConsistency ()

const double &	chisqValue () const

const double &	nDOF () const

Public Member Functions inherited from Consistency
	Consistency ()

	Consistency (double consistency, double likelihood=0.)

	Consistency (const Consistency &rhs)

virtual	~Consistency ()

Consistency &	operator= (const Consistency &rhs)

bool	operator== (const Consistency &rhs) const

bool	operator< (const Consistency &rhs) const

bool	operator> (const Consistency &rhs) const

double	significanceLevel () const

double	likelihood () const

double	consistency () const

ConsistentStatus	status () const

void	setStatus (ConsistentStatus s)

ConsistentSign	sign () const

void	setSign (ConsistentSign s)

virtual const ConsistencySet *	genealogy () const

virtual void	print (std::ostream &) const

Protected Attributes
double	_chisq

double	_nDof

Protected Attributes inherited from Consistency
ConsistentStatus	_stat

ConsistentSign	_sign

double	_value

double	_likelihood

Additional Inherited Members
Public Types inherited from Consistency
enum	ConsistentStatus { OK =0 , noMeasure , underFlow , unPhysical }

enum	ConsistentSign { left =-1 , unknown =0 , right =1 }

Static Public Member Functions inherited from Consistency
static const Consistency &	badMeasurement ()

Detailed Description

Definition at line 27 of file ChisqConsistency.h.

Constructor & Destructor Documentation

◆ ChisqConsistency() [1/4]

ChisqConsistency::ChisqConsistency ( )

Definition at line 38 of file ChisqConsistency.cxx.

                                   :
  _chisq(-1.0), _nDof(0)
{}

◆ ChisqConsistency() [2/4]

ChisqConsistency::ChisqConsistency	(	double	chisq,
		double	nDof
	)

Definition at line 42 of file ChisqConsistency.cxx.

                                                            :
   _chisq(chisq), _nDof(nDof)
{
  double z2 = 0.5*_chisq;
  double n2 = 0.5*_nDof;
  
  if (n2<=0 || z2<0) { 
    std::cout << "ErrMsg(warning)" << " Got unphysical values: chisq = " << chisq 
            << " #dof = " << nDof << std::endl;
    _value=0;
    _likelihood=0;
    setStatus(Consistency::unPhysical);
    return;
  }
  setStatus(Consistency::OK);
  
// given that n2>0 && z2>=0, gammq will NOT abort
  _value = NumRecipes::gammq(n2,z2);
  
  if (_chisq==0) { 
    _likelihood=1; 
  } else {
    double loglike=(n2-1)*log(z2)-z2-NumRecipes::gammln(n2);
    if ( loglike < DBL_MIN_EXP ) {
      _likelihood = 0; 
      setStatus(Consistency::underFlow);
    } else {
      _likelihood = 0.5*exp(loglike);
    }
  }
}

◆ ChisqConsistency() [3/4]

ChisqConsistency::ChisqConsistency	(	unsigned	nDof,
		double	consistency
	)

Definition at line 76 of file ChisqConsistency.cxx.

                                                             :
  _nDof(nDof)
{
  if(prob >= 0.0|| prob <= 1.0 || nDof < 0)
    _value = prob;
  else {
    std::cout << "ErrMsg(warning)" << " Got unphysical values: prob = " << prob 
            << " #dof = " << nDof << std::endl;
    _value=0;
    _likelihood=0;
    setStatus(Consistency::unPhysical);
    return;
  }
  setStatus(Consistency::OK);
  if(prob != 1.0){
// use the cernlib function to get chisq.  Note the funny convention on prob!!
    float value = 1.0-float(_value);
    int ndof = nDof;
    if(value < 1.0)
      _chisq = chisin_(value,ndof);
    else
      _chisq = log(double(FLT_MAX));
// use the same algorithm as above to get loglikelihood
    double z2 = 0.5*_chisq;
    double n2 = 0.5*_nDof;
    if (_chisq==0) { 
      _likelihood=1; 
    } else {
      double loglike=(n2-1)*log(z2)-z2-NumRecipes::gammln(n2);
      if ( loglike < DBL_MIN_EXP ) {
        _likelihood = 0; 
        setStatus(Consistency::underFlow);
      } else {
    _likelihood = 0.5*exp(loglike);
      }
    }
  }
}

◆ ChisqConsistency() [4/4]

ChisqConsistency::ChisqConsistency ( const ChisqConsistency & other )

Definition at line 116 of file ChisqConsistency.cxx.

                                                                :
  Consistency(other), _chisq(other._chisq), _nDof(other._nDof)
{}

◆ ~ChisqConsistency()

virtual ChisqConsistency::~ChisqConsistency ( )

inlinevirtual

Definition at line 44 of file ChisqConsistency.h.

44{}

Member Function Documentation

◆ chisqValue()

const double & ChisqConsistency::chisqValue ( ) const

inline

Definition at line 46 of file ChisqConsistency.h.

46{ return _chisq; }

◆ nDOF()

const double & ChisqConsistency::nDOF ( ) const

inline

Definition at line 47 of file ChisqConsistency.h.

47{ return _nDof; }

◆ operator=()

ChisqConsistency & ChisqConsistency::operator= ( const ChisqConsistency & other )

Definition at line 121 of file ChisqConsistency.cxx.

                                                          {
  if(this != &other){
    Consistency::operator =(other);
    _chisq = other._chisq;
    _nDof = other._nDof;
  }
  return *this;
}

Member Data Documentation

◆ _chisq

double ChisqConsistency::_chisq

protected

Definition at line 50 of file ChisqConsistency.h.

Referenced by ChisqConsistency(), chisqValue(), and operator=().

◆ _nDof

double ChisqConsistency::_nDof