Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
1
2 {
3 gROOT->Reset();
4
5 // Draw histogram fill by Geant4 TestBruce simulation
6 TFile f("./Au3.local.root");
7 TH1D* h1d = (TH1D*) f.Get("4");
8 h1d->SetTitle("Fluence distribution of 13 MeV e- in Au (93.7 mg/cm2)");
9 h1d->GetXaxis()->SetTitle("r (mm)");
10 h1d->GetYaxis()->SetTitle("Fluence");
11 h1d->SetStats(kFALSE); // Eliminate statistics box
12 h1d->Draw("HIST");
13
14 /* data
15 * Bruce et al.
16 */
17
18 ifstream in;
19 in.open("../data/Au3.13MeV.ascii");
20
21 TMarker *pt;
22 Double_t x, y;
23 // First indicate number of data
24 int nbdata = 0;
25 in >> nbdata;
26 for ( int i = 0 ; i < nbdata ; i++ ) {
27 in >> x >> y ;
28 if (!in.good()) break;
29 pt = new TMarker(x,y,22); // 22 for triangle TMatker
30 pt->SetMarkerColor(kRed);
31 pt->Draw();
32 }
33 in.close();
34
35 // Print the histograms legend
36 TLegend* legend = new TLegend(0.6,0.55,0.8,0.68);
37 legend->AddEntry(h1d,"msc95","l");
38 legend->AddEntry(pt,"Faddegon data","P");
39 legend->Draw();
40 }
41