LArSoft  v07_13_02
Liquid Argon Software toolkit - http://larsoft.org/
hanson.C
Go to the documentation of this file.
1 
2 {
3  gROOT->Reset();
4 
5  // Draw histograms fill by Geant4 TestEm11 simulation
6  TFile f1("./19um.opt0.root");
7  TH1D* h1 = (TH1D*) f1.Get("12");
8  h1->SetTitle("Angular distribution of 15.7 MeV e- after 19um Au foil");
9  //h1->SetTitle("Angular distribution of 15.7 MeV e- after 9um Au foil");
10  h1->GetXaxis()->SetTitle("theta (deg)");
11  h1->GetYaxis()->SetTitle("dN/dOmega");
12  h1->SetStats(kFALSE); // Eliminate statistics box
13  h1->SetLineColor(kBlack);
14  h1->Draw("HIST");
15 
16  TFile f2("./19um.opt3.root");
17  TH1D* h2 = (TH1D*) f2.Get("12");
18  h2->SetStats(kFALSE); // Eliminate statistics box
19  h2->SetLineColor(kBlue);
20  h2->Draw("SAME HIST");
21 
22  TFile f3("./19um.local.root");
23  TH1D* h3 = (TH1D*) f3.Get("12");
24  h3->SetStats(kFALSE); // Eliminate statistics box
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 634.
32 */
33 
34  ifstream in;
35  in.open("./19um.ascii");
36  //in.open("./9um.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,22); // 22 for triangle TMatker
47  pt->SetMarkerColor(kRed);
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.68);
54  legend->AddEntry(h1,"opt0 ","l");
55  legend->AddEntry(h2,"opt3 ","l");
56  legend->AddEntry(h3,"local","l");
57  legend->AddEntry(pt,"Hanson data","P");
58  legend->Draw();
59 }
ifstream in
Definition: hanson.C:34
Double_t x
Definition: hanson.C:39
TH1D * h2
Definition: hanson.C:17
Double_t y
Definition: hanson.C:39
TFile f3("./19um.local.root")
TLegend * legend
Definition: hanson.C:53
TMarker * pt
Definition: hanson.C:38
TFile f2("./19um.opt3.root")
TH1D * h3
Definition: hanson.C:23
int nbdata
Definition: hanson.C:41
TFile f1("./19um.opt0.root")
TH1D * h1
Definition: hanson.C:7