Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
examples/extended/hadronic/Hadr04/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 70755 2013-06-05 12:17:48Z ihrivnac $
30 //
31 
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 #include "DetectorConstruction.hh"
36 #include "DetectorMessenger.hh"
37 #include "G4Material.hh"
38 #include "G4NistManager.hh"
39 
40 #include "G4Box.hh"
41 #include "G4LogicalVolume.hh"
42 #include "G4PVPlacement.hh"
43 
44 #include "G4GeometryManager.hh"
45 #include "G4PhysicalVolumeStore.hh"
46 #include "G4LogicalVolumeStore.hh"
47 #include "G4SolidStore.hh"
48 #include "G4RunManager.hh"
49 
50 #include "G4UnitsTable.hh"
51 #include "G4SystemOfUnits.hh"
52 
53 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
54 
57  fPBox(0), fLBox(0), fMaterial(0), fDetectorMessenger(0)
58 {
59  fBoxSize = 1*m;
61  SetMaterial("Water_ts");
63 }
64 
65 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
66 
68 { delete fDetectorMessenger;}
69 
70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
71 
73 {
74  return ConstructVolumes();
75 }
76 
77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
78 
80 {
81  // specific element name for thermal neutronHP
82  // (see G4ParticleHPThermalScatteringNames.cc)
83 
84  G4int ncomponents, natoms;
85 
86  // pressurized water
87  G4Element* H = new G4Element("TS_H_of_Water" ,"H" , 1., 1.0079*g/mole);
88  G4Element* O = new G4Element("Oxygen" ,"O" , 8., 16.00*g/mole);
89  G4Material* H2O =
90  new G4Material("Water_ts", 1.000*g/cm3, ncomponents=2,
91  kStateLiquid, 593*kelvin, 150*bar);
92  H2O->AddElement(H, natoms=2);
93  H2O->AddElement(O, natoms=1);
95 
96  // heavy water
97  G4Isotope* H2 = new G4Isotope("H2",1,2);
98  G4Element* D = new G4Element("TS_D_of_Heavy_Water", "D", 1);
99  D->AddIsotope(H2, 100*perCent);
100  G4Material* D2O = new G4Material("HeavyWater", 1.11*g/cm3, ncomponents=2,
101  kStateLiquid, 293.15*kelvin, 1*atmosphere);
102  D2O->AddElement(D, natoms=2);
103  D2O->AddElement(O, natoms=1);
104 
105  // graphite
106  G4Isotope* C12 = new G4Isotope("C12", 6, 12);
107  G4Element* C = new G4Element("TS_C_of_Graphite","C", ncomponents=1);
108  C->AddIsotope(C12, 100.*perCent);
109  G4Material* graphite =
110  new G4Material("graphite", 2.27*g/cm3, ncomponents=1,
111  kStateSolid, 293*kelvin, 1*atmosphere);
112  graphite->AddElement(C, natoms=1);
113 
115 }
116 
117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
118 
120  G4String symbol, G4double density, G4int Z, G4int A)
121 {
122  // define a material from an isotope
123  //
124  G4int ncomponents;
125  G4double abundance, massfraction;
126 
127  G4Isotope* isotope = new G4Isotope(symbol, Z, A);
128 
129  G4Element* element = new G4Element(name, symbol, ncomponents=1);
130  element->AddIsotope(isotope, abundance= 100.*perCent);
131 
132  G4Material* material = new G4Material(name, density, ncomponents=1);
133  material->AddElement(element, massfraction=100.*perCent);
134 
135  return material;
136 }
137 
138 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
139 
141 {
142  // Cleanup old geometry
147 
148  G4Box*
149  sBox = new G4Box("Container", //its name
150  fBoxSize/2,fBoxSize/2,fBoxSize/2); //its dimensions
151 
152  fLBox = new G4LogicalVolume(sBox, //its shape
153  fMaterial, //its material
154  fMaterial->GetName()); //its name
155 
156  fPBox = new G4PVPlacement(0, //no rotation
157  G4ThreeVector(), //at (0,0,0)
158  fLBox, //its logical volume
159  fMaterial->GetName(), //its name
160  0, //its mother volume
161  false, //no boolean operation
162  0); //copy number
163 
164  PrintParameters();
165 
166  //always return the root volume
167  //
168  return fPBox;
169 }
170 
171 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
172 
174 {
175  G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
176  << " of " << fMaterial->GetName()
177  << "\n \n" << fMaterial << G4endl;
178 }
179 
180 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
181 
182 void DetectorConstruction::SetMaterial(G4String materialChoice)
183 {
184  // search the material by its name
185  G4Material* pttoMaterial =
186  G4NistManager::Instance()->FindOrBuildMaterial(materialChoice);
187 
188  if (pttoMaterial) {
189  if(fMaterial != pttoMaterial) {
190  fMaterial = pttoMaterial;
191  if(fLBox) { fLBox->SetMaterial(pttoMaterial); }
193  }
194  } else {
195  G4cout << "\n--> warning from DetectorConstruction::SetMaterial : "
196  << materialChoice << " not found" << G4endl;
197  }
198 }
199 
200 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
201 
203 {
204  fBoxSize = value;
206 }
207 
208 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
209 
static constexpr double kelvin
Definition: G4SIunits.hh:281
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:80
const XML_Char * name
Definition: expat.h:151
static constexpr double atmosphere
Definition: G4SIunits.hh:237
static void Clean()
Definition: G4SolidStore.cc:73
CLHEP::Hep3Vector G4ThreeVector
#define G4endl
Definition: G4ios.hh:61
double D(double temp)
static constexpr double perCent
Definition: G4SIunits.hh:332
void SetMeanExcitationEnergy(G4double value)
const G4String & GetName() const
Definition: G4Material.hh:179
void AddIsotope(G4Isotope *isotope, G4double RelativeAbundance)
Definition: G4Element.cc:152
Float_t Z
static constexpr double g
Definition: G4SIunits.hh:183
static constexpr double m
Definition: G4SIunits.hh:129
double G4double
Definition: G4Types.hh:76
const XML_Char int const XML_Char * value
Definition: expat.h:331
static G4GeometryManager * GetInstance()
Definition: G4Box.hh:64
G4Material * FindOrBuildMaterial(const G4String &name, G4bool isotopes=true, G4bool warning=false)
double A(double temperature)
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()
int G4int
Definition: G4Types.hh:78
void ReinitializeGeometry(G4bool destroyFirst=false, G4bool prop=true)
double C(double temp)
void AddElement(G4Element *element, G4int nAtoms)
Definition: G4Material.cc:368
G4IonisParamMat * GetIonisation() const
Definition: G4Material.hh:227
void PhysicsHasBeenModified()
G4GLOB_DLL std::ostream G4cout
void SetMaterial(G4Material *pMaterial)
static constexpr double mole
Definition: G4SIunits.hh:286
static constexpr double bar
Definition: G4SIunits.hh:236
static constexpr double cm3
Definition: G4SIunits.hh:121
static G4LogicalVolumeStore * GetInstance()
static G4NistManager * Instance()
G4Material * MaterialWithSingleIsotope(G4String, G4String, G4double, G4int, G4int)