#if !defined(__CINT__) || defined(__MAKECINT__)

#include <cmath>
#include <iostream>

using namespace std;

#include "TMath.h"

#include "MyBayesUL.h"

#endif

void
ul()
{
   MyBayesUL myul;

   printf("PDG 2004, Statistics page 22"
	  " in http://pdg.lbl.gov/2004/reviews/statrpp.pdf\n");
   printf("Table 32.3, nu_up\n");
   printf("\n\n\n");

   printf("CERNLIB CHISIN test\n");

   printf("------------------------\n");
   printf("n\tUP(90\%)\tUP(95\%)\n");
   printf("------------------------\n");
   for (Int_t nobs=0; nobs<=10; ++nobs) {
      Double_t cl90 = myul.GetUL(nobs, 0, 0.90);
      Double_t cl95 = myul.GetUL(nobs, 0, 0.95);
      printf("%d\t%.2f\t%.2f\n",
	     nobs, cl90, cl95);
   }
   printf("------------------------\n");
   printf("\n");
   printf("\n");
   printf("\n");


   printf("TRolke::Chi2Percentile() test\n");

   printf("------------------------\n");
   printf("n\tUP(90\%)\tUP(95\%)\n");
   printf("------------------------\n");
   for (Int_t nobs=0; nobs<=10; ++nobs) {
      Double_t cl90 = myul.GetUL_TRolke(nobs, 0.90);
      Double_t cl95 = myul.GetUL_TRolke(nobs, 0.95);
      printf("%d\t%.2f\t%.2f\n",
	     nobs, cl90, cl95);
   }
   printf("------------------------\n");

}