Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/electromagnetic/TestEm5/pixe.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("pixe.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 pixe in the pad 1
 21   c1->cd(1);
 22   gPad->SetLogy(1);
 23   gPad->SetGridy(1);
 24   TH1D* hist1 = (TH1D*)f.Get("h55");
 25   hist1->SetLineColor(kBlue); 
 26   hist1->Draw("HIST");
 27   
 28   // Draw deexcitation in the pad 2
 29   c1->cd(2);
 30   gPad->SetLogy(1);
 31   gPad->SetGridy(0);  
 32   TH1D* hist2 = (TH1D*)f.Get("h51");
 33   hist2->SetLineColor(kBlue);   
 34   hist2->Draw("HIST");
 35   
 36   // Draw total in the pad 3
 37   c1->cd(3);
 38   gPad->SetLogy(1);
 39   gPad->SetGridy(1);
 40   TH1D* hist3 = (TH1D*)f.Get("h3");
 41   hist3->SetLineColor(kRed); 
 42   hist3->Draw("HIST");
 43 }  
 44