Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/basic/B4/macros/plotHisto.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 // ROOT macro file for plotting example B4 histograms
  2 //
  3 // Can be run from ROOT session:
  4 // root[0] .x plotHisto.C
  5 
  6 {
  7   gROOT->Reset();
  8   gROOT->SetStyle("Plain");
  9 
 10   // Draw histos filled by Geant4 simulation
 11   //
 12 
 13   // Open file filled by Geant4 simulation
 14   TFile f("B4.root");
 15 
 16   // Create a canvas and divide it into 2x2 pads
 17   TCanvas* c1 = new TCanvas("c1", "", 20, 20, 1000, 1000);
 18   c1->Divide(2,2);
 19 
 20   // Draw Eabs histogram in the pad 1
 21   c1->cd(1);
 22   TH1D* hist1 = (TH1D*)f.Get("Eabs");
 23   hist1->Draw("HIST");
 24 
 25   // Draw Labs histogram in the pad 2
 26   c1->cd(2);
 27   TH1D* hist2 = (TH1D*)f.Get("Labs");
 28   hist2->Draw("HIST");
 29 
 30   // Draw Egap histogram in the pad 3
 31 
 32   TH1D* hist3 = (TH1D*)f.Get("Egap");
 33   c1->cd(3);
 34   // set logarithmic scale for y  
 35   //gPad->SetLogy(1);
 36   hist3->Draw("HIST");
 37 
 38   // Draw Lgap histogram in the pad 4
 39   c1->cd(4);
 40   // set logarithmic scale for y  
 41   //gPad->SetLogy(1);  
 42   TH1D* hist4 = (TH1D*)f.Get("Lgap");
 43   hist4->Draw("HIST");
 44 }
 45