Geant4  9.6.p02
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DetectorConstruction.cc
Go to the documentation of this file.
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 //
30 // $Id$
31 //
32 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
33 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
34 
35 #include "DetectorConstruction.hh"
36 #include "DetectorMessenger.hh"
37 
38 #include "G4Material.hh"
39 #include "G4Box.hh"
40 #include "G4LogicalVolume.hh"
41 #include "G4PVPlacement.hh"
42 
43 #include "G4GeometryManager.hh"
44 #include "G4PhysicalVolumeStore.hh"
45 #include "G4LogicalVolumeStore.hh"
46 #include "G4SolidStore.hh"
47 
48 #include "G4UnitsTable.hh"
49 #include "G4PhysicalConstants.hh"
50 #include "G4SystemOfUnits.hh"
51 
52 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
53 
54 DetectorConstruction::DetectorConstruction()
55 :fPBox(0), fLBox(0), fMaterial(0)
56 {
57  fBoxSize = 1*m;
58  DefineMaterials();
59  SetMaterial("Iron");
60  fDetectorMessenger = new DetectorMessenger(this);
61 }
62 
63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
64 
65 DetectorConstruction::~DetectorConstruction()
66 { delete fDetectorMessenger;}
67 
68 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
69 
70 G4VPhysicalVolume* DetectorConstruction::Construct()
71 {
72  return ConstructVolumes();
73 }
74 
75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
76 
77 void DetectorConstruction::DefineMaterials()
78 {
79  //
80  // define Elements
81  //
82  G4double z,a;
83 
84  G4Element* C = new G4Element("Carbon" ,"C" , z= 6., a= 12.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* Ca = new G4Element("Calcium" ,"Ca", z=20., a= 40.08*g/mole);
88  G4Element* H = new G4Element("Hydrogen" ,"H" , z= 1., a= 1.01*g/mole);
89  G4Element* I = new G4Element("Iodine" ,"I" , z= 53., a= 126.90*g/mole);
90 
91  //
92  // define materials
93  //
94  G4double density;
95  G4int ncomponents, natoms;
96  G4double fractionmass;
97 
98  new G4Material("galactic", z=1., a= 1.01*g/mole, universe_mean_density,
99  kStateGas, 2.73*kelvin, 3.e-18*pascal);
100 
101  G4Material* Air =
102  new G4Material("Air", density= 1.290*mg/cm3, ncomponents=2);
103  Air->AddElement(N, fractionmass=70.*perCent);
104  Air->AddElement(O, fractionmass=30.*perCent);
105 
106  G4Material* CaCO3 =
107  new G4Material("CaCO3", density= 2.80*g/cm3, ncomponents=3);
108  CaCO3->AddElement(Ca, natoms= 1);
109  CaCO3->AddElement(C , natoms= 1);
110  CaCO3->AddElement(O , natoms= 3);
111 
112  new G4Material("Carbon", z= 6., a= 12.01*g/mole, density= 2.265*g/cm3);
113  new G4Material("Iron" , z=26., a= 55.85*g/mole, density= 7.870*g/cm3);
114  new G4Material("Tin" , z=50., a= 118.7*g/mole, density= 7.310*g/cm3);
115 
116  G4Material* HI =
117  new G4Material("HI", density= 2.8*g/cm3, ncomponents=2); //-35oC
118  HI->AddElement(H, natoms= 1);
119  HI->AddElement(I, natoms= 1);
120 
121  G4cout << *(G4Material::GetMaterialTable()) << G4endl;
122 }
123 
124 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
125 
126 G4VPhysicalVolume* DetectorConstruction::ConstructVolumes()
127 {
128  // Cleanup old geometry
129  G4GeometryManager::GetInstance()->OpenGeometry();
130  G4PhysicalVolumeStore::GetInstance()->Clean();
131  G4LogicalVolumeStore::GetInstance()->Clean();
132  G4SolidStore::GetInstance()->Clean();
133 
134  G4Box*
135  sBox = new G4Box("Container", //its name
136  fBoxSize/2,fBoxSize/2,fBoxSize/2); //its dimensions
137 
138  fLBox = new G4LogicalVolume(sBox, //its shape
139  fMaterial, //its material
140  fMaterial->GetName()); //its name
141 
142  fPBox = new G4PVPlacement(0, //no rotation
143  G4ThreeVector(), //at (0,0,0)
144  fLBox, //its logical volume
145  fMaterial->GetName(), //its name
146  0, //its mother volume
147  false, //no boolean operation
148  0); //copy number
149 
150  PrintParameters();
151 
152  //always return the root volume
153  //
154  return fPBox;
155 }
156 
157 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
158 
159 void DetectorConstruction::PrintParameters()
160 {
161  G4cout << "\n The Box is " << G4BestUnit(fBoxSize,"Length")
162  << " of " << fMaterial->GetName() << G4endl;
163 }
164 
165 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
166 
167 void DetectorConstruction::SetMaterial(G4String materialChoice)
168 {
169  // search the material by its name
170  G4Material* pttoMaterial = G4Material::GetMaterial(materialChoice);
171  if (pttoMaterial) fMaterial = pttoMaterial;
172 }
173 
174 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
175 
176 void DetectorConstruction::SetSize(G4double value)
177 {
178  fBoxSize = value;
179 }
180 
181 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
182 
183 #include "G4RunManager.hh"
184 
185 void DetectorConstruction::UpdateGeometry()
186 {
187  G4RunManager::GetRunManager()->DefineWorldVolume(ConstructVolumes());
188 }
189 
190 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Generated on Sat May 25 2013 14:31:55 for Geant4 by   doxygen 1.8.4