{
  delete gRandom;
  gRandom = new TRandom(123456); // Setting the seed
  TH2F *g2D = new TH2F("g2D", "Gauss 2D", 200, 0, 20, 200, 0, 20) ;
  for(int i=0; i<100000; ++i) {
    double a=gRandom->Gaus(10, 2);
    double b=gRandom->Gaus(7, 1) ;
    g2D->Fill(a,b) ;
  }
  c1=new TCanvas("c1", "Gauss 2D", 600, 400);
  c1->cd();
  g2D->Draw() ;
}