Geant4 Cross Reference |
1 // ------------------------------------------- 1 2 // ------------------------------------------- 3 // 4 // ******************************************* 5 // To execute this macro under ROOT after your 6 // 1 - launch ROOT (usually type 'root' at y 7 // 2 - type '.X plot.C' at the ROOT session 8 // ******************************************* 9 10 void plot() 11 { 12 gROOT->Reset(); 13 gStyle->SetPalette(1); 14 gROOT->SetStyle("Plain"); 15 16 TCanvas* c1 = new TCanvas ("c1","",20,20,100 17 c1->Divide(2,1); 18 19 // Uncomment if merging should be done 20 //system ("rm -rf dna.root"); 21 //system ("hadd dna.root dna_*.root"); 22 23 TFile* f = new TFile("dna.root"); 24 25 TNtuple* ntuple; 26 ntuple = (TNtuple*)f->Get("dna"); 27 bool rowWise = true; 28 TBranch* eventBranch = ntuple->FindBranch("r 29 if ( ! eventBranch ) rowWise = false; 30 // std::cout << "rowWise: " << rowWise << s 31 32 // Canvas tab 1 33 c1->cd(1); 34 gStyle->SetOptStat(000000); 35 36 // All 37 ntuple->SetFillStyle(1001); 38 ntuple->SetFillColor(2); 39 ntuple->Draw("flagProcess","","B"); 40 41 // Excitation 42 ntuple->SetFillStyle(1001); 43 ntuple->SetFillColor(3); 44 ntuple->Draw("flagProcess","flagProcess==12| 45 46 // Elastic 47 ntuple->SetFillStyle(1001); 48 ntuple->SetFillColor(4); 49 ntuple->Draw("flagProcess","flagProcess==11| 50 51 // Ionisation 52 ntuple->SetFillStyle(1001); 53 ntuple->SetFillColor(5); 54 ntuple->Draw("flagProcess","flagProcess==13| 55 56 // Charge decrease 57 //ntuple->SetFillStyle(1001); 58 //ntuple->SetFillColor(6); 59 //ntuple->Draw("flagProcess","flagProcess==2 60 61 // Charge increase 62 //ntuple->SetFillStyle(1001); 63 //ntuple->SetFillColor(7); 64 //ntuple->Draw("flagProcess","flagProcess==3 65 66 gPad->SetLogy(); 67 68 // Canvas tab 2 69 c1->cd(2); 70 71 // Electrons 72 ntuple->SetMarkerColor(2); 73 ntuple->SetMarkerStyle(20); 74 ntuple->SetMarkerSize(.2); 75 ntuple->Draw("x:y:z","flagParticle==1"); 76 77 // Protons, hydrogen 78 ntuple->SetMarkerColor(4); 79 ntuple->Draw("x:y:z","flagParticle==2 || fla 80 } 81