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 // This example is provided by the Geant4-DNA collaboration 27 // Any report or published results obtained using the Geant4-DNA software 28 // shall cite the following Geant4-DNA collaboration publications: 29 // Med. Phys. 45 (2018) e722-e739 30 // Phys. Med. 31 (2015) 861-874 31 // Med. Phys. 37 (2010) 4692-4708 32 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178 33 // 34 // The Geant4-DNA web site is available at http://geant4-dna.org 35 // 36 /// \file medical/dna/svalue/src/DetectorConstruction.cc 37 /// \brief Implementation of the DetectorConstruction class 38 39 #include "DetectorConstruction.hh" 40 41 #include "DetectorMessenger.hh" 42 43 #include "G4NistManager.hh" 44 #include "G4PVPlacement.hh" 45 #include "G4PhysicalConstants.hh" 46 #include "G4RunManager.hh" 47 #include "G4Sphere.hh" 48 #include "G4UnitsTable.hh" 49 #include "G4UserLimits.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 53 DetectorConstruction::DetectorConstruction() : G4VUserDetectorConstruction() 54 { 55 // Default tracking cut 56 57 fTrackingCut = 7.4 * CLHEP::eV; 58 59 // Default parameter values 60 61 fWorldRadius = 10 * CLHEP::m; 62 fCytoThickness = 20 * CLHEP::nm; 63 fNuclRadius = 10 * CLHEP::nm; 64 65 DefineMaterials(); 66 67 // Create commands for interactive definition of the detector 68 69 fDetectorMessenger = new DetectorMessenger(this); 70 } 71 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 73 74 DetectorConstruction::~DetectorConstruction() 75 { 76 delete fDetectorMessenger; 77 } 78 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 80 81 void DetectorConstruction::DefineMaterials() 82 { 83 G4NistManager* man = G4NistManager::Instance(); 84 85 fWorldMaterial = man->FindOrBuildMaterial("G4_WATER"); 86 fCytoMaterial = fNuclMaterial = man->FindOrBuildMaterial("G4_WATER"); 87 } 88 89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 90 91 G4VPhysicalVolume* DetectorConstruction::Construct() 92 { 93 if (fWorld) return fWorld; 94 95 // Spherical world 96 97 G4Sphere* sWorld = new G4Sphere("World", 0., 1000 * fNuclRadius, 0., twopi, 0., pi); 98 99 fLogicalWorld = new G4LogicalVolume(sWorld, fWorldMaterial, "World"); 100 101 fWorld = new G4PVPlacement(0, G4ThreeVector(), fLogicalWorld, "World", 0, false, 0); 102 103 // Spherical nucleus 104 105 G4Sphere* sNucl = new G4Sphere("Nucl", 0., fNuclRadius, 0., twopi, 0., pi); 106 107 fLogicalNucl = new G4LogicalVolume(sNucl, fNuclMaterial, "Nucl"); 108 109 fNucl = new G4PVPlacement(0, G4ThreeVector(), "Nucl", fLogicalNucl, fWorld, false, 0); 110 111 // Spherical shell for cytoplasm 112 113 G4Sphere* sCyto = 114 new G4Sphere("Cyto", fNuclRadius, fNuclRadius + fCytoThickness, 0., twopi, 0., pi); 115 116 fLogicalCyto = new G4LogicalVolume(sCyto, fCytoMaterial, "Cyto"); 117 118 fCyto = new G4PVPlacement(0, G4ThreeVector(), "Cyto", fLogicalCyto, fWorld, false, 0); 119 // 120 121 PrintParameters(); 122 123 // Tracking cut 124 125 fLogicalNucl->SetUserLimits(new G4UserLimits(DBL_MAX, DBL_MAX, DBL_MAX, fTrackingCut)); 126 fLogicalCyto->SetUserLimits(new G4UserLimits(DBL_MAX, DBL_MAX, DBL_MAX, fTrackingCut)); 127 fLogicalWorld->SetUserLimits(new G4UserLimits(DBL_MAX, DBL_MAX, DBL_MAX, fTrackingCut)); 128 129 // 130 return fWorld; 131 } 132 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 135 void DetectorConstruction::PrintParameters() const 136 { 137 G4cout << "\n---------------------------------------------------------\n"; 138 G4cout << "---> The tracking cut is set to " << G4BestUnit(fTrackingCut, "Energy") << G4endl; 139 G4cout << "---> The World is a sphere of " << G4BestUnit(1000 * fNuclRadius, "Length") 140 << "radius of " << fWorldMaterial->GetName() << G4endl; 141 G4cout << "---> The Nucleus is a sphere of " << G4BestUnit(fNuclRadius, "Length") << "radius of " 142 << fWorldMaterial->GetName() << " of mass " << G4BestUnit(GetNuclMass(), "Mass") << G4endl; 143 G4cout << "---> The Cytoplasm is a spherical shell of thickness " 144 << G4BestUnit(fCytoThickness, "Length") << "of " << fWorldMaterial->GetName() 145 << " of mass " << G4BestUnit(GetCytoMass(), "Mass") << G4endl; 146 G4cout << "\n---------------------------------------------------------\n"; 147 } 148 149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 150 151 void DetectorConstruction::SetTrackingCut(G4double value) 152 { 153 fTrackingCut = value; 154 } 155 156 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 157 158 void DetectorConstruction::SetNuclRadius(G4double value) 159 { 160 fNuclRadius = value; 161 } 162 163 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 164 165 void DetectorConstruction::SetCytoThickness(G4double value) 166 { 167 fCytoThickness = value; 168 } 169 170 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 171 172 void DetectorConstruction::SetWorldMaterial(const G4String& materialChoice) 173 { 174 // Search the material by its name 175 G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 176 177 if (pttoMaterial && pttoMaterial != fWorldMaterial) { 178 fWorldMaterial = pttoMaterial; 179 if (fLogicalWorld) fLogicalWorld->SetMaterial(pttoMaterial); 180 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 181 } 182 } 183 184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 185 186 void DetectorConstruction::SetCytoMaterial(const G4String& materialChoice) 187 { 188 // Search the material by its name 189 G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 190 191 if (pttoMaterial && pttoMaterial != fCytoMaterial) { 192 fCytoMaterial = pttoMaterial; 193 if (fLogicalCyto) fLogicalCyto->SetMaterial(pttoMaterial); 194 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 195 } 196 } 197 198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 199 200 void DetectorConstruction::SetNuclMaterial(const G4String& materialChoice) 201 { 202 // Search the material by its name 203 G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 204 205 if (pttoMaterial && pttoMaterial != fNuclMaterial) { 206 fNuclMaterial = pttoMaterial; 207 if (fLogicalNucl) fLogicalNucl->SetMaterial(pttoMaterial); 208 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 209 } 210 } 211