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");
12 h1->SetLineColor(kBlack);
20 in.open(
"./EGSnrc/Ta_1000keV_EGSnrc.ascii");
27 for (
int i = 0 ; i <
nbdata ; i++ ) {
29 if (!in.good())
break;
30 pt =
new TMarker(x,y,22);
31 pt->SetMarkerColor(kRed);
37 TLegend*
legend =
new TLegend(0.6,0.5,0.8,0.70);
38 legend->AddEntry(h1,
"ref10 ",
"l");
39 legend->AddEntry(pt,
"EGSnrc",
"P");