void NormalCoverage()
{
  double rnd1,rnd2;
  int Nthrows=1000000;
  int nblue;
  int ncov=0;
  int ncov2=0;
  double meas;
  for (int it=0;it!=Nthrows;it++) {
    meas=gRandom->Gaus(10.,1);
    if (abs(meas-10.)<1.) ncov++;
  }
  double epsilon = (float) ncov/Nthrows;

  // Use Binomial errors
  cout << "Fraction of throws not covered " << 100.*epsilon << " +/-" << 100.*sqrt(epsilon*(1-epsilon)/Nthrows) << " %" << endl;

}