R__LOAD_LIBRARY (libMathMore)
R__LOAD_LIBRARY (runglauber_v3.2_C)

int nucloop () 
{
  gStyle->SetOptStat (0);
  gStyle->SetNdivisions ( 505 , "X");
  gStyle->SetNdivisions ( 508 , "Y");
  gStyle->SetLabelSize (0.055, "XY");

  TGlauberMC* gmc = new TGlauberMC ("Pb", "Pb", 25.);    // sigma_NN = 25 mb

  TH1F* hNucColl = new TH1F ("hnuccoll", "", 13, -0.5, 12.5);

  TObjArray*    nuclArray ;
  TGlauNucleon* nuc ;

  for (int iEvent = 0; iEvent < 100 ; iEvent++) 
  {
    gmc->NextEvent ( -1 );                               // Random bpar [fm]

    nuclArray = gmc->GetNucleons ();
    
    for (int iNuc = 0; iNuc < nuclArray->GetEntries() ; iNuc++) 
    {
      nuc = (TGlauNucleon*) nuclArray->At (iNuc) ;
      hNucColl->Fill ( nuc->GetNColl() );
    }
  }

  TCanvas* c1 = new TCanvas ("c1", "", 800, 600);
  c1->SetLogy ();

  hNucColl->Draw ();

  return 0;
}