6 TFile
f1(
"./Ta.1000keV.opt0.root");
7 TH1D*
h1 = (TH1D*)
f1.Get(
"8");
8 h1->SetTitle(
"Depth dose distribution of 1000 keV e- in Ta");
9 h1->GetXaxis()->SetTitle(
"Edep (Mev.cm2/g) along x/r0 x/r0");
10 h1->GetYaxis()->SetTitle(
"MeV*cm2/g");
15 TFile
f2(
"./Ta.1000keV.opt2.root");
16 TH1D*
h2 = (TH1D*)
f2.Get(
"8");
19 h2->Draw(
"SAME HIST");
27 in.open(
"./data/Ta_1000keV.ascii");
34 for (
int i = 0 ; i <
nbdata ; i++ ) {
36 if (!in.good())
break;
37 pt =
new TMarker(x,y,22);
38 pt->SetMarkerColor(kRed);
44 TLegend*
legend =
new TLegend(0.6,0.55,0.8,0.68);
45 legend->AddEntry(h1,
"ref10-opt0 ",
"l");
46 legend->AddEntry(h2,
"ref10-opt2 ",
"l");
47 legend->AddEntry(pt,
"Sandia data",
"P");