Geant4 Cross Reference |
1 { 1 2 3 //Open File where data has been stored! 4 TFile* f = new TFile("radioprotection_NEW.root 5 TDirectory* dir = (TDirectory*)f->Get("radiopr 6 7 TTree * ntuple1 = (TTree*)dir->Get("101"); 8 TTree * ntuple2 = (TTree*)dir->Get("102"); 9 TTree * ntuple3 = (TTree*)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, mi 23 TH1F* edep1DistributionPrimary = new TH1F("h0" 24 numberOfBinsX, Xmin, Xmax); 25 26 //the type of histogram, how many variable, mi 27 TH1F* edep1Distribution = new TH1F("h1", "Ener 28 numberOfBinsX, Xmin, Xmax); 29 30 31 TH2F* edep2DDistribution = new TH2F("h2", "Ene 32 numberOfBinsX, Xmin, Xmax, // 33 numberOfBinsY, Ymin, Ymax); //Ed 34 35 36 TH3F* edep3DDistribution = new TH3F("h3", "3De 37 numberOfBinsX, Xmin, Xmin, // ed 38 numberOfBinsY, Ymin, Ymax, // A 39 numberOfBinsZ, Zmin, Zmax); // Z 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 47 //ntuple3.Draw("Z:edep>>h2", "", ""); 48 //ntuple3.Draw("Z:A:edep>>h3", "", ""); 49 } 50