| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
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.14/
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 DicomBeamCompensator::DicomBeamCompensator(DRTCompensatorSequence::Item bcompItem)
35 {
36 OFString fstr;
37 Sint32 fint;
38 Float64 ffloat;
39 OFVector<Float64> fvfloat;
40 OFCondition cond;
41 G4cout << " DicomBeamCompensator::DicomBeamCompensator " << G4endl;
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 center of adjacent rows, and the second value
51 // (column spacing) is the spacing between the center of adjacent columns.
52 cond = bcompItem.getCompensatorPixelSpacing(fvfloat);
53 theCompensatorPixelSpacing = fvfloat;
54
55 cond = bcompItem.getCompensatorPosition(fvfloat);
56 theCompensatorPosition = fvfloat;
57
58 cond = bcompItem.getCompensatorTransmissionData(fvfloat);
59 if (cond.good()) theCompensatorTransmissionData = fvfloat;
60 cond = bcompItem.getCompensatorThicknessData(fvfloat);
61 if (cond.good()) theCompensatorThicknessData = fvfloat;
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.getSourceToCompensatorDistance(fvfloat);
69 if (cond.good()) theSourceToCompensatorDistance = fvfloat;
70 cond = bcompItem.getSourceToCompensatorTrayDistance(ffloat);
71 theSourceToCompensatorTrayDistance = ffloat;
72
73 cond = bcompItem.getCompensatorDescription(fstr);
74 cond = bcompItem.getCompensatorDivergence(fstr);
75 cond = bcompItem.getCompensatorID(fstr);
76 cond = bcompItem.getCompensatorMountingPosition(fstr);
77 cond = bcompItem.getAccessoryCode(fstr);
78 }
79
80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
81 void DicomBeamCompensator::Print(std::ostream&) {}
82
83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 void DicomBeamCompensator::DumpToFile(std::ofstream& fout)
85 {
86 std::string name = ":P COMP_" + G4UIcommand::ConvertToString(theCompensatorNumber) + "_";
87 fout << name << "PixelSpacing_1 " << theCompensatorPixelSpacing[0] << G4endl;
88 fout << name << "PixelSpacing_2 " << theCompensatorPixelSpacing[1] << G4endl;
89
90 fout << name << "POSX " << theCompensatorPosition[0] << G4endl;
91 fout << name << "POSY " << theCompensatorPosition[1] << G4endl;
92
93 fout << name << "SourceToCompensatorTrayDistance " << theSourceToCompensatorTrayDistance
94 << G4endl;
95
96 for (size_t ii = 0; ii < theSourceToCompensatorDistance.size(); ii++) {
97 int iCol = ii % theCompensatorColumns;
98 int iRow = ii / theCompensatorColumns;
99 fout << name
100 << "SourceToCompensatorDistance_" + G4UIcommand::ConvertToString(iRow) + "_"
101 + G4UIcommand::ConvertToString(iCol)
102 << " " << theSourceToCompensatorDistance[ii] << G4endl;
103 }
104
105 /* for( size_t ii = 0; ii < theCompensatorTransmissionData.size(); ii++ ) {
106 int iCol = ii%theCompensatorColumns;
107 int iRow = ii/theCompensatorColumns;
108 fout << name << "Transmission_" +G4UIcommand::ConvertToString(iRow) + "_"
109 +G4UIcommand::ConvertToString(iCol) << " " << theCompensatorTransmissionData[ii] << G4endl;
110 }
111
112 for( size_t ii = 0; ii < theCompensatorThicknessData.size(); ii++ ) {
113 int iCol = ii%theCompensatorColumns;
114 int iRow = ii/theCompensatorColumns;
115 fout << name << "Thickness_" +G4UIcommand::ConvertToString(iRow) + "_"
116 +G4UIcommand::ConvertToString(iCol) << " " << theCompensatorThicknessData[ii] << G4endl;
117 }
118 */
119 }
120