Version:
[ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]
Diff markup
1 { 1 {
2 gROOT -> Reset(); 2 gROOT -> Reset();
3 TFile f("human_phantom.root"); 3 TFile f("human_phantom.root");
4 4
5 TDirectory* dir = (TDirectory*)f.Get("human_ph 5 TDirectory* dir = (TDirectory*)f.Get("human_phantom_ntuple");
6 TTree* ntuple = (TTree*)dir->Get("1"); 6 TTree* ntuple = (TTree*)dir->Get("1");
7 ntuple -> Print(); 7 ntuple -> Print();
8 8
9 // Print the content of the ntuple 9 // Print the content of the ntuple
10 Int_t nevent = Int_t(ntuple->GetEntries()); 10 Int_t nevent = Int_t(ntuple->GetEntries());
11 11
12 Double_t xx; 12 Double_t xx;
13 Double_t edep; 13 Double_t edep;
14 ntuple->GetBranch("organID")->SetAddress(&xx); 14 ntuple->GetBranch("organID")->SetAddress(&xx);
15 ntuple->GetBranch("edep")->SetAddress(&edep); 15 ntuple->GetBranch("edep")->SetAddress(&edep);
16 16
17 for ( Int_t i=0; i<nevent; i++ ) { 17 for ( Int_t i=0; i<nevent; i++ ) {
18 ntuple->GetEvent(i); 18 ntuple->GetEvent(i);
19 cout << "organ ID, edep (MeV): " 19 cout << "organ ID, edep (MeV): "
20 << xx << ", " << edep << endl; 20 << xx << ", " << edep << endl;
21 } 21 }
22 22
23 23
24 } 24 }
25 25