Geant4  v4-10.4-release
 모두 클래스 네임스페이스들 파일들 함수 변수 타입정의 열거형 타입 열거형 멤버 Friends 매크로 그룹들 페이지들
examples/extended/electromagnetic/TestEm16/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 84365 2014-10-14 12:43:52Z gcosmo $
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 #include "G4UserLimits.hh"
43 #include "G4PropagatorInField.hh"
44 
45 #include "G4GeometryManager.hh"
46 #include "G4PhysicalVolumeStore.hh"
47 #include "G4LogicalVolumeStore.hh"
48 #include "G4SolidStore.hh"
49 #include "G4RunManager.hh"
50 
51 #include "G4UnitsTable.hh"
52 #include "G4SystemOfUnits.hh"
53 
54 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
55 
58  fLBox(0), fBox(0), fBoxSize(500*m), fMaterial(0),
59  fUserLimits(0), fDetectorMessenger(0)
60 {
62  SetMaterial("Iron");
63 
64  // create UserLimits
65  fUserLimits = new G4UserLimits();
66 
67  // create commands for interactive definition of the detector
69 }
70 
71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72 
74 { delete fDetectorMessenger;}
75 
76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
77 
79 {
80  return ConstructVolumes();
81 }
82 
83 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 
86 {
87  G4double a, z, density;
88 
89  new G4Material("Beryllium", z= 4., a= 9.012182*g/mole, density= 1.848*g/cm3);
90  new G4Material("Carbon", z= 6., a= 12.011*g/mole, density= 2.265*g/cm3);
91  new G4Material("Iron", z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
92 
93  // define a vacuum with a restgas pressure typical for accelerators
94  G4double const Torr = atmosphere/760.; // 1 Torr
95  G4double pressure = 10e-9*Torr, temperature = 296.150*kelvin; // 23 Celsius
96  new G4Material("Vacuum", z=7., a=14.01*g/mole, density= 1.516784e-11*kg/m3,
97  kStateGas, temperature, pressure);
98 
100 }
101 
102 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
103 
105 {
110 
111  G4Box*
112  sBox = new G4Box("Container", //its name
113  fBoxSize/2,fBoxSize/2,fBoxSize/2); //its dimensions
114 
115  fLBox = new G4LogicalVolume(sBox, //its shape
116  fMaterial, //its material
117  fMaterial->GetName()); //its name
118 
120 
121  fBox = new G4PVPlacement(0, //no rotation
122  G4ThreeVector(), //at (0,0,0)
123  fLBox, //its logical volume
124  fMaterial->GetName(), //its name
125  0, //its mother volume
126  false, //no boolean operation
127  0); //copy number
128 
129  PrintParameters();
130 
131  //
132  //always return the root volume
133  //
134  return fBox;
135 }
136 
137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
138 
140 {
141  G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
142  << " of " << fMaterial->GetName() << G4endl;
143 }
144 
145 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
146 
147 void DetectorConstruction::SetMaterial(G4String materialChoice)
148 {
149  // search the material by its name
150  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
151  if (pttoMaterial) {
152  fMaterial = pttoMaterial;
153  if ( fLBox ) fLBox->SetMaterial(fMaterial);
154  }
156 }
157 
158 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
159 
161 {
162  fBoxSize = value;
164 }
165 
166 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
167 
169 #include "G4AutoDelete.hh"
170 
172 {
173  if ( fFieldMessenger.Get() == 0 ) {
174  // Create global magnetic field messenger.
175  // Uniform magnetic field is then created automatically if
176  // the field value is not zero.
177  G4ThreeVector fieldValue = G4ThreeVector();
179  new G4GlobalMagFieldMessenger(fieldValue);
180  //msg->SetVerboseLevel(1);
182  fFieldMessenger.Put( msg );
183 
184  }
185 }
186 
187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
188 
190 {
191  // set the maximum allowed step size
192  //
193  if (val <= DBL_MIN)
194  { G4cout << "\n --->warning from SetMaxStepSize: maxStep "
195  << val << " out of range. Command refused" << G4endl;
196  return;
197  }
199 }
200 
201 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
202 
204 {
205  // set the maximum length of tracking step
206  //
207  if (val <= DBL_MIN)
208  { G4cout << "\n --->warning from SetMaxStepLength: maxStep "
209  << val << " out of range. Command refused" << G4endl;
210  return;
211  }
212  G4TransportationManager* tmanager =
214  tmanager->GetPropagatorInField()
216 }
217 
218 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
static constexpr double kelvin
Definition: G4SIunits.hh:281
static G4RunManager * GetRunManager()
Definition: G4RunManager.cc:80
static constexpr double atmosphere
Definition: G4SIunits.hh:237
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
void SetLargestAcceptableStep(G4double newBigDist)
static constexpr double kg
Definition: G4SIunits.hh:182
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
Definition: G4Material.cc:608
void Register(T *inst)
Definition: G4AutoDelete.hh:65
const G4String & GetName() const
Definition: G4Material.hh:179
static constexpr double g
Definition: G4SIunits.hh:183
static constexpr double m
Definition: G4SIunits.hh:129
static constexpr double m3
Definition: G4SIunits.hh:131
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
value_type & Get() const
Definition: G4Cache.hh:314
void OpenGeometry(G4VPhysicalVolume *vol=0)
void Put(const value_type &val) const
Definition: G4Cache.hh:318
void SetUserLimits(G4UserLimits *pULimits)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
static G4TransportationManager * GetTransportationManager()
static G4SolidStore * GetInstance()
static G4PhysicalVolumeStore * GetInstance()
void ReinitializeGeometry(G4bool destroyFirst=false, G4bool prop=true)
#define DBL_MIN
Definition: templates.hh:75
void PhysicsHasBeenModified()
G4Cache< G4GlobalMagFieldMessenger * > fFieldMessenger
G4GLOB_DLL std::ostream G4cout
void SetMaterial(G4Material *pMaterial)
G4PropagatorInField * GetPropagatorInField() const
static constexpr double mole
Definition: G4SIunits.hh:286
virtual void SetMaxAllowedStep(G4double ustepMax)
static constexpr double cm3
Definition: G4SIunits.hh:121
static G4LogicalVolumeStore * GetInstance()