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 medical/fanoCavity/src/DetectorConst 26 /// \file medical/fanoCavity/src/DetectorConstruction.cc 27 /// \brief Implementation of the DetectorConst 27 /// \brief Implementation of the DetectorConstruction class 28 // 28 // 29 29 30 // 30 // 31 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 32 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 33 33 34 #include "DetectorConstruction.hh" 34 #include "DetectorConstruction.hh" 35 << 36 #include "DetectorMessenger.hh" 35 #include "DetectorMessenger.hh" 37 36 38 #include "G4GeometryManager.hh" << 39 #include "G4LogicalVolume.hh" << 40 #include "G4LogicalVolumeStore.hh" << 41 #include "G4Material.hh" 37 #include "G4Material.hh" 42 #include "G4NistManager.hh" 38 #include "G4NistManager.hh" >> 39 #include "G4Tubs.hh" >> 40 #include "G4LogicalVolume.hh" >> 41 #include "G4VPhysicalVolume.hh" 43 #include "G4PVPlacement.hh" 42 #include "G4PVPlacement.hh" 44 #include "G4PhysicalConstants.hh" << 43 >> 44 #include "G4GeometryManager.hh" 45 #include "G4PhysicalVolumeStore.hh" 45 #include "G4PhysicalVolumeStore.hh" >> 46 #include "G4LogicalVolumeStore.hh" 46 #include "G4SolidStore.hh" 47 #include "G4SolidStore.hh" 47 #include "G4SystemOfUnits.hh" << 48 48 #include "G4Tubs.hh" << 49 #include "G4UnitsTable.hh" 49 #include "G4UnitsTable.hh" 50 #include "G4VPhysicalVolume.hh" << 50 #include "G4PhysicalConstants.hh" >> 51 #include "G4SystemOfUnits.hh" 51 52 52 //....oooOO0OOooo........oooOO0OOooo........oo 53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 53 54 54 DetectorConstruction::DetectorConstruction() : << 55 DetectorConstruction::DetectorConstruction() >> 56 : fWall(nullptr), fCavity(nullptr) 55 { 57 { 56 // default parameter values 58 // default parameter values 57 fCavityThickness = 2 * mm; << 59 fCavityThickness = 2*mm; 58 fCavityRadius = 1 * cm; << 60 fCavityRadius = 1*cm; 59 << 61 60 fWallThickness = 5 * mm; << 62 fWallThickness = 5*mm; 61 << 63 62 DefineMaterials(); 64 DefineMaterials(); 63 SetWallMaterial("G4_WATER"); 65 SetWallMaterial("G4_WATER"); 64 SetCavityMaterial("g4Water_gas"); 66 SetCavityMaterial("g4Water_gas"); 65 << 67 66 // create commands for interactive definitio << 68 // create commands for interactive definition of the detector 67 fDetectorMessenger = new DetectorMessenger(t 69 fDetectorMessenger = new DetectorMessenger(this); 68 } 70 } 69 71 70 //....oooOO0OOooo........oooOO0OOooo........oo 72 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 71 73 72 DetectorConstruction::~DetectorConstruction() 74 DetectorConstruction::~DetectorConstruction() 73 { << 75 { 74 delete fDetectorMessenger; 76 delete fDetectorMessenger; 75 } 77 } 76 78 77 //....oooOO0OOooo........oooOO0OOooo........oo 79 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 78 80 79 void DetectorConstruction::DefineMaterials() 81 void DetectorConstruction::DefineMaterials() 80 { 82 { 81 G4double z, a; << 83 G4double z,a; 82 << 84 83 G4Element* H = new G4Element("Hydrogen", "H" << 85 G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole); 84 G4Element* O = new G4Element("Oxygen", "O", << 86 G4Element* O = new G4Element("Oxygen" ,"O" , z= 8., a= 16.00*g/mole); 85 87 86 G4Material* H2O = new G4Material("Water", 1. << 88 G4Material* H2O = >> 89 new G4Material("Water", 1.0*g/cm3, 2); 87 H2O->AddElement(H, 2); 90 H2O->AddElement(H, 2); 88 H2O->AddElement(O, 1); 91 H2O->AddElement(O, 1); 89 H2O->GetIonisation()->SetMeanExcitationEnerg << 92 H2O->GetIonisation()->SetMeanExcitationEnergy(78.0*eV); 90 << 93 91 G4Material* gas = new G4Material("Water_gas" << 94 G4Material* gas = >> 95 new G4Material("Water_gas", 1.0*mg/cm3, 2); 92 gas->AddElement(H, 2); 96 gas->AddElement(H, 2); 93 gas->AddElement(O, 1); 97 gas->AddElement(O, 1); 94 gas->GetIonisation()->SetMeanExcitationEnerg << 98 gas->GetIonisation()->SetMeanExcitationEnergy(78.0*eV); 95 << 99 96 new G4Material("Graphite", 6, 12.01 * g / mo << 100 new G4Material("Graphite", 6, 12.01*g/mole, 2.265*g/cm3); 97 new G4Material("Graphite_gas", 6, 12.01 * g << 101 new G4Material("Graphite_gas", 6, 12.01*g/mole, 2.265*mg/cm3); 98 << 102 99 new G4Material("Aluminium", 13, 26.98 * g / << 103 new G4Material("Aluminium", 13, 26.98*g/mole, 2.700*g/cm3); 100 new G4Material("Aluminium_gas", 13, 26.98 * << 104 new G4Material("Aluminium_gas", 13, 26.98*g/mole, 2.700*mg/cm3); 101 << 105 102 // alternatively, use G4 data base 106 // alternatively, use G4 data base 103 // 107 // 104 G4NistManager* nist = G4NistManager::Instanc 108 G4NistManager* nist = G4NistManager::Instance(); 105 109 106 nist->FindOrBuildMaterial("G4_WATER"); 110 nist->FindOrBuildMaterial("G4_WATER"); 107 nist->BuildMaterialWithNewDensity("g4Water_g << 111 nist->BuildMaterialWithNewDensity("g4Water_gas","G4_WATER", 1.0*mg/cm3); 108 112 109 // printout 113 // printout 110 G4cout << *(G4Material::GetMaterialTable()) 114 G4cout << *(G4Material::GetMaterialTable()) << G4endl; 111 } 115 } 112 116 113 //....oooOO0OOooo........oooOO0OOooo........oo 117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 114 << 118 115 G4VPhysicalVolume* DetectorConstruction::Const 119 G4VPhysicalVolume* DetectorConstruction::Construct() 116 { 120 { 117 if (fWall) { << 121 if(fWall) { return fWall; } 118 return fWall; << 119 } << 120 122 121 // Chamber 123 // Chamber 122 // 124 // 123 fTotalThickness = fCavityThickness + 2 * fWa << 125 fTotalThickness = fCavityThickness + 2*fWallThickness; 124 fWallRadius = fCavityRadius + fWallThickness << 126 fWallRadius = fCavityRadius + fWallThickness; 125 << 127 126 G4Tubs* sChamber = new G4Tubs("Chamber", // << 128 G4Tubs* 127 0., fWallRadiu << 129 sChamber = new G4Tubs("Chamber", //name 128 << 130 0.,fWallRadius,0.5*fTotalThickness,0.,twopi); //size 129 G4LogicalVolume* lChamber = new G4LogicalVol << 131 130 << 132 G4LogicalVolume* 131 << 133 lChamber = new G4LogicalVolume(sChamber, //solid 132 << 134 fWallMaterial, //material 133 fWall = new G4PVPlacement(0, // no rotation << 135 "Chamber"); //name 134 G4ThreeVector(), << 136 135 lChamber, // logi << 137 fWall = new G4PVPlacement(0, //no rotation 136 "Wall", // name << 138 G4ThreeVector(), //at (0,0,0) 137 0, // mother vol << 139 lChamber, //logical volume 138 false, // no bool << 140 "Wall", //name 139 0); // copy numbe << 141 0, //mother volume >> 142 false, //no boolean operation >> 143 0); //copy number 140 144 141 // Cavity 145 // Cavity 142 // << 146 // 143 G4Tubs* sCavity = new G4Tubs("Cavity", 0., f << 147 G4Tubs* 144 << 148 sCavity = new G4Tubs("Cavity", 145 G4LogicalVolume* lCavity = new G4LogicalVolu << 149 0.,fCavityRadius,0.5*fCavityThickness,0.,twopi); 146 << 150 147 << 151 G4LogicalVolume* 148 << 152 lCavity = new G4LogicalVolume(sCavity, //shape 149 fCavity = new G4PVPlacement(0, // no rotati << 153 fCavityMaterial, //material 150 G4ThreeVector(), << 154 "Cavity"); //name 151 lCavity, // log << 155 152 "Cavity", // na << 156 fCavity = new G4PVPlacement(0, //no rotation 153 lChamber, // mo << 157 G4ThreeVector(), //at (0,0,0) 154 false, // no bo << 158 lCavity, //logical volume 155 1); // copy num << 159 "Cavity", //name 156 << 160 lChamber, //mother volume >> 161 false, //no boolean operation >> 162 1); //copy number >> 163 157 PrintParameters(); 164 PrintParameters(); 158 << 165 159 // << 160 // always return the root volume << 161 // 166 // >> 167 //always return the root volume >> 168 // 162 return fWall; 169 return fWall; 163 } 170 } 164 171 165 //....oooOO0OOooo........oooOO0OOooo........oo 172 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 166 173 167 void DetectorConstruction::PrintParameters() 174 void DetectorConstruction::PrintParameters() 168 { 175 { 169 G4cout << "\n------------------------------- 176 G4cout << "\n---------------------------------------------------------\n"; 170 G4cout << "---> The Wall is " << G4BestUnit( << 177 G4cout << "---> The Wall is " << G4BestUnit(fWallThickness,"Length") 171 << fWallMaterial->GetName() << " ( " << 178 << " of " << fWallMaterial->GetName() << " ( " 172 << G4BestUnit(fWallMaterial->GetDensi << 179 << G4BestUnit(fWallMaterial->GetDensity(),"Volumic Mass") << " )\n"; 173 G4cout << " The Cavity is " << G4BestUni << 180 G4cout << " The Cavity is " << G4BestUnit(fCavityThickness,"Length") 174 << fCavityMaterial->GetName() << " ( << 181 << " of " << fCavityMaterial->GetName() << " ( " 175 << G4BestUnit(fCavityMaterial->GetDen << 182 << G4BestUnit(fCavityMaterial->GetDensity(),"Volumic Mass") << " )"; 176 G4cout << "\n------------------------------- 183 G4cout << "\n---------------------------------------------------------\n"; 177 G4cout << G4endl; 184 G4cout << G4endl; 178 } 185 } 179 186 180 //....oooOO0OOooo........oooOO0OOooo........oo 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 181 188 182 void DetectorConstruction::SetWallThickness(G4 189 void DetectorConstruction::SetWallThickness(G4double value) 183 { 190 { 184 fWallThickness = value; 191 fWallThickness = value; 185 } 192 } 186 << 193 187 //....oooOO0OOooo........oooOO0OOooo........oo 194 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 188 195 189 void DetectorConstruction::SetWallMaterial(con 196 void DetectorConstruction::SetWallMaterial(const G4String& materialChoice) 190 { 197 { 191 // search the material by its name << 198 // search the material by its name 192 G4Material* pttoMaterial = G4Material::GetMa << 199 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); 193 if (pttoMaterial) fWallMaterial = pttoMateri 200 if (pttoMaterial) fWallMaterial = pttoMaterial; 194 } 201 } 195 202 196 //....oooOO0OOooo........oooOO0OOooo........oo 203 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 197 204 198 void DetectorConstruction::SetCavityThickness( 205 void DetectorConstruction::SetCavityThickness(G4double value) 199 { 206 { 200 fCavityThickness = value; 207 fCavityThickness = value; 201 } 208 } 202 << 209 203 //....oooOO0OOooo........oooOO0OOooo........oo 210 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 204 211 205 void DetectorConstruction::SetCavityRadius(G4d 212 void DetectorConstruction::SetCavityRadius(G4double value) 206 { 213 { 207 fCavityRadius = value; << 214 fCavityRadius = value; 208 } << 215 } 209 216 210 //....oooOO0OOooo........oooOO0OOooo........oo 217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 211 218 212 void DetectorConstruction::SetCavityMaterial(c 219 void DetectorConstruction::SetCavityMaterial(const G4String& materialChoice) 213 { 220 { 214 // search the material by its name 221 // search the material by its name 215 G4Material* pttoMaterial = G4Material::GetMa << 222 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice); 216 if (pttoMaterial) fCavityMaterial = pttoMate 223 if (pttoMaterial) fCavityMaterial = pttoMaterial; 217 } 224 } 218 225 219 //....oooOO0OOooo........oooOO0OOooo........oo 226 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 220 227