| Geant4 Cross Reference |
1 { 1
2 gROOT->Reset();
3
4 // Draw histograms fill by Geant4 TestEm5 s
5 TFile f1("./local.root");
6 TH1D* h1 = (TH1D*) f1.Get("h1");
7 h1->SetTitle("Energy deposited in 530 um of
8 h1->GetXaxis()->SetTitle("Ekine (MeV)");
9 h1->GetYaxis()->SetTitle("nb / MeV");
10 h1->SetStats(kFALSE); // Eliminate statist
11 h1->SetLineColor(kBlue);
12 h1->Draw("HIST");
13
14 TFile f2("./opt4.root");
15 TH1D* h2 = (TH1D*) f2.Get("h1");
16 h2->SetStats(kFALSE); // Eliminate statist
17 h2->SetLineColor(kRed);
18 h2->Draw("SAME HIST");
19 /*
20 * e- 1 MeV in Silicon 530 um
21 * M.J. Berger et al. NIM 69 (p.181) 1969
22 * distribution of energy deposition
23 * (from 110 keV to 1.03 MeV by bin of 10 keV -
24 */
25
26 ifstream in;
27 in.open("530um.ascii");
28
29 TMarker *pt;
30 Double_t x, y;
31 // First indicate number of data
32 int nbdata = 0;
33 in >> nbdata;
34 for ( int i = 0 ; i < nbdata ; i++ ) {
35 in >> x >> y ;
36 if (!in.good()) break;
37 pt = new TMarker(x,y,32); // 32 for open
38 pt->SetMarkerColor(kGreen);
39 pt->Draw();
40 }
41 in.close();
42
43 // Print the histograms legend
44 TLegend *legend = new TLegend(0.6,0.6,0.8,0.
45 legend->AddEntry(h1,"local","l");
46 legend->AddEntry(h2,"opt4 ","l");
47 legend->AddEntry(pt,"Berger data","P");
48 legend->Draw();
49 }
50