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