Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/exp_microdosimetry/LoadPlotNtuple.C

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ] 

  1 {
  2 
  3 //Open File where data has been stored!
  4 TFile* f = new TFile("radioprotection_NEW.root");
  5 TDirectory* dir = (TDirectory*)f->Get("radioprotection_ntuple"); 
  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, min and max values plus size of bins
 23 TH1F* edep1DistributionPrimary = 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 }     
 50