source: trunk/examples/extended/electromagnetic/TestEm17/src/DetectorConstruction.cc@ 1230

Last change on this file since 1230 was 1230, checked in by garnier, 16 years ago

update to geant4.9.3
File size: 6.6 KB
Line 
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//
26//
27// $Id: DetectorConstruction.cc,v 1.3 2007/05/12 16:33:53 vnivanch Exp $
28// GEANT4 tag $Name: geant4-09-03-cand-01 $
29//
30//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
31//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
32
33#include "DetectorConstruction.hh"
34#include "DetectorMessenger.hh"
35
36#include "G4Material.hh"
37#include "G4Box.hh"
38#include "G4LogicalVolume.hh"
39#include "G4PVPlacement.hh"
40
41#include "G4GeometryManager.hh"
42#include "G4PhysicalVolumeStore.hh"
43#include "G4LogicalVolumeStore.hh"
44#include "G4SolidStore.hh"
45
46#include "G4UnitsTable.hh"
47
48//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49
50DetectorConstruction::DetectorConstruction()
51:pBox(0), lBox(0), aMaterial(0)
52{
53 BoxSize = 1*m;
54 DefineMaterials();
55 SetMaterial("Iron");
56 detectorMessenger = new DetectorMessenger(this);
57}
58
59//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60
61DetectorConstruction::~DetectorConstruction()
62{ delete detectorMessenger;}
63
64//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
65
66G4VPhysicalVolume* DetectorConstruction::Construct()
67{
68 return ConstructVolumes();
69}
70
71//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
72
73void DetectorConstruction::DefineMaterials()
74{
75 //
76 // define Elements
77 //
78 G4double z,a;
79
80 G4Element* C = new G4Element("Carbon" ,"C" , z= 6., a= 12.01*g/mole);
81 G4Element* N = new G4Element("Nitrogen" ,"N" , z= 7., a= 14.01*g/mole);
82 G4Element* O = new G4Element("Oxygen" ,"O" , z= 8., a= 16.00*g/mole);
83 G4Element* Ca = new G4Element("Calcium" ,"Ca", z=20., a= 40.08*g/mole);
84 G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole);
85 G4Element* I = new G4Element("Iodine" ,"I" , z= 53., a= 126.90*g/mole);
86
87 //
88 // define materials
89 //
90 G4double density;
91 G4int ncomponents, natoms;
92 G4double fractionmass;
93
94 new G4Material("galactic", z=1., a= 1.01*g/mole, universe_mean_density,
95 kStateGas, 2.73*kelvin, 3.e-18*pascal);
96
97 G4Material* Air =
98 new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
99 Air->AddElement(N, fractionmass=70.*perCent);
100 Air->AddElement(O, fractionmass=30.*perCent);
101
102 G4Material* CaCO3 =
103 new G4Material("CaCO3", density= 2.80*g/cm3, ncomponents=3);
104 CaCO3->AddElement(Ca, natoms= 1);
105 CaCO3->AddElement(C , natoms= 1);
106 CaCO3->AddElement(O , natoms= 3);
107
108 new G4Material("Carbon", z= 6., a= 12.01*g/mole, density= 2.265*g/cm3);
109 new G4Material("Iron" , z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
110 new G4Material("Tin" , z=50., a= 118.7*g/mole, density= 7.310*g/cm3);
111
112 G4Material* HI =
113 new G4Material("HI", density= 2.8*g/cm3, ncomponents=2); //-35oC
114 HI->AddElement(H, natoms= 1);
115 HI->AddElement(I, natoms= 1);
116
117 G4cout << *(G4Material::GetMaterialTable()) << G4endl;
118}
119
120//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
121
122G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
123{
124 // Cleanup old geometry
125 G4GeometryManager::GetInstance()->OpenGeometry();
126 G4PhysicalVolumeStore::GetInstance()->Clean();
127 G4LogicalVolumeStore::GetInstance()->Clean();
128 G4SolidStore::GetInstance()->Clean();
129
130 G4Box*
131 sBox = new G4Box("Container", //its name
132 BoxSize/2,BoxSize/2,BoxSize/2); //its dimensions
133
134 lBox = new G4LogicalVolume(sBox, //its shape
135 aMaterial, //its material
136 aMaterial->GetName()); //its name
137
138 pBox = new G4PVPlacement(0, //no rotation
139 G4ThreeVector(), //at (0,0,0)
140 lBox, //its logical volume
141 aMaterial->GetName(), //its name
142 0, //its mother volume
143 false, //no boolean operation
144 0); //copy number
145
146 PrintParameters();
147
148 //always return the root volume
149 //
150 return pBox;
151}
152
153//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
154
155void DetectorConstruction::PrintParameters()
156{
157 G4cout << "\n The Box is " << G4BestUnit(BoxSize,"Length")
158 << " of " << aMaterial->GetName() << G4endl;
159}
160
161//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
162
163void DetectorConstruction::SetMaterial(G4String materialChoice)
164{
165 // search the material by its name
166 G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
167 if (pttoMaterial) aMaterial = pttoMaterial;
168}
169
170//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
171
172void DetectorConstruction::SetSize(G4double value)
173{
174 BoxSize = value;
175}
176
177//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
178
179#include "G4RunManager.hh"
180
181void DetectorConstruction::UpdateGeometry()
182{
183 G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
184}
185
186//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.