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 // 27 /// \file DicomRegularDetectorConstruction.cc 28 /// \brief Implementation of the DicomRegularDetectorConstruction clas 29 // 30 // History: 31 // Pedro Arce 32 // 33 //******************************************************* 34 35 #include "DicomRegularDetectorConstruction.hh" 36 37 #include "DicomPhantomParameterisationColour.hh" 38 39 #include "G4Box.hh" 40 #include "G4Colour.hh" 41 #include "G4Element.hh" 42 #include "G4LogicalVolume.hh" 43 #include "G4Material.hh" 44 #include "G4PVParameterised.hh" 45 #include "G4PVPlacement.hh" 46 #include "G4VPhysicalVolume.hh" 47 #include "G4VisAttributes.hh" 48 #include "G4ios.hh" 49 #include "globals.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 DicomRegularDetectorConstruction::DicomRegularDetectorConstruction() : DicomDetectorConstruction() 53 {} 54 55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 56 DicomRegularDetectorConstruction::~DicomRegularDetectorConstruction() {} 57 58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 59 void DicomRegularDetectorConstruction::ConstructPhantom() 60 { 61 #ifdef G4VERBOSE 62 G4cout << "DicomRegularDetectorConstruction::ConstructPhantom " << G4endl; 63 #endif 64 65 //----- Create parameterisation 66 DicomPhantomParameterisationColour* param = new DicomPhantomParameterisationColour(); 67 68 //----- Set voxel dimensions 69 param->SetVoxelDimensions(fVoxelHalfDimX, fVoxelHalfDimY, fVoxelHalfDimZ); 70 71 //----- Set number of voxels 72 param->SetNoVoxels(fNoVoxelsX, fNoVoxelsY, fNoVoxelsZ); 73 74 //----- Set list of materials 75 param->SetMaterials(fMaterials); 76 77 //----- Set list of material indices: for each voxel it is a number that 78 // correspond to the index of its material in the vector of materials 79 // defined above 80 param->SetMaterialIndices(fMateIDs); 81 82 //----- Define voxel logical volume 83 G4Box* voxel_solid = new G4Box("Voxel", fVoxelHalfDimX, fVoxelHalfDimY, fVoxelHalfDimZ); 84 G4LogicalVolume* voxel_logic = 85 new G4LogicalVolume(voxel_solid, fMaterials[0], "VoxelLogical", 0, 0, 0); 86 // material is not relevant, it will be changed by the 87 // ComputeMaterial method of the parameterisation 88 89 voxel_logic->SetVisAttributes(new G4VisAttributes(G4VisAttributes::GetInvisible())); 90 91 //--- Assign the fContainer volume of the parameterisation 92 param->BuildContainerSolid(fContainer_phys); 93 94 //--- Assure yourself that the voxels are completely filling the 95 // fContainer volume 96 param->CheckVoxelsFillContainer(fContainer_solid->GetXHalfLength(), 97 fContainer_solid->GetYHalfLength(), 98 fContainer_solid->GetZHalfLength()); 99 100 //----- The G4PVParameterised object that uses the created parameterisation 101 // should be placed in the fContainer logical volume 102 G4PVParameterised* phantom_phys = new G4PVParameterised( 103 "phantom", voxel_logic, fContainer_logic, kXAxis, fNoVoxelsX * fNoVoxelsY * fNoVoxelsZ, param); 104 // if axis is set as kUndefined instead of kXAxis, GEANT4 will 105 // do an smart voxel optimisation 106 // (not needed if G4RegularNavigation is used) 107 108 //----- Set this physical volume as having a regular structure of type 1, 109 // so that G4RegularNavigation is used 110 phantom_phys->SetRegularStructureId(1); // if not set, G4VoxelNavigation 111 // will be used instead 112 113 SetScorer(voxel_logic); 114 } 115