Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/medical/dna/pdb4dna/analysis.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 // -------------------------------------------------------------------
  3 
  4 // This macro requires the pdb4dna_output.root file generated from PDB4DNA example
  5 
  6 {
  7   gROOT->Reset();
  8  
  9   gStyle->SetOptStat("em");
 10 
 11   TCanvas *c1;
 12   TPad *pad1, *pad2, *pad3;
 13   c1 = new TCanvas("c1","PDB DNA outputs",200,10,700,780);
 14   c1->SetFillColor(0);
 15 
 16   pad1 = new TPad("pad1","pad1",0.02,0.52,0.98,0.98,21);
 17   pad2 = new TPad("pad2","pad2",0.02,0.02,0.48,0.48,21);
 18   pad3 = new TPad("pad3","pad3",0.52,0.02,0.98,0.48,21);
 19 
 20   pad1->SetFillColor(0);
 21   pad1->Draw();
 22   pad2->SetFillColor(0);
 23   pad2->Draw();
 24   pad3->SetFillColor(0);
 25   pad3->Draw();
 26 
 27 
 28   TFile f("pdb4dna_output.root");
 29 
 30   // Draw histograms
 31 
 32   TH1D* hist1 = (TH1D*)f.Get("1");
 33   pad1->cd();
 34   hist1->Draw("HIST");
 35 
 36   TH1D* hist2 = (TH1D*)f.Get("2");
 37   pad2->cd();
 38   hist2->Draw("HIST");
 39 
 40   TH1D* hist3 = (TH1D*)f.Get("3");
 41   pad3->cd();
 42   hist3->Draw("HIST");
 43 
 44   c1->Modified();
 45   c1->Update();
 46 
 47 
 48   // Read stats to get global quantities
 49 
 50   double* pdbStats=new double[4];
 51 
 52   hist1->GetStats(pdbStats);
 53   cout << "-> Edep in the target : " << pdbStats[2]/1E6 << " MeV" << endl;
 54 
 55   hist2->GetStats(pdbStats);
 56   cout << "-> Number of SSB : " << pdbStats[2] << endl;
 57 
 58   hist3->GetStats(pdbStats);
 59   cout << "-> Number of DSB : " << pdbStats[2] << endl;
 60 }
 61