{
// Create some histograms, a profile histogram and an ntuple
  hpx    = new TH1F("hpx","This is the px distribution",100,1,10);


// Fill histograms randomly
  Float_t px, py, pz;
  for (Int_t i = 0; i < 25000; i++) {
     gRandom->Rannor(px,py);
     hpx->Fill(px);
  }
hpx->Draw();
c1->SetLogy();
c1->Update();
}