dna/microdosimetry/plot.C

0001 // -------------------------------------------------------------------
0002 // -------------------------------------------------------------------
0003 //
0004 // *********************************************************************
0005 // To execute this macro under ROOT after your simulation ended,
0006 //   1 - launch ROOT (usually type 'root' at your machine's prompt)
0007 //   2 - type '.X plot.C' at the ROOT session prompt
0008 // *********************************************************************
0009
0010 void plot()
0011 {
0012   gROOT->Reset();
0013   gStyle->SetPalette(1);
0014   gROOT->SetStyle("Plain");
0015
0016   TCanvas* c1 = new TCanvas ("c1","",20,20,1000,500);
0017   c1->Divide(2,1);
0018
0019   // Uncomment if merging should be done
0020   //system ("rm -rf dna.root");
0021   //system ("hadd dna.root dna_*.root");
0022
0023   TFile* f = new TFile("dna.root");
0024
0025   TNtuple* ntuple;
0026   ntuple = (TNtuple*)f->Get("dna");
0027   bool rowWise = true;
0028   TBranch* eventBranch = ntuple->FindBranch("row_wise_branch");
0029   if ( ! eventBranch ) rowWise = false;
0030   // std::cout <<  "rowWise: " << rowWise << std::endl;
0031
0032   // Canvas tab 1
0033   c1->cd(1);
0034   gStyle->SetOptStat(000000);
0035
0036   // All
0037   ntuple->SetFillStyle(1001);
0038   ntuple->SetFillColor(2);
0039   ntuple->Draw("flagProcess","","B");
0040
0041   // Excitation
0042   ntuple->SetFillStyle(1001);
0043   ntuple->SetFillColor(3);
0044   ntuple->Draw("flagProcess","flagProcess==12||flagProcess==15||flagProcess==22||flagProcess==32||flagProcess==42||flagProcess==52||flagProcess==62","Bsame");
0045
0046   // Elastic
0047   ntuple->SetFillStyle(1001);
0048   ntuple->SetFillColor(4);
0049   ntuple->Draw("flagProcess","flagProcess==11||flagProcess==21||flagProcess==31||flagProcess==41||flagProcess==51||flagProcess==61||flagProcess==110||flagProcess==210||flagProcess==410||flagProcess==510||flagProcess==710||flagProcess==120||flagProcess==220||flagProcess==420||flagProcess==520||flagProcess==720","Bsame");
0050
0051   // Ionisation
0052   ntuple->SetFillStyle(1001);
0053   ntuple->SetFillColor(5);
0054   ntuple->Draw("flagProcess","flagProcess==13||flagProcess==23||flagProcess==33||flagProcess==43||flagProcess==53||flagProcess==63||flagProcess==73||flagProcess==130||flagProcess==230||flagProcess==430||flagProcess==530||flagProcess==730","Bsame");
0055
0056   // Charge decrease
0057   //ntuple->SetFillStyle(1001);
0058   //ntuple->SetFillColor(6);
0059   //ntuple->Draw("flagProcess","flagProcess==24||flagProcess==44||flagProcess==54","Bsame");
0060
0061   // Charge increase
0062   //ntuple->SetFillStyle(1001);
0063   //ntuple->SetFillColor(7);
0064   //ntuple->Draw("flagProcess","flagProcess==35||flagProcess==55||flagProcess==65","Bsame");
0065
0066   gPad->SetLogy();
0067
0068   // Canvas tab 2
0069   c1->cd(2);
0070
0071   // Electrons
0072   ntuple->SetMarkerColor(2);
0073   ntuple->SetMarkerStyle(20);
0074   ntuple->SetMarkerSize(.2);
0075   ntuple->Draw("x:y:z","flagParticle==1");
0076
0077   // Protons, hydrogen
0078   ntuple->SetMarkerColor(4);
0079   ntuple->Draw("x:y:z","flagParticle==2 || flagParticle==3 ","same");
0080 }