BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
ChisqConsistency Class Reference

#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
 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
 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
 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
 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 ,
  OK = 0 , noMeasure , underFlow , unPhysical ,
  OK = 0 , noMeasure , underFlow , unPhysical
}
enum  ConsistentSign {
  left = -1 , unknown = 0 , right = 1 , left = -1 ,
  unknown = 0 , right = 1 , left = -1 , unknown = 0 ,
  right = 1
}
enum  ConsistentStatus {
  OK = 0 , noMeasure , underFlow , unPhysical ,
  OK = 0 , noMeasure , underFlow , unPhysical ,
  OK = 0 , noMeasure , underFlow , unPhysical
}
enum  ConsistentSign {
  left = -1 , unknown = 0 , right = 1 , left = -1 ,
  unknown = 0 , right = 1 , left = -1 , unknown = 0 ,
  right = 1
}
enum  ConsistentStatus {
  OK = 0 , noMeasure , underFlow , unPhysical ,
  OK = 0 , noMeasure , underFlow , unPhysical ,
  OK = 0 , noMeasure , underFlow , unPhysical
}
enum  ConsistentSign {
  left = -1 , unknown = 0 , right = 1 , left = -1 ,
  unknown = 0 , right = 1 , left = -1 , unknown = 0 ,
  right = 1
}
Static Public Member Functions inherited from Consistency
static const Consistency & badMeasurement ()
static const Consistency & badMeasurement ()
static const Consistency & badMeasurement ()

Detailed Description

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

Constructor & Destructor Documentation

◆ ChisqConsistency() [1/12]

ChisqConsistency::ChisqConsistency ( )

Definition at line 38 of file ChisqConsistency.cxx.

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

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

◆ ChisqConsistency() [2/12]

ChisqConsistency::ChisqConsistency ( double chisq,
double nDof )

Definition at line 40 of file ChisqConsistency.cxx.

41 : _chisq( chisq ), _nDof( nDof ) {
42 double z2 = 0.5 * _chisq;
43 double n2 = 0.5 * _nDof;
44
45 if ( n2 <= 0 || z2 < 0 )
46 {
47 std::cout << "ErrMsg(warning)"
48 << " Got unphysical values: chisq = " << chisq << " #dof = " << nDof
49 << std::endl;
50 _value = 0;
51 _likelihood = 0;
52 setStatus( Consistency::unPhysical );
53 return;
54 }
55 setStatus( Consistency::OK );
56
57 // given that n2>0 && z2>=0, gammq will NOT abort
58 _value = NumRecipes::gammq( n2, z2 );
59
60 if ( _chisq == 0 ) { _likelihood = 1; }
61 else
62 {
63 double loglike = ( n2 - 1 ) * log( z2 ) - z2 - NumRecipes::gammln( n2 );
64 if ( loglike < DBL_MIN_EXP )
65 {
66 _likelihood = 0;
67 setStatus( Consistency::underFlow );
68 }
69 else { _likelihood = 0.5 * exp( loglike ); }
70 }
71}
exp
EvtComplex exp(const EvtComplex &c)
Definition EvtComplex.hh:209
n2
int n2
Definition SD0Tag.cxx:59
NumRecipes::gammq
static double gammq(double a, double x)
Definition NumRecipes.cxx:72
NumRecipes::gammln
static double gammln(double x)
Definition NumRecipes.cxx:40

◆ ChisqConsistency() [3/12]

ChisqConsistency::ChisqConsistency ( unsigned nDof,
double consistency )

Definition at line 73 of file ChisqConsistency.cxx.

73 : _nDof( nDof ) {
74 if ( prob >= 0.0 || prob <= 1.0 || nDof < 0 ) _value = prob;
75 else
76 {
77 std::cout << "ErrMsg(warning)"
78 << " Got unphysical values: prob = " << prob << " #dof = " << nDof << std::endl;
79 _value = 0;
80 _likelihood = 0;
81 setStatus( Consistency::unPhysical );
82 return;
83 }
84 setStatus( Consistency::OK );
85 if ( prob != 1.0 )
86 {
87 // use the cernlib function to get chisq. Note the funny convention on prob!!
88 float value = 1.0 - float( _value );
89 int ndof = nDof;
90 if ( value < 1.0 ) _chisq = chisin_( value, ndof );
91 else _chisq = log( double( FLT_MAX ) );
92 // use the same algorithm as above to get loglikelihood
93 double z2 = 0.5 * _chisq;
94 double n2 = 0.5 * _nDof;
95 if ( _chisq == 0 ) { _likelihood = 1; }
96 else
97 {
98 double loglike = ( n2 - 1 ) * log( z2 ) - z2 - NumRecipes::gammln( n2 );
99 if ( loglike < DBL_MIN_EXP )
100 {
101 _likelihood = 0;
102 setStatus( Consistency::underFlow );
103 }
104 else { _likelihood = 0.5 * exp( loglike ); }
105 }
106 }
107}
chisin_
float chisin_(const float &, const int &)

◆ ChisqConsistency() [4/12]

ChisqConsistency::ChisqConsistency ( const ChisqConsistency & other)

Definition at line 109 of file ChisqConsistency.cxx.

110 : Consistency( other ), _chisq( other._chisq ), _nDof( other._nDof ) {}
EvtCyclic3::other
Index other(Index i, Index j)
Definition EvtCyclic3.cc:116

◆ ~ChisqConsistency() [1/3]

virtual ChisqConsistency::~ChisqConsistency ( )
inlinevirtual

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

42{}

◆ ChisqConsistency() [5/12]

ChisqConsistency::ChisqConsistency ( )

◆ ChisqConsistency() [6/12]

ChisqConsistency::ChisqConsistency ( double chisq,
double nDof )

◆ ChisqConsistency() [7/12]

ChisqConsistency::ChisqConsistency ( unsigned nDof,
double consistency )

◆ ChisqConsistency() [8/12]

ChisqConsistency::ChisqConsistency ( const ChisqConsistency & )

◆ ~ChisqConsistency() [2/3]

virtual ChisqConsistency::~ChisqConsistency ( )
inlinevirtual

Definition at line 42 of file InstallArea/x86_64-el9-gcc13-opt/include/ProbTools/ChisqConsistency.h.

42{}

◆ ChisqConsistency() [9/12]

ChisqConsistency::ChisqConsistency ( )

◆ ChisqConsistency() [10/12]

ChisqConsistency::ChisqConsistency ( double chisq,
double nDof )

◆ ChisqConsistency() [11/12]

ChisqConsistency::ChisqConsistency ( unsigned nDof,
double consistency )

◆ ChisqConsistency() [12/12]

ChisqConsistency::ChisqConsistency ( const ChisqConsistency & )

◆ ~ChisqConsistency() [3/3]

virtual ChisqConsistency::~ChisqConsistency ( )
inlinevirtual

Definition at line 42 of file Reconstruction/MdcPatRec/ProbTools/include/ProbTools/ChisqConsistency.h.

42{}

Member Function Documentation

◆ chisqValue() [1/3]

const double & ChisqConsistency::chisqValue ( ) const
inline

Definition at line 44 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

44{ return _chisq; }

◆ chisqValue() [2/3]

const double & ChisqConsistency::chisqValue ( ) const
inline

Definition at line 44 of file InstallArea/x86_64-el9-gcc13-opt/include/ProbTools/ChisqConsistency.h.

44{ return _chisq; }

◆ chisqValue() [3/3]

const double & ChisqConsistency::chisqValue ( ) const
inline

Definition at line 44 of file Reconstruction/MdcPatRec/ProbTools/include/ProbTools/ChisqConsistency.h.

44{ return _chisq; }

◆ nDOF() [1/3]

const double & ChisqConsistency::nDOF ( ) const
inline

Definition at line 45 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

45{ return _nDof; }

◆ nDOF() [2/3]

const double & ChisqConsistency::nDOF ( ) const
inline

Definition at line 45 of file InstallArea/x86_64-el9-gcc13-opt/include/ProbTools/ChisqConsistency.h.

45{ return _nDof; }

◆ nDOF() [3/3]

const double & ChisqConsistency::nDOF ( ) const
inline

Definition at line 45 of file Reconstruction/MdcPatRec/ProbTools/include/ProbTools/ChisqConsistency.h.

45{ return _nDof; }

◆ operator=() [1/3]

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

Definition at line 112 of file ChisqConsistency.cxx.

112 {
113 if ( this != &other )
114 {
115 Consistency::operator=( other );
116 _chisq = other._chisq;
117 _nDof = other._nDof;
118 }
119 return *this;
120}
Consistency::operator=
Consistency & operator=(const Consistency &rhs)
Definition Consistency.cxx:65

◆ operator=() [2/3]

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

◆ operator=() [3/3]

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

Member Data Documentation

◆ _chisq

double ChisqConsistency::_chisq
protected

Definition at line 48 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

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

◆ _nDof

double ChisqConsistency::_nDof
protected

Definition at line 49 of file InstallArea/x86_64-el9-gcc13-dbg/include/ProbTools/ChisqConsistency.h.

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

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