Gabriele Cosmo
50 lines
1.9 KiB
C++
50 lines
1.9 KiB
C++
{
|
|
|
|
//Open File where data has been stored!
|
|
TFile* f = new TFile("radioprotection_NEW.root");
|
|
TDirectory* dir = (TDirectory*)f->Get("radioprotection_ntuple");
|
|
|
|
TTree * ntuple1 = (TTree*)dir->Get("101");
|
|
TTree * ntuple2 = (TTree*)dir->Get("102");
|
|
TTree * ntuple3 = (TTree*)dir->Get("103");
|
|
|
|
int numberOfBinsX = 500;
|
|
int numberOfBinsY = 500;
|
|
int numberOfBinsZ = 500;
|
|
|
|
int Xmin = 0;
|
|
int Xmax = 1000;
|
|
int Ymin = 0;
|
|
int Ymax = 1000;
|
|
int Zmin = 0;
|
|
int Zmax = 1000;
|
|
|
|
//the type of histogram, how many variable, min and max values plus size of bins
|
|
TH1F* edep1DistributionPrimary = new TH1F("h0", "Primary Particle Energy Spectrum; Energy (Mev);Frequency",
|
|
numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
|
|
|
|
//the type of histogram, how many variable, min and max values plus size of bins
|
|
TH1F* edep1Distribution = new TH1F("h1", "Energy deposition; Edep (kev);Frequency",
|
|
numberOfBinsX, Xmin, Xmax); //Edep // binning, xmin, xmax, along x direction
|
|
|
|
|
|
TH2F* edep2DDistribution = new TH2F("h2", "Energy deposited by Ions; Z; Edep (keV)",
|
|
numberOfBinsX, Xmin, Xmax, // Z // binning, xmin, xmax, along x direction
|
|
numberOfBinsY, Ymin, Ymax); //Edep // binning, xmin, xmax, along y direction
|
|
|
|
|
|
TH3F* edep3DDistribution = new TH3F("h3", "3Dedep; Edep (keV) ; A; Z",
|
|
numberOfBinsX, Xmin, Xmin, // edep // binning, xmin, xmax, along x direction
|
|
numberOfBinsY, Ymin, Ymax, // A //binning, xmin, xmax, along y direction
|
|
numberOfBinsZ, Zmin, Zmax); // Z //binning, xmin, xmax, along z direction
|
|
|
|
|
|
//Plot Primary Energy of Incident Particle
|
|
ntuple1->Draw("Ek>>h0","","");
|
|
//Plot Energy Deposition within SV
|
|
//ntuple2.Draw("edep>>h1", "", "");
|
|
//Plot 2D/3D Histogram of energy with particle type using A and Z
|
|
//ntuple3.Draw("Z:edep>>h2", "", "");
|
|
//ntuple3.Draw("Z:A:edep>>h3", "", "");
|
|
}
|