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 polarisation/Pol01/src/DetectorConst 26 /// \file polarisation/Pol01/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConst 27 /// \brief Implementation of the DetectorConstruction class 28 // 28 // 29 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 32 33 #include "DetectorConstruction.hh" 33 #include "DetectorConstruction.hh" 34 << 35 #include "DetectorMessenger.hh" 34 #include "DetectorMessenger.hh" 36 35 >> 36 #include "G4Material.hh" 37 #include "G4Box.hh" 37 #include "G4Box.hh" 38 #include "G4GeometryManager.hh" << 39 #include "G4LogicalVolume.hh" 38 #include "G4LogicalVolume.hh" 40 #include "G4LogicalVolumeStore.hh" << 41 #include "G4Material.hh" << 42 #include "G4NistManager.hh" << 43 #include "G4PVPlacement.hh" 39 #include "G4PVPlacement.hh" 44 #include "G4PhysicalVolumeStore.hh" << 40 45 #include "G4PolarizationManager.hh" << 46 #include "G4RunManager.hh" 41 #include "G4RunManager.hh" >> 42 #include "G4GeometryManager.hh" >> 43 #include "G4PhysicalVolumeStore.hh" >> 44 #include "G4LogicalVolumeStore.hh" 47 #include "G4SolidStore.hh" 45 #include "G4SolidStore.hh" 48 #include "G4SystemOfUnits.hh" << 46 49 #include "G4UnitsTable.hh" 47 #include "G4UnitsTable.hh" 50 48 >> 49 #include "G4PolarizationManager.hh" >> 50 #include "G4NistManager.hh" >> 51 #include "G4SystemOfUnits.hh" >> 52 51 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 52 54 53 DetectorConstruction::DetectorConstruction() 55 DetectorConstruction::DetectorConstruction() 54 : G4VUserDetectorConstruction(), fWorld(0), << 56 : G4VUserDetectorConstruction(), >> 57 fWorld(0), fBox(0), fTargetMaterial(0), fWorldMaterial(0) 55 { 58 { 56 fBoxSizeXY = 50 * mm; << 59 fBoxSizeXY = 50*mm; 57 fBoxSizeZ = 5 * mm; << 60 fBoxSizeZ = 5*mm; 58 fWorldSize = 1. * m; << 61 fWorldSize = 1.*m; 59 ConstructMaterials(); 62 ConstructMaterials(); 60 fMessenger = new DetectorMessenger(this); 63 fMessenger = new DetectorMessenger(this); 61 } 64 } 62 65 63 //....oooOO0OOooo........oooOO0OOooo........oo 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 64 67 65 DetectorConstruction::~DetectorConstruction() 68 DetectorConstruction::~DetectorConstruction() 66 { 69 { 67 delete fMessenger; 70 delete fMessenger; 68 } 71 } 69 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 70 73 71 void DetectorConstruction::ConstructMaterials( 74 void DetectorConstruction::ConstructMaterials() 72 { 75 { 73 auto nistManager = G4NistManager::Instance() 76 auto nistManager = G4NistManager::Instance(); 74 77 75 fTargetMaterial = nistManager->FindOrBuildMa 78 fTargetMaterial = nistManager->FindOrBuildMaterial("G4_Fe"); 76 if (fTargetMaterial == nullptr) { << 79 if(fTargetMaterial == nullptr) { 77 G4cerr << "### ERROR - Material: <" 80 G4cerr << "### ERROR - Material: <" 78 << "G4_Fe" << 81 << "G4_Fe" << "> not found" << G4endl; 79 << "> not found" << G4endl; << 80 } 82 } 81 83 82 fWorldMaterial = nistManager->FindOrBuildMat 84 fWorldMaterial = nistManager->FindOrBuildMaterial("G4_Galactic"); 83 if (fWorldMaterial == nullptr) { << 85 if(fWorldMaterial == nullptr) { 84 G4cerr << "### ERROR - Material: <" 86 G4cerr << "### ERROR - Material: <" 85 << "G4_Galactic" << 87 << "G4_Galactic" << "> not found" << G4endl; 86 << "> not found" << G4endl; << 87 } 88 } 88 } 89 } 89 90 90 //....oooOO0OOooo........oooOO0OOooo........oo 91 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 91 92 92 G4LogicalVolume* lBox; 93 G4LogicalVolume* lBox; 93 94 94 G4VPhysicalVolume* DetectorConstruction::Const 95 G4VPhysicalVolume* DetectorConstruction::Construct() 95 { 96 { 96 // G4PolarizationManager::GetInstance()->Cle 97 // G4PolarizationManager::GetInstance()->Clean(); 97 98 98 // World 99 // World 99 // 100 // 100 G4Box* sWorld = new G4Box("World", // name << 101 G4Box* 101 fWorldSize / 2, fW << 102 sWorld = new G4Box("World", //name 102 << 103 fWorldSize/2,fWorldSize/2,fWorldSize/2); //dimensions 103 G4LogicalVolume* lWorld = new G4LogicalVolum << 104 104 << 105 G4LogicalVolume* 105 << 106 lWorld = new G4LogicalVolume(sWorld, //shape 106 << 107 fWorldMaterial, //material 107 fWorld = new G4PVPlacement(0, // no rotatio << 108 "World"); //name 108 G4ThreeVector(), << 109 109 lWorld, // logic << 110 fWorld = new G4PVPlacement(0, //no rotation 110 "World", // name << 111 G4ThreeVector(), //at (0,0,0) 111 0, // mother vol << 112 lWorld, //logical volume 112 false, // no boo << 113 "World", //name 113 0); // copy numb << 114 0, //mother volume >> 115 false, //no boolean operation >> 116 0); //copy number 114 117 115 // Box 118 // Box 116 // 119 // 117 G4Box* sBox = new G4Box("Container", // its << 120 G4Box* 118 fBoxSizeXY / 2., fBo << 121 sBox = new G4Box("Container", //its name 119 << 122 fBoxSizeXY/2.,fBoxSizeXY/2.,fBoxSizeZ/2.);//its dimensions 120 // G4LogicalVolume* << 123 121 lBox = new G4LogicalVolume(sBox, // its sha << 124 //G4LogicalVolume* 122 fTargetMaterial, << 125 lBox = new G4LogicalVolume(sBox, //its shape 123 "theBox"); // it << 126 fTargetMaterial, //its material 124 << 127 "theBox"); //its name 125 fBox = new G4PVPlacement(0, // no rotation << 128 126 G4ThreeVector(), / << 129 fBox = new G4PVPlacement(0, //no rotation 127 lBox, // its logic << 130 G4ThreeVector(), //at (0,0,0) 128 fTargetMaterial->Ge << 131 lBox, //its logical volume 129 lWorld, // its mot << 132 fTargetMaterial->GetName(), //its name 130 false, // no boole << 133 lWorld, //its mother volume 131 0); // copy number << 134 false, //no boolean operation >> 135 0); //copy number 132 136 133 PrintParameters(); 137 PrintParameters(); 134 138 135 // always return the root volume << 139 //always return the root volume 136 // 140 // 137 return fWorld; 141 return fWorld; 138 } 142 } 139 143 140 //....oooOO0OOooo........oooOO0OOooo........oo 144 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 141 145 142 void DetectorConstruction::ConstructSDandField 146 void DetectorConstruction::ConstructSDandField() 143 { 147 { 144 // register logical Volume in PolarizationMa 148 // register logical Volume in PolarizationManager with zero polarization 145 G4PolarizationManager* polMgr = G4Polarizati << 149 G4PolarizationManager * polMgr = G4PolarizationManager::GetInstance(); 146 polMgr->SetVolumePolarization(lBox, G4ThreeV << 150 polMgr->SetVolumePolarization(lBox,G4ThreeVector(0.,0.,0.)); 147 } 151 } 148 152 149 //....oooOO0OOooo........oooOO0OOooo........oo 153 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 150 154 151 void DetectorConstruction::PrintParameters() 155 void DetectorConstruction::PrintParameters() 152 { 156 { 153 G4cout << "\n The Box is " << G4BestUnit(fBo << 157 G4cout << "\n The Box is " << G4BestUnit(fBoxSizeXY,"Length") 154 << G4BestUnit(fBoxSizeXY, "Length") < << 158 << " x " << G4BestUnit(fBoxSizeXY,"Length") 155 << fTargetMaterial->GetName() << G4en << 159 << " x " << G4BestUnit(fBoxSizeZ,"Length") >> 160 << " of " << fTargetMaterial->GetName() << G4endl; 156 } 161 } 157 162 158 //....oooOO0OOooo........oooOO0OOooo........oo 163 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 159 164 160 void DetectorConstruction::SetTargetMaterial(G 165 void DetectorConstruction::SetTargetMaterial(G4String materialChoice) 161 { 166 { 162 // search the material by its name 167 // search the material by its name 163 G4Material* mat = G4NistManager::Instance()- << 168 G4Material* mat = >> 169 G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 164 if (mat != fTargetMaterial) { 170 if (mat != fTargetMaterial) { 165 if (mat) { << 171 if(mat) { 166 fTargetMaterial = mat; 172 fTargetMaterial = mat; 167 G4RunManager::GetRunManager()->PhysicsHa 173 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 168 UpdateGeometry(); 174 UpdateGeometry(); 169 } << 175 } else { 170 else { << 176 G4cout << "### Warning! Target material: <" 171 G4cout << "### Warning! Target material << 177 << materialChoice << "> not found" << G4endl; 172 } 178 } 173 } 179 } 174 } 180 } 175 181 176 //....oooOO0OOooo........oooOO0OOooo........oo 182 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 177 183 178 void DetectorConstruction::SetWorldMaterial(G4 184 void DetectorConstruction::SetWorldMaterial(G4String materialChoice) 179 { 185 { 180 // search the material by its name 186 // search the material by its name 181 G4Material* mat = G4NistManager::Instance()- << 187 G4Material* mat = >> 188 G4NistManager::Instance()->FindOrBuildMaterial(materialChoice); 182 if (mat != fWorldMaterial) { 189 if (mat != fWorldMaterial) { 183 if (mat) { << 190 if(mat) { 184 fWorldMaterial = mat; 191 fWorldMaterial = mat; 185 G4RunManager::GetRunManager()->PhysicsHa 192 G4RunManager::GetRunManager()->PhysicsHasBeenModified(); 186 UpdateGeometry(); 193 UpdateGeometry(); 187 } << 194 } else { 188 else { << 195 G4cout << "### Warning! World material: <" 189 G4cout << "### Warning! World material: << 196 << materialChoice << "> not found" << G4endl; 190 } 197 } 191 } 198 } 192 } 199 } 193 200 194 //....oooOO0OOooo........oooOO0OOooo........oo 201 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 195 202 196 void DetectorConstruction::SetSizeXY(G4double 203 void DetectorConstruction::SetSizeXY(G4double value) 197 { 204 { 198 fBoxSizeXY = value; 205 fBoxSizeXY = value; 199 if (fWorldSize < fBoxSizeXY) fWorldSize = 1. << 206 if (fWorldSize<fBoxSizeXY) fWorldSize = 1.2*fBoxSizeXY; 200 UpdateGeometry(); 207 UpdateGeometry(); 201 } 208 } 202 209 203 void DetectorConstruction::SetSizeZ(G4double v 210 void DetectorConstruction::SetSizeZ(G4double value) 204 { 211 { 205 fBoxSizeZ = value; 212 fBoxSizeZ = value; 206 if (fWorldSize < fBoxSizeZ) fWorldSize = 1.2 << 213 if (fWorldSize<fBoxSizeZ) fWorldSize = 1.2*fBoxSizeZ; 207 UpdateGeometry(); 214 UpdateGeometry(); 208 } 215 } 209 216 210 //....oooOO0OOooo........oooOO0OOooo........oo 217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 211 218 212 #include "G4RunManager.hh" 219 #include "G4RunManager.hh" 213 220 214 void DetectorConstruction::UpdateGeometry() 221 void DetectorConstruction::UpdateGeometry() 215 { 222 { 216 // if (fWorld) 223 // if (fWorld) 217 // G4RunManager::GetRunManager()->DefineWo 224 // G4RunManager::GetRunManager()->DefineWorldVolume(Construct()); 218 G4RunManager::GetRunManager()->GeometryHasBe 225 G4RunManager::GetRunManager()->GeometryHasBeenModified(); 219 } 226 } 220 227 221 //....oooOO0OOooo........oooOO0OOooo........oo 228 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 222 229