DifNumber.h File Reference

#include "CLHEP/Matrix/Matrix.h"
#include "CLHEP/Matrix/SymMatrix.h"
#include "CLHEP/Matrix/Vector.h"
#include "MdcRecoUtil/Code.h"
#include <assert.h>
#include <math.h>
#include <iosfwd>
#include "MdcRecoUtil/DifNumber.icc"

Go to the source code of this file.

Classes
class	DifNumber

Macros
#define	MATRIX_BOUND_CHECK

Functions
DifNumber	solveQuad (const DifNumber &a, const DifNumber &b, const DifNumber &c, int pref, Code &code)

Macro Definition Documentation

MATRIX_BOUND_CHECK

#define MATRIX_BOUND_CHECK

Definition at line 26 of file Reconstruction/MdcPatRec/MdcRecoUtil/include/MdcRecoUtil/DifNumber.h.

Function Documentation

solveQuad()

DifNumber solveQuad	(	const DifNumber &	a,
		const DifNumber &	b,
		const DifNumber &	c,
		int	pref,
		Code &	code )

Definition at line 91 of file DifNumber.cxx.

{
  DifNumber descr = b * b - 4.0 * a * c;
  if ( descr < 0.0 )
  { // solution not real
    code.setFail( 1341 );
    return DifNumber( 0.0 );
  }
  if ( a.number() == 0.0 )
  {
    if ( b.number() == 0.0 )
    {
      code.setFail( 1342 );
      return DifNumber( 0.0 );
    }
    code.setSuccess( 40 );
    return -c / b + a * c / pow( b, 3 );
  }
  code.setSuccess( 40 );
  descr       = sqrt( descr );
  DifNumber s = -b;
 
  if ( pref == +1 )
  { // positive solution
    s += descr;
  }
  else if ( pref == -1 )
  { // negative solution
    s -= descr;
  }
  else if ( pref == 0 )
  { // smallest solution
    if ( s > 0.0 ) { s -= descr; }
    else { s += descr; }
  }
  else
  { // illegal prefrence
    code.setFail( 1343 );
    return DifNumber( 0.0 );
  }
  s /= 2.0 * a;
  return s;
}