3 TFile
f(
"human_phantom.root");
5 TDirectory*
dir = (TDirectory*)
f.Get(
"human_phantom_ntuple");
6 TTree*
ntuple = (TTree*)dir->Get(
"1");
10 Int_t
nevent = Int_t(ntuple->GetEntries());
14 ntuple->GetBranch(
"organID")->SetAddress(&xx);
15 ntuple->GetBranch(
"edep")->SetAddress(&edep);
19 cout <<
"organ ID, edep (MeV): " 20 << xx <<
", " << edep << endl;
ntupleExperimental Reset()
TFile f("brachytherapy.root")