Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
examples/extended/electromagnetic/TestEm15/src/DetectorConstruction.cc
이 파일의 문서화 페이지로 가기
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 //
28 //
29 // $Id: DetectorConstruction.cc 68315 2013-03-21 18:26:24Z maire $
30 //
31 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 
34 #include "DetectorConstruction.hh"
35 #include "DetectorMessenger.hh"
36 
37 #include "G4Material.hh"
38 #include "G4Box.hh"
39 #include "G4LogicalVolume.hh"
40 #include "G4PVPlacement.hh"
41 
42 #include "G4GeometryManager.hh"
43 #include "G4PhysicalVolumeStore.hh"
44 #include "G4LogicalVolumeStore.hh"
45 #include "G4SolidStore.hh"
46 
47 #include "G4UnitsTable.hh"
48 #include "G4PhysicalConstants.hh"
49 #include "G4SystemOfUnits.hh"
50 
51 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 
54 :G4VUserDetectorConstruction(),fWorld(0), fBox(0), fMaterial(0),fWMaterial(0),
55  fDetectorMessenger(0)
56 {
57  fBoxSize = 100*m; fWorldSize = 1.2*fBoxSize;
59  SetMaterial("Water");
61 }
62 
63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
64 
66 { delete fDetectorMessenger;}
67 
68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69 
71 {
72  return ConstructVolumes();
73 }
74 
75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
76 
78 {
79  //
80  // define Elements
81  //
82  G4double z,a;
83 
84  G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole);
85  G4Element* N = new G4Element("Nitrogen" ,"N" , z= 7., a= 14.01*g/mole);
86  G4Element* O = new G4Element("Oxygen" ,"O" , z= 8., a= 16.00*g/mole);
87  G4Element* Na = new G4Element("Sodium" ,"Na", z=11., a= 22.99*g/mole);
88  G4Element* Ge = new G4Element("Germanium","Ge", z=32., a= 72.59*g/mole);
89  G4Element* I = new G4Element("Iodine" ,"I" , z=53., a= 126.90*g/mole);
90  G4Element* Bi = new G4Element("Bismuth" ,"Bi", z=83., a= 208.98*g/mole);
91 
92  //
93  // define materials
94  //
95  G4double density;
96  G4int ncomponents, natoms;
97  G4double fractionmass;
98 
99  G4Material* Air =
100  new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
101  Air->AddElement(N, fractionmass=70.*perCent);
102  Air->AddElement(O, fractionmass=30.*perCent);
103 
104  G4Material* H2l =
105  new G4Material("H2liquid", density= 70.8*mg/cm3, ncomponents=1);
106  H2l->AddElement(H, fractionmass=1.);
107 
108  G4Material* H2O =
109  new G4Material("Water", density= 1.000*g/cm3, ncomponents=2);
110  H2O->AddElement(H, natoms=2);
111  H2O->AddElement(O, natoms=1);
112  H2O->SetChemicalFormula("H_2O");
114 
115  new G4Material("liquidArgon", z=18., a= 39.95*g/mole, density= 1.390*g/cm3);
116 
117  new G4Material("Carbon" , z=6., a= 12.01*g/mole, density= 2.267*g/cm3);
118 
119  new G4Material("Aluminium" , z=13., a= 26.98*g/mole, density= 2.700*g/cm3);
120 
121  new G4Material("Silicon" , z=14., a= 28.09*g/mole, density= 2.330*g/cm3);
122 
123  new G4Material("Germanium" , z=32., a= 72.61*g/mole, density= 5.323*g/cm3);
124 
125  G4Material* NaI =
126  new G4Material("NaI", density= 3.67*g/cm3, ncomponents=2);
127  NaI->AddElement(Na, natoms=1);
128  NaI->AddElement(I , natoms=1);
130 
131  G4Material* Iod =
132  new G4Material("Iodine", density= 4.93*g/cm3, ncomponents=1);
133  Iod->AddElement(I , natoms=1);
134 
135  G4Material* BGO =
136  new G4Material("BGO", density= 7.10*g/cm3, ncomponents=3);
137  BGO->AddElement(O , natoms=12);
138  BGO->AddElement(Ge, natoms= 3);
139  BGO->AddElement(Bi, natoms= 4);
140 
141  new G4Material("Iron" , z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
142 
143  new G4Material("Tungsten" , z=74., a=183.85*g/mole, density= 19.30*g/cm3);
144 
145  new G4Material("Lead" , z=82., a=207.19*g/mole, density= 11.35*g/cm3);
146 
147  new G4Material("Uranium" , z=92., a=238.03*g/mole, density= 18.95*g/cm3);
148 
149  density = universe_mean_density; //from PhysicalConstants.h
150  G4double pressure = 3.e-18*pascal;
151  G4double temperature = 2.73*kelvin;
152  G4Material* vacuum =
153  new G4Material("Galactic",z= 1,a= 1.008*g/mole,density,
154  kStateGas,temperature,pressure);
155 
157 
158  //default material
159  fWMaterial = vacuum;
160 }
161 
162 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
163 
165 {
166  // Cleanup old geometry
171 
172  // World
173  //
174  G4Box*
175  sWorld = new G4Box("World", //name
176  fWorldSize/2,fWorldSize/2,fWorldSize/2); //dimensions
177 
179  lWorld = new G4LogicalVolume(sWorld, //shape
180  fWMaterial, //material
181  "World"); //name
182 
183  fWorld = new G4PVPlacement(0, //no rotation
184  G4ThreeVector(), //at (0,0,0)
185  lWorld, //logical volume
186  "World", //name
187  0, //mother volume
188  false, //no boolean operation
189  0); //copy number
190 
191  // Box
192  //
193  G4Box*
194  sBox = new G4Box("Container", //its name
195  fBoxSize/2,fBoxSize/2,fBoxSize/2); //its dimensions
196 
198  lBox = new G4LogicalVolume(sBox, //its shape
199  fMaterial, //its material
200  fMaterial->GetName()); //its name
201 
202  fBox = new G4PVPlacement(0, //no rotation
203  G4ThreeVector(), //at (0,0,0)
204  lBox, //its logical volume
205  fMaterial->GetName(), //its name
206  lWorld, //its mother volume
207  false, //no boolean operation
208  0); //copy number
209 
210  PrintParameters();
211 
212  //always return the root volume
213  //
214  return fWorld;
215 }
216 
217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
218 
220 {
221  G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
222  << " of " << fMaterial->GetName() << G4endl;
223 }
224 
225 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
226 
227 void DetectorConstruction::SetMaterial(G4String materialChoice)
228 {
229  // search the material by its name
230  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
231  if (pttoMaterial) {
232  fMaterial = pttoMaterial;
233  UpdateGeometry();
234  } else {
235  G4cout << "\n--> warning from DetectorConstruction::SetMaterial : "
236  << materialChoice << " not found" << G4endl;
237  }
238 }
239 
240 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
241 
243 {
245  UpdateGeometry();
246 }
247 
248 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
249 
250 #include "G4RunManager.hh"
251 
253 {
254  if (fWorld)
256 }
257 
258 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static constexpr double kelvin
Definition: G4SIunits.hh:281
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:80
static void Clean()
Definition: G4SolidStore.cc:73
CLHEP::Hep3Vector G4ThreeVector
std::vector< ExP01TrackerHit * > a
Definition: ExP01Classes.hh:33
static G4MaterialTable * GetMaterialTable()
Definition: G4Material.cc:593
#define G4endl
Definition: G4ios.hh:61
Double_t z
static constexpr double perCent
Definition: G4SIunits.hh:332
void SetMeanExcitationEnergy(G4double value)
virtual void DefineWorldVolume(G4VPhysicalVolume *worldVol, G4bool topologyIsChanged=true)
static constexpr double universe_mean_density
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
Definition: G4Material.cc:608
const G4String & GetName() const
Definition: G4Material.hh:179
static constexpr double g
Definition: G4SIunits.hh:183
static constexpr double m
Definition: G4SIunits.hh:129
double G4double
Definition: G4Types.hh:76
**D E S C R I P T I O N
const XML_Char int const XML_Char * value
Definition: expat.h:331
static constexpr double mg
Definition: G4SIunits.hh:184
static G4GeometryManager * GetInstance()
Definition: G4Box.hh:64
static constexpr double eV
Definition: G4SIunits.hh:215
void OpenGeometry(G4VPhysicalVolume *vol=0)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
static G4SolidStore * GetInstance()
static G4PhysicalVolumeStore * GetInstance()
void SetChemicalFormula(const G4String &chF)
Definition: G4Material.hh:174
int G4int
Definition: G4Types.hh:78
void AddElement(G4Element *element, G4int nAtoms)
Definition: G4Material.cc:368
G4IonisParamMat * GetIonisation() const
Definition: G4Material.hh:227
G4GLOB_DLL std::ostream G4cout
static constexpr double mole
Definition: G4SIunits.hh:286
#define pascal
static constexpr double cm3
Definition: G4SIunits.hh:121
static G4LogicalVolumeStore * GetInstance()