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 /// \file polarisation/Pol01/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConstruction class 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 33 #include "DetectorConstruction.hh" 34 35 #include "DetectorMessenger.hh" 36 37 #include "G4Box.hh" 38 #include "G4GeometryManager.hh" 39 #include "G4LogicalVolume.hh" 40 #include "G4LogicalVolumeStore.hh" 41 #include "G4Material.hh" 42 #include "G4NistManager.hh" 43 #include "G4PVPlacement.hh" 44 #include "G4PhysicalVolumeStore.hh" 45 #include "G4PolarizationManager.hh" 46 #include "G4RunManager.hh" 47 #include "G4SolidStore.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "G4UnitsTable.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 53 DetectorConstruction::DetectorConstruction() 54 : G4VUserDetectorConstruction(), fWorld(0), fBox(0), fTargetMaterial(0), fWorldMaterial(0) 55 { 56 fBoxSizeXY = 50 * mm; 57 fBoxSizeZ = 5 * mm; 58 fWorldSize = 1. * m; 59 ConstructMaterials(); 60 fMessenger = new DetectorMessenger(this); 61 } 62 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 65 DetectorConstruction::~DetectorConstruction() 66 { 67 delete fMessenger; 68 } 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 70 71 void DetectorConstruction::ConstructMaterials() 72 { 73 auto nistManager = G4NistManager::Instance(); 74 75 fTargetMaterial = nistManager->FindOrBuildMaterial("G4_Fe"); 76 if (fTargetMaterial == nullptr) { 77 G4cerr << "### ERROR - Material: <" 78 << "G4_Fe" 79 << "> not found" << G4endl; 80 } 81 82 fWorldMaterial = nistManager->FindOrBuildMaterial("G4_Galactic"); 83 if (fWorldMaterial == nullptr) { 84 G4cerr << "### ERROR - Material: <" 85 << "G4_Galactic" 86 << "> not found" << G4endl; 87 } 88 } 89 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 91 92 G4LogicalVolume* lBox; 93 94 G4VPhysicalVolume* DetectorConstruction::Construct() 95 { 96 // G4PolarizationManager::GetInstance()->Clean(); 97 98 // World 99 // 100 G4Box* sWorld = new G4Box("World", // name 101 fWorldSize / 2, fWorldSize / 2, fWorldSize / 2); // dimensions 102 103 G4LogicalVolume* lWorld = new G4LogicalVolume(sWorld, // shape 104 fWorldMaterial, // material 105 "World"); // name 106 107 fWorld = new G4PVPlacement(0, // no rotation 108 G4ThreeVector(), // at (0,0,0) 109 lWorld, // logical volume 110 "World", // name 111 0, // mother volume 112 false, // no boolean operation 113 0); // copy number 114 115 // Box 116 // 117 G4Box* sBox = new G4Box("Container", // its name 118 fBoxSizeXY / 2., fBoxSizeXY / 2., fBoxSizeZ / 2.); // its dimensions 119 120 // G4LogicalVolume* 121 lBox = new G4LogicalVolume(sBox, // its shape 122 fTargetMaterial, // its material 123 "theBox"); // its name 124 125 fBox = new G4PVPlacement(0, // no rotation 126 G4ThreeVector(), // at (0,0,0) 127 lBox, // its logical volume 128 fTargetMaterial->GetName(), // its name 129 lWorld, // its mother volume 130 false, // no boolean operation 131 0); // copy number 132 133 PrintParameters(); 134 135 // always return the root volume 136 // 137 return fWorld; 138 } 139 140 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 141 142 void DetectorConstruction::ConstructSDandField() 143 { 144 // register logical Volume in PolarizationManager with zero polarization 145 G4PolarizationManager* polMgr = G4PolarizationManager::GetInstance(); 146 polMgr->SetVolumePolarization(lBox, G4ThreeVector(0., 0., 0.)); 147 } 148 149 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 150 151 void DetectorConstruction::PrintParameters() 152 { 153 G4cout << "\n The Box is " << G4BestUnit(fBoxSizeXY, "Length") << " x " 154 << G4BestUnit(fBoxSizeXY, "Length") << " x " << G4BestUnit(fBoxSizeZ, "Length") << " of " 155 << fTargetMaterial->GetName() << G4endl; 156 } 157 158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 159 160 void DetectorConstruction::SetTargetMaterial(G4String materialChoice) 161 { 162 // search the material by its name 163 G4Material* mat = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 164 if (mat != fTargetMaterial) { 165 if (mat) { 166 fTargetMaterial = mat; 167 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 168 UpdateGeometry(); 169 } 170 else { 171 G4cout << "### Warning! Target material: <" << materialChoice << "> not found" << G4endl; 172 } 173 } 174 } 175 176 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 177 178 void DetectorConstruction::SetWorldMaterial(G4String materialChoice) 179 { 180 // search the material by its name 181 G4Material* mat = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 182 if (mat != fWorldMaterial) { 183 if (mat) { 184 fWorldMaterial = mat; 185 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 186 UpdateGeometry(); 187 } 188 else { 189 G4cout << "### Warning! World material: <" << materialChoice << "> not found" << G4endl; 190 } 191 } 192 } 193 194 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 195 196 void DetectorConstruction::SetSizeXY(G4double value) 197 { 198 fBoxSizeXY = value; 199 if (fWorldSize < fBoxSizeXY) fWorldSize = 1.2 * fBoxSizeXY; 200 UpdateGeometry(); 201 } 202 203 void DetectorConstruction::SetSizeZ(G4double value) 204 { 205 fBoxSizeZ = value; 206 if (fWorldSize < fBoxSizeZ) fWorldSize = 1.2 * fBoxSizeZ; 207 UpdateGeometry(); 208 } 209 210 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 211 212 #include "G4RunManager.hh" 213 214 void DetectorConstruction::UpdateGeometry() 215 { 216 // if (fWorld) 217 // G4RunManager::GetRunManager()->DefineWorldVolume(Construct()); 218 G4RunManager::GetRunManager()->GeometryHasBeenModified(); 219 } 220 221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 222