dna/pdb4dna/analysis.C

0001 // -------------------------------------------------------------------
0002 // -------------------------------------------------------------------
0003
0004 // This macro requires the pdb4dna_output.root file generated from PDB4DNA example
0005
0006 {
0007   gROOT->Reset();
0008
0009   gStyle->SetOptStat("em");
0010
0011   TCanvas *c1;
0012   TPad *pad1, *pad2, *pad3;
0013   c1 = new TCanvas("c1","PDB DNA outputs",200,10,700,780);
0014   c1->SetFillColor(0);
0015
0016   pad1 = new TPad("pad1","pad1",0.02,0.52,0.98,0.98,21);
0017   pad2 = new TPad("pad2","pad2",0.02,0.02,0.48,0.48,21);
0018   pad3 = new TPad("pad3","pad3",0.52,0.02,0.98,0.48,21);
0019
0020   pad1->SetFillColor(0);
0021   pad1->Draw();
0022   pad2->SetFillColor(0);
0023   pad2->Draw();
0024   pad3->SetFillColor(0);
0025   pad3->Draw();
0026
0027
0028   TFile f("pdb4dna_output.root");
0029
0030   // Draw histograms
0031
0032   TH1D* hist1 = (TH1D*)f.Get("1");
0033   pad1->cd();
0034   hist1->Draw("HIST");
0035
0036   TH1D* hist2 = (TH1D*)f.Get("2");
0037   pad2->cd();
0038   hist2->Draw("HIST");
0039
0040   TH1D* hist3 = (TH1D*)f.Get("3");
0041   pad3->cd();
0042   hist3->Draw("HIST");
0043
0044   c1->Modified();
0045   c1->Update();
0046
0047
0048   // Read stats to get global quantities
0049
0050   double* pdbStats=new double[4];
0051
0052   hist1->GetStats(pdbStats);
0053   cout << "-> Edep in the target : " << pdbStats[2]/1E6 << " MeV" << endl;
0054
0055   hist2->GetStats(pdbStats);
0056   cout << "-> Number of SSB : " << pdbStats[2] << endl;
0057
0058   hist3->GetStats(pdbStats);
0059   cout << "-> Number of DSB : " << pdbStats[2] << endl;
0060 }