double Zech(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;
  //  for (s=-Nbackground/2;s<Nobserved*2; s++) {
  for (s=-Nbackground/2;s<Nobserved*2;) {
        int k=0;
        double poisson=0,gauss=0,zech=0;
        poisson=ROOT::Math::poisson_cdf(Nobserved,Nbackground+s);
        zech=(ROOT::Math::poisson_cdf(Nobserved,Nbackground+s))/(ROOT::Math::poisson_cdf(Nobserved,Nbackground));
        //printf("s= %f CDF Poisson = %f CDF Zech = %f \n",s,poisson,zech);
        if (zech<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 );
  return s;
}