#include <TMath.h>

void random()
{
  const Int_t nEv = 1000;
  const Int_t xMin = 0;
  const Int_t xMax = 6;
  
  Double_t parBino = 0.3;
  
  Int_t parPois = 3;

  Double_t parGaus[3] = {1.0,3.0,0.8}; 
  
  gRandom->SetSeed(0);

  TH1F* hBino = new TH1F("hBino","Binomiale",100,0,xMax);
  TH1F* hPoisson = new TH1F("hPoisson","Poisson",100,0,xMax);
  TH1F* hGauss = new TH1F("hGauss","Gauss",100,xMin,xMax);
  
  for (Int_t i = 0; i<nEv; i++)
    { 
      hBino->Fill(gRandom->Binomial(xMax,parBino));
      hPoisson->Fill(gRandom->Poisson(parPois));
      hGauss->Fill(gRandom->Gaus(parGaus[1],parGaus[2]));
    }
  
  TCanvas* c2 = new TCanvas("c2","Binomiale",220,50,700,500);
  hBino->Draw();
  TCanvas* c3 = new TCanvas("c3","Poisson",240,70,700,500);
  hPoisson->Draw();
  TCanvas* c4 = new TCanvas("c4","Gauss",260,90,700,500);
  hGauss->Draw();
  
}