| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 //
27 /// \file medical/DICOM/src/DicomPhantomZSlice
28 /// \brief Implementation of the DicomPhantomZ
29 //
30 //
31 // The code was written by :
32 // * Jonathan Madsen : jonathan.madsen@ce
33 //
34 // Texas A&M University
35 // 3133 TAMU, Zachry Building
36 // College Station, TX 77843, USA
37 //
38 //********************************************
39
40 #include "DicomPhantomZSliceMerged.hh"
41
42 #include "G4SystemOfUnits.hh"
43 #include "G4UnitsTable.hh"
44
45 //....oooOO0OOooo........oooOO0OOooo........oo
46
47 DicomPhantomZSliceMerged::DicomPhantomZSliceMe
48
49 //....oooOO0OOooo........oooOO0OOooo........oo
50
51 DicomPhantomZSliceMerged::~DicomPhantomZSliceM
52 {
53 fSlices.clear();
54 }
55
56 //....oooOO0OOooo........oooOO0OOooo........oo
57
58 void DicomPhantomZSliceMerged::CheckSlices()
59 {
60 G4cout << "\nDicomPhantomZSliceMerged::Check
61
62 G4bool uniformSliceThickness = true;
63
64 if (fSlices.size() > 1) {
65 if (fSlices.size() == 2) {
66 DicomPhantomZSliceHeader* one = fSlices.
67 DicomPhantomZSliceHeader* two = fSlices.
68
69 G4double real_distance = (two->GetSliceL
70
71 if (one->GetMaxZ() != two->GetMinZ()) {
72 one->SetMaxZ(one->GetSliceLocation() +
73 two->SetMinZ(two->GetSliceLocation() -
74 // one->SetMinZ(one->GetSliceLocation(
75 // two->SetMaxZ(two->GetSliceLocation(
76 if (uniformSliceThickness) {
77 one->SetMinZ(one->GetSliceLocation()
78 two->SetMaxZ(two->GetSliceLocation()
79 }
80 }
81 }
82 else {
83 auto ite0 = fSlices.begin();
84 auto ite1 = fSlices.begin();
85 auto ite2 = fSlices.begin();
86 ++ite1;
87 ++ite2;
88 ++ite2;
89
90 for (; ite2 != fSlices.end(); ++ite0, ++
91 DicomPhantomZSliceHeader* prev = ite0-
92 DicomPhantomZSliceHeader* slice = ite1
93 DicomPhantomZSliceHeader* next = ite2-
94 G4double real_max_distance = (next->Ge
95 G4double real_min_distance = (slice->G
96 G4double real_distance = real_max_dist
97 G4double stated_distance = slice->GetM
98 if (real_distance != stated_distance)
99 uintmax_t sliceNum = std::distance(f
100 G4cout << "\tDicomPhantomZSliceMerge
101 Slice Distance Error in sl
102 << sliceNum << "]: Real Dista
103 << " mm, Stated Distance = "
104 slice->SetMinZ(slice->GetSliceLocati
105 slice->SetMaxZ(slice->GetSliceLocati
106
107 if (ite0 == fSlices.begin()) {
108 prev->SetMaxZ(slice->GetMinZ());
109 // Using below would make all slic
110 // prev->SetMinZ(prev->GetSliceLoc
111 if (uniformSliceThickness) {
112 prev->SetMinZ(prev->GetSliceLoca
113 }
114 }
115 if (static_cast<unsigned int>(std::d
116 {
117 next->SetMinZ(slice->GetMaxZ());
118 // Using below would make all slic
119 // next->SetMaxZ(next->GetSliceLoc
120 if (uniformSliceThickness) {
121 next->SetMaxZ(next->GetSliceLoca
122 }
123 }
124 }
125 }
126 }
127 }
128 G4cout << G4endl;
129
130 for (auto ite = fSlices.cbegin(); ite != fSl
131 ite->second->DumpToFile();
132 }
133 }
134 //....oooOO0OOooo........oooOO0OOooo........oo
135