// Calcola l'UL Poissoniano
//per run root> .L PoissonUL.C
//PoissonUL(N)


//int n = 0;  // number of observed events
//double alpha = 1. - 0.9; // 0.9 is the CL 
//ROOT::Math::gamma_quantile_c( alpha, n+1, 1);

//For varying n you can reproduce the numbers in the table 40.3 of https://pdg.lbl.gov/2021/reviews/rpp2020-rev-statistics.pdf

 


void PoissonUL(int n, double alpha=0.10){

	double value =   ROOT::Math::gamma_quantile_c(alpha,n+1,1);
       std::cout << value << std::endl;

}