6 TFile
f1(
"./Ta_1000keV_opt3.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");
28 in.open(
"./data/Ta_1000keV.ascii");
35 for (
int i = 0 ; i <
nbdata ; i++ ) {
37 if (!in.good())
break;
38 pt =
new TMarker(x,y,22);
39 pt->SetMarkerColor(kRed);
45 TLegend*
legend =
new TLegend(0.6,0.55,0.8,0.68);
46 legend->AddEntry(h1,
"ref10-opt3 ",
"l");
48 legend->AddEntry(pt,
"Sandia data",
"P");
TFile f1("./Ta_1000keV_opt3.root")