Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
LoadPlotNtuple.C
이 파일의 문서화 페이지로 가기
1 {
2 
3 //Open File where data has been stored!
4 TFile f("radioprotection_NEW.root");
5 TDirectory* dir = f.Get("radioprotection_ntuple");
6 
7 TNtuple * ntuple1 = (TNtuple*)dir->Get("101");
8 TNtuple * ntuple2 = (TNtuple*)dir->Get("102");
9 TNtuple * ntuple3 = (TNtuple*)dir->Get("103");
10 
11 int numberOfBinsX = 500;
12 int numberOfBinsY = 500;
13 int numberOfBinsZ = 500;
14 
15 int Xmin = 0;
16 int Xmax = 1000;
17 int Ymin = 0;
18 int Ymax = 1000;
19 int Zmin = 0;
20 int Zmax = 1000;
21 
22 //the type of histogram, how many variable, min and max values plus size of bins
23 TH1F* edep1Distribution = new TH1F("h0", "Primary Particle Energy Spectrum; Energy (Mev);Frequency",
24  numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
25 
26 //the type of histogram, how many variable, min and max values plus size of bins
27 TH1F* edep1Distribution = new TH1F("h1", "Energy deposition; Edep (kev);Frequency",
28  numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
29 
30 
31 TH2F* edep2DDistribution = new TH2F("h2", "Energy deposited by Ions; Z; Edep (keV)",
32  numberOfBinsX, Xmin, Xmax, // Z // binning, xmin, xmax, along x direction
33  numberOfBinsY, Ymin, Ymax); //Edep // binning, xmin, xmax, along y direction
34 
35 
36 TH3F* edep3DDistribution = new TH3F("h3", "3Dedep; Edep (keV) ; A; Z",
37  numberOfBinsX, Xmin, Xmin, // edep // binning, xmin, xmax, along x direction
38  numberOfBinsY, Ymin, Ymax, // A //binning, xmin, xmax, along y direction
39  numberOfBinsZ, Zmin, Zmax); // Z //binning, xmin, xmax, along z direction
40 
41 
42 //Plot Primary Energy of Incident Particle
43 ntuple1.Draw("Ek>>h0","","");
44 //Plot Energy Deposition within SV
45 //ntuple2.Draw("edep>>h1", "", "");
46 //Plot 2D/3D Histogram of energy with particle type using A and Z
47 //ntuple3.Draw("Z:edep>>h2", "", "");
48 //ntuple3.Draw("Z:A:edep>>h3", "", "");
49 }
int Zmin
TH1F * edep1Distribution
TH2F * edep2DDistribution
int numberOfBinsX
TNtuple * ntuple2
int Ymin
int Ymax
int Zmax
int Xmax
TH3F * edep3DDistribution
TNtuple * ntuple1
TNtuple * ntuple3
int numberOfBinsZ
TDirectory * dir
int numberOfBinsY
int Xmin