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