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 electromagnetic/TestEm13/src/Detecto 27 /// \brief Implementation of the DetectorConst 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 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 "G4RunManager.hh" 46 #include "G4SolidStore.hh" 47 #include "G4SystemOfUnits.hh" 48 #include "G4UnitsTable.hh" 49 50 //....oooOO0OOooo........oooOO0OOooo........oo 51 52 DetectorConstruction::DetectorConstruction() 53 { 54 fBoxSize = 1 * cm; 55 DefineMaterials(); 56 SetMaterial("Water"); 57 fDetectorMessenger = new DetectorMessenger(t 58 } 59 60 //....oooOO0OOooo........oooOO0OOooo........oo 61 62 DetectorConstruction::~DetectorConstruction() 63 { 64 delete fDetectorMessenger; 65 } 66 67 //....oooOO0OOooo........oooOO0OOooo........oo 68 69 G4VPhysicalVolume* DetectorConstruction::Const 70 { 71 return ConstructVolumes(); 72 } 73 74 //....oooOO0OOooo........oooOO0OOooo........oo 75 76 void DetectorConstruction::DefineMaterials() 77 { 78 // 79 // define Elements 80 // 81 G4double z, a; 82 83 G4Element* H = new G4Element("Hydrogen", "H" 84 G4Element* N = new G4Element("Nitrogen", "N" 85 G4Element* O = new G4Element("Oxygen", "O", 86 G4Element* Na = new G4Element("Sodium", "Na" 87 G4Element* Ge = new G4Element("Germanium", " 88 G4Element* I = new G4Element("Iodine", "I", 89 G4Element* Bi = new G4Element("Bismuth", "Bi 90 91 // 92 // define materials 93 // 94 G4double density; 95 G4int ncomponents, natoms; 96 G4double fractionmass; 97 98 G4Material* Air = new G4Material("Air", dens 99 Air->AddElement(N, fractionmass = 70. * perC 100 Air->AddElement(O, fractionmass = 30. * perC 101 102 G4Material* H2l = new G4Material("H2liquid", 103 H2l->AddElement(H, fractionmass = 1.); 104 105 G4Material* H2O = new G4Material("Water", de 106 H2O->AddElement(H, natoms = 2); 107 H2O->AddElement(O, natoms = 1); 108 H2O->SetChemicalFormula("H_2O"); 109 H2O->GetIonisation()->SetMeanExcitationEnerg 110 111 new G4Material("liquidArgon", z = 18., a = 3 112 113 new G4Material("Carbon", z = 6., a = 12.01 * 114 115 new G4Material("Aluminium", z = 13., a = 26. 116 117 new G4Material("Silicon", z = 14., a = 28.09 118 119 new G4Material("Germanium", z = 32., a = 72. 120 121 G4Material* NaI = new G4Material("NaI", dens 122 NaI->AddElement(Na, natoms = 1); 123 NaI->AddElement(I, natoms = 1); 124 NaI->GetIonisation()->SetMeanExcitationEnerg 125 126 G4Material* Iod = new G4Material("Iodine", d 127 Iod->AddElement(I, natoms = 1); 128 129 G4Material* BGO = new G4Material("BGO", dens 130 BGO->AddElement(O, natoms = 12); 131 BGO->AddElement(Ge, natoms = 3); 132 BGO->AddElement(Bi, natoms = 4); 133 134 new G4Material("Iron", z = 26., a = 55.85 * 135 136 new G4Material("Tungsten", z = 74., a = 183. 137 138 new G4Material("Lead", z = 82., a = 207.19 * 139 140 new G4Material("Uranium", z = 92., a = 238.0 141 142 /// G4cout << *(G4Material::GetMaterialTable 143 } 144 145 //....oooOO0OOooo........oooOO0OOooo........oo 146 147 G4VPhysicalVolume* DetectorConstruction::Const 148 { 149 // Cleanup old geometry 150 G4GeometryManager::GetInstance()->OpenGeomet 151 G4PhysicalVolumeStore::GetInstance()->Clean( 152 G4LogicalVolumeStore::GetInstance()->Clean() 153 G4SolidStore::GetInstance()->Clean(); 154 155 G4Box* sBox = new G4Box("Container", // its 156 fBoxSize / 2, fBoxSi 157 158 fLBox = new G4LogicalVolume(sBox, // its sh 159 fMaterial, // i 160 fMaterial->GetNa 161 162 fPBox = new G4PVPlacement(0, // no rotation 163 G4ThreeVector(), 164 fLBox, // its log 165 fMaterial->GetName 166 0, // its mother 167 false, // no bool 168 0); // copy numbe 169 170 PrintParameters(); 171 172 // always return the root volume 173 // 174 return fPBox; 175 } 176 177 //....oooOO0OOooo........oooOO0OOooo........oo 178 179 void DetectorConstruction::PrintParameters() 180 { 181 G4cout << "\n The Box is " << G4BestUnit(fBo 182 << G4endl; 183 G4cout << "\n" << fMaterial << G4endl; 184 } 185 186 //....oooOO0OOooo........oooOO0OOooo........oo 187 188 void DetectorConstruction::SetMaterial(G4Strin 189 { 190 // search the material by its name, or build 191 G4Material* pttoMaterial = G4NistManager::In 192 193 if (pttoMaterial) { 194 fMaterial = pttoMaterial; 195 G4RunManager::GetRunManager()->PhysicsHasB 196 } 197 else { 198 G4cout << "\n--> warning from DetectorCons 199 << " not found" << G4endl; 200 } 201 } 202 203 //....oooOO0OOooo........oooOO0OOooo........oo 204 205 void DetectorConstruction::SetSize(G4double va 206 { 207 fBoxSize = value; 208 G4RunManager::GetRunManager()->ReinitializeG 209 } 210 211 //....oooOO0OOooo........oooOO0OOooo........oo 212