{
   gROOT->Reset();

   // Draw histograms fill by Geant4 TestEm11 simulation
   TFile f1("./19um-local.root");
   TH1D* h1 = (TH1D*) f1.Get("h12");
   h1->SetTitle("Angular distribution of 15.7 MeV e- after 19um Au foil");
   //h1->SetTitle("Angular distribution of 15.7 MeV e- after 9um Au foil");   
   h1->GetXaxis()->SetTitle("theta (deg)");
   h1->GetYaxis()->SetTitle("dN/dOmega");
   h1->SetStats(kFALSE);  // Eliminate statistics box
   h1->SetLineColor(kOrange);
   h1->Draw("HIST");
   
   TFile f2("./19um-opt3.root");
   TH1D* h2 = (TH1D*) f2.Get("h12");
   h2->SetStats(kFALSE);  // Eliminate statistics box
   h2->SetLineColor(kBlue);
   h2->Draw("SAME HIST");
   
   TFile f3("./19um-opt4.root");
   TH1D* h3 = (TH1D*) f3.Get("h12");
   h3->SetStats(kFALSE);  // Eliminate statistics box
   h3->SetLineColor(kGreen);
   h3->Draw("SAME HIST"); 
     
/* data
* angle distribution of  15.7 MeV electrons
* transmitted through thin gold foils.
* A.O.Hanson et al. Phys.Rev.84 (1951) page 634.
*/

   ifstream in;
   in.open("./19um-wwols.ascii");
   //in.open("./9um-wwols.ascii");
   
   TMarker *pt;
   Double_t x, y;
   // First indicate number of data
   int nbdata = 0;
   in >> nbdata;
   for ( int i = 0 ; i < nbdata ; i++ ) {
      in >> x >> y ;
      if (!in.good()) break;
      pt = new TMarker(x,y,23); // 23 for triangle-down TMarker
      pt->SetMarkerColor(kMagenta);
      pt->Draw();
   }
   in.close();

   // Print the histograms legend
   TLegend* legend = new TLegend(0.6,0.5,0.8,0.68);
   legend->AddEntry(h1,"local ","l");
   legend->AddEntry(h2,"opt3 ","l");
   legend->AddEntry(h3,"opt4","l");      
   legend->AddEntry(pt,"Hanson data","P");
   legend->Draw();
}