LArSoft
v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
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
(
"./opt3.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
36
Float_t
x
,
y
;
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
}
h1f
TH1F * h1f
Definition:
tramu.C:34
h1
TH1D * h1
Definition:
tramu.C:6
x_max
float x_max
Definition:
tramu.C:33
y
Float_t y
Definition:
tramu.C:36
in
ifstream in
Definition:
tramu.C:27
legend
TLegend * legend
Definition:
tramu.C:49
x_min
float x_min
Definition:
tramu.C:32
nb_bins
int nb_bins
Definition:
tramu.C:31
x
Float_t x
Definition:
tramu.C:36
f1
TFile f1("./opt3.root")
geant4.10.06.p01
examples
extended
electromagnetic
TestEm5
tramu
tramu.C
Generated on Thu May 2 2024 20:59:32 for LArSoft by
1.8.11