Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/microprox/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 after your simulation ended,
  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 // Author: Sebastien Incerti
  7 // Date: March 2nd, 2019
  8 // The Geant4-DNA collaboration
  9 // *********************************************************************
 10 
 11 void SetLeafAddress(TNtuple* ntuple, const char* name, void* address);
 12 
 13 void plot()
 14 {
 15 gROOT->Reset();
 16 gStyle->SetPalette(1);
 17 gROOT->SetStyle("Plain");
 18 gStyle->SetOptStat(00000);
 19 
 20 //********************
 21 Int_t nbRadius = 101;
 22 //********************
 23 
 24 TCanvas *c1 = new TCanvas ("c1","",60,60,800,800);
 25 Int_t mycolor;
 26 
 27 TFile f("t.root");
 28 mycolor=4;
 29 
 30 TNtuple* ntuple;
 31 ntuple = (TNtuple*)f.Get("t");
 32 
 33 bool rowWise = true;
 34 TBranch* eventBranch = ntuple->FindBranch("row_wise_branch");
 35 if ( ! eventBranch ) rowWise = false;
 36 
 37 Double_t radius1,nofHits,nbEdep,edep,radius2,Einc;
 38 Int_t noRadius;
 39 
 40 if ( ! rowWise ) {
 41   ntuple->SetBranchAddress("radius1",&radius1);
 42   ntuple->SetBranchAddress("noRadius",&noRadius);
 43   ntuple->SetBranchAddress("nbHits",&nofHits);
 44   ntuple->SetBranchAddress("nbScoredHits",&nbEdep);
 45   ntuple->SetBranchAddress("edep",&edep);
 46   ntuple->SetBranchAddress("radius2",&radius2);
 47   ntuple->SetBranchAddress("Einc",&Einc);
 48   }
 49 else {
 50   SetLeafAddress(ntuple, "radius1",&radius1);
 51   SetLeafAddress(ntuple, "noRadius",&noRadius);
 52   SetLeafAddress(ntuple, "nbHits",&nofHits);
 53   SetLeafAddress(ntuple, "nbScoredHits",&nbEdep);
 54   SetLeafAddress(ntuple, "edep",&edep);
 55   SetLeafAddress(ntuple, "radius2",&radius2);
 56   SetLeafAddress(ntuple, "Einc",&Einc);
 57 }
 58 
 59 Int_t nentries = (Int_t)ntuple->GetEntries();
 60 
 61 //
 62 
 63 Double_t t[1000]; // 1000 is the max number of radius values
 64 Double_t population[1000];
 65 Double_t myRad[1000];
 66 
 67 for (Int_t i=0; i<1000; i++)
 68 {
 69  t[i]=0;
 70  population[i]=0;
 71  myRad[i]=0;
 72 }
 73 
 74 Int_t event = 0;
 75 
 76 for (Int_t i=0; i<nentries; i++)
 77 {
 78   ntuple->GetEntry(i);
 79   t[noRadius] = t[noRadius] + edep;
 80   population[noRadius]=population[noRadius]+1;
 81   myRad[noRadius] = radius1;
 82 }
 83 
 84 // Mean
 85 
 86 for (Int_t j=1; j<nbRadius; j++)
 87 {
 88   t[j] = t[j]/population[j];
 89   t[j] = t[j]/(myRad[j+1]-myRad[j]);
 90   //cout << j << " " << myRad[j] << " " << myRad[j+1]
 91   // << " " << t[j] << " " << population[j] << endl;
 92 }
 93 
 94 //
 95 
 96  c1->cd(1);
 97 
 98  TGraph* gr1 =new TGraph(nbRadius,myRad,t);
 99  gr1->SetMarkerColor(2);
100  gr1->SetMarkerStyle(20);
101  gr1->SetMarkerSize(1);
102  gr1->SetLineColor(2);
103  gr1->SetTitle("");
104  gr1->GetXaxis()->SetLimits(0.1,100);
105  gr1->GetYaxis()->SetLimits(0.,22);
106  gr1->GetXaxis()->SetLabelSize(0.025);
107  gr1->GetYaxis()->SetLabelSize(0.025);
108  gr1->GetXaxis()->SetTitleSize(0.035);
109  gr1->GetYaxis()->SetTitleSize(0.035);
110  gr1->GetXaxis()->SetTitleOffset(1.4);
111  gr1->GetYaxis()->SetTitleOffset(1.4);
112  gr1->GetXaxis()->SetTitle("r (nm)");
113  gr1->GetYaxis()->SetTitle("t (eV/nm)");
114  gr1->Draw("");
115  gPad->SetLogx();
116 
117 }
118 
119 void SetLeafAddress(TNtuple* ntuple, const char* name, void* address) {
120   TLeaf* leaf = ntuple->FindLeaf(name);
121   if ( ! leaf ) {
122     std::cerr << "Error in <SetLeafAddress>: unknown leaf --> " << name << std::endl;
123     return;
124   }
125   leaf->SetAddress(address);
126 }
127