JMathSupportkit.cc File Reference

Test application of consistency between ROOT and C++ or Jpp functions. More...

#include <string>
#include <iostream>
#include <iomanip>
#include <cmath>
#include "TMath.h"
#include "Math/PdfFuncMathCore.h"
#include "Math/ProbFuncMathCore.h"
#include "JMath/JMathSupportkit.hh"
#include "Jeep/JParser.hh"
#include "Jeep/JMessage.hh"
#include "Jeep/JPrint.hh"

Go to the source code of this file.

Macros
#define	PRINT(A)

Functions
int	main (int argc, char **argv)

Detailed Description

Test application of consistency between ROOT and C++ or Jpp functions.

Author

mdejong

Definition in file JMathSupportkit.cc.

Macro Definition Documentation

PRINT

#define PRINT

(

A

)

Value:

#A << " = " << SCIENTIFIC(12,3) << A << ';'

Function Documentation

main()

int main	(	int	argc,
		char **	argv )

Definition at line 23 of file JMathSupportkit.cc.

{
  using namespace std;
  using namespace JPP;
 
  double       precision;
  int          debug;
 
  try {
 
    JParser<> zap;
 
    zap['p'] = make_field(precision)    = 1.0e-10;
    zap['d'] = make_field(debug)        = 3;
 
    zap(argc, argv);
  }
  catch(const exception &error) {
    FATAL(error.what() << endl);
  }
 
 
  for (int n = 1; n != 12; ++n) {
    for (int m = 1; m <= n; ++m) {
 
      DEBUG(setw(2) << n << ' ' << setw(2) << m << ' ' << TMath::Binomial(n,m) << ' ' << binomial(n,m) << endl);
 
      ASSERT(TMath::Binomial(n,m) == binomial(n,m), "Test of binomial function");
    }
  }
 
  for (const double x : { -1.0e10, -10.0, -1.0, 0.0, +1.0, +10.0, +1.0e10 } ) {
 
    DEBUG("TMath::Erf  (" << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << TMath::Erf(x)   << ' ' <<
          "erf("          << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << erf(x)          << endl);
 
    ASSERT(fabs(TMath::Erf(x) - erf(x))      <= precision, "Test of erf function");
  }
 
  for (const double x : { -1.0e10, -10.0, -1.0, 0.0, +1.0, +10.0, +1.0e10 } ) {
 
    DEBUG("TMath::Erfc (" << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << TMath::Erfc (x) << ' ' << 
          "erfc("         << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << erfc(x)         << endl);
 
    ASSERT(fabs(TMath::Erfc(x) - erfc(x))    <= precision, "Test of erfc function");
  }
 
  for (const double x : { 1.0e-10, +1.0, +10.0, 100.0 } ) {
 
    DEBUG("TMath::Gamma(" << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << TMath::Gamma(x) << ' ' <<
          "tgamma("       << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << tgamma(x)       << endl);
 
    ASSERT(fabs(TMath::Gamma(x) - tgamma(x)) <= precision * tgamma(x), "Test of Gamma function");
  }
 
  for (const double a : { 1.0e-10, +1.0, +10.0, +1.0e10 } ) {
    for (const double x : { 1.0e-10, +1.0, +10.0, +1.0e10 } ) {
 
      DEBUG("TMath::Gamma(" << SCIENTIFIC(12,5) << a << "," << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << TMath::Gamma(a, x) << ' ' << 
            "Gamma("        << SCIENTIFIC(12,5) << a << "," << SCIENTIFIC(12,5) << x << ") " << SCIENTIFIC(12,5) << Gamma(a, x)        << endl);
 
      ASSERT(fabs(TMath::Gamma(a,x) - Gamma(a,x)) <= precision, "Test of Gamma function");
    }
  }
 
  for (int n : { 0, 1, 10, 100 }) {
    for (const double mu : { 1.0e-10, +1.0, +10.0, +1.0e10 } ) {
 
      DEBUG("ROOT::Math::poisson_pdf(" << n << "," << SCIENTIFIC(12,5) << mu << ") " << SCIENTIFIC(12,5) << ROOT::Math::poisson_pdf(n, mu) << ' ' << 
            "poisson("                 << n << "," << SCIENTIFIC(12,5) << mu << ") " << SCIENTIFIC(12,5) << poisson(n, mu)                 << endl);
 
      ASSERT(fabs(ROOT::Math::poisson_cdf(n, mu) - Poisson(n,mu)) <= precision, "Test of Poisson cumulative density function");
    }
  }
 
  for (int n : { 0, 1, 10, 100 }) {
    for (const double mu : { 1.0e-10, +1.0, +10.0, +1.0e10 } ) {
 
      DEBUG("ROOT::Math::poisson_cdf(" << n << "," << SCIENTIFIC(12,5) << mu << ") " << SCIENTIFIC(12,5) << ROOT::Math::poisson_cdf(n, mu) << ' ' << 
            "Poisson("                 << n << "," << SCIENTIFIC(12,5) << mu << ") " << SCIENTIFIC(12,5) << Poisson(n, mu)                 << endl);
 
      ASSERT(fabs(ROOT::Math::poisson_cdf(n, mu) - Poisson(n,mu)) <= precision, "Test of Poisson cumulative density function");
    }
  }
 
  for (double P : { 1.0e-3, 1.0e-4, 1.0e-5 }) {
    for (double x : { 1.0e-1, 1.0e0, 1.0e1 }) {
 
      DEBUG("getNs(" << FIXED(7,3) << x << "," << FIXED(9,5) << P << ") = " << setw(3)    << getNs(x,P)     << ' ' <<
            "getFs(" << FIXED(7,3) << x << "," << FIXED(9,5) << P << ") = " << FIXED(7,3) << getFs(x,P,1.0) << endl);
 
      ASSERT(fabs(getFs(x,P,1.0) - getNs(x,P)) <= precision, "Test of signal events."); 
    }
  }
  
  for (double a : { -1.0, +1.0, -1.0e-99, +1.0e-99, 0.0 }) {
    for (double b : { -1.0, +1.0, -1.0e-99, +1.0e-99, 0.0 }) {
      for (double c : { -1.0, +1.0, -1.0e-99, +1.0e-99, 0.0 }) {
 
        auto f1 = [a,b,c](const double x) { return a*x*x + b*x + c; };
 
        const roots_type rs = roots(a, b, c, precision);
 
        const double x1 = rs.x1;
        const double x2 = rs.x2;
 
#define PRINT(A) #A << " = " << SCIENTIFIC(12,3) << A << ';'
 
        if (rs.roots == roots_type::DOUBLE_t) {
          ASSERT(fabs(f1(x1)) <= precision, PRINT(a) << ' ' << PRINT(b) << ' ' << PRINT(c) << ' ' << PRINT(x1) << " -> " << PRINT(f1(x1)));
          ASSERT(fabs(f1(x1)) <= precision, PRINT(a) << ' ' << PRINT(b) << ' ' << PRINT(c) << ' ' << PRINT(x2) << " -> " << PRINT(f1(x2)));
        }
 
        if (rs.roots == roots_type::SINGLE_t) {
          ASSERT(fabs(f1(x1)) <= precision, PRINT(a) << ' ' << PRINT(b) << ' ' << PRINT(c) << ' ' << PRINT(x1) << " -> " << PRINT(f1(x1)));
        }
      }
    }
  }
 
  return 0;
}