Geant4  10.01.p01
DicomDetectorConstruction.hh
Go to the documentation of this file.
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 // $Id: DicomDetectorConstruction.hh 74809 2013-10-22 09:49:26Z gcosmo $
27 //
30 //
31 
32 #ifndef DicomDetectorConstruction_h
33 #define DicomDetectorConstruction_h 1
34 
35 #include "globals.hh"
38 
39 #include <set>
40 
41 class G4Material;
42 class G4Box;
43 class G4LogicalVolume;
45 
46 //*******************************************************
56 //*******************************************************
57 
59 {
60 public:
61 
64 
65  virtual G4VPhysicalVolume* Construct();
66  // trigger the construction of the geometry
67 
68 protected:
70  // create the original materials
71 
72  void ReadPhantomData();
73  // read the DICOM files describing the phantom
74 
75  void ReadPhantomDataFile(const G4String& fname);
76  // read one of the DICOM files describing the phantom (usually one per Z slice).
77  // Build a DicomPhantomZSliceHeader for each file
78 
79  void MergeZSliceHeaders();
80  // merge the slice headers of all the files
81 
83  float density, G4String newMateName );
84  // build a new material if the density of the voxel is different to the other voxels
85 
87  virtual void ConstructPhantom() = 0;
88  // construct the phantom volumes.
89  // This method should be implemented for each of the derived classes
90 
91  void SetScorer(G4LogicalVolume* voxel_logic);
92 
93  virtual void ConstructSDandField();
94 
95 protected:
97 
98  // World ...
102 
106 
107  G4int fNoFiles; // number of DICOM files
108  std::vector<G4Material*> fOriginalMaterials; // list of original materials
109  std::vector<G4Material*> fMaterials;
110  // list of new materials created to distinguish different density
111  // voxels that have the same original materials
112  size_t* fMateIDs; // index of material of each voxel
113  //unsigned int* fMateIDs; // index of material of each voxel
114 
115  std::map<G4int,G4double> fDensityDiffs;
116  // Density difference to distinguish material for each original material (by index)
117 
118  std::vector<DicomPhantomZSliceHeader*> fZSliceHeaders;
119  // list of z slice header (one per DICOM files)
121  // z slice header resulted from merging all z slice headers
122 
125 
127 
128  std::set<G4LogicalVolume*> scorers;
129 
131 };
132 
133 #endif
134 
void SetScorer(G4LogicalVolume *voxel_logic)
Definition: G4Box.hh:64
std::vector< DicomPhantomZSliceHeader * > fZSliceHeaders
std::set< G4LogicalVolume * > scorers
virtual G4VPhysicalVolume * Construct()
int G4int
Definition: G4Types.hh:78
DicomPhantomZSliceHeader * fZSliceHeaderMerged
G4double density
Definition: TRTMaterials.hh:39
G4Material * BuildMaterialWithChangingDensity(const G4Material *origMate, float density, G4String newMateName)
bool G4bool
Definition: G4Types.hh:79
DicomPhantomZSliceHeader class.
DicomPhantomZSliceMerged * mergedSlices
Dicom detector construction.
Definition of the DicomPhantomZSliceHeader class.
double G4double
Definition: G4Types.hh:76
std::vector< G4Material * > fMaterials
std::vector< G4Material * > fOriginalMaterials
void ReadPhantomDataFile(const G4String &fname)
virtual void ConstructPhantom()=0
std::map< G4int, G4double > fDensityDiffs