Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
tramu.C
Go to the documentation of this file.
1 {
2  gROOT->Reset();
3 
4  // Draw histograms fill by Geant4 TestEm5 simulation
5  TFile f1("./local.root");
6  TH1D* h1 = (TH1D*) f1.Get("10");
7  h1->SetTitle("1 TeV muon in 3 m iron : kinetic energy at exit (GeV)");
8  h1->GetXaxis()->SetTitle("Ekine (GeV)");
9  h1->GetYaxis()->SetTitle("nb/GeV");
10  h1->SetStats(kFALSE); // Eliminate statistics box
11  h1->SetLineColor(kBlue);
12  h1->Draw("HIST");
13 /*
14  TFile f2("./local.root");
15  TH1D* h2 = (TH1D*) f2.Get("10");
16  h2->SetStats(kFALSE); // Eliminate statistics box
17  h2->SetLineColor(kRed);
18  h2->Draw("SAME HIST");
19 */
20 /*
21 * muon 1 TeV/c in 3 m Iron
22 * Particle Data Group. Physics Letters B 592 (2004) page 251
23 * distribution of the muon kinetic energy
24 * (from 950 GeV to 1000 GeV by bin of 0.5 GeV --> 100 bins)
25 */
26 
27  ifstream in;
28  in.open("mars14.ascii");
29 
30  // Create a new histogram with mars14.acsii values
31  int nb_bins = 100;
32  float x_min = 950.;
33  float x_max = 1000.;
34  TH1F* h1f = new TH1F("h1f","",nb_bins,x_min,x_max);
35 
37  while (1) {
38  in >> x >> y ;
39  if (!in.good()) break;
40  h1f->Fill(x,y);
41  }
42  in.close();
43 
44  // Draw histogram fill by mars14.acsii values
45  h1f->SetLineColor(kRed);
46  h1f->Draw("SAME");
47 
48  // Print the histograms legend
49  TLegend* legend = new TLegend(0.2,0.55,0.45,0.70);
50  legend->AddEntry(h1,"local (Urban90)","l");
51  //legend->AddEntry(h2,"local (msc90)","l");
52  legend->AddEntry(h1f,"Mars14 simul ","L");
53  legend->Draw();
54 
55 
56 }