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 electromagnetic/TestEm17/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConstruction class 28 // 29 // 30 // 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 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 "G4PhysicalConstants.hh" 46 #include "G4PhysicalVolumeStore.hh" 47 #include "G4RunManager.hh" 48 #include "G4SolidStore.hh" 49 #include "G4SystemOfUnits.hh" 50 #include "G4UnitsTable.hh" 51 52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 54 DetectorConstruction::DetectorConstruction() 55 : G4VUserDetectorConstruction(), fPBox(0), fLBox(0), fMaterial(0), fDetectorMessenger(0) 56 { 57 fBoxSize = 1 * m; 58 DefineMaterials(); 59 SetMaterial("G4_Fe"); 60 fDetectorMessenger = new DetectorMessenger(this); 61 } 62 63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 65 DetectorConstruction::~DetectorConstruction() 66 { 67 delete fDetectorMessenger; 68 } 69 70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 71 72 G4VPhysicalVolume* DetectorConstruction::Construct() 73 { 74 return ConstructVolumes(); 75 } 76 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 78 79 void DetectorConstruction::DefineMaterials() 80 { 81 // 82 // define Elements 83 // 84 G4double z, a; 85 86 G4Element* C = new G4Element("Carbon", "C", z = 6., a = 12.01 * g / mole); 87 G4Element* N = new G4Element("Nitrogen", "N", z = 7., a = 14.01 * g / mole); 88 G4Element* O = new G4Element("Oxygen", "O", z = 8., a = 16.00 * g / mole); 89 G4Element* Ca = new G4Element("Calcium", "Ca", z = 20., a = 40.08 * g / mole); 90 G4Element* H = new G4Element("Hydrogen", "H", z = 1., a = 1.01 * g / mole); 91 G4Element* I = new G4Element("Iodine", "I", z = 53., a = 126.90 * g / mole); 92 93 // 94 // define materials 95 // 96 G4double density; 97 G4int ncomponents, natoms; 98 G4double fractionmass; 99 100 new G4Material("galactic", z = 1., a = 1.01 * g / mole, universe_mean_density, kStateGas, 101 2.73 * kelvin, 3.e-18 * pascal); 102 103 G4Material* Air = new G4Material("Air", density = 1.290 * mg / cm3, ncomponents = 2); 104 Air->AddElement(N, fractionmass = 70. * perCent); 105 Air->AddElement(O, fractionmass = 30. * perCent); 106 107 G4Material* CaCO3 = new G4Material("CaCO3", density = 2.80 * g / cm3, ncomponents = 3); 108 CaCO3->AddElement(Ca, natoms = 1); 109 CaCO3->AddElement(C, natoms = 1); 110 CaCO3->AddElement(O, natoms = 3); 111 112 new G4Material("Carbon", z = 6., a = 12.01 * g / mole, density = 2.265 * g / cm3); 113 new G4Material("Iron", z = 26., a = 55.85 * g / mole, density = 7.870 * g / cm3); 114 new G4Material("Tin", z = 50., a = 118.7 * g / mole, density = 7.310 * g / cm3); 115 116 G4Material* HI = new G4Material("HI", density = 2.8 * g / cm3, ncomponents = 2); //-35oC 117 HI->AddElement(H, natoms = 1); 118 HI->AddElement(I, natoms = 1); 119 120 G4cout << *(G4Material::GetMaterialTable()) << G4endl; 121 } 122 123 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 124 125 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes() 126 { 127 if (fPBox) { 128 return fPBox; 129 } 130 131 G4Box* sBox = new G4Box("Container", // its name 132 fBoxSize / 2, fBoxSize / 2, fBoxSize / 2); // its dimensions 133 134 fLBox = new G4LogicalVolume(sBox, // its shape 135 fMaterial, // its material 136 fMaterial->GetName()); // its name 137 138 fPBox = new G4PVPlacement(0, // no rotation 139 G4ThreeVector(), // at (0,0,0) 140 fLBox, // its logical volume 141 fMaterial->GetName(), // its name 142 0, // its mother volume 143 false, // no boolean operation 144 0); // copy number 145 146 PrintParameters(); 147 148 // always return the root volume 149 // 150 return fPBox; 151 } 152 153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 154 155 void DetectorConstruction::PrintParameters() 156 { 157 G4cout << "\n The Box is " << G4BestUnit(fBoxSize, "Length") << " of " << fMaterial->GetName() 158 << G4endl; 159 } 160 161 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 162 163 void DetectorConstruction::SetMaterial(const G4String& materialChoice) 164 { 165 // search the material by its name 166 G4Material* pttoMaterial = G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 167 168 if (pttoMaterial && fMaterial != pttoMaterial) { 169 fMaterial = pttoMaterial; 170 if (fLBox) { 171 fLBox->SetMaterial(fMaterial); 172 } 173 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 174 } 175 } 176 177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 178 179 void DetectorConstruction::SetSize(G4double value) 180 { 181 fBoxSize = value; 182 } 183 184 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 185 186 void DetectorConstruction::UpdateGeometry() 187 { 188 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes()); 189 } 190 191 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 192