Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 /// \file DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConst 28 // 29 // 30 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 34 #include "DetectorConstruction.hh" 35 36 #include "DetectorMessenger.hh" 37 38 #include "G4Box.hh" 39 #include "G4GeometryManager.hh" 40 #include "G4LogicalVolume.hh" 41 #include "G4LogicalVolumeStore.hh" 42 #include "G4Material.hh" 43 #include "G4NistManager.hh" 44 #include "G4PVPlacement.hh" 45 #include "G4PhysicalVolumeStore.hh" 46 #include "G4RunManager.hh" 47 #include "G4SolidStore.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "G4UnitsTable.hh" 50 51 //....oooOO0OOooo........oooOO0OOooo........oo 52 53 DetectorConstruction::DetectorConstruction() 54 { 55 fBoxSize = 10 * m; 56 DefineMaterials(); 57 SetMaterial("Molybdenum98"); 58 fDetectorMessenger = new DetectorMessenger(t 59 } 60 61 //....oooOO0OOooo........oooOO0OOooo........oo 62 63 DetectorConstruction::~DetectorConstruction() 64 { 65 delete fDetectorMessenger; 66 } 67 68 //....oooOO0OOooo........oooOO0OOooo........oo 69 70 G4VPhysicalVolume* DetectorConstruction::Const 71 { 72 return ConstructVolumes(); 73 } 74 75 //....oooOO0OOooo........oooOO0OOooo........oo 76 77 void DetectorConstruction::DefineMaterials() 78 { 79 // define a Material from isotopes 80 // 81 MaterialWithSingleIsotope("Molybdenum98", "M 82 83 // NE213 84 G4Element* H = new G4Element("Hydrogen", "H" 85 G4Element* C = new G4Element("Carbon", "C", 86 G4Material* ne213 = new G4Material("NE213", 87 ne213->AddElement(H, 9.2 * perCent); 88 ne213->AddElement(C, 90.8 * perCent); 89 90 G4Material* hydrogen = new G4Material("hydro 91 hydrogen->AddElement(H, 1); 92 93 G4Material* carbon = new G4Material("carbon" 94 carbon->AddElement(C, 1); 95 96 G4Material* plastic = new G4Material("plasti 97 plastic->AddElement(H, 1); 98 plastic->AddElement(C, 1); 99 100 // or use G4-NIST materials data base 101 // 102 G4NistManager* man = G4NistManager::Instance 103 man->FindOrBuildMaterial("G4_B"); 104 105 G4Element* O = man->FindOrBuildElement("O"); 106 G4Element* Hf = man->FindOrBuildElement("Hf" 107 108 G4Material* HfO2 = new G4Material("HfO2", 9. 109 HfO2->AddElement(Hf, 1); 110 HfO2->AddElement(O, 2); 111 112 /// G4cout << *(G4Material::GetMaterialTable 113 } 114 115 //....oooOO0OOooo........oooOO0OOooo........oo 116 117 G4Material* DetectorConstruction::MaterialWith 118 119 { 120 // define a material from an isotope 121 // 122 G4int ncomponents; 123 G4double abundance, massfraction; 124 125 G4Isotope* isotope = new G4Isotope(symbol, Z 126 127 G4Element* element = new G4Element(name, sym 128 element->AddIsotope(isotope, abundance = 100 129 130 G4Material* material = new G4Material(name, 131 material->AddElement(element, massfraction = 132 133 return material; 134 } 135 136 //....oooOO0OOooo........oooOO0OOooo........oo 137 138 G4VPhysicalVolume* DetectorConstruction::Const 139 { 140 // Cleanup old geometry 141 G4GeometryManager::GetInstance()->OpenGeomet 142 G4PhysicalVolumeStore::GetInstance()->Clean( 143 G4LogicalVolumeStore::GetInstance()->Clean() 144 G4SolidStore::GetInstance()->Clean(); 145 146 G4Box* sBox = new G4Box("Container", // its 147 fBoxSize / 2, fBoxSi 148 149 fLBox = new G4LogicalVolume(sBox, // its sh 150 fMaterial, // i 151 fMaterial->GetNa 152 153 fPBox = new G4PVPlacement(0, // no rotation 154 G4ThreeVector(), 155 fLBox, // its log 156 fMaterial->GetName 157 0, // its mother 158 false, // no bool 159 0); // copy numbe 160 161 PrintParameters(); 162 163 // always return the root volume 164 // 165 return fPBox; 166 } 167 168 //....oooOO0OOooo........oooOO0OOooo........oo 169 170 void DetectorConstruction::PrintParameters() 171 { 172 G4cout << "\n The Box is " << G4BestUnit(fBo 173 << "\n \n" 174 << fMaterial << G4endl; 175 } 176 177 //....oooOO0OOooo........oooOO0OOooo........oo 178 179 void DetectorConstruction::SetMaterial(G4Strin 180 { 181 // search the material by its name 182 G4Material* pttoMaterial = G4NistManager::In 183 184 if (pttoMaterial) { 185 if (fMaterial != pttoMaterial) { 186 fMaterial = pttoMaterial; 187 if (fLBox) { 188 fLBox->SetMaterial(pttoMaterial); 189 } 190 G4RunManager::GetRunManager()->PhysicsHa 191 } 192 } 193 else { 194 G4cout << "\n--> warning from DetectorCons 195 << " not found" << G4endl; 196 } 197 } 198 199 //....oooOO0OOooo........oooOO0OOooo........oo 200 201 void DetectorConstruction::SetSize(G4double va 202 { 203 fBoxSize = value; 204 G4RunManager::GetRunManager()->ReinitializeG 205 } 206 207 //....oooOO0OOooo........oooOO0OOooo........oo 208