+ //
+ // Draw Binary Cross-Section
+ //
+ TCanvas *c6 = new TCanvas("c6","Binary Cross-Section",400,10,600,700);
+ c6->cd();
+ Double_t max=fgWSbinary->GetMaximum(0,fgBMax)*1.01;
+ TH2F *h2f=new TH2F("h2fbinary","Differential Binary Cross-Section: #sigma^{hard}_{NN} dT_{AB}/db [fm]",2,0,fgBMax,2,0,max);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("b [fm]");
+ h2f->GetYaxis()->SetTitle("d#sigma^{hard}_{AB}/db [fm]");
+ h2f->Draw();
+ fgWSbinary->Draw("same");
+ TLegend *l1a = new TLegend(0.50,0.8,.90,0.9);
+ l1a->SetFillStyle(0);
+ l1a->SetBorderSize(0);
+ Char_t label[100];
+ sprintf(label,"#sigma_{NN}^{hard} = %.1f mbarn",fSigmaHard);
+ l1a->AddEntry(fgWSb,label,"");
+ l1a->Draw();
+ c6->Update();
+}
+
+void AliFastGlauber::DrawN() const
+{
+ //
+ // Draw Binaries per event (Ncoll)
+ //
+ TCanvas *c7 = new TCanvas("c7","Binaries per event",400,10,600,700);
+ c7->cd();
+ Double_t max=fgWSN->GetMaximum(0,fgBMax)*1.01;
+ TH2F *h2f=new TH2F("h2fhardcols","Number of hard collisions: T_{AB} #sigma^{hard}_{NN}/#sigma_{AB}^{geo}",2,0,fgBMax,2,0,max);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("b [fm]");
+ h2f->GetYaxis()->SetTitle("N_{coll}");
+ h2f->Draw();
+ fgWSN->Draw("same");
+ TLegend *l1a = new TLegend(0.50,0.75,.90,0.9);
+ l1a->SetFillStyle(0);
+ l1a->SetBorderSize(0);
+ Char_t label[100];
+ sprintf(label,"#sigma^{hard}_{NN} = %.1f mbarn",fSigmaHard);
+ l1a->AddEntry(fgWSN,label,"");
+ sprintf(label,"#sigma^{inel.}_{NN} = %.1f mbarn",fSigmaNN);
+ l1a->AddEntry(fgWSN,label,"");
+ l1a->Draw();
+ c7->Update();
+}
+
+void AliFastGlauber::DrawKernel(Double_t b) const
+{
+ //
+ // Draw Kernel
+ //
+ TCanvas *c8 = new TCanvas("c8","Kernel",400,10,600,700);
+ c8->cd();
+ fgWStarfi->SetParameter(0, b);
+ TH2F *h2f=new TH2F("h2fkernel","Kernel of Overlap function: d^{2}T_{AB}/dr/d#phi [fm^{-3}]",2,0,fgBMax,2,0,TMath::Pi());
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("r [fm]");
+ h2f->GetYaxis()->SetTitle("#phi [rad]");
+ h2f->Draw();
+ fgWStarfi->Draw("same");
+ TLegend *l1a = new TLegend(0.65,0.8,.90,0.9);
+ l1a->SetFillStyle(0);
+ l1a->SetBorderSize(0);
+ Char_t label[100];
+ sprintf(label,"b = %.1f fm",b);
+ l1a->AddEntry(fgWStarfi,label,"");
+ l1a->Draw();
+ c8->Update();
+}
+
+void AliFastGlauber::DrawAlmond(Double_t b) const
+{
+ //
+ // Draw Interaction Almond
+ //
+ TCanvas *c9 = new TCanvas("c9","Almond",400,10,600,700);
+ c9->cd();
+ fgWAlmond->SetParameter(0, b);
+ TH2F *h2f=new TH2F("h2falmond","Interaction Almond [fm^{-4}]",2,0,fgBMax,2,0,fgBMax);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("x [fm]");
+ h2f->GetYaxis()->SetTitle("y [fm]");
+ h2f->Draw();
+ fgWAlmond->Draw("same");
+ TLegend *l1a = new TLegend(0.65,0.8,.90,0.9);
+ l1a->SetFillStyle(0);
+ l1a->SetBorderSize(0);
+ Char_t label[100];
+ sprintf(label,"b = %.1f fm",b);
+ l1a->AddEntry(fgWAlmond,label,"");
+ l1a->Draw();
+ c9->Update();
+}
+
+void AliFastGlauber::DrawEnergyDensity() const
+{
+ //
+ // Draw energy density
+ //
+ TCanvas *c10 = new TCanvas("c10","Energy Density",400, 10, 600, 700);
+ c10->cd();
+ fgWEnergyDensity->SetMinimum(0.);
+ Double_t max=fgWEnergyDensity->GetMaximum(0,fgWEnergyDensity->GetParameter(0))*1.01;
+ TH2F *h2f=new TH2F("h2fenergydens","Energy density",2,0,fgBMax,2,0,max);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("b [fm]");
+ h2f->GetYaxis()->SetTitle("fm^{-4}");
+ h2f->Draw();
+ fgWEnergyDensity->Draw("same");
+ c10->Update();
+}
+
+void AliFastGlauber::DrawPathLength0(Double_t b, Int_t iopt) const
+{
+ //
+ // Draw Path Length
+ //
+ TCanvas *c11 = new TCanvas("c11","Path Length",400,10,600,700);
+ c11->cd();
+ fgWPathLength0->SetParameter(0, b);
+ fgWPathLength0->SetParameter(1, Double_t(iopt));
+ fgWPathLength0->SetMinimum(0.);
+ fgWPathLength0->SetMaximum(10.);
+ TH2F *h2f=new TH2F("h2fpathlength0","Path length",2,-TMath::Pi(), TMath::Pi(),2,0,10.);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("#phi [rad]");
+ h2f->GetYaxis()->SetTitle("l [fm]");
+ h2f->Draw();
+ fgWPathLength0->Draw("same");
+}
+
+void AliFastGlauber::DrawPathLength(Double_t b , Int_t ni, Int_t iopt) const
+{
+ //
+ // Draw Path Length
+ //
+ TCanvas *c12 = new TCanvas("c12","Path Length",400,10,600,700);
+ c12->cd();
+ fgWAlmond->SetParameter(0, b);
+ fgWPathLength->SetParameter(0, b);
+ fgWPathLength->SetParameter(1, Double_t (ni));
+ fgWPathLength->SetParameter(2, Double_t (iopt));
+ fgWPathLength->SetMinimum(0.);
+ fgWPathLength->SetMaximum(10.);
+ TH2F *h2f=new TH2F("h2fpathlength","Path length",2,-TMath::Pi(), TMath::Pi(),2,0,10.);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("#phi [rad]");
+ h2f->GetYaxis()->SetTitle("l [fm]");
+ h2f->Draw();
+ fgWPathLength->Draw("same");
+}
+
+void AliFastGlauber::DrawIntRadius(Double_t b) const
+{
+ //
+ // Draw Interaction Radius
+ //
+ TCanvas *c13 = new TCanvas("c13","Interaction Radius",400,10,600,700);
+ c13->cd();
+ fgWIntRadius->SetParameter(0, b);
+ fgWIntRadius->SetMinimum(0);
+ Double_t max=fgWIntRadius->GetMaximum(0,fgBMax)*1.01;
+ TH2F *h2f=new TH2F("h2fintradius","Interaction Density",2,0.,fgBMax,2,0,max);
+ h2f->SetStats(0);
+ h2f->GetXaxis()->SetTitle("r [fm]");
+ h2f->GetYaxis()->SetTitle("[fm^{-3}]");
+ h2f->Draw();
+ fgWIntRadius->Draw("same");