| Geant4 Cross Reference |
1 { 1 {
2 2
3 //Open File where data has been stored! 3 //Open File where data has been stored!
4 TFile* f = new TFile("radioprotection_NEW.root 4 TFile* f = new TFile("radioprotection_NEW.root");
5 TDirectory* dir = (TDirectory*)f->Get("radiopr 5 TDirectory* dir = (TDirectory*)f->Get("radioprotection_ntuple");
6 6
7 TTree * ntuple1 = (TTree*)dir->Get("101"); 7 TTree * ntuple1 = (TTree*)dir->Get("101");
8 TTree * ntuple2 = (TTree*)dir->Get("102"); 8 TTree * ntuple2 = (TTree*)dir->Get("102");
9 TTree * ntuple3 = (TTree*)dir->Get("103"); 9 TTree * ntuple3 = (TTree*)dir->Get("103");
10 10
11 int numberOfBinsX = 500; 11 int numberOfBinsX = 500;
12 int numberOfBinsY = 500; 12 int numberOfBinsY = 500;
13 int numberOfBinsZ = 500; 13 int numberOfBinsZ = 500;
14 14
15 int Xmin = 0; 15 int Xmin = 0;
16 int Xmax = 1000; 16 int Xmax = 1000;
17 int Ymin = 0; 17 int Ymin = 0;
18 int Ymax = 1000; 18 int Ymax = 1000;
19 int Zmin = 0; 19 int Zmin = 0;
20 int Zmax = 1000; 20 int Zmax = 1000;
21 21
22 //the type of histogram, how many variable, mi 22 //the type of histogram, how many variable, min and max values plus size of bins
23 TH1F* edep1DistributionPrimary = new TH1F("h0" 23 TH1F* edep1DistributionPrimary = new TH1F("h0", "Primary Particle Energy Spectrum; Energy (Mev);Frequency",
24 numberOfBinsX, Xmin, Xmax); 24 numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
25 25
26 //the type of histogram, how many variable, mi 26 //the type of histogram, how many variable, min and max values plus size of bins
27 TH1F* edep1Distribution = new TH1F("h1", "Ener 27 TH1F* edep1Distribution = new TH1F("h1", "Energy deposition; Edep (kev);Frequency",
28 numberOfBinsX, Xmin, Xmax); 28 numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
29 29
30 30
31 TH2F* edep2DDistribution = new TH2F("h2", "Ene 31 TH2F* edep2DDistribution = new TH2F("h2", "Energy deposited by Ions; Z; Edep (keV)",
32 numberOfBinsX, Xmin, Xmax, // 32 numberOfBinsX, Xmin, Xmax, // Z // binning, xmin, xmax, along x direction
33 numberOfBinsY, Ymin, Ymax); //Ed 33 numberOfBinsY, Ymin, Ymax); //Edep // binning, xmin, xmax, along y direction
34 34
35 35
36 TH3F* edep3DDistribution = new TH3F("h3", "3De 36 TH3F* edep3DDistribution = new TH3F("h3", "3Dedep; Edep (keV) ; A; Z",
37 numberOfBinsX, Xmin, Xmin, // ed 37 numberOfBinsX, Xmin, Xmin, // edep // binning, xmin, xmax, along x direction
38 numberOfBinsY, Ymin, Ymax, // A 38 numberOfBinsY, Ymin, Ymax, // A //binning, xmin, xmax, along y direction
39 numberOfBinsZ, Zmin, Zmax); // Z 39 numberOfBinsZ, Zmin, Zmax); // Z //binning, xmin, xmax, along z direction
40 40
41 41
42 //Plot Primary Energy of Incident Particle 42 //Plot Primary Energy of Incident Particle
43 ntuple1->Draw("Ek>>h0","",""); 43 ntuple1->Draw("Ek>>h0","","");
44 //Plot Energy Deposition within SV 44 //Plot Energy Deposition within SV
45 //ntuple2.Draw("edep>>h1", "", ""); 45 //ntuple2.Draw("edep>>h1", "", "");
46 //Plot 2D/3D Histogram of energy with particle 46 //Plot 2D/3D Histogram of energy with particle type using A and Z
47 //ntuple3.Draw("Z:edep>>h2", "", ""); 47 //ntuple3.Draw("Z:edep>>h2", "", "");
48 //ntuple3.Draw("Z:A:edep>>h3", "", ""); 48 //ntuple3.Draw("Z:A:edep>>h3", "", "");
49 } 49 }
50 50