| Geant4 Cross Reference |
1 1
2 {
3 gROOT->Reset();
4
5 // Draw histograms fill by Geant4 TestEm11
6 TFile f1("./19um-local.root");
7 TH1D* h1 = (TH1D*) f1.Get("h12");
8 h1->SetTitle("Angular distribution of 15.7
9 //h1->SetTitle("Angular distribution of 15.
10 h1->GetXaxis()->SetTitle("theta (deg)");
11 h1->GetYaxis()->SetTitle("dN/dOmega");
12 h1->SetStats(kFALSE); // Eliminate statist
13 h1->SetLineColor(kOrange);
14 h1->Draw("HIST");
15
16 TFile f2("./19um-opt3.root");
17 TH1D* h2 = (TH1D*) f2.Get("h12");
18 h2->SetStats(kFALSE); // Eliminate statist
19 h2->SetLineColor(kBlue);
20 h2->Draw("SAME HIST");
21
22 TFile f3("./19um-opt4.root");
23 TH1D* h3 = (TH1D*) f3.Get("h12");
24 h3->SetStats(kFALSE); // Eliminate statist
25 h3->SetLineColor(kGreen);
26 h3->Draw("SAME HIST");
27
28 /* data
29 * angle distribution of 15.7 MeV electrons
30 * transmitted through thin gold foils.
31 * A.O.Hanson et al. Phys.Rev.84 (1951) page 63
32 */
33
34 ifstream in;
35 in.open("./19um-wwols.ascii");
36 //in.open("./9um-wwols.ascii");
37
38 TMarker *pt;
39 Double_t x, y;
40 // First indicate number of data
41 int nbdata = 0;
42 in >> nbdata;
43 for ( int i = 0 ; i < nbdata ; i++ ) {
44 in >> x >> y ;
45 if (!in.good()) break;
46 pt = new TMarker(x,y,23); // 23 for tria
47 pt->SetMarkerColor(kMagenta);
48 pt->Draw();
49 }
50 in.close();
51
52 // Print the histograms legend
53 TLegend* legend = new TLegend(0.6,0.5,0.8,0
54 legend->AddEntry(h1,"local ","l");
55 legend->AddEntry(h2,"opt3 ","l");
56 legend->AddEntry(h3,"opt4","l");
57 legend->AddEntry(pt,"Hanson data","P");
58 legend->Draw();
59 }
60