Geant4 Cross Reference |
1 // ******************************************* 1 // ******************************************************************** 2 // * License and Disclaimer 2 // * License and Disclaimer * 3 // * 3 // * * 4 // * The Geant4 software is copyright of th 4 // * The Geant4 software is copyright of the Copyright Holders of * 5 // * the Geant4 Collaboration. It is provided 5 // * the Geant4 Collaboration. It is provided under the terms and * 6 // * conditions of the Geant4 Software License 6 // * conditions of the Geant4 Software License, included in the file * 7 // * LICENSE and available at http://cern.ch/ 7 // * LICENSE and available at http://cern.ch/geant4/license . These * 8 // * include a list of copyright holders. 8 // * include a list of copyright holders. * 9 // * 9 // * * 10 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 11 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 12 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 13 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 14 // * use. Please see the license in the file 14 // * use. Please see the license in the file LICENSE and URL above * 15 // * for the full disclaimer and the limitatio 15 // * for the full disclaimer and the limitation of liability. * 16 // * 16 // * * 17 // * This code implementation is the result 17 // * This code implementation is the result of the scientific and * 18 // * technical work of the GEANT4 collaboratio 18 // * technical work of the GEANT4 collaboration. * 19 // * By using, copying, modifying or distri 19 // * By using, copying, modifying or distributing the software (or * 20 // * any work based on the software) you ag 20 // * any work based on the software) you agree to acknowledge its * 21 // * use in resulting scientific publicati 21 // * use in resulting scientific publications, and indicate your * 22 // * acceptance of all terms of the Geant4 Sof 22 // * acceptance of all terms of the Geant4 Software license. * 23 // ******************************************* 23 // ******************************************************************** 24 24 25 //GEANT4 - Depth-of-Interaction enabled Positr 25 //GEANT4 - Depth-of-Interaction enabled Positron emission tomography (PET) advanced example 26 26 27 //Contributors 27 //Contributors 28 28 29 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 29 // Abdella M. Ahmed (1, 2), Andrew Chacon (1, 2), Harley Rutherford (1, 2), 30 // Hideaki Tashima (3), Go Akamatsu (3), Akram 30 // Hideaki Tashima (3), Go Akamatsu (3), Akram Mohammadi (3), Eiji Yoshida (3), Taiga Yamaya (3) 31 // Susanna Guatelli (2), and Mitra Safavi-Naei 31 // Susanna Guatelli (2), and Mitra Safavi-Naeini (1, 2) 32 32 33 // (1) Australian Nuclear Science and Technolo 33 // (1) Australian Nuclear Science and Technology Organisation, Australia 34 // (2) University of Wollongong, Australia 34 // (2) University of Wollongong, Australia 35 // (3) National Institute of Radiological Scie 35 // (3) National Institute of Radiological Sciences, Japan 36 36 37 37 38 38 39 #ifndef doiPETDetectorConstruction_h 39 #ifndef doiPETDetectorConstruction_h 40 #define doiPETDetectorConstruction_h 1 40 #define doiPETDetectorConstruction_h 1 41 41 42 #include "G4VUserDetectorConstruction.hh" 42 #include "G4VUserDetectorConstruction.hh" 43 #include "globals.hh" 43 #include "globals.hh" 44 #include "G4UnitsTable.hh" 44 #include "G4UnitsTable.hh" 45 #include "G4Element.hh" 45 #include "G4Element.hh" 46 #include "doiPETGlobalParameters.hh" 46 #include "doiPETGlobalParameters.hh" 47 47 48 class G4Material; 48 class G4Material; 49 class G4VPhysicalVolume; 49 class G4VPhysicalVolume; 50 class G4LogicalVolume; 50 class G4LogicalVolume; 51 // 51 // 52 class doiPETAnalysis; 52 class doiPETAnalysis; 53 class doiPETDetectorConstructionMessenger; 53 class doiPETDetectorConstructionMessenger; 54 54 55 /// Detector construction class to define mate 55 /// Detector construction class to define materials and geometry. 56 /// 56 /// 57 /// Crystals are positioned in Ring, with an a 57 /// Crystals are positioned in Ring, with an appropriate rotation matrix. 58 /// Several copies of Ring are placed in the f 58 /// Several copies of Ring are placed in the full detector. 59 59 60 class doiPETDetectorConstruction : public G4VU 60 class doiPETDetectorConstruction : public G4VUserDetectorConstruction 61 { 61 { 62 public: 62 public: 63 doiPETDetectorConstruction(); << 63 doiPETDetectorConstruction(); 64 virtual ~doiPETDetectorConstruction(); << 64 virtual ~doiPETDetectorConstruction(); 65 65 66 public: 66 public: 67 virtual G4VPhysicalVolume* Construct(); << 67 virtual G4VPhysicalVolume* Construct(); 68 68 69 void ConstructPhantom(G4LogicalVolume*); << 69 void ConstructPhantom(G4LogicalVolume*); 70 void ChangePhantom(G4String choice); << 70 void ChangePhantom(G4String choice); 71 void SetPhantomPosition(G4ThreeVector); << 71 void SetPhantomPosition(G4ThreeVector); 72 void SetPhantomRadius(G4double); << 72 void SetPhantomRadius(G4double); 73 void SetPhantomLength(G4double); << 73 void SetPhantomLength(G4double); 74 void SetNumberOfSleeves(G4int); << 75 74 76 75 77 private: 76 private: 78 void DefineMaterials(); << 77 void DefineMaterials(); 79 doiPETDetectorConstructionMessenger* fDetect << 78 doiPETDetectorConstructionMessenger* fDetectorMessenger; 80 doiPETAnalysis* pAnalysis; << 79 doiPETAnalysis* pAnalysis; 81 80 82 G4LogicalVolume* phantom_logicalV; << 81 G4LogicalVolume* phantom_logicalV; 83 G4VPhysicalVolume* phantom_physicalV; << 82 G4VPhysicalVolume* phantom_physicalV; 84 //G4LogicalVolume* gelatin_logicalV; << 83 G4LogicalVolume* gelatin_logicalV; 85 //G4VPhysicalVolume* gelatin_physicalV; << 84 G4VPhysicalVolume* gelatin_physicalV; 86 85 87 //world << 86 //wolrd 88 G4LogicalVolume* world_logicalV; << 87 G4LogicalVolume* world_logicalV; 89 G4VPhysicalVolume* world_physicalV; << 88 G4VPhysicalVolume* world_physicalV; 90 89 91 //detector block << 90 //detector block 92 G4LogicalVolume* blockDetector_logicalV; << 91 G4LogicalVolume* blockDetector_logicalV; 93 G4VPhysicalVolume* blockDetector_physicalV; << 92 G4VPhysicalVolume* blockDetector_physicalV; 94 93 95 //air volume to fill the detector block << 94 //air volume to fill the detector block 96 G4LogicalVolume* airBox_logicalV; << 95 G4LogicalVolume* airBox_logicalV; 97 G4VPhysicalVolume* airBox_physicalV; << 96 G4VPhysicalVolume* airBox_physicalV; 98 97 99 //crystals << 98 //crystals 100 G4LogicalVolume* crystal_logicalV; << 99 G4LogicalVolume* crystal_logicalV; 101 G4VPhysicalVolume* crystal_physicalV; << 100 G4VPhysicalVolume* crystal_physicalV; 102 101 103 102 104 //water << 103 //water 105 G4LogicalVolume* water_logicalV; << 104 G4LogicalVolume* water_logicalV; 106 G4VPhysicalVolume* water_physicalV; << 105 G4VPhysicalVolume* water_physicalV; 107 106 108 //lung << 107 //lung 109 G4LogicalVolume* lung_logicalV; << 108 G4LogicalVolume* lung_logicalV; 110 G4VPhysicalVolume* lung_physicalV; << 109 G4VPhysicalVolume* lung_physicalV; 111 110 112 //test line phantom << 111 //test line phantom 113 G4LogicalVolume* test_logicalV; << 112 G4LogicalVolume* test_logicalV; 114 G4VPhysicalVolume* test_physicalV; << 113 G4VPhysicalVolume* test_physicalV; 115 114 116 G4LogicalVolume* lung_logicalV_PMMA; << 115 G4LogicalVolume* lung_logicalV_PMMA; 117 G4VPhysicalVolume* lung_physicalVPMMA; << 116 G4VPhysicalVolume* lung_physicalVPMMA; 118 117 119 //cold regions << 118 //cold regions 120 //G4LogicalVolume* coldRegion_logicalV; << 119 G4LogicalVolume* coldRegion_logicalV; 121 //G4VPhysicalVolume* coldRegion_physicalV; << 120 G4VPhysicalVolume* coldRegion_physicalV; 122 121 123 // << 122 // 124 //Surrounding PMMA for hot sphere << 123 //Surrounding PMMA for hot sphere 125 G4LogicalVolume* hotSpherePMMA_logicalV; << 124 G4LogicalVolume* hotSpherePMMA_logicalV; 126 G4VPhysicalVolume* hotSpherePMMA_physicalV; << 125 G4VPhysicalVolume* hotSpherePMMA_physicalV; 127 126 128 //hot water phantom (activity is distributed << 127 //hot water phantom (activity is distributed) 129 G4LogicalVolume* hotSphereWater_logicalV; << 128 G4LogicalVolume* hotSphereWater_logicalV; 130 G4VPhysicalVolume* hotSphereWater_physicalV; << 129 G4VPhysicalVolume* hotSphereWater_physicalV; 131 130 132 //surrounding PMMA cold sphere << 131 //surrounding PMMA cold sphere 133 G4LogicalVolume* coldSpherePMMA_logicalV; << 132 G4LogicalVolume* coldSpherePMMA_logicalV; 134 G4VPhysicalVolume* coldSpherePMMA_physicalV; << 133 G4VPhysicalVolume* coldSpherePMMA_physicalV; 135 134 136 //cold Water phantom in the cold PMMA sphere << 135 //cold Water phantom in the cold PMMA sphere 137 G4LogicalVolume* coldSphereWater_logicalV; << 136 G4LogicalVolume* coldSphereWater_logicalV; 138 G4VPhysicalVolume* coldSphereWater_physicalV << 137 G4VPhysicalVolume* coldSphereWater_physicalV; 139 138 140 //fillable polyethylene phantom for sensitiv << 139 //fillable polyethylene phantom for sensitivity 141 G4LogicalVolume* phantomPE_logicalV; << 140 G4LogicalVolume* phantomPE_logicalV; 142 G4VPhysicalVolume* phantomPE_physicalV; << 141 G4VPhysicalVolume* phantomPE_physicalV; 143 142 144 //Image quality phantom for small animal NEM << 143 //Image quality phantom for small animal NEMA NU-4 145 G4LogicalVolume* waterPhantom_logicalV; << 144 G4LogicalVolume* waterPhantom_logicalV; 146 G4VPhysicalVolume* WaterPhantom_physicalV; << 145 G4VPhysicalVolume* WaterPhantom_physicalV; 147 146 148 G4LogicalVolume* rod_phantom_logicalV; << 147 G4LogicalVolume* rod_phantom_logicalV; 149 G4VPhysicalVolume* rod_phantom_physicalV; << 148 G4VPhysicalVolume* rod_phantom_physicalV; 150 149 151 G4LogicalVolume* chamberPMMA_logicalV; << 150 G4LogicalVolume* chamberPMMA_logicalV; 152 G4VPhysicalVolume* chamberPMMA_physicalV; << 151 G4VPhysicalVolume* chamberPMMA_physicalV; 153 152 154 // << 153 // 155 G4LogicalVolume* chamberWater_logicalV; << 154 G4LogicalVolume* chamberWater_logicalV; 156 G4VPhysicalVolume* chamberWater_physicalV; << 155 G4VPhysicalVolume* chamberWater_physicalV; 157 156 158 // << 157 // 159 G4LogicalVolume* chamberAir_logicalV; << 158 G4LogicalVolume* chamberAir_logicalV; 160 G4VPhysicalVolume* chamberAir_physicalV; << 159 G4VPhysicalVolume* chamberAir_physicalV; 161 160 162 161 163 //Dimension of the sphere << 162 //Dimension of the sphere 164 G4double spherePositionX, spherePositionY; / << 163 G4double spherePositionX, spherePositionY, spherePositionZ; 165 G4double sphereDiameter; << 164 G4double sphereDiameter; 166 G4double distanceFromCenter; << 165 G4double distanceFromCenter; 167 G4int numberOfSpheres; << 166 G4int numberOfSpheres; 168 G4double sphereWallThickness; << 167 G4double sphereWallThickness; 169 G4double zOffsetSpherePhantom; << 168 G4double zOffsetSpherePhantom; 170 169 171 G4String PhantomType; << 170 G4String PhantomType; 172 171 173 //materials << 172 //materials 174 G4Material* air; << 173 G4Material* air; 175 G4Material* pmma; << 174 G4Material* pmma; 176 G4Material* water; << 175 G4Material* water; 177 G4Material* polyethylene; << 176 G4Material* polyethylene; 178 G4Material* polyethylene_NEMA; << 177 G4Material* polyethylene_NEMA; 179 //G4Material* inflatedLung; << 178 G4Material* inflatedLung; 180 G4Material* polystyrene; << 179 G4Material* polystyrene; 181 G4Material* Aluminum; << 180 G4Material* Aluminum; 182 181 183 //elements for GSO << 182 //elements for GSO 184 G4Element* O; << 183 G4Element* O; 185 G4Element* Si; << 184 G4Element* Si; 186 G4Element* Gd; << 185 G4Element* Gd; 187 G4Material* GSO; << 186 G4Material* GSO; 188 187 189 G4Material* crystalMaterial; << 188 G4Material* crystalMaterial; 190 //G4Material* phantomMaterial; << 189 //G4Material* phantomMaterial; 191 190 192 G4bool fCheckOverlaps; << 191 G4bool fCheckOverlaps; 193 G4bool isotopes; << 192 G4bool isotopes; 194 193 195 194 196 //size of world << 195 //size of world 197 G4double worldSizeX; << 196 G4double worldSizeX; 198 G4double worldSizeY; << 197 G4double worldSizeY; 199 G4double worldSizeZ; << 198 G4double worldSizeZ; 200 199 201 200 202 //The following is moved to doiPETGlobalPara << 201 //The following is moved to doiPETGlobalParameters.hh 203 //G4int numberOfCrystal_DOI; << 202 //G4int numberOfCrystal_DOI; 204 //G4int numberOfCrystal_tangential; << 203 //G4int numberOfCrystal_tangential; 205 //G4int numberOfCrystal_axial; << 204 //G4int numberOfCrystal_axial; 206 205 207 //// << 206 //// 208 //G4double sizeOfCrystal_DOI; << 207 //G4double sizeOfCrystal_DOI; 209 //G4double sizeOfCrystal_tangential; << 208 //G4double sizeOfCrystal_tangential; 210 //G4double sizeOfCrystal_axial; << 209 //G4double sizeOfCrystal_axial; 211 210 212 //// << 211 //// 213 //G4double crystalGap_DOI; << 212 //G4double crystalGap_DOI; 214 //G4double crystalGap_tangential; << 213 //G4double crystalGap_tangential; 215 //G4double crystalGap_axial; << 214 //G4double crystalGap_axial; 216 215 217 G4double sizeOfAirBox_DOI; << 216 G4double sizeOfAirBox_DOI; 218 G4double sizeOfAirBox_axial; << 217 G4double sizeOfAirBox_axial; 219 G4double sizeOfAirBox_tangential; << 218 G4double sizeOfAirBox_tangential; 220 219 221 220 222 G4double sizeOfBlockDetector_DOI; << 221 G4double sizeOfBlockDetector_DOI; 223 G4double sizeOfBlockDetector_axial; << 222 G4double sizeOfBlockDetector_axial; 224 G4double sizeOfBlockDetector_tangential; << 223 G4double sizeOfBlockDetector_tangential; 225 224 226 //G4double AluminumCoverThickness; << 225 //G4double AluminumCoverThickness; 227 226 228 227 229 //G4int numberOfPETDetector; << 228 //G4int numberOfPETDetector; 230 //G4int numberOfRings; << 229 //G4int numberOfRings; 231 230 232 231 233 //G4double scannerRadius; << 232 //G4double scannerRadius; 234 G4double thetaDetector; //The azimuthal angl << 233 G4double thetaDetector; //The azimuthal angle for arranging the detector in the PET ring 235 //G4double ringGap; << 234 //G4double ringGap; 236 G4int blockIndex; << 235 G4int blockIndex; 237 //G4int AlCase_Index; << 236 G4int AlCase_Index; 238 G4int crystalIndex; << 237 G4int crystalIndex; 239 238 240 //detector position << 239 //detector position 241 G4double detectorPositionX; << 240 G4double detectorPositionX; 242 G4double detectorPositionY; << 241 G4double detectorPositionY; 243 G4double detectorPositionZ; << 242 G4double detectorPositionZ; 244 243 245 //crystal position << 244 //crystal position 246 G4double crystalPositionX; << 245 G4double crystalPositionX; 247 G4double crystalPositionY; << 246 G4double crystalPositionY; 248 G4double crystalPositionZ; << 247 G4double crystalPositionZ; 249 248 250 249 251 250 252 G4ThreeVector phantomPosition; << 251 G4ThreeVector phantomPosition; 253 252 254 // << 253 // 255 G4double phantomRadius; << 254 G4double phantomRadius; 256 G4double phantomLength; << 255 G4double phantomLength; 257 256 258 //Phantom dimension for rectangular box (pla << 257 //Phantom dimension for rectangular box (placed for therapy study) 259 //G4double phantomSizeX, phantomSizeY, phant << 258 G4double phantomSizeX, phantomSizeY, phantomSizeZ; 260 259 261 //the following is to make the body phantom << 260 //the following is to make the body phantom 262 G4double yOffsetBodyPhantom; << 261 G4double yOffsetBodyPhantom; 263 G4double zOffsetBodyPhantom; << 262 G4double zOffsetBodyPhantom; 264 G4double lengthOfBodyPhantom; //Interior len << 263 G4double lengthOfBodyPhantom; //Interior length ( = 180m mm) + wallthickness (= 3mm) 265 G4double radiusOfBodyPhantom; << 264 G4double radiusOfBodyPhantom; 266 G4double wallThicknessOfBodyPhantom; << 265 G4double wallThicknessOfBodyPhantom; 267 G4double radiusOfLungPhantom; << 266 G4double radiusOfLungPhantom; 268 267 269 //Test phantom defnition. The phantom has th << 268 //Test phantom defnition. The phantom has the same as that that of NECR phantom except 270 G4double hieghtOfTestPhantom; << 269 G4double hieghtOfTestPhantom; 271 G4double diameterOfTestPhantom; << 270 G4double diameterOfTestPhantom; 272 271 273 //To the cylindrical phantom to make the bod << 272 //To the cylindrical phantom to make the body phantom 274 G4double radiusOfSmallcyl; << 273 G4double radiusOfSmallcyl; 275 G4double boxWidth; << 274 G4double boxWidth; 276 G4double boxHeight; << 275 G4double boxHeight; 277 276 278 277 279 //Image quality phantom for small animals << 278 //Image quality phantom for small animals 280 G4double waterPhantomRadius; << 279 G4double waterPhantomRadius; 281 G4double waterPhantomLength; << 280 G4double waterPhantomLength; 282 281 283 G4double rodPhantomLength; << 282 G4double rodPhantomLength; 284 G4double rodDiameter; << 283 G4double rodDiameter; 285 G4int numberOfRods; << 284 G4int numberOfRods; 286 285 287 //Declare position for the rod phantoms << 286 //Declare position for the rod phantoms 288 G4double rodPositionX, rodPositionY, rodPosi << 287 G4double rodPositionX, rodPositionY, rodPositionZ; >> 288 >> 289 //Declare position for cold region chanmbers >> 290 G4double chamberPositionX, chamberPositionY, chamberPositionZ; >> 291 G4double chamberPhantomLength; >> 292 G4double chamberDiameter; >> 293 G4double wallThicknessOfChamber; 289 294 290 //Declare position for cold region chanmbers << 291 G4double chamberPositionX, chamberPositionY, << 292 G4double chamberPhantomLength; << 293 G4double chamberDiameter; << 294 G4double wallThicknessOfChamber; << 295 295 296 //Declare the number of concentric (Al) slee << 297 G4int numOfSleeves; << 298 }; 296 }; 299 297 300 ////////////////////////////////////////////// 298 ///////////////////////////////////////////////////////////////////////////////////////// 301 299 302 #endif 300 #endif 303 301