LArSoft  v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
plot.C
Go to the documentation of this file.
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 }
Double_t y
Definition: plot.C:276
TFile f("microbeam.root")
Double_t x
Definition: plot.C:276
TLegend * legend
Definition: plot.C:36
TMarker * pt
Definition: plot.C:21
ifstream in
Definition: plot.C:18
int nbdata
Definition: plot.C:24
TH1D * h1d
Definition: plot.C:7