{
  delete gRandom;
  gRandom = new TRandom(123456); // Setting the seed
  TH1F *gauss = new TH1F("gauss", "Gauss Distribution", 200, 0, 20);
  TH1F *landau = new TH1F("landau", "landau Distribution", 200, 0, 20);
  for(int i=0; i<100000; ++i) {
    double a=gRandom->Gaus(10,2);
    gauss->Fill(a) ;
    double b=gRandom->Landau(12, 0.5) ;
    landau->Fill(b) ;
  }
  c1=new TCanvas("c1", "Gauss - Landau", 900, 400);
  c1->Divide(2,1);
  c1->cd(1);
  gauss->Draw() ;
  c1->cd(2);
  landau->Draw() ;
}