LArSoft
v09_90_00
Liquid Argon Software toolkit - https://larsoft.org/
Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
egs4.C
Go to the documentation of this file.
1
{
2
gROOT->Reset();
3
4
// Draw histogram filled by Geant4 TestEm2 simulation
5
TFile
f
(
"./93ref0.root"
);
6
TH1D*
h1d
= (TH1D*)
f
.Get(
"4"
);
7
h1d->SetTitle(
"30 GeV e- on 20 X0 Fe : energy dep, longit profil"
);
8
h1d->GetXaxis()->SetTitle(
"depth (X0)"
);
9
h1d->GetYaxis()->SetTitle(
"(100/E0) (dE/dt)"
);
10
h1d->SetStats(kFALSE);
// Eliminate statistics box
11
h1d->Draw(
"HIST"
);
12
13
// pdg.ascii came from egs4 simulation
14
ifstream
in
;
15
in.open(
"pdg.ascii"
);
16
17
// Create a new histogramm which egs4.acsii values
18
int
nb_bins
= 40;
19
float
x_min
= 0;
20
float
x_max
= 20;
21
TH1F*
h1f
=
new
TH1F(
"h1f"
,
""
,nb_bins,x_min,x_max);
22
23
Float_t
x
,
y
;
24
while
(1) {
25
in >> x >>
y
;
26
if
(!in.good())
break
;
27
h1f->Fill(x,y);
28
}
29
in.close();
30
31
// Draw histogram fill by egs4.acsii values
32
h1f->SetLineColor(2);
33
h1f->Draw(
"SAME"
);
34
35
// Print the histograms legend
36
TLegend *
legend
=
new
TLegend(0.65,0.55,0.85,0.68);
37
legend->AddEntry(h1d,
"93ref0"
,
"l"
);
38
legend->AddEntry(h1f,
"EGS4"
,
"L"
);
39
legend->Draw();
40
}
nb_bins
int nb_bins
Definition:
egs4.C:18
h1d
TH1D * h1d
Definition:
egs4.C:6
x_min
float x_min
Definition:
egs4.C:19
x_max
float x_max
Definition:
egs4.C:20
y
Float_t y
Definition:
egs4.C:23
legend
TLegend * legend
Definition:
egs4.C:36
h1f
TH1F * h1f
Definition:
egs4.C:21
f
TFile f("./93ref0.root")
x
Float_t x
Definition:
egs4.C:23
in
ifstream in
Definition:
egs4.C:14
geant4.10.06.p01
examples
extended
electromagnetic
TestEm2
egs4
egs4.C
Generated on Thu May 2 2024 20:59:32 for LArSoft by
1.8.11