void LimitiGauss(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 s;
  double val=1.-CL;
  // devo assumere una varianza per il processo
  double sigma=sqrt(Nbackground);
  for (s=-Nbackground/2;s<Nobserved*2;) {
        int k=0;
        double gauss=0;
        gauss=ROOT::Math::gaussian_cdf(Nobserved,sigma,Nbackground+s);
        printf("s= %f  CDF Gauss = %f \n",s, gauss);
        if (gauss<val) break;
	s+=0.01;
  }
  printf("Al %f di C.L. con Nosservati = %f e Nfondo = %f il limite sul segnale e' = %f \n", CL, Nobserved,Nbackground,s );

}