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 // $Id: DetectorConstruction.cc,v 1.2 2006/06/29 16:43:00 gunter Exp $ 27 /// \brief Implementation of the DetectorConst << 27 // GEANT4 tag $Name: geant4-09-01 $ 28 // << 29 // 28 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 31 33 #include "DetectorConstruction.hh" 32 #include "DetectorConstruction.hh" 34 << 35 #include "DetectorMessenger.hh" 33 #include "DetectorMessenger.hh" 36 34 37 #include "G4GeometryManager.hh" << 38 #include "G4LogicalVolume.hh" << 39 #include "G4LogicalVolumeStore.hh" << 40 #include "G4NistManager.hh" 35 #include "G4NistManager.hh" >> 36 #include "G4Sphere.hh" >> 37 #include "G4LogicalVolume.hh" >> 38 #include "G4VPhysicalVolume.hh" 41 #include "G4PVPlacement.hh" 39 #include "G4PVPlacement.hh" 42 #include "G4PVReplica.hh" 40 #include "G4PVReplica.hh" 43 #include "G4PhysicalConstants.hh" << 41 #include "G4UniformMagField.hh" >> 42 >> 43 #include "G4GeometryManager.hh" 44 #include "G4PhysicalVolumeStore.hh" 44 #include "G4PhysicalVolumeStore.hh" 45 #include "G4RunManager.hh" << 45 #include "G4LogicalVolumeStore.hh" 46 #include "G4SolidStore.hh" 46 #include "G4SolidStore.hh" 47 #include "G4Sphere.hh" << 47 48 #include "G4SystemOfUnits.hh" << 49 #include "G4UnitsTable.hh" 48 #include "G4UnitsTable.hh" 50 #include "G4VPhysicalVolume.hh" << 51 49 52 //....oooOO0OOooo........oooOO0OOooo........oo 50 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 51 54 DetectorConstruction::DetectorConstruction() 52 DetectorConstruction::DetectorConstruction() 55 { 53 { 56 // default parameter values 54 // default parameter values 57 fAbsorRadius = 3 * cm; << 55 absorRadius = 3*cm; 58 fNbOfLayers = 1; << 56 nbOfLayers = 1; 59 << 57 >> 58 absorMaterial = 0; >> 59 magField = 0; >> 60 pAbsor = 0; >> 61 60 DefineMaterials(); 62 DefineMaterials(); 61 SetMaterial("G4_WATER"); 63 SetMaterial("G4_WATER"); 62 64 63 // create commands for interactive definitio << 65 // create commands for interactive definition of the detector 64 fDetectorMessenger = new DetectorMessenger(t << 66 detectorMessenger = new DetectorMessenger(this); 65 } 67 } 66 68 67 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 68 70 69 DetectorConstruction::~DetectorConstruction() 71 DetectorConstruction::~DetectorConstruction() >> 72 { delete detectorMessenger;} >> 73 >> 74 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 75 >> 76 G4VPhysicalVolume* DetectorConstruction::Construct() 70 { 77 { 71 delete fDetectorMessenger; << 78 return ConstructVolumes(); 72 } 79 } 73 80 74 //....oooOO0OOooo........oooOO0OOooo........oo 81 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 75 82 76 void DetectorConstruction::DefineMaterials() 83 void DetectorConstruction::DefineMaterials() 77 { << 84 { 78 G4NistManager* man = G4NistManager::Instance 85 G4NistManager* man = G4NistManager::Instance(); 79 << 86 80 man->FindOrBuildMaterial("G4_Al"); << 87 G4bool isotopes = false; 81 man->FindOrBuildMaterial("G4_Si"); << 88 82 man->FindOrBuildMaterial("G4_Fe"); << 89 man->FindOrBuildMaterial("G4_Al", isotopes); 83 man->FindOrBuildMaterial("G4_Ge"); << 90 man->FindOrBuildMaterial("G4_Si", isotopes); 84 man->FindOrBuildMaterial("G4_Gd"); << 91 man->FindOrBuildMaterial("G4_Fe", isotopes); 85 man->FindOrBuildMaterial("G4_W"); << 92 man->FindOrBuildMaterial("G4_Ge", isotopes); 86 man->FindOrBuildMaterial("G4_Pb"); << 93 man->FindOrBuildMaterial("G4_W" , isotopes); 87 << 94 man->FindOrBuildMaterial("G4_Pb", isotopes); 88 man->FindOrBuildMaterial("G4_AIR"); << 95 89 man->FindOrBuildMaterial("G4_WATER"); << 96 man->FindOrBuildMaterial("G4_AIR" , isotopes); 90 man->FindOrBuildMaterial("G4_ALUMINUM_OXIDE" << 97 man->FindOrBuildMaterial("G4_WATER", isotopes); 91 << 98 92 /// G4cout << *(G4Material::GetMaterialTable << 99 G4cout << *(G4Material::GetMaterialTable()) << G4endl; 93 } 100 } 94 101 95 //....oooOO0OOooo........oooOO0OOooo........oo 102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 96 << 103 97 G4VPhysicalVolume* DetectorConstruction::Const << 104 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes() 98 { << 105 { >> 106 G4GeometryManager::GetInstance()->OpenGeometry(); >> 107 G4PhysicalVolumeStore::GetInstance()->Clean(); >> 108 G4LogicalVolumeStore::GetInstance()->Clean(); >> 109 G4SolidStore::GetInstance()->Clean(); >> 110 99 // Absorber 111 // Absorber 100 // 112 // 101 G4Sphere* sAbsor = new G4Sphere("Absorber", << 113 G4Sphere* 102 0., fAbsorRa << 114 sAbsor = new G4Sphere("Absorber", //name 103 << 115 0., absorRadius, 0., twopi, 0., pi); //size 104 fSpheres.push_back(sAbsor); << 116 105 << 117 G4LogicalVolume* 106 G4LogicalVolume* lAbsor = new G4LogicalVolum << 118 lAbsor = new G4LogicalVolume(sAbsor, //solid 107 << 119 absorMaterial, //material 108 << 120 "Absorber"); //name 109 fLVolumes.push_back(lAbsor); << 121 110 << 122 pAbsor = new G4PVPlacement(0, //no rotation 111 fAbsor = new G4PVPlacement(0, // no rotatio << 123 G4ThreeVector(), //at (0,0,0) 112 G4ThreeVector(), << 124 lAbsor, //logical volume 113 lAbsor, // logic << 125 "Absorber", //name 114 "Absorber", // n << 126 0, //mother volume 115 0, // mother vo << 127 false, //no boolean operation 116 false, // no boo << 128 0); //copy number 117 0); // copy numb << 118 129 119 // Layers 130 // Layers 120 // 131 // 121 fLayerThickness = fAbsorRadius / fNbOfLayers << 132 layerThickness = absorRadius/nbOfLayers; 122 << 133 123 for (G4int i = 1; i <= fNbOfLayers; i++) { << 134 for (G4int i=1; i<=nbOfLayers; i++) { 124 G4Sphere* sLayer = << 135 G4Sphere* 125 new G4Sphere("Layer", (i - 1) * fLayerTh << 136 sLayer = new G4Sphere("Layer", (i-1)*layerThickness, i*layerThickness, 126 << 137 0., twopi, 0., pi); 127 fSpheres.push_back(sLayer); << 138 128 << 139 G4LogicalVolume* 129 G4LogicalVolume* lLayer = new G4LogicalVol << 140 lLayer = new G4LogicalVolume(sLayer, //shape 130 << 141 absorMaterial, //material 131 << 142 "Layer"); //name 132 fLVolumes.push_back(lLayer); << 143 133 << 144 new G4PVPlacement(0, //no rotation 134 new G4PVPlacement(0, // no rotation << 145 G4ThreeVector(), //at (0,0,0) 135 G4ThreeVector(), // at << 146 lLayer, //logical volume 136 lLayer, // logical volu << 147 "Layer", //name 137 "Layer", // name << 148 lAbsor, //mother volume 138 lAbsor, // mother volu << 149 false, //no boolean operation 139 false, // no boolean op << 150 i); //copy number 140 i); // copy number << 151 141 } << 152 } 142 153 143 PrintParameters(); 154 PrintParameters(); 144 << 155 145 // << 146 // always return the root volume << 147 // 156 // 148 return fAbsor; << 157 //always return the root volume >> 158 // >> 159 return pAbsor; 149 } 160 } 150 161 151 //....oooOO0OOooo........oooOO0OOooo........oo 162 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 152 163 153 void DetectorConstruction::PrintParameters() 164 void DetectorConstruction::PrintParameters() 154 { 165 { 155 G4cout << "\n------------------------------- 166 G4cout << "\n---------------------------------------------------------\n"; 156 G4cout << "---> The Absorber is a sphere of << 167 G4cout << "---> The Absorber is a sphere of " 157 << " radius of " << fAbsorMaterial->G << 168 << G4BestUnit(absorRadius,"Length") << " radius of " 158 << " slices of " << G4BestUnit(fLayer << 169 << absorMaterial->GetName() << " divided in " << nbOfLayers 159 << fAbsorMaterial << G4endl; << 170 << " slices of " << G4BestUnit(layerThickness,"Length") << G4endl; 160 G4cout << "\n------------------------------- 171 G4cout << "\n---------------------------------------------------------\n"; 161 } 172 } 162 173 163 //....oooOO0OOooo........oooOO0OOooo........oo 174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 164 175 165 void DetectorConstruction::SetRadius(G4double 176 void DetectorConstruction::SetRadius(G4double value) 166 { 177 { 167 // geometry was already constructed - scale << 178 absorRadius = value; 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; << 182 } 179 } 183 180 184 //....oooOO0OOooo........oooOO0OOooo........oo 181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 185 182 186 void DetectorConstruction::SetMaterial(G4Strin 183 void DetectorConstruction::SetMaterial(G4String materialChoice) 187 { 184 { 188 // search the material by its name << 185 // search the material by its name 189 G4Material* pttoMaterial = G4Material::GetMa << 186 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); >> 187 if (pttoMaterial) absorMaterial = pttoMaterial; >> 188 } 190 189 191 if (pttoMaterial && pttoMaterial != fAbsorMa << 190 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 192 fAbsorMaterial = pttoMaterial; << 191 193 G4RunManager::GetRunManager()->PhysicsHasB << 192 void DetectorConstruction::SetNbOfLayers(G4int value) 194 << 193 { 195 // geometry was already constructed - only << 194 nbOfLayers = value; 196 if (fAbsor) { << 195 } 197 for (auto lv : fLVolumes) { << 196 198 lv->SetMaterial(fAbsorMaterial); << 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 199 } << 198 >> 199 #include "G4FieldManager.hh" >> 200 #include "G4TransportationManager.hh" >> 201 >> 202 void DetectorConstruction::SetMagField(G4double fieldValue) >> 203 { >> 204 //apply a global uniform magnetic field along Z axis >> 205 G4FieldManager* fieldMgr >> 206 = G4TransportationManager::GetTransportationManager()->GetFieldManager(); >> 207 >> 208 if (magField) delete magField; //delete the existing magn field >> 209 >> 210 if (fieldValue!=0.) // create a new one if non nul >> 211 { >> 212 magField = new G4UniformMagField(G4ThreeVector(0.,0.,fieldValue)); >> 213 fieldMgr->SetDetectorField(magField); >> 214 fieldMgr->CreateChordFinder(magField); >> 215 } >> 216 else >> 217 { >> 218 magField = 0; >> 219 fieldMgr->SetDetectorField(magField); 200 } 220 } 201 } << 202 } 221 } 203 222 204 //....oooOO0OOooo........oooOO0OOooo........oo 223 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 205 224 206 void DetectorConstruction::SetNbOfLayers(G4int << 225 #include "G4RunManager.hh" >> 226 >> 227 void DetectorConstruction::UpdateGeometry() 207 { 228 { 208 fNbOfLayers = value; << 229 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes()); 209 } 230 } 210 231 211 //....oooOO0OOooo........oooOO0OOooo........oo 232 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 212 233