Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer 3 // * License and Disclaimer * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. 9 // * include a list of copyright holders. * 10 // * 10 // * * 11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. * 17 // * 17 // * * 18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************* 24 // ******************************************************************** 25 // 25 // 26 /// \file electromagnetic/TestEm12/src/Detecto 26 /// \file electromagnetic/TestEm12/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConst 27 /// \brief Implementation of the DetectorConstruction class 28 // 28 // >> 29 // $Id: DetectorConstruction.cc 78723 2014-01-20 10:32:17Z gcosmo $ 29 // 30 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 33 33 #include "DetectorConstruction.hh" 34 #include "DetectorConstruction.hh" 34 << 35 #include "DetectorMessenger.hh" 35 #include "DetectorMessenger.hh" 36 36 37 #include "G4GeometryManager.hh" << 38 #include "G4LogicalVolume.hh" << 39 #include "G4LogicalVolumeStore.hh" << 40 #include "G4NistManager.hh" 37 #include "G4NistManager.hh" >> 38 #include "G4Sphere.hh" >> 39 #include "G4LogicalVolume.hh" >> 40 #include "G4VPhysicalVolume.hh" 41 #include "G4PVPlacement.hh" 41 #include "G4PVPlacement.hh" 42 #include "G4PVReplica.hh" 42 #include "G4PVReplica.hh" 43 #include "G4PhysicalConstants.hh" << 43 >> 44 #include "G4GeometryManager.hh" 44 #include "G4PhysicalVolumeStore.hh" 45 #include "G4PhysicalVolumeStore.hh" 45 #include "G4RunManager.hh" << 46 #include "G4LogicalVolumeStore.hh" 46 #include "G4SolidStore.hh" 47 #include "G4SolidStore.hh" 47 #include "G4Sphere.hh" << 48 #include "G4RunManager.hh" 48 #include "G4SystemOfUnits.hh" << 49 49 #include "G4UnitsTable.hh" 50 #include "G4UnitsTable.hh" 50 #include "G4VPhysicalVolume.hh" << 51 #include "G4PhysicalConstants.hh" >> 52 #include "G4SystemOfUnits.hh" 51 53 52 //....oooOO0OOooo........oooOO0OOooo........oo 54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 55 54 DetectorConstruction::DetectorConstruction() 56 DetectorConstruction::DetectorConstruction() >> 57 : G4VUserDetectorConstruction(), >> 58 fAbsorMaterial(0), >> 59 fAbsor(0), >> 60 fDetectorMessenger(0) 55 { 61 { 56 // default parameter values 62 // default parameter values 57 fAbsorRadius = 3 * cm; << 63 fAbsorRadius = 3*cm; 58 fNbOfLayers = 1; 64 fNbOfLayers = 1; 59 << 65 60 DefineMaterials(); 66 DefineMaterials(); 61 SetMaterial("G4_WATER"); 67 SetMaterial("G4_WATER"); 62 68 63 // create commands for interactive definitio << 69 // create commands for interactive definition of the detector 64 fDetectorMessenger = new DetectorMessenger(t 70 fDetectorMessenger = new DetectorMessenger(this); 65 } 71 } 66 72 67 //....oooOO0OOooo........oooOO0OOooo........oo 73 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 74 69 DetectorConstruction::~DetectorConstruction() 75 DetectorConstruction::~DetectorConstruction() >> 76 { delete fDetectorMessenger;} >> 77 >> 78 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 79 >> 80 G4VPhysicalVolume* DetectorConstruction::Construct() 70 { 81 { 71 delete fDetectorMessenger; << 82 return ConstructVolumes(); 72 } 83 } 73 84 74 //....oooOO0OOooo........oooOO0OOooo........oo 85 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 75 86 76 void DetectorConstruction::DefineMaterials() 87 void DetectorConstruction::DefineMaterials() 77 { << 88 { 78 G4NistManager* man = G4NistManager::Instance 89 G4NistManager* man = G4NistManager::Instance(); 79 << 90 80 man->FindOrBuildMaterial("G4_Al"); 91 man->FindOrBuildMaterial("G4_Al"); 81 man->FindOrBuildMaterial("G4_Si"); 92 man->FindOrBuildMaterial("G4_Si"); 82 man->FindOrBuildMaterial("G4_Fe"); 93 man->FindOrBuildMaterial("G4_Fe"); 83 man->FindOrBuildMaterial("G4_Ge"); 94 man->FindOrBuildMaterial("G4_Ge"); 84 man->FindOrBuildMaterial("G4_Gd"); << 85 man->FindOrBuildMaterial("G4_W"); 95 man->FindOrBuildMaterial("G4_W"); 86 man->FindOrBuildMaterial("G4_Pb"); 96 man->FindOrBuildMaterial("G4_Pb"); 87 << 97 88 man->FindOrBuildMaterial("G4_AIR"); 98 man->FindOrBuildMaterial("G4_AIR"); 89 man->FindOrBuildMaterial("G4_WATER"); 99 man->FindOrBuildMaterial("G4_WATER"); 90 man->FindOrBuildMaterial("G4_ALUMINUM_OXIDE" << 100 91 << 101 ///G4cout << *(G4Material::GetMaterialTable()) << G4endl; 92 /// G4cout << *(G4Material::GetMaterialTable << 93 } 102 } 94 103 95 //....oooOO0OOooo........oooOO0OOooo........oo 104 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 96 << 105 97 G4VPhysicalVolume* DetectorConstruction::Const << 106 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes() 98 { 107 { >> 108 G4GeometryManager::GetInstance()->OpenGeometry(); >> 109 G4PhysicalVolumeStore::GetInstance()->Clean(); >> 110 G4LogicalVolumeStore::GetInstance()->Clean(); >> 111 G4SolidStore::GetInstance()->Clean(); >> 112 99 // Absorber 113 // Absorber 100 // 114 // 101 G4Sphere* sAbsor = new G4Sphere("Absorber", << 115 G4Sphere* 102 0., fAbsorRa << 116 sAbsor = new G4Sphere("Absorber", //name 103 << 117 0., fAbsorRadius, 0., twopi, 0., pi); //size 104 fSpheres.push_back(sAbsor); << 118 105 << 119 G4LogicalVolume* 106 G4LogicalVolume* lAbsor = new G4LogicalVolum << 120 lAbsor = new G4LogicalVolume(sAbsor, //solid 107 << 121 fAbsorMaterial, //material 108 << 122 "Absorber"); //name 109 fLVolumes.push_back(lAbsor); << 123 110 << 124 fAbsor = new G4PVPlacement(0, //no rotation 111 fAbsor = new G4PVPlacement(0, // no rotatio << 125 G4ThreeVector(), //at (0,0,0) 112 G4ThreeVector(), << 126 lAbsor, //logical volume 113 lAbsor, // logic << 127 "Absorber", //name 114 "Absorber", // n << 128 0, //mother volume 115 0, // mother vo << 129 false, //no boolean operation 116 false, // no boo << 130 0); //copy number 117 0); // copy numb << 118 131 119 // Layers 132 // Layers 120 // 133 // 121 fLayerThickness = fAbsorRadius / fNbOfLayers << 134 fLayerThickness = fAbsorRadius/fNbOfLayers; 122 << 135 123 for (G4int i = 1; i <= fNbOfLayers; i++) { << 136 for (G4int i=1; i<=fNbOfLayers; i++) { 124 G4Sphere* sLayer = << 137 G4Sphere* 125 new G4Sphere("Layer", (i - 1) * fLayerTh << 138 sLayer = new G4Sphere("Layer", (i-1)*fLayerThickness, i*fLayerThickness, 126 << 139 0., twopi, 0., pi); 127 fSpheres.push_back(sLayer); << 140 128 << 141 G4LogicalVolume* 129 G4LogicalVolume* lLayer = new G4LogicalVol << 142 lLayer = new G4LogicalVolume(sLayer, //shape 130 << 143 fAbsorMaterial, //material 131 << 144 "Layer"); //name 132 fLVolumes.push_back(lLayer); << 145 133 << 146 new G4PVPlacement(0, //no rotation 134 new G4PVPlacement(0, // no rotation << 147 G4ThreeVector(), //at (0,0,0) 135 G4ThreeVector(), // at << 148 lLayer, //logical volume 136 lLayer, // logical volu << 149 "Layer", //name 137 "Layer", // name << 150 lAbsor, //mother volume 138 lAbsor, // mother volu << 151 false, //no boolean operation 139 false, // no boolean op << 152 i); //copy number 140 i); // copy number << 153 141 } << 154 } 142 155 143 PrintParameters(); 156 PrintParameters(); 144 << 157 145 // << 146 // always return the root volume << 147 // 158 // >> 159 //always return the root volume >> 160 // 148 return fAbsor; 161 return fAbsor; 149 } 162 } 150 163 151 //....oooOO0OOooo........oooOO0OOooo........oo 164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 152 165 153 void DetectorConstruction::PrintParameters() 166 void DetectorConstruction::PrintParameters() 154 { 167 { 155 G4cout << "\n------------------------------- 168 G4cout << "\n---------------------------------------------------------\n"; 156 G4cout << "---> The Absorber is a sphere of << 169 G4cout << "---> The Absorber is a sphere of " 157 << " radius of " << fAbsorMaterial->G << 170 << G4BestUnit(fAbsorRadius,"Length") << " radius of " 158 << " slices of " << G4BestUnit(fLayer << 171 << fAbsorMaterial->GetName() << " divided in " << fNbOfLayers 159 << fAbsorMaterial << G4endl; << 172 << " slices of " << G4BestUnit(fLayerThickness,"Length") >> 173 << "\n \n" << fAbsorMaterial << G4endl; 160 G4cout << "\n------------------------------- 174 G4cout << "\n---------------------------------------------------------\n"; 161 } 175 } 162 176 163 //....oooOO0OOooo........oooOO0OOooo........oo 177 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 164 178 165 void DetectorConstruction::SetRadius(G4double 179 void DetectorConstruction::SetRadius(G4double value) 166 { 180 { 167 // geometry was already constructed - scale << 168 if (fAbsor) { << 169 G4double scale = value / fAbsorRadius; << 170 for (auto solid : fSpheres) { << 171 if (scale > 1.0) { << 172 solid->SetOuterRadius(solid->GetOuterR << 173 solid->SetInnerRadius(solid->GetInnerR << 174 } << 175 else { << 176 solid->SetInnerRadius(solid->GetInnerR << 177 solid->SetOuterRadius(solid->GetOuterR << 178 } << 179 } << 180 } << 181 fAbsorRadius = value; 181 fAbsorRadius = value; >> 182 G4RunManager::GetRunManager()->ReinitializeGeometry(); 182 } 183 } 183 184 184 //....oooOO0OOooo........oooOO0OOooo........oo 185 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 185 186 186 void DetectorConstruction::SetMaterial(G4Strin 187 void DetectorConstruction::SetMaterial(G4String materialChoice) 187 { 188 { 188 // search the material by its name << 189 // search the material by its name 189 G4Material* pttoMaterial = G4Material::GetMa << 190 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); 190 << 191 if (pttoMaterial) fAbsorMaterial = pttoMaterial; 191 if (pttoMaterial && pttoMaterial != fAbsorMa << 192 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 192 fAbsorMaterial = pttoMaterial; << 193 G4RunManager::GetRunManager()->PhysicsHasB << 194 << 195 // geometry was already constructed - only << 196 if (fAbsor) { << 197 for (auto lv : fLVolumes) { << 198 lv->SetMaterial(fAbsorMaterial); << 199 } << 200 } << 201 } << 202 } 193 } 203 194 204 //....oooOO0OOooo........oooOO0OOooo........oo 195 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 205 196 206 void DetectorConstruction::SetNbOfLayers(G4int 197 void DetectorConstruction::SetNbOfLayers(G4int value) 207 { 198 { 208 fNbOfLayers = value; 199 fNbOfLayers = value; >> 200 G4RunManager::GetRunManager()->ReinitializeGeometry(); 209 } 201 } 210 << 202 211 //....oooOO0OOooo........oooOO0OOooo........oo 203 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 212 204