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/TestEm16/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 "G4GDMLParser.hh" 39 #include "G4GeometryManager.hh" 40 #include "G4LogicalVolume.hh" 41 #include "G4LogicalVolumeStore.hh" 42 #include "G4Material.hh" 43 #include "G4PVPlacement.hh" 44 #include "G4PhysicalVolumeStore.hh" 45 #include "G4PropagatorInField.hh" 46 #include "G4RunManager.hh" 47 #include "G4SolidStore.hh" 48 #include "G4SystemOfUnits.hh" 49 #include "G4TransportationManager.hh" 50 #include "G4UnitsTable.hh" 51 #include "G4UserLimits.hh" 52 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 54 55 DetectorConstruction::DetectorConstruction() : G4VUserDetectorConstruction(), fBoxSize(500 * m) 56 { 57 DefineMaterials(); 58 SetMaterial("Iron"); 59 60 // create UserLimits 61 fUserLimits = new G4UserLimits(); 62 63 // create commands for interactive definition of the detector 64 fDetectorMessenger = new DetectorMessenger(this); 65 } 66 67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 69 G4VPhysicalVolume* DetectorConstruction::Construct() 70 { 71 return ConstructVolumes(); 72 } 73 74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 75 76 void DetectorConstruction::DefineMaterials() 77 { 78 G4double a, z, density; 79 80 new G4Material("Beryllium", z = 4., a = 9.012182 * g / mole, density = 1.848 * g / cm3); 81 new G4Material("Carbon", z = 6., a = 12.011 * g / mole, density = 2.265 * g / cm3); 82 new G4Material("Iron", z = 26., a = 55.85 * g / mole, density = 7.870 * g / cm3); 83 84 // define a vacuum with a restgas pressure typical for accelerators 85 G4double const Torr = atmosphere / 760.; // 1 Torr 86 G4double pressure = 10e-9 * Torr, 87 temperature = 296.150 * kelvin; // 23 Celsius 88 new G4Material("Vacuum", z = 7., a = 14.01 * g / mole, density = 1.516784e-11 * kg / m3, 89 kStateGas, temperature, pressure); 90 G4cout << *(G4Material::GetMaterialTable()) << G4endl; // print the materials defined 91 } 92 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 94 95 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes() 96 { 97 G4GeometryManager::GetInstance()->OpenGeometry(); 98 G4PhysicalVolumeStore::GetInstance()->Clean(); 99 G4LogicalVolumeStore::GetInstance()->Clean(); 100 G4SolidStore::GetInstance()->Clean(); 101 102 if (fGeomFileName == G4String()) { 103 auto sBox = new G4Box("Container", fBoxSize / 2, fBoxSize / 2, fBoxSize / 2); 104 fLBox = new G4LogicalVolume(sBox, fMaterial, fMaterial->GetName()); 105 fLBox->SetUserLimits(fUserLimits); 106 fBox = new G4PVPlacement(0, G4ThreeVector(), fLBox, fMaterial->GetName(), 0, false, 0); 107 PrintParameters(); 108 } 109 else { 110 G4GDMLParser parser; 111 std::size_t nmat = G4Material::GetNumberOfMaterials(); 112 parser.Read(fGeomFileName.c_str()); // volumes and materials 113 if (G4Material::GetNumberOfMaterials() > nmat) { 114 const std::vector<G4Material*> MatPtrVec = *G4Material::GetMaterialTable(); 115 if (G4Material::GetNumberOfMaterials() > nmat) 116 G4cout << "Materials defined by " << fGeomFileName << " :" << G4endl; 117 for (std::size_t imat = nmat; imat < MatPtrVec.size(); ++imat) 118 G4cout << MatPtrVec[imat] << G4endl; 119 } 120 fBox = parser.GetWorldVolume(); 121 } 122 return fBox; 123 } 124 125 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 126 127 void DetectorConstruction::PrintParameters() 128 { 129 G4cout << "\n The Box is " << G4BestUnit(fBoxSize, "Length") << " of " << fMaterial->GetName() 130 << G4endl; 131 } 132 133 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 134 135 void DetectorConstruction::SetMaterial(G4String materialChoice) 136 { 137 // search the material by its name 138 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); 139 if (pttoMaterial) { 140 fMaterial = pttoMaterial; 141 if (fLBox) fLBox->SetMaterial(fMaterial); 142 } 143 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 144 } 145 146 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 147 148 void DetectorConstruction::SetSize(G4double value) 149 { 150 fBoxSize = value; 151 G4RunManager::GetRunManager()->ReinitializeGeometry(); 152 } 153 154 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 155 156 #include "G4AutoDelete.hh" 157 #include "G4GlobalMagFieldMessenger.hh" 158 159 void DetectorConstruction::ConstructSDandField() 160 { 161 if (fFieldMessenger.Get() == 0) { 162 // Create global magnetic field messenger. 163 // Uniform magnetic field is then created automatically if 164 // the field value is not zero. 165 G4ThreeVector fieldValue = G4ThreeVector(); 166 auto msg = new G4GlobalMagFieldMessenger(fieldValue); 167 // msg->SetVerboseLevel(1); 168 G4AutoDelete::Register(msg); 169 fFieldMessenger.Put(msg); 170 } 171 } 172 173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 174 175 void DetectorConstruction::SetMaxStepSize(G4double val) 176 { 177 // set the maximum allowed step size 178 // 179 if (val <= DBL_MIN) { 180 G4cout << "\n --->warning from SetMaxStepSize: maxStep " << val 181 << " out of range. Command refused" << G4endl; 182 return; 183 } 184 fUserLimits->SetMaxAllowedStep(val); 185 } 186 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 188 189 void DetectorConstruction::SetMaxStepLength(G4double val) 190 { 191 // set the maximum length of tracking step 192 // 193 if (val <= DBL_MIN) { 194 G4cout << "\n --->warning from SetMaxStepLength: maxStep " << val 195 << " out of range. Command refused" << G4endl; 196 return; 197 } 198 G4TransportationManager* tmanager = G4TransportationManager::GetTransportationManager(); 199 tmanager->GetPropagatorInField()->SetLargestAcceptableStep(val); 200 } 201 202 void DetectorConstruction::SetGeomFileName(G4String GeomFileName) 203 { 204 fGeomFileName = GeomFileName; 205 } 206 207 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 208