| 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 #include "DicomBeamCompensator.hh"
27
28 #include "dcmtk/dcmrt/seq/drtcos.h"
29
30 #include "G4UIcommand.hh"
31
32 // DOC at https://www.dabsoft.ch/dicom/3/C.8.8
33 //....oooOO0OOooo........oooOO0OOooo........oo
34 DicomBeamCompensator::DicomBeamCompensator(DRT
35 {
36 OFString fstr;
37 Sint32 fint;
38 Float64 ffloat;
39 OFVector<Float64> fvfloat;
40 OFCondition cond;
41 G4cout << " DicomBeamCompensator::DicomBeamC
42 cond = bcompItem.getCompensatorNumber(fint);
43 theCompensatorNumber = fint;
44 G4cout << " Number " << fint << G4endl;
45
46 cond = bcompItem.getCompensatorColumns(fint)
47 theCompensatorColumns = fint;
48 cond = bcompItem.getCompensatorRows(fint);
49 theCompensatorRows = fint;
50 // first value is the spacing between the c
51 // (column spacing) is the spacing between t
52 cond = bcompItem.getCompensatorPixelSpacing(
53 theCompensatorPixelSpacing = fvfloat;
54
55 cond = bcompItem.getCompensatorPosition(fvfl
56 theCompensatorPosition = fvfloat;
57
58 cond = bcompItem.getCompensatorTransmissionD
59 if (cond.good()) theCompensatorTransmissionD
60 cond = bcompItem.getCompensatorThicknessData
61 if (cond.good()) theCompensatorThicknessData
62
63 cond = bcompItem.getCompensatorTrayID(fstr);
64 cond = bcompItem.getCompensatorType(fstr);
65
66 cond = bcompItem.getMaterialID(fstr);
67 if (cond.good()) theMaterialID = fstr;
68 cond = bcompItem.getSourceToCompensatorDista
69 if (cond.good()) theSourceToCompensatorDista
70 cond = bcompItem.getSourceToCompensatorTrayD
71 theSourceToCompensatorTrayDistance = ffloat;
72
73 cond = bcompItem.getCompensatorDescription(f
74 cond = bcompItem.getCompensatorDivergence(fs
75 cond = bcompItem.getCompensatorID(fstr);
76 cond = bcompItem.getCompensatorMountingPosit
77 cond = bcompItem.getAccessoryCode(fstr);
78 }
79
80 //....oooOO0OOooo........oooOO0OOooo........oo
81 void DicomBeamCompensator::Print(std::ostream&
82
83 //....oooOO0OOooo........oooOO0OOooo........oo
84 void DicomBeamCompensator::DumpToFile(std::ofs
85 {
86 std::string name = ":P COMP_" + G4UIcommand:
87 fout << name << "PixelSpacing_1 " << theComp
88 fout << name << "PixelSpacing_2 " << theComp
89
90 fout << name << "POSX " << theCompensatorPos
91 fout << name << "POSY " << theCompensatorPos
92
93 fout << name << "SourceToCompensatorTrayDist
94 << G4endl;
95
96 for (size_t ii = 0; ii < theSourceToCompensa
97 int iCol = ii % theCompensatorColumns;
98 int iRow = ii / theCompensatorColumns;
99 fout << name
100 << "SourceToCompensatorDistance_" + G
101 + G4UIcommand::ConvertToString(i
102 << " " << theSourceToCompensatorDista
103 }
104
105 /* for( size_t ii = 0; ii < theCompensatorT
106 int iCol = ii%theCompensatorColumns;
107 int iRow = ii/theCompensatorColumns;
108 fout << name << "Transmission_" +G4UIcomma
109 +G4UIcommand::ConvertToString(iCol) <<
110 }
111
112 for( size_t ii = 0; ii < theCompensatorThick
113 int iCol = ii%theCompensatorColumns;
114 int iRow = ii/theCompensatorColumns;
115 fout << name << "Thickness_" +G4UIcommand:
116 +G4UIcommand::ConvertToString(iCol) << " " <
117 }
118 */
119 }
120