Geant4
v4-10.4-release
메인 페이지
관련된 페이지
모듈
네임스페이스
클래스
파일들
파일 목록
파일 멤버
모두
클래스
네임스페이스들
파일들
함수
변수
타입정의
열거형 타입
열거형 멤버
Friends
매크로
그룹들
페이지들
examples
extended
electromagnetic
TestEm12
berger
TestEm12/berger/berger.C
이 파일의 문서화 페이지로 가기
1
{
2
gROOT->Reset();
3
4
/*
5
* Mainegra et al. Med. Phys. 32, 685-699 (2005)
6
* Seltzer, Appl. Radiat. Isot. 42 (1991) page 917
7
* L. Ferrer et al. Cancer Bio. Rad. 22-1 (2007)
8
*/
9
10
ifstream
in
;
11
in.open(
"./EGSnrc/100keV-DPK.ascii"
);
12
//in.open("./Etran/100keV-DPK.ascii");
13
14
// Create a new histogram
15
double
x_min
= 0;
16
double
x_max
= 1.2;
17
18
int
nbdata
= 0;
19
in >>
nbdata
;
20
21
TH1D*
h1d
=
new
TH1D(
"h1d"
,
""
,nbdata,x_min,x_max);
22
23
Double_t
x
,
y
;
24
for
(
int
i = 0 ; i <
nbdata
; i++ ) {
25
in >> x >>
y
;
26
if
(!in.good())
break
;
27
h1d->Fill(x,y);
28
}
29
in.close();
30
31
// Draw histogram
32
h1d->SetTitle(
"Dose point kernel : energy deposition profile, e- 100 keV"
);
33
h1d->GetXaxis()->SetTitle(
"d(E/E0)/d(r/r0) along r/r0 r/r0"
);
34
h1d->GetYaxis()->SetTitle(
"DPK"
);
35
h1d->SetStats(kFALSE);
// Eliminate statistics box
36
h1d->SetLineColor(kRed);
37
h1d->Draw(
"L"
);
38
39
40
// Draw histograms fill by Geant4 TestEm12 simulation
41
TFile
f1
(
"./100keV_local.root"
);
42
TH1D*
h1
= (TH1D*)
f1
.Get(
"8"
);
43
h1->SetLineColor(kBlack);
44
h1->Draw(
"SAME HIST"
);
45
/*
46
TFile f2("./10keV.liver.root");
47
TH1D* h2 = (TH1D*) f2.Get("8");
48
h2->SetLineColor(kBlue);
49
h2->Draw("SAME HIST");
50
51
TFile f3("./10keV.penel.root");
52
TH1D* h3 = (TH1D*) f3.Get("8");
53
h3->SetLineColor(kGreen);
54
h3->Draw("SAME HIST");
55
*/
56
// Print the histograms legend
57
TLegend *
legend
=
new
TLegend(0.7,0.6,0.86,0.78);
58
legend->AddEntry(h1,
"local "
,
"l"
);
61
legend->AddEntry(h1d,
"EGSnrc"
,
"l"
);
62
legend->Draw();
63
64
}
x
Float_t x
Definition:
compare.C:6
G4INCL::Double_t
G4double Double_t
Definition:
G4INCLEventInfo.hh:63
nbdata
int nbdata
Definition:
egs.C:28
y
Float_t y
Definition:
compare.C:6
h1d
TH1D * h1d
Definition:
TestEm12/berger/berger.C:21
x_min
double x_min
Definition:
TestEm12/berger/berger.C:15
f1
Float_t f1
Definition:
comparison_ascii.C:56
legend
TLegend * legend
Definition:
egs.C:40
x_max
double x_max
Definition:
TestEm12/berger/berger.C:16
in
ifstream in
Definition:
comparison.C:7
h1
TH1F * h1
Definition:
advanced/microbeam/plot.C:42
다음에 의해 생성됨 :
1.8.5