| Geant4 Cross Reference |
1 { 1 {
2 gROOT -> Reset(); 2 gROOT -> Reset();
3 TFile f("brachytherapy.root"); 3 TFile f("brachytherapy.root");
>> 4
>> 5 ntuple -> Print();
>> 6
>> 7 Int_t index;
>> 8 Double_t xx;
>> 9 Double_t yy;
>> 10 Double_t zz;
>> 11 Double_t edep;
>> 12 ntuple->GetBranch("xx")->SetAddress(&xx);
>> 13 ntuple->GetBranch("yy")->SetAddress(&yy);
>> 14 ntuple->GetBranch("zz")->SetAddress(&zz);
>> 15 ntuple->GetBranch("edep")->SetAddress(&edep);
>> 16
>> 17 // Print the content of the ntuple
>> 18 /*Int_t nevent = Int_t(ntuple->GetEntries());
>> 19
>> 20 for ( Int_t i=0; i<nevent; i++ ) {
>> 21 ntuple->GetEvent(i);
>> 22 cout << "xx, yy, zz, edep: "
>> 23 << xx << ", " << yy << ", " << zz << ", " << edep << endl;
>> 24 }
>> 25 */
>> 26
>> 27 // The phantom is 30 cm wide along x, y, z
>> 28 // the voxel size is 1 mm. The number of voxels is 300 along x, y, z
>> 29
>> 30 // Plot the energy deposition in the phantom in 3D
4 TCanvas* c1 = new TCanvas("c1", " "); 31 TCanvas* c1 = new TCanvas("c1", " ");
5 h20->Draw(""); << 32
>> 33 TH3F* edepDDistribution3D = new TH3F("h30", "3Dedepxyz",
>> 34 300, -150, 150, // binning, xmin, xmax, along x direction
>> 35 300, -150, 150, // binning, xmin, xmax, along y direction
>> 36 300, -150, 150);// binning, xmin, xmax, along z direction
>> 37
>> 38 gStyle->SetPalette(1);
>> 39
>> 40 ntuple.Draw("xx:yy:zz:edep>>h30", "", "colz");
6 } 41 }
7 42