double LimitiPoisson(double Nobserved=14,double Nbackground=12, double CL=0.9)
{
  // Calcolare l'upper limit nel caso di
  // n eventi rivelati contro un fondo atteso di b eventi

  int i,j;
  double val=1.-CL;
  double s;
  double poisson=0,gauss;
  for (s=-Nbackground/2;s<Nobserved*2;) {
        int k=0;
        poisson=ROOT::Math::poisson_cdf(Nobserved,Nbackground+s);
        //printf("s= %f CDF Poisson = %f \n",s,poisson);
        if (poisson<val) break;
	s+=0.01;
  }
  printf("Al %f di C.L. con Nosservati = %f e Nfondo = %f il limite sul segnale e' = %f cum. = %f \n", CL, Nobserved,Nbackground,s,poisson );
  return s;
}