10 gStyle->SetOptStat(
"em");
14 c1 =
new TCanvas(
"c1",
"PDB DNA outputs",200,10,700,780);
17 pad1 =
new TPad(
"pad1",
"pad1",0.02,0.52,0.98,0.98,21);
18 pad2 =
new TPad(
"pad2",
"pad2",0.02,0.02,0.48,0.48,21);
19 pad3 =
new TPad(
"pad3",
"pad3",0.52,0.02,0.98,0.48,21);
21 pad1->SetFillColor(0);
23 pad2->SetFillColor(0);
25 pad3->SetFillColor(0);
29 TFile
f(
"pdb4dna_output.root");
53 hist1->GetStats(pdbStats);
54 cout <<
"-> Edep in the target : " << pdbStats[2]/1E6 <<
" MeV" << endl;
56 hist2->GetStats(pdbStats);
57 cout <<
"-> Number of SSB : " << pdbStats[2] << endl;
59 hist3->GetStats(pdbStats);
60 cout <<
"-> Number of DSB : " << pdbStats[2] << endl;