Geant4_10
macro.C
Go to the documentation of this file.
1 {
2 gROOT -> Reset();
3 TFile f("human_phantom.root");
4 
5 TDirectory* dir = (TDirectory*)f.Get("human_phantom_ntuple");
6 TTree* ntuple = (TTree*)dir->Get("1");
7 ntuple -> Print();
8 
9 // Print the content of the ntuple
10 Int_t nevent = Int_t(ntuple->GetEntries());
11 
14 ntuple->GetBranch("organID")->SetAddress(&xx);
15 ntuple->GetBranch("edep")->SetAddress(&edep);
16 
17 for ( Int_t i=0; i<nevent; i++ ) {
18  ntuple->GetEvent(i);
19  cout << "organ ID, edep (MeV): "
20  << xx << ", " << edep << endl;
21  }
22 
23 
24 }
Double_t xx
Definition: macro.C:10
TTree * ntuple
Definition: macro.C:6
Double_t edep
Definition: macro.C:13
TFile f
Definition: plotHisto.C:6
G4int Int_t
ntupleExperimental Reset()
G4double Double_t
TDirectory * dir
Definition: macro.C:5
Int_t nevent
Definition: macro.C:10
void Print(G4Element &ele)
Definition: pyG4Element.cc:56