Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/microelectronics/plot.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 // To execute this macro under ROOT, 
  3 //   1 - launch ROOT (usually type 'root' at your machine's prompt)
  4 //   2 - type '.X plot.C' at the ROOT session prompt
  5 // *********************************************************************
  6 {
  7 gROOT->Reset();
  8 gStyle->SetPalette(1);
  9 gROOT->SetStyle("Plain");
 10 Double_t scale;
 11   
 12 c1 = new TCanvas ("c1","",20,20,1000,500);
 13 c1->Divide(2,1);
 14 
 15 //system ("rm -rf microelectronics.root");
 16 //system ("hadd microelectronics.root microelectronics_*.root");
 17 
 18 TFile f("microelectronics.root"); 
 19 
 20 TNtuple* ntuple;
 21 ntuple = (TNtuple*)f.Get("microelectronics"); 
 22      
 23 c1->cd(1);
 24   gStyle->SetOptStat(000000);
 25   
 26   // All
 27   ntuple->Draw("flagProcess","","B");
 28   ntuple->SetFillColor(2);
 29   
 30   // Elastic
 31   ntuple->Draw("flagProcess","flagProcess==11","same");
 32   ntuple->SetFillColor(3); 
 33   
 34   // Ionisation
 35   ntuple->Draw("flagProcess","flagProcess==12||flagProcess==14||flagProcess==15||flagProcess==16||flagProcess==17","same");
 36   ntuple->SetFillColor(4);
 37   
 38   gPad->SetLogy();
 39 
 40 c1->cd(2);
 41 
 42   // Electrons
 43   ntuple->SetMarkerColor(2);
 44   ntuple->Draw("x:y:z/1000","flagParticle==1");
 45 
 46   // Protons
 47   ntuple->SetMarkerColor(4);
 48   ntuple->SetMarkerSize(4);
 49   ntuple->Draw("x:y:z/1000","flagParticle==2","same");
 50 
 51   // Ions
 52   ntuple->SetMarkerColor(3);
 53   ntuple->SetMarkerSize(3);
 54   ntuple->Draw("x:y:z/1000","flagParticle==3","same");
 55   
 56 }
 57